问题描述
rm(list = ls())
(可选)设置工作目录
setwd("/Users/i/DropBox /")
mydata<-read.csv("total.csv")
加载相关库
library(tidyverse)
library(apollo)
library(haven)
library(dplyr)
必填步骤
apollo_initialise()
----------------------------------------------- ---------------------------------
步骤 1-2 - 设置 Apollo 控件apollo_control = list(modelName ="Model_ORL",modelDescr ="ORL_final ",indivID ="X")
tbl <- mydata
database <- as.data.frame(tbl)
读取SPSS文件
tbl <- database %>% mutate(tripsC = case_when(
tripslevel == 1 ~ "one",tripslevel == 2 ~ "two",tripslevel == 3 ~ "three"),trips_factor = factor(tripsC,ordered = TRUE,levels = c("one","two","three")))
apollo_beta = c(b_temp = 0,b_humidity =0,b_windspeed =0,b_precipitation =0,b_gasprice =0,tau_1 = 1,tau_2 = 2)
apollo_fixed = c()
----------------------------------------------- ---------------------------------
第 4 步 - 验证数据
apollo_inputs = apollo_validateInputs()
----------------------------------------------- ---------------------------------
第 5 步 - 定义阿波罗概率
apollo_probabilities=function(apollo_beta,apollo_inputs,functionality="estimate"){
在函数退出后附加输入并分离
apollo_attach(apollo_beta,apollo_inputs)
on.exit(apollo_detach(apollo_beta,apollo_inputs))
在这里定义效用方程 *************************************
没有常数的效用/倾向方程
V = b_temp * Temp + b_humidity * Humidity + b_windspeed * Windspeed+ b_precipitation*Precipitation+ b_gasprice*Gasprice
P = list()
定义 ORL 模型组件的设置
ol_settings = list(
outcomeOrdered = tripslevel,V = V,tau = c(tau_1,tau_2),coding = c(0,1 ) )
----------------------------------------------- ------------------------------
使用 MNL 模型计算概率
P[['model']] = apollo_ol(ol_settings,functionality)
准备并返回函数的输出
P = apollo_prepareProb(P,functionality)
return(P)}
可选:在模型估计之前计算LL
apollo_llCalc(apollo_beta,apollo_probabilities,apollo_inputs)
----------------------------------------------- ---------------------------------
第 6 步 - 估计模型
model = apollo_estimate(apollo_beta,apollo_fixed,apollo_inputs)
- apollo_preprocess 中的错误(输入 = ol_settings,modelType,function,: 模型组件“model”的“outcomeOrdered”参数需要是一个标量或向量,“数据库”中的每个观察值都有一个条目
- 我一直在寻找相同的错误消息,但找不到解决方案。
解决方法
暂无找到可以解决该程序问题的有效方法,小编努力寻找整理中!
如果你已经找到好的解决方法,欢迎将解决方案带上本链接一起发送给小编。
小编邮箱:dio#foxmail.com (将#修改为@)