R中拟合模型的时候出现下面的情况,请问大神们怎么破

> Hawaii=read.table("Hawaii.txt",header = TRUE) > Birds=c(Hawaii$Stilt.Oahu,Hawaii$Stilt.Maui, + Hawaii$Coot.Oahu,Hawaii$Coot.Maui) > Time=rep(Hawaii$Year,4) > Rain=rep(Hawaii$Rainfall,4) > ID=factor(rep(c("Stilt.Oahu","Stilt.Maui","Coot.Oahu","Coot.Maui"), + each=length(Hawaii$Year))) > library(lattice) > xyplot(Birds~Time|ID,col=1) > library(mgcv) > BM1=gamm(Birds~Rain+ID+ + s(Time,by=as.numeric(ID=="Stilt.Oahu"))+ + s(Time,by=as.numeric(ID=="Stilt.Maui"))+ + s(Time,by=as.numeric(ID=="Coot.Oahu"))+ + s(Time,by=as.numeric(ID=="Coot.Maui")), + na.action=na.omit, + weights = varIdent(form = ~1|ID)) Error in MEestimate(lmeSt, grps) : Singularity in backsolve at level 0, block 1

> Hawaii=read.table("Hawaii.txt",header = TRUE) > Birds=c(Hawaii$Stilt.Oahu,Hawaii$Stilt.Maui, Hawaii$Coot.Oahu,Hawaii$Coot.Maui) > Time=rep(Hawaii$Year,4) > Rain=rep(Hawaii$Rainfall,4) > ID=factor(rep(c("Stilt.Oahu","Stilt.Maui","Coot.Oahu","Coot.Maui"), each=length(Hawaii$Year))) > library(lattice) > xyplot(Birds~Time|ID,col=1) > library(mgcv) > BM1=gamm(Birds~Rain+ID+ s(Time,by=as.numeric(ID=="Stilt.Oahu"))+ s(Time,by=as.numeric(ID=="Stilt.Maui"))+ s(Time,by=as.numeric(ID=="Coot.Oahu"))+ s(Time,by=as.numeric(ID=="Coot.Maui")), na.action=na.omit, weights = varIdent(form = ~1|ID)) 出现这样的错误 Error in MEestimate(lmeSt, grps) : Singularity in backsolve at level 0, block 1

你好,我是有问必答小助手。为了技术专家团更好地为您解答问题,烦请您补充下(1)问题背景详情,(2)您想解决的具体问题,(3)问题相关代码图片或者报错信息。便于技术专家团更好地理解问题,并给出解决方案。

您可以点击问题下方的【编辑】,进行补充修改问题。