有2个R程序，运行出现点问题，请帮忙解决，具体加Q：865120364

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

错误于`[.data.frame`(obs, , all.cols) : undefined columns selected

Getting baseline for SF6 ...
错误于axis(side, at = z, labels = labels, lwd = lwd, lwd.ticks = lwd.ticks,  :
  所有的位置值都是无限的

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：R程序 Undefined Baseline Selected getting 程序

这是用什么程序做的？

challenge24 发表于 2012-3-28 12:53
这是用什么程序做的？

R软件！

challenge24 发表于 2012-3-28 12:53
这是用什么程序做的？

在别人机子上能够运行，我这里出现点错误，应该不难 加Q说吧

：####################################################################
# baseline estimation
####################################################################
for (ii in 1:length(para)){
cat("Getting baseline for", para[ii], "...\n")
para.flag = paste(para[ii], "flag", sep="")

if (para.flag %in% names(obs)){
msk = which(!is.na(obs[[para[ii]]]) & obs[[para.flag]] %in% flag)
} else {
msk = which(!is.na(obs[[para[ii]]]))
}
if (length(msk)==0){
obs[[para[ii]]] = NULL
} else {
tmp = summary(factor(round(diff(obs$dtm[msk]),5)), maxsum=2)[1]
# main time step between measurements in days
main.dt = as.numeric(names(tmp))
NoXP = round(width/main.dt)
delta = diff(range(obs$dtm[msk]))/200

fit = rfbaseline(obs$dtm[msk], obs[[para[ii]]][msk], NoXP=NoXP, Scale=Scale, maxit=maxit, delta=delta)

obs[[paste(para[ii], "_bg", sep="")]] = approx(fit$x, fit$fit, xout=obs$dtm)$y
attr(obs[[paste(para[ii], "_bg", sep="")]], "units") = fit$sigma

# add non-random (scale uncertainty)
idx = which(para.table$Parameter==para[ii])
if (paste(id, "_Medusa.u.scale", sep="") %in% names(para.table)){
obs[[paste("unr", para[ii], sep=".")]] = obs[[para[ii]]]*para.table[[paste(id, "_Medusa.u.scale", sep="")]][idx]/100.
} else {
obs[[paste("unr", para[ii], sep=".")]] = obs[[para[ii]]]*para.table[[paste(id, "_GCMD.u.scale", sep="")]][idx]/100.
}

fn.png = file.path(obs.dir, "plots", paste(out.id[1], '.', para[ii], '_time.series.png', sep=''))
bitmap(fn.png, 'png16m', res=300, width=7, height=4, taa=2, gaa=2, pointsize=9)
plot.baseline.fit(obs, para[ii], sigma=fit$sigma, dtm.col="dtm")
dev.off()

}
}

Getting baseline for SF6 ...
错误于axis(side, at = z, labels = labels, lwd = lwd, lwd.ticks = lwd.ticks, :
所有的位置值都是无限的

#sort columns, remove undesired columns
> all.cols = c("dtm", paste(rep(c("", "ur.", "unr.", ""), length(para)), rep(sort(para), each=4), rep(c("", "", "", "_bg"), length(para)), sep=""))
> obs = obs[, all.cols]
错误于[.data.frame(obs, , all.cols) : undefined columns selected
> obs.agg = obs.agg[, all.cols]
错误于[.data.frame(obs.agg, , all.cols) : undefined columns selected