data<-read.csv("D:/课题相关/连续变量插补.csv")
DBIL<-subset(data,select = c(ID,survtime,DBIL,TBIL))
print(str(DBIL))
DBIL.bds <- matrix(c(0,3.9,87.7), nrow = 1, ncol = 3)
a.DBIL<-amelia(DBIL,m=5,ts="survtime",cs="ID",lags="DBIL",leads="DBIL",polytime=2,bounds = DBIL.bds, resample=1000)
运行上述代码出现报错:Amelia Error Code: 50 The bounds argument is a three-column matrix.
在Amelia函数中指定的所有列名都存在于我的数据集中:
且同样一段代码运行另一变量时就能成功:
data<-read.csv("D:/课题相关/连续变量插补.csv")
ALT<-subset(data,select = c(ID,survtime,TBIL,ALT))
ALT.bds <- matrix(c(2,22,500), nrow = 1, ncol = 3)
a.ALT<-amelia(ALT,m=5,ts="survtime",cs="ID",lags="ALT",leads="ALT",polytime=2,bounds = ALT.bds, resample=1000)
summary(a.ALT)
tscsPlot(a.ALT,cs=c(1,5,6,9),main="Combined Time and Lags/Leads",var="ALT")
write.csv(a.ALT$imputations[[1]], file = "ALT1.csv")
请问有无大神可以指导一下这个问题,谢谢!