将 for 循环的输出保存到数据框

Question

我正在尝试将函数应用于每行的 netcdf 数据帧。这是脚本的阅读部分。

library(ncdf4)
library(car)
library(forecast)
ncname <- "plus1"             
ncfname <- paste(ncname,".nc", sep="")
dname <- "precip"
ncin <- nc_open(ncfname)
lon <- ncvar_get(ncin,"lon")
nlon <- dim(lon)
lat <- ncvar_get(ncin,"lat",verbose=F)
nlat <- dim(lat)
t <- ncvar_get(ncin,"time")
tunits <- ncatt_get(ncin,"time","units")
nt <- dim(t)
tmp_array <- ncvar_get(ncin,dname)
nc_close(ncin)
tmp_vec_long <- as.vector(tmp_array)
tmp_mat <- matrix(tmp_vec_long, nrow=nlon*nlat, ncol=nt)
lonlat <- as.matrix(expand.grid(lon,lat))
tmp_df02 <- data.frame(cbind(lonlat,tmp_mat))
tmp_df03<-na.omit(tmp_df02)

tmp_df03$lamb<-apply(tmp_df02[3:365],1,function(x)BoxCox.lambda(x,method="loglik",lower=-2,upper=2))

我想使用 tmp_df03$lamb 列中的值并应用于数据框的每一行。例如（x应该是数据框的名字）：

trans<-yjPower(x,tmp_df03$lamb,jacobian.adjusted=TRUE)

如果我把它放在一个 for 循环中：

d<-matrix(ncol=363,nrow=899)
for (i in 1:899){
trans[i]<-yjPower(as.numeric(tmp_df03[i,][3:365]),tmp_df03$lamb[i],jacobian.adjusted=TRUE)
d<-rbind(d,trans) 
}

我遇到了以下警告信息。

警告消息：

1: In trans[i] <- yjPower(as.numeric(tmp_df03[i, ][3:365]), tmp_df03$lamb[i],  :

要替换的项目数不是替换长度的倍数

问题：

如何组合 for 循环的输出并创建结果的新数据框？
谁能建议解决这个问题的方法？

我将不胜感激任何帮助。非常感谢，

Answer 1

使用 sapply() - 我不太了解 yjPower() ，所以只是给你一个提示！

vec <- sapply(1:nrow(x), function(i) yjPower(x[i,], tmp_df03$lamb[i], jacobian.adjusted = TRUE))

所以上面的代码将给出一个向量：将其作为 data.frame，在这个向量上使用 as.data.frame()。

如果你想加入这个 "x" 然后使用 cbind(x, as.data.frame(vec))

将 for 循环的输出保存到数据框

Saving the output of a for loop to a data frame

r

for-loop

dataframe

netcdf