gmestan:用Stan估算银河系质量 源码

gmestan 该存储库提供了一些示例,这些示例使用Stan来创建多级模型,以使用来自示踪剂的运动学数据来估计银河系的质量。 提供了两个数据集: :球状星团数据集从 :从矮星系的数据集 斯坦代码在 。 有几种可用的模型,每种变体都有两个单独的文件来处理GC数据集和DG数据集。 从最简单到最复杂的顺序: 不包括测量不确定性两款车型中给出和 。 匹配模型 ,用于再现的结果和Slizewski等。 人在准备,是在和分别。 这些模型假定两个适当运动(RA和DEC)之间是独立的,并且结果仅使用完整的数据得出。 两种型号这对于正确的动作一体化协方差矩阵中给出和 。 如果给定的相关性为零,则可以简化为先前的模型。 两个模型中给出其处理为可适当运动丢失的数据(那里有没有可用适当的运动和/或线的视向速度),并纳入协方差矩阵和 。 如果所有对象都具有完整的数据,这将简化为先前的模型。 P