r - R 包 gmnl 和 Rchoice 中可能存在错误？-6ren

r - R 包 gmnl 和 Rchoice 中可能存在错误？

转载作者：行者123 更新时间：2023-12-04 15:53:51

我想估计 gmnl 中的二进制随机参数 logit，同时考虑未观察到和观察到的异质性。当 Gender 是一个具有 2 个级别 (0,1) 的因素时，它会导致以下错误。当 Gender 以与数值变量完全相同的结构进入 gmnl 函数时，公式运行正常。有没有人知道一个错误，其中没有因子变量(例如虚拟变量)可以在 gmnl 函数中输入公式作为社会人口特征，而只能输入数字？

类似的R包Rchoice也是如此。

有没有可能没有错误，但公式有其他我无法发现的错误？

MIXL_socdem_IX <- gmnl(formula = route~ cost + total_time |1 |0 |Gender-1, data = dataIX_MIXL, model = "mixl", ranp=c(cost="ln", total_time="n"), mvar = list(cost=c("Gender"), total_time=c("Gender")), panel=T, R=200, haltons = list("primes"= c(3,5), drop= rep(10,2)), method="bfgs", print.init = T)

错误是:

Error in gmnl(formula = route ~ cost + total_time | 1 | 0 | Gender - 1, : The following variables are not specified in the formula: Gender

谢谢R 版本 3.5.1 (2018-07-02) -- “羽毛喷雾”

      ind     alt route  cost total_time chid Gender Freq_2.4w
1.A     1 Egnatia FALSE -11.0          9    1      0         0
1.B     1     PEO  TRUE -12.0          9    1      0         0
2.A     1 Egnatia FALSE -11.0          9    2      0         0
2.B     1     PEO  TRUE -12.0         10    2      0         0
3.A     1 Egnatia  TRUE -10.5          9    3      0         0
3.B     1     PEO FALSE -12.0          9    3      0         0
4.A     1 Egnatia  TRUE -10.5          9    4      0         0
4.B     1     PEO FALSE -12.0         10    4      0         0
5.A     1 Egnatia  TRUE -10.0          9    5      0         0
5.B     1     PEO FALSE -12.0          9    5      0         0
6.A     1 Egnatia FALSE -11.0          9    6      0         0
6.B     1     PEO  TRUE -12.0         11    6      0         0
7.A     1 Egnatia  TRUE -10.0          9    7      0         0
7.B     1     PEO FALSE -12.0         10    7      0         0
8.A     1 Egnatia  TRUE -10.5          9    8      0         0
8.B     1     PEO FALSE -12.0         11    8      0         0
9.A     1 Egnatia  TRUE -10.0          9    9      0         0
9.B     1     PEO FALSE -12.0         11    9      0         0
10.A    2 Egnatia FALSE -13.5          5   10      0         1
10.B    2     PEO  TRUE  -8.0          6   10      0         1
11.A    2 Egnatia FALSE -13.5          5   11      0         1
11.B    2     PEO  TRUE  -8.0          7   11      0         1
12.A    2 Egnatia FALSE -13.0          5   12      0         1
12.B    2     PEO  TRUE  -8.0          6   12      0         1
13.A    2 Egnatia FALSE -13.0          5   13      0         1
13.B    2     PEO  TRUE  -8.0          7   13      0         1
14.A    2 Egnatia FALSE -12.5          5   14      0         1
14.B    2     PEO  TRUE  -8.0          6   14      0         1
15.A    2 Egnatia FALSE -13.5          5   15      0         1
15.B    2     PEO  TRUE  -8.0          8   15      0         1
16.A    2 Egnatia FALSE -12.5          5   16      0         1
16.B    2     PEO  TRUE  -8.0          7   16      0         1
17.A    2 Egnatia FALSE -13.0          5   17      0         1
17.B    2     PEO  TRUE  -8.0          8   17      0         1
18.A    2 Egnatia FALSE -12.5          5   18      0         1
18.B    2     PEO  TRUE  -8.0          8   18      0         1
19.A    3 Egnatia FALSE -13.0         12   19      0         0
19.B    3     PEO  TRUE -10.0          9   19      0         0
20.A    3 Egnatia FALSE -13.0         12   20      0         0
20.B    3     PEO  TRUE -10.0         10   20      0         0
21.A    3 Egnatia  TRUE -12.5         12   21      0         0
21.B    3     PEO FALSE -10.0          9   21      0         0
22.A    3 Egnatia  TRUE -12.5         12   22      0         0
22.B    3     PEO FALSE -10.0         10   22      0         0
23.A    3 Egnatia FALSE -12.0         12   23      0         0
23.B    3     PEO  TRUE -10.0          9   23      0         0
24.A    3 Egnatia FALSE -13.0         12   24      0         0
24.B    3     PEO  TRUE -10.0         11   24      0         0
25.A    3 Egnatia  TRUE -12.0         12   25      0         0
25.B    3     PEO FALSE -10.0         10   25      0         0
26.A    3 Egnatia  TRUE -12.5         12   26      0         0
26.B    3     PEO FALSE -10.0         11   26      0         0
27.A    3 Egnatia  TRUE -12.0         12   27      0         0
27.B    3     PEO FALSE -10.0         11   27      0         0
28.A    4 Egnatia FALSE -16.0          9   28      0         1
28.B    4     PEO  TRUE -10.0          9   28      0         1
29.A    4 Egnatia FALSE -16.0          9   29      0         1
29.B    4     PEO  TRUE -10.0         10   29      0         1
30.A    4 Egnatia FALSE -15.5          9   30      0         1
30.B    4     PEO  TRUE -10.0          9   30      0         1
31.A    4 Egnatia FALSE -15.5          9   31      0         1
31.B    4     PEO  TRUE -10.0         10   31      0         1
610.A  69 Egnatia  TRUE -12.0          9  610      0         0
610.B  69     PEO FALSE -11.0         10  610      0         0
611.A  70 Egnatia FALSE -11.0          9  611      1         0
611.B  70     PEO  TRUE -10.0         10  611      1         0
612.A  70 Egnatia FALSE -11.0          9  612      1         0
612.B  70     PEO  TRUE -10.0         11  612      1         0
613.A  70 Egnatia  TRUE -10.5          9  613      1         0
613.B  70     PEO FALSE -10.0         10  613      1         0
614.A  70 Egnatia  TRUE -10.5          9  614      1         0
614.B  70     PEO FALSE -10.0         11  614      1         0
615.A  70 Egnatia  TRUE -10.0          9  615      1         0
615.B  70     PEO FALSE -10.0         10  615      1         0
616.A  70 Egnatia FALSE -11.0          9  616      1         0
616.B  70     PEO  TRUE -10.0         12  616      1         0
617.A  70 Egnatia  TRUE -10.0          9  617      1         0
617.B  70     PEO FALSE -10.0         11  617      1         0
618.A  70 Egnatia  TRUE -10.5          9  618      1         0
618.B  70     PEO FALSE -10.0         12  618      1         0
619.A  70 Egnatia  TRUE -10.0          9  619      1         0
619.B  70     PEO FALSE -10.0         12  619      1         0
620.A  71 Egnatia FALSE -11.0          9  620      0         0
620.B  71     PEO  TRUE  -8.0          9  620      0         0
621.A  71 Egnatia FALSE -11.0          9  621      0         0
621.B  71     PEO  TRUE  -8.0         10  621      0         0
622.A  71 Egnatia FALSE -10.5          9  622      0         0
622.B  71     PEO  TRUE  -8.0          9  622      0         0
623.A  71 Egnatia FALSE -10.5          9  623      0         0
623.B  71     PEO  TRUE  -8.0         10  623      0         0
624.A  71 Egnatia FALSE -10.0          9  624      0         0
624.B  71     PEO  TRUE  -8.0          9  624      0         0
625.A  71 Egnatia FALSE -11.0          9  625      0         0
625.B  71     PEO  TRUE  -8.0         11  625      0         0
626.A  71 Egnatia FALSE -10.0          9  626      0         0
626.B  71     PEO  TRUE  -8.0         10  626      0         0
627.A  71 Egnatia  TRUE -10.5          9  627      0         0
627.B  71     PEO FALSE  -8.0         11  627      0         0
628.A  71 Egnatia  TRUE -10.0          9  628      0         0
628.B  71     PEO FALSE  -8.0         11  628      0         0

structure(list(ind = c(1L, 1L, 1L, 1L, 1L, 1L), route = c(FALSE, 
TRUE, FALSE, TRUE, TRUE, FALSE), Freq_2.4w = c(0L, 0L, 0L, 0L, 
0L, 0L), Gender = c(0L, 0L, 0L, 0L, 0L, 0L), alt = c("Egnatia", 
"PEO", "Egnatia", "PEO", "Egnatia", "PEO"), cost = c(-11, -12, 
-11, -12, -10.5, -12), total_time = c(9, 9, 9, 10, 9, 9), chid = c(1L, 
1L, 2L, 2L, 3L, 3L)), row.names = c("1.A", "1.B", "2.A", "2.B", 
"3.A", "3.B"), class = c("mlogit.data", "data.frame"))

最佳答案

因为 Gender 是一个因子变量。 R 会自动将 Gender=0 设置为基线，并使用名为 Gender1 的变量估计 Gender=1 的效果。

诚然还没有测试过，但这应该可以工作:

MIXL_socdem_IX <- gmnl(formula = route~ cost + total_time |1 |0 |Gender-1, data = dataIX_MIXL, model = "mixl", ranp=c(cost="ln", total_time="n"), mvar = list(cost=c("Gender1"), total_time=c("Gender1")), panel=T, R=200, haltons = list("primes"= c(3,5), drop= rep(10,2)), method="bfgs", print.init = T)

对于看到这篇文章的其他人，例如，如果您有四个类别的性别编码为 0、1、2、3，您将设置:

mvar = list(cost=c("Gender1","Gender2","Gender3"),total_time=("Gender1","Gender2","Gender3"))

关于r - R 包 gmnl 和 Rchoice 中可能存在错误？，我们在Stack Overflow上找到一个类似的问题： https://stackoverflow.com/questions/52720535/

文章推荐： protocol-buffers - protoc 在哪里以及如何安装它？

文章推荐： scala - Play 2.1 在启动时崩溃时如何进行错误记录

文章推荐： google-maps - 谷歌地图。指定地点的交通信息

javascript - 我需要将文本放在一个中，它位于一个 Div 中，该 Div 位于另一个 Div 中，该 Div 位于另一个 Div 中
我需要将文本放在中在一个 Div 中，在另一个 Div 中，在另一个 Div 中。所以这是它的样子: #document Change PIN
html - 两个背景图像。一个在 HTML 中，一个在 BODY 中。在 Firefox 中，主体图像未呈现
奇怪的事情发生了。我有一个基本的 html 代码。 html，头部， body 。(因为我收到了一些反对票，这里是完整的代码) 这是我的CSS: html { backgroun
ios - 将图像从 asset.xcassets 加载到 imageArray 中，并将其动态加载到 UIImageView 中，该 UIImageView 存在于 UICollectionView 中 - swift
我正在尝试将 Assets 中的一组图像加载到 UICollectionview 中存在的 ImageView 中，但每当我运行应用程序时它都会显示错误。而且也没有显示图像。我在ViewDidLoa
linux - 在 BASH 中，我需要根据 perl 脚本的输出更改一些环境变量。在 tcsh 中，我可以使用别名 eval 组合。不能在 bash 中
我需要根据带参数的 perl 脚本的输出更改一些环境变量。在 tcsh 中，我可以使用别名命令来评估 perl 脚本的输出。 tcsh: alias setsdk 'eval `/localhome/
asp.net - Windows 身份验证适用于 IIS，但不适用于 Kestrel/Microsoft.AspNetCore.Authentication.Negotiate(不在 Chrome 中，有时在 Edge 中，始终在 IE 中)？
我使用 Windows 身份验证创建了一个新的 Blazor(服务器端)应用程序，并使用 IIS Express 运行它。它将显示一条消息“Hello Domain\User!”来自右上方的以下 Ra
java - java 中 Kotlin 中的等价物是什么？
这是我的方法 void login(Event event);我想知道 Kotlin 中应该如何最佳答案在 Kotlin 中通配符运算符是 * 。它指示编译器它是未知的，但一旦知道，就不会有其他类
express - 在 Jade 中，为什么有时我可以按原样使用变量而有时必须将它们包含在#{......} 中？
看下面的代码 for story in book if story.title.length < 140 - var story
c - C 中 strstr() 中 for 循环的错误使用
我正在尝试用 C 语言学习字符串处理。我写了一个程序，它存储了一些音乐轨道，并帮助用户检查他/她想到的歌曲是否存在于存储的轨道中。这是通过要求用户输入一串字符来完成的。然后程序使用 strstr()
c - * 在 sscanf 中，* 在 [] 中
我正在学习 sscanf 并遇到如下格式字符串: sscanf("%[^:]:%[^*=]%*[*=]%n",a,b,&c); 我理解 %[^:] 部分意味着扫描直到遇到 ':' 并将其分配给 a。:
python - 在 Python (2.7.3) 中，如果 str(x) 中的任何字符在 str(y) 中(或 str(y) 在 str(x) 中)，我如何编写一个函数来回答？
def char_check(x,y): if (str(x) in y or x.find(y) > -1) or (str(y) in x or y.find(x) > -1):
ansible - 在 Ansible 中，如何将一行移动到一个 block 中？
我有一种情况，我想将文本文件中的现有行包含到一个新 block 中。 line 1 line 2 line in block line 3 line 4 应该变成 line 1 line 2 line
Django 调试工具栏显示在根 URL 中，但不显示在应用程序 URL 中
我有一个新项目，我正在尝试设置 Django 调试工具栏。首先，我尝试了快速设置，它只涉及将 'debug_toolbar' 添加到我的已安装应用程序列表中。有了这个，当我转到我的根 URL 时，调试
r - 在 R 中，Matlab 中 @ 函数句柄的等价物是什么？
在 Matlab 中，如果我有一个函数 f，例如签名是 f(a,b,c)，我可以创建一个只有一个变量 b 的函数，它将使用固定的 a=a1 和 c=c1 调用 f: g = @(b) f(a1, b,
swiftui - SwiftUI 中 ScrollView 中 VStack 元素中的神秘间距或填充
我不明白为什么 ForEach 中的元素之间有多余的垂直间距在 VStack 里面在 ScrollView 里面使用 GeometryReader 时渲染自定义水平分隔线。 Scrol
cookies - 什么应该存储在 session 中，什么应该存储在 cookie 中？
我想知道，是否有关于何时使用 session 和 cookie 的指南或最佳实践？什么应该和什么不应该存储在其中？谢谢! 最佳答案这些文档很好地了解了 session cookie 的安全问题以及
python - Python 中 matplotlib 中 3d 直方图的奇怪行为
我在 scipy/numpy 中有一个 Nx3 矩阵，我想用它制作一个 3 维条形图，其中 X 轴和 Y 轴由矩阵的第一列和第二列的值、高度确定每个条形的是矩阵中的第三列，条形的数量由 N 确定。
c - c 中 sem_init(...) 中 value 参数的不同用法
假设我用两种不同的方式初始化信号量 sem_init(&randomsem,0,1) sem_init(&randomsem,0,0) 现在， sem_wait(&randomsem) 在这两种情况下
c - 实际值存储在 pstr 中，但是该值如何存储在数组 "WORD"中
我怀疑该值如何存储在“WORD”中，因为 PStr 包含实际输出。？既然Pstr中存储的是小写到大写的字母，那么在printf中如何将其给出为“WORD”。有人可以吗？解释一下？ #include
javascript - 数组索引选择像在 numpy 中，但在 javascript 中
我有一个 3x3 数组: var my_array = [[0,1,2], [3,4,5], [6,7,8]]; 并想获得它的第一个 2
javascript - 在 Javascript 中，如何检测浏览器窗口何时在 View 中？
我意识到您可以使用如下方式轻松检查焦点: var hasFocus = true; $(window).blur(function(){ hasFocus = false; }); $(win

行者123

个人简介

我是一名优秀的程序员,十分优秀！

作者热门文章

滴滴打车优惠券免费领取

全站热门文章

首页

博学

6Ren·AI

商城

r - R 包 gmnl 和 Rchoice 中可能存在错误？