在sequoia包中运行以下重复检查和亲子关系分配时:
ParOUT <- sequoia(GenoM = Geno, LifeHistData = LH, Err=0.005, Module="ped",
UseAge = "no",
quiet = FALSE, Plot = TRUE)
字符串
我得到以下输出和error(最后两行):
Warning: There are 2111 SNPs scored for <5% of individuals, these will be excluded
Warning: There are 1 monomorphic (fixed) SNPs, these will be excluded
Warning: In addition, there are 52874 SNPs scored for <50% of individuals
Warning: *********** There are 8 individuals scored for <5% of SNPs, these WILL BE IGNORED ***********
After exclusion, There are 40 individuals and 185843 SNPs.
Ageprior: Flat 0/1, overlapping generations, MaxAgeParent = 4,4
Error in if (!is.double(MAF) || any(MAF < 0 | MAF > 1)) { :
missing value where TRUE/FALSE needed
型
什么是**if(!is.double(MAF)||任何(MAF < 0| MAF > 1)){}**检查我的数据中是否可能缺少此内容?
我尝试了上面的方法,并将Module=同时作为“ped”和“par”运行,结果没有什么不同。
1条答案
按热度按时间mec1mxoz1#
这是软件包中的一个错误:如果SNP对零个个体进行评分,则次要等位基因频率(MAF)为,这会破坏检查所有MAF是否在0和1之间。它应该在CRAN上的当前版本中工作。