w*****1 发帖数: 473 | 1 Hi,
I am trying to use lme function to find correlation between phenotype and
genotype. My code is as follows:
library(nlme)
> pheno<-read.table('length2.dat',head=T,as.is=T)
> unique.class<-unique(pheno$cnvregion)
for(x in unique.class){
+{subclass<-pheno[pheno$cnvregion==x,]
+ templocus <-subclass$cnvclass
+ fm1 <- lme(pheno~templocus, data=subclass, random=~1|famid, na.action=na.
omit)
+ temp.p <- summary(fm1)$tTable[dim(summary(fm1)$tTable)[1], dim(summary(fm1
)$tTable)[2]]
+ cat(temp.p, "\n", file="p.txt")}
Then here is a error like this:
Error in MEEM(object, conLin, control$niterEM) :
Singularity in backsolve at level 0, block 1
Is there anyone who can tell me what the error mean and how to fix it? I can
give you baozi for it. Thanks a lot! |
|