Question

我正在使用带有R的<{1}}包。

使用以下方法导入数据后

mlogit

并致电：

t <-read.csv('junk.csv',header=TRUE, sep=",", dec=".")

我收到以下错误：

x <- mlogit.data(t,choice="D",shape="long",id.var="key",alt.var="altkey")

任何想法如何解决？

我的数据在csv文件中以以下格式存在：

[junk.csv]

Error in `row.names<-.data.frame`(`*tmp*`, value = c("1.1", "1.2", "1.3",  : 
  duplicate 'row.names' are not allowed
In addition: Warning message:
non-unique values when setting 'row.names': ‘1.1’, ‘1.2’, ‘1.3’, ‘1.4’, ‘1.5’, ‘1.6’

Answer 1

我对mlogit的体验是，对于那些不完全的数据，它并不是很宽容。

在您的情况下，我注意到第一个受访者有6个选择，而第二个受访者有7个选择。如果您将数据格式化为每个响应者具有相同数量的备选方案，mlogit.data函数将起作用：

dat <- read.table(sep=",",text="
key,altkey,A,B,C,D
201005131,1, 2.6,118.17,117,0
201005131,2,1.4,117.11,115,0
201005131,3,1.1,117.38,122,1
201005131,4,24.6,,122,0
201005131,5,48.6,91.90,122,0
201005131,6,59.8,,122,0
201005132,1,20.2,118.23,113,0
201005132,2,2.5,123.67,120,1
201005132,3,7.4,116.30,120,0
201005132,4,2.8,118.86,120,0
201005132,5,6.9,124.72,120,0
201005132,6,2.5,123.81,120,0
201005132,7,8.5,119.23,115,0
", header=TRUE)

对所有数据运行mlogit会重现错误：

> mlogit.data(dat, choice="D", shape="long", id.var="key", alt.var="altkey")
Error in `row.names<-.data.frame`(`*tmp*`, value = c("1.1", "1.2", "1.3",  : 
  duplicate 'row.names' are not allowed
In addition: Warning message:
non-unique values when setting 'row.names': '1.1', '1.2', '1.3', '1.4', '1.5', '1.6'

但是，删除第13行（即第7个备选方案）可以起作用：

> mlogit.data(dat[-13, ], choice="D", shape="long", id.var="key", alt.var="altkey")
          key altkey    A      B   C     D
1.1 201005131      1  2.6 118.17 117 FALSE
1.2 201005131      2  1.4 117.11 115 FALSE
1.3 201005131      3  1.1 117.38 122  TRUE
1.4 201005131      4 24.6     NA 122 FALSE
1.5 201005131      5 48.6  91.90 122 FALSE
1.6 201005131      6 59.8     NA 122 FALSE
2.1 201005132      1 20.2 118.23 113 FALSE
2.2 201005132      2  2.5 123.67 120  TRUE
2.3 201005132      3  7.4 116.30 120 FALSE
2.4 201005132      4  2.8 118.86 120 FALSE
2.5 201005132      5  6.9 124.72 120 FALSE
2.6 201005132      6  2.5 123.81 120 FALSE

当然，这不是很令人满意，因为它会破坏一些数据。更好的解决方案是以mlogit()期望的格式构建数据，然后直接调用mlogit()：

dat$key <- factor(as.numeric(as.factor(dat$key)))
dat$altkey <- as.factor(dat$altkey)
dat$D <- as.logical(dat$D)
row.names(dat) <- paste(dat$key, dat$altkey, sep = ".")

现在数据如下所示：

    key altkey    A      B   C     D
1.1   1      1  2.6 118.17 117 FALSE
1.2   1      2  1.4 117.11 115 FALSE
1.3   1      3  1.1 117.38 122  TRUE
1.4   1      4 24.6     NA 122 FALSE
1.5   1      5 48.6  91.90 122 FALSE
1.6   1      6 59.8     NA 122 FALSE
2.1   2      1 20.2 118.23 113 FALSE
2.2   2      2  2.5 123.67 120  TRUE
2.3   2      3  7.4 116.30 120 FALSE
2.4   2      4  2.8 118.86 120 FALSE
2.5   2      5  6.9 124.72 120 FALSE
2.6   2      6  2.5 123.81 120 FALSE
2.7   2      7  8.5 119.23 115 FALSE

您可以直接致电mlogit()：

mlogit(D ~ A + B + C, dat, 
       chid.var = "key", 
       alt.var = "altkey", 
       choice = "D", 
       shape = "long")

结果：

Call:
mlogit(formula = D ~ A + B + C, data = dat, chid.var = "key",     alt.var = "altkey", choice = "D", shape = "long", method = "nr",     print.level = 0)

Coefficients:
2:(intercept)  3:(intercept)  4:(intercept)  5:(intercept)  6:(intercept)  
      10.7774         4.8129         5.2257       -17.2522        -7.7364  
7:(intercept)              A              B              C  
      10.0389         1.6010         2.7156         2.9888

Answer 2

JohnP，

我认为你要找的答案是maxabet。它说：

“如果读者尝试使用R来计算第6章中介绍的MLR，可能会出现以下错误消息：

`row.names＆lt; - .data.frame ...中的错误（之后的一些数据）不允许重复'row.names'

此错误是当前版本的R（2011年4月发布）的问题。只有在每场比赛中“行”（在我们的例子中为“跑步者”）的数量相同时，MLR才会运行。因此，其中一种方法是制造一些“假马”，这样每个种族都有14个跑步者。读者可能要么等待R来解决这个问题，要么尝试某种方式来下载早期版本（我使用的是2.9.1），而不是制作虚拟跑步者。 “

所以我猜你可以使用较旧版本的R或更好的版本，用虚拟值填充现有数据以完成每场比赛。我要做的一件事是为每个种族添加一个新的“字段大小”变量。通过投掷假人，但在新列中具有正确数量的跑步者应该照顾任何场大小差异。

显然，你得到的答案将与你正在阅读的书有所不同，但这将是一个良好的开端，以便更好地设置障碍：）

我应该如何格式化R mlogit包的数据？

2 个答案: