通過循環中的數據幀行追加值

-2

我正在運行具有各種初始值的模型，並且我試圖按行向循環中的數據框添加值（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.

來源

2017-09-26 Kay

請註明maxLik屬於庫 –

1.做到這一點全部在列表中，然後do.call（rbind，mylist），2.我寧願用RUN-INDEX創建一個新列，而不是自己的行。在錯誤函數中，我會記錄失敗的運行。無論如何，你可以擺脫它。 –

也請提一下f是什麼？ –

我運行一個例子，通過改變函數f

library("maxLik") 

t <- rexp(100, 2) 
loglik <- function(theta) log(theta) - theta*t 


df <- data.frame(alpha = numeric(), rho = numeric(), lambda = numeric()) 

for (j in 1:20){ 
    tryCatch({ 
    ml <- maxLik(loglik, start = c(alpha = runif(1, 0, 2), rho = runif(1, 0, 4), 
           lambda = runif(1, 0, 10)), method = "NM") 

    df <- rbind(df, data.frame(alpha = ml$estimate[1], 
           rho = ml$estimate[2], 
           lambda = ml$estimate[3])) 
    # I tried to append values for each column   
    }, error = function(e) {NA})} 

> row.names(df) <- NULL 
> head(df) 
    alpha  rho lambda 
1 2.368739 2.322220 2.007375 
2 2.367607 2.322328 2.007093 
3 2.368324 2.322105 2.007597 
4 2.368515 2.322072 2.007334 
5 2.368269 2.322071 2.007142 
6 2.367998 2.322438 2.007391

來源

2017-09-26 09:45:12

謝謝，數據框仍然沒有從代碼中返回任何東西。我注意到的一件事是你的代碼中沒有索引'j'，所以這不是迭代的。 – Kay

謝謝，這是'f'問題？ – Kay

可能是，請檢查您的f函數，否則數據框將如上所示追加 –

通過循環中的數據幀行追加值

回答

相關問題