我正在运行具有各种初始值的模型,并且我试图按行向循环中的数据框添加值(3个估计值)。我给循环内的估计器赋值,但我不记得它们产生一个数据帧。通过循环中的数据帧行追加值
我的代码:f
是估算的模型。三个参数:模型中的alpha,rho和lambda。我想输出这3个值。发生
library("maxLik")
f <- function(param) {
alpha <- param[1]
rho <- param[2]
lambda <- param[3]
u <- 0.5 * (dataset$v_50_1)^alpha - 0.5 * lambda * (dataset$v_50_2)^alpha
p <- 1/(1 + exp(-rho * u))
logl <- sum(dataset$gamble * log(p) + (1 - dataset$gamble) * log(1 - p))
}
df <- data.frame(alpha = numeric(), rho = numeric(), lambda = numeric())
for (j in 1:20) {
tryCatch({
ml <- maxLik(f, start = c(alpha = runif(1, 0, 2), rho = runif(1, 0, 4), lambda = runif(1,
0, 10)), method = "NM")
alpha[j] <- ml$estimate[1]
rho[j] <- ml$estimate[2]
lambda[j] <- ml$estimate[3]
}, error = function(e) {NA})
}
output <- data.frame(alpha, rho, lambda)
错误:
Error in data.frame(alpha, rho, lambda) : object 'alpha' not found
预计输出
alpha rho lambda
0.4 1 2 # estimators append by row.
0.6 1.1 3 # each row has estimators that are estimated
0.7 1.5 4 # by one set of initial values, there are 20
# rows, as the estimation loops for 20 times.
请注明maxLik属于库 –
1.做到这一点全部在列表中,然后do.call(rbind,mylist),2.我宁愿用RUN-INDEX创建一个新列,而不是自己的行。在错误函数中,我会记录失败的运行。无论如何,你可以摆脱它。 –
也请提一下f是什么? –