numpy.ma.cov-allow_masked-您可以使其双向运行吗?

问题描述

numpy文档对ma.cov函数说以下话:

allow_maskedbool 可选

如果True,则将屏蔽值成对传播:如果在x中屏蔽一个值,则在y中屏蔽了相应的值。

我想知道您是否可以同时使用它(即,如果某个值在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 (将#修改为@)