数据帧的csv在这里。
基本上,列模式表示用户可以进入机场的交通选择,选择代表用户的选择,所有其他列描述影响选择的变量。我在下面发布一些代码。该mFormula功能只是分离了变量分为3份分离“|” 即
替代相关w泛型系数| 个别特定的w替代特定系数| 替代的依赖w替代特定系数
我尝试了两个带有mFormula函数的公式,其中一个公式输出系数(f1),但没有为效用函数建模变量(f)。
library(mlogit)
path=""
a<-read.table(path,sep=",", header=TRUE)
Ta<-mlogit.data(a,shape="long",alt.var="Mode",choice="Choice")
f <- mFormula(Choice ~ CostPersonal+CostTrainBus+CostCarpool+CostTaxi+CostSharedTaxi+CostCS+CostSCS+DistSYR+DistNYC+DistALB+DistEWR+SafteyPersonal+SafteyTrainBus+SafteyCarpool+SafteyTaxi+SafteySharedTaxi+SafteyCS+SafteySCS+PrivacyPersonal+PrivacyTrainBus+PrivacyCarpool+PrivacyTaxi+PrivacySharedTaxi+PrivacyCS+PrivacySCS+ConviniencePersonal+ConvinienceTrainBus+ConvinienceCarpool+ConvinienceTaxi+ ConvinienceSharedTaxi+ConvinienceCS+ConvinienceSCS+payUticaSYR+payUticaALB+flightsYear+daysBreturn+incomeB+familySize+ageB|Airport+TimePersonal+TimeTrainBus+TimeCarpool+TimeTaxi+TimeSharedTaxi+TimeCS+TimeSCS+luggage+ownVehicle+shuttleMaxPersonsB)
f1<- mFormula(Choice ~ 1| +0+CostPersonal+CostTrainBus+CostCarpool+CostTaxi+CostSharedTaxi+CostCS+CostSCS+Airport+DistSYR+DistNYC+DistALB+SafteyPersonal+SafteyTrainBus+SafteyCarpool+SafteyTaxi+SafteySharedTaxi+PrivacyPersonal+PrivacyTrainBus+PrivacyCarpool+PrivacyTaxi+PrivacySharedTaxi+ConviniencePersonal+ConvinienceTrainBus+ConvinienceCarpool+ConvinienceTaxi+ConvinienceSharedTaxi+TimePersonal+TimeTrainBus+TimeCarpool+TimeTaxi+TimeSharedTaxi+flightsYear+daysBreturn+luggage+ownVehicle+incomeB+familySize+ageB+shuttleMaxPersonsB+payUticaSYR+payUticaALB)
ml.a<-mlogit(f,Ta)#does not output coefficients
ml.a<-mlogit(f1,Ta)#does output coefficients
summary(ml.a)`
The error is:
Error in solve.default(H, g[!fixed]) :
Lapack routine dgesv: system is exactly singular: U[3,3] = 0
虽然带有f1的mlogit确实用a输出系数
Likelihood: -2.2039e-07
显然,我不明白应该和不应该包括哪些变量,以及哪些变量之间存在共同关系。
发布于 2019-04-08 09:38:23
random是每个选项选择的唯一数字在1-100之间的列
f3<-mFormula(Choice ~ random)
ml.a<-mlogit(f3,Ta)
summary(ml.a)
结果
Error in solve.default(H, g[!fixed]) :
Lapack routine dgesv: system is exactly singular: U[4,4] = 0
https://stackoverflow.com/questions/-100006544
复制相似问题