首页 > 解决方案 > Apollo R:apollo_estimate 中的错误(apollo_beta,apollo_fixed,apollo_probabilities,:未进行估计

问题描述

当使用 Apollo 包估计二元选择模型时,我收到错误:“apollo_estimate(apollo_beta, apollo_fixed, apollo_probabilities, : 不进行估计!”我检查了我的代码并查看了一些论坛以尝试找出我的错误在哪里是,但我似乎找不到它。我在下面发布了代码,希望有人能帮助我。提前谢谢你!

-伊里尼

'''

rm(list=ls())
install.packages("apollo")
library(apollo)
database <- read.csv2("Dutch_daten.csv")
str(database)
database$cost <- as.numeric(as.character(database$cost))
database$ivtt <- as.numeric(as.character(database$ivtt))
database$rpovt <- as.numeric(as.character(database$rpovt))

### Initialise code
apollo_initialise()

### Set core controls
apollo_control= list(
  modelName="ApolloNetherlands_Exc",
  modelDescr="BinaryChoiceModel",
  indivID = "idx.alt"
)

### Define model parameters for specification 0

apollo_beta <- c(b0_rail =0,
                 b1_cost = 0,
                 b2_time = 0)
apollo_fixed = c()

### GROUP AND VALIDATE INPUTS
apollo_inputs = apollo_validateInputs()

### DEFINE MODEL AND LIKELIHOOD FUNCTION
apollo_probabilities=function(apollo_beta, apollo_inputs, functionality="estimate"){
  
  
  apollo_attach(apollo_beta, apollo_inputs)
  on.exit(apollo_detach(apollo_beta, apollo_inputs))
  
  
  P = list()
  
  
  V <- list()
  
  V[['rail']] <- b0_rail + b1_cost*cost+ b2_time*(ivtt+rpovt)
  V[['car']] <- b1_cost*cost+ b2_time*(ivtt+rpovt)
  
  mnl_settings = list(
    alternatives = c(car=1, rail=2),
    avail = list(car=1, rail=1),
    choiceVar = idx.alt,
    V = V
  )
  
  P[["model"]]=apollo_mnl(mnl_settings, functionality)
  
  
  P = apollo_panelProd(P, apollo_inputs, functionality)
  
  P = apollo_prepareProb(P, apollo_inputs, functionality)
  return(P)
}

# ################################################################# #
#### MODEL ESTIMATION                                            ####
# ################################################################# #

model = apollo_estimate(apollo_beta, apollo_fixed, apollo_probabilities, apollo_inputs)

标签: rmodelbinaryapolloestimation

解决方案


推荐阅读