如何解决numpy.ma.cov-allow_masked-您可以使其双向运行吗?
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 举报,一经查实,本站将立刻删除。