在计算两个trait之间的genetic correaltion 的时候，如果是用ldsc的脚本计算的话，有很大情况会出现NA值，一开始也很疑惑，后来搜集了一些资料，对该情况进行了解释。
（1）出现NA的情况，第一种情况是，两种traits，一种或者两种的trait的遗传力太低，导致在计算两个trait之间的方差，基于方差计算rg的时候，出现NA的现象
（2） 第二种情况是h2的SE太高。解决方法，可以检查与基线LD评分重叠的SNPs数量是否足够高(整个基因组至少需要200k SNPs，这个可以直接在log文件中查看，如果不够200k的SNP，程序会warning，check一下是否有这个警告)。但是如果重叠是高,还是这样一个低h2,你可以试着重新计算LD分数使用千人基因组3期的基线LDSC提供的注释生成一组新的基线和只使用单核苷酸多态性数据计算相应的文件。这可能有助于稍微提高功率。如果这不起作用，那么很可能是你的GWAS没有足够的power。
问题描述 在计算两个trait之间的genetic correaltion 的时候，如果是用ldsc的脚本计算的话...
问题 在R中，如何用非NA值替换向量或因素型变量中的NA值？ 指南 下列代码将会告诉我们如何填补向量中的空缺，如果...
有时候处理Excel数据时会出现很多NA值，比如将多行数据合并的时候，导入R里面就会发现NA值，对于NA值的填充有...
1、缺失数据 缺失值 NA，在 R 中，NA 代表缺失值，NA是不可用，not available 的简称，用来存...
NA值处理 数据NA值查看 NA值处理 直接删掉含NA的那行 用某值将其替换掉，如0 通过探索相关性将NA值填补上...
缺失信息 R中用NA表示缺失信息。往往很多数据中包含NA值，需要去除。 查看数据集中是否含义缺失值(is.na()...
删除缺失值 R中缺失值使用NA表示。NA是R中的三个逻辑值之一，另外两个是TRUE和FALSE。NULL表示无效对...
2.7NA与NULL值 NA：缺失值向量中有一个缺失值NA，会导致mean()等函数计算出的结果均为NA。但可将n...
举例1：向量类型判断缺失值is.na和缺失值的填补which (x<-c(1,2,3,NA)) is.na(x) ...
在R中缺失值是使用“NA”表示的。可以通过函数is.na()来判断一个对象是否是缺失值。在处理数据时，经常就会遇到...
本文标题：Genetic correaltion 中出现NA值的问题
本文链接：https://www.haomeiwen.com/subject/ajhgmltx.html
网友评论