R语言amelia函数报错

用户在使用Amelia函数进行连续变量插补时,遇到DBIL变量的bounds参数报错,指出bounds矩阵不正确。对比ALT变量的成功案例,问题可能出在DBIL的bounds矩阵格式上,需要检查其列数和范围是否符合要求。
摘要由CSDN通过智能技术生成
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")

请问有无大神可以指导一下这个问题,谢谢!

评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值