Question

我已经查看了一些已在此处发布的答案，但却找不到能够帮助解决我的特殊情况的答案。我是R的新人，所以请耐心等待。

    rm(list = ls()) 
dat=read.csv('C:\\Users\\Casandra\\Downloads\\roaches.csv',as.is=T) 
ndata=nrow(dat) 
param=read.csv('C:\\Users\\Casandra\\Downloads\\roaches+posterior+after+burnin.csv',as.is=T) 

param.med=apply(param,2,mean) 

yhat=param.med[1]+param.med[2]*dat$x 

png('pred distribution.png') 
plot(dat$x,dat$y,xlim=c(0,1),ylim=c(0,1.5),xlab='Covariate x',col='grey',ylab='') 
lines(dat$x,yhat,col='red') 
nsim=nrow(param) 
yhat1=yhat2=matrix(NA,ndata,3) 
for (i in 1:ndata){ 
 media=param[,1]+param[,2]*dat$x[i] 
 yhat1[i,]=quantile(media,c(0.025,0.5,0.975)) 
 tmp=rnorm(nsim,mean=media,sd=sqrt(param[,3])) 
 yhat2[i,]=quantile(tmp,c(0.025,0.5,0.975)) 
} 
lines(dat$x,yhat1[,1],col='orange',lty=2) 
lines(dat$x,yhat1[,3],col='orange',lty=2) 
lines(dat$x,yhat2[,1],col='grey',lty=2) 
lines(dat$x,yhat2[,3],col='grey',lty=2) 
dev.off()

以上是我的代码，这是我收到的错误

Error in [{data.frame {1}}

(param, , 3) : undefined columns selected数据集有2列，dat数据集也有2列。

有人能指出我正确的方向来修复我的代码吗？

数据集如下：

dat：https://drive.google.com/file/d/0B9MIgQ2O0SHnakRzZ2p2bDhIZ2c/edit?usp=sharing

param：https://drive.google.com/file/d/0B9MIgQ2O0SHnTTVMU0E2TTRDR2M/edit?usp=sharing

选择了未定义的列

0 个答案: