因此,我刚开始使用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?
任何建议或解释将不胜感激!