问题描述
我正在尝试运行保存for循环结果的代码。每次运行for循环时,我都希望将Landa和dlanda的值保存在一个名为best的数据帧中,其中两列分别用于landa和dlanda,但我不知道如何创建此数据帧(称为optimize)在循环内。这是我的代码:
Mean_Precipitation_July$mean_precip_monthYear
landa <- seq(0,1,by=0.1)
landa
(i <- 0)
(x <- 1)
for(i in landa) {
T1 <- -(Mean_Precipitation_July$mean_precip_monthYear)^(landa)
dlanda <- (mean(T1)-median(T1))/var(T1)
optimal[x,1] <- landa
optimal[x,2] <- dlanda
i <- i + 0.1
x<- x + 1
}
解决方法
这是一种方法。我使用了伪造的测试数据集,因为OP尚未发布一个。
测试数据向量被命名为data_vec
。
set.seed(2020)
data_vec <- runif(20,2,5)
landa <- seq(0,1,by = 0.1)
optimal <- matrix(nrow = length(landa),ncol = 2)
optimal <- as.data.frame(optimal)
names(optimal) <- c("landa","dlanda")
for(i in seq_along(landa)) {
T1 <- -(data_vec)^landa[i]
dlanda <- (mean(T1) - median(T1))/var(T1)
optimal$landa[i] <- landa[i]
optimal$dlanda[i] <- dlanda
}
optimal
# landa dlanda
#1 0.0 NaN
#2 0.1 7.0957224
#3 0.2 3.0064453
#4 0.3 1.6915265
#5 0.4 1.0659035
#6 0.5 0.7129269
#7 0.6 0.4939882
#8 0.7 0.3498969
#9 0.8 0.2512264
#10 0.9 0.1817638
#11 1.0 0.1318914