SEQUOIA(R)-运行父级分配步骤sequoia()时出错

o0lyfsai  于 12个月前  发布在  其他
关注(0)|答案(1)|浏览(100)

在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”运行,结果没有什么不同。

mec1mxoz

mec1mxoz1#

这是软件包中的一个错误:如果SNP对零个个体进行评分,则次要等位基因频率(MAF)为,这会破坏检查所有MAF是否在0和1之间。它应该在CRAN上的当前版本中工作。

相关问题