MEEM中的错误（object，conLin，control $ niterEM）：在0级，第1区的backsolve中的奇点

Question

我正在尝试在nlme包中使用线性混合效果模型。我的数据是：

x=1:100
y.1=rbeta(x,2,3)
y.2=rbeta(x,2,3)
y.3=rbeta(x,2,3)
y.4=rbeta(x,2,3)
y.5=rbeta(x,2,3)
x=t(rbind(rep(x,5)))
X=cbind(rep(1,length(x)),x,x,x,x,x)
colnames(X)=c("intercept","X1","X2","X3","X4","X5")
y=rbind(y.1,y.2,y.3,y.4,y.5)
y=t(y)
y=as.data.frame(y)
y=data.frame(stack(y[1:ncol(y)]))
y=as.matrix(y[1])


n=500    
dummy <- rep(1, n)
dimnames(y)[[1]] <- 1:n
tendon <- groupedData(y~X|dummy,data=data.frame(X, y))
colnames(tendon)[7]="y"
fit=lme(y~-1+X,data=tendon,method = "REML")

然后出现上述错误。任何人都可以对此有所了解。