问题描述
numpy文档对ma.cov函数说以下话:
allow_maskedbool ,可选
我想知道您是否可以同时使用它(即,如果某个值在y中被屏蔽,则相应的值在x中被屏蔽)。
我有一个二维numpy数组,其中包含多个总体中一组基因座的等位基因频率(初始化为site_freq = np.zeros((num_Loci,num_pops),dtype=np.float64)
)
我想计算群体中等位基因频率的方差(仅使用频率> 0的基因座)和群体之间等位基因频率的协方差(仅使用等位基因在两个群体中均分离的基因座)。
masked_freq=ma.masked_equal(site_freq,0.0)
ma.var(masked_freq,axis=0)
似乎可以正常工作,但是
ma.cov(masked_freq,rowvar=False,allow_masked=True)
我不太确定。如果对于一个总体(称为x
),某个基因座没有被掩盖并包含在方差计算中,那么我如何确定如果该基因座在另一种群中被掩盖,则该基因座将不包含在协方差计算中(称为y
)。
解决方法
暂无找到可以解决该程序问题的有效方法,小编努力寻找整理中!
如果你已经找到好的解决方法,欢迎将解决方案带上本链接一起发送给小编。
小编邮箱:dio#foxmail.com (将#修改为@)