2

我有一个看起来像这样的数据。我下面的代码只是计算一些值并将输出向量绑定到原始数​​据帧。

options(width=200)

args<-commandArgs(trailingOnly=FALSE)
dat <- read.table("http://dpaste.com/89376/plain/",fill=T);

problist <- c();

for (lmer in 1:10) {
   meanl <- lmer;
   stdevl <- (0.17*sqrt(lmer));
   err_prob <- pnorm(dat$V4,mean=meanl, sd=stdevl);
   problist <- cbind(problist,err_prob);
}

dat <- cbind(dat,problist)
#print(dat,row.names=F, column.names=F,justify=left)

# Why this breaks?
write(dat, file="output.txt", sep="\t",append=F);

关于上述问题，我有几个问题：

  1. 但是为什么上面的'write()'函数会给出这个错误。有没有办法解决它？

    cat(list(...)、file、sep、fill、labels、append) 中的错误：参数 1（类型 'list'）不能由 'cat' 处理调用：write -> cat 执行停止

  2. 数据框中绑定向量的名称被添加为所有 10 个新列的“errprob”。有没有办法将它们命名为“errprob1”、“errprob2”等？

4

2 回答 2

4

首先，不需要分号，R 知道行尾是一个中断。


for (lmer in 1:10){
    meanl <- lmer
    stdevl <- (0.17*sqrt(lmer))
    err_prob <- pnorm(dat$V4,mean=meanl, sd=stdevl)
    problist <- cbind(problist,err_prob)
}
colnames(problist)<-paste("errorprob",1:10,sep="")
dat <- cbind(dat,problist)
write.table(dat, file="output.txt", sep="\t",append=F)

  1. 我相信您正在寻找 write.table 功能

  2. 使用 colnames 函数

于 2009-09-04T05:02:16.743 回答
2

  1. 您可以使用 write.table() 而不是 write() 来使用上面指定的参数。后者适用于打印矩阵（但可能需要指定 ncol 或转置输入矩阵），但前者更通用，我将它用于矩阵和数据帧。

  2. 你可以更换

    err_prob <- pnorm(dat$V4,mean=meanl, sd=stdevl)

    problist <- cbind(problist,err_prob)

和

assign(sprintf("err_prob%d",lmer),pnorm(dat$V4,mean=meanl, sd=stdevl))
problist <- eval(parse(text=sprintf("cbind(problist,err_prob%d)", lmer)))

最后一行将字符串解析为表达式，然后对其求值。你也可以做

colnames(problist) <- sprintf("err_prob%d",1:10)

后验的

于 2009-09-04T04:50:48.967 回答