Apollo R: ошибка в apollo_estimate (apollo_beta, apollo_fixed, apollo_probabilities,: оценка не выполняется - PullRequest
Новая
       88

Apollo R: ошибка в apollo_estimate (apollo_beta, apollo_fixed, apollo_probabilities,: оценка не выполняется

1 голос
/ 03 августа 2020

при использовании пакета Apollo для оценки модели бинарного выбора я получаю сообщение об ошибке: «apollo_estimate (apollo_beta, apollo_fixed, apollo_probabilities,: Не продолжать оценку!» Я проверил свой код и просмотрел несколько форумов, чтобы попытаться найти где моя ошибка, но я не могу ее найти. Я разместил приведенный ниже код, надеясь, что кто-то может мне помочь. Заранее спасибо !!

-Irini

'' '

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)
...