用lme's运行分段SEM

时间:2019-11-26 21:54:09

标签: r nlme piecewise semplot

因此,我刚开始使用psem,但遇到了一个奇怪的问题。每当我在刚运行的summary()上运行psem函数时,它都会加载大约60%的代码,然后出现以下错误:r in solve.default(estimates[dimE[1] - (p:1), dimE[2] - (p:1), drop = FALSE]) : system is computationally singular: reciprocal condition number = 4.92235e-17 ...关于正在发生的事情的任何想法这里?每当我自己运行lme所用的psem时,我就没有问题,只是当我将它们插入psem时才发生此错误。甚至更陌生的是,从lmes的输入中删除一个或两个psem将使其成功运行...请参见下面的示例:

lme(A~1+2+3+4,random=~1|5,data=EX) Works fine 
lme(B~1+2+3+4,random=~1|5,data=EX) Works fine 
lme(C~A+B,random=~1|5,data=EX) Works fine 
lme(D~A+B,random=~1|5,data=EX) Works fine
lme(E~C+D,random=~1|5,data=EX) Works fine 

Trialpsem<-psem( 
lme(A~1+2+3+4,random=~1|5,data=EX)
lme(B~1+2+3+4,random=~1|5,data=EX)
lme(C~A+B,random=~1|5,data=EX)
lme(D~A+B,random=~1|5,data=EX)
lme(E~C+D,random=~1|5,data=EX)

This receives error code 

Trialpsem<-psem( 
lme(C~A+B,random=~1|5,data=EX)
lme(D~A+B,random=~1|5,data=EX)
lme(E~C+D,random=~1|5,data=EX)

This works?

任何建议或解释将不胜感激!

0 个答案:

没有答案