当使用分类独立变量时,mlogit导致“系统是计算上奇异的”错误

时间:2015-04-11 01:23:39

标签: r

我在R Studio中使用R x64 3.1.1并且一直在尝试使用mlogit开发多项logit模型时遇到困难,而mlogit包含相当大的数据集,包含分类和数字IV。该模型使用来自卫星图像(自变量)的土地覆盖数据来预测土地利用(因变量)。这些数据适用于北卡罗来纳州的地块。

当我为完整型号运行mlogit时,我收到:  " solve.default中的错误(H,g [!fixed]):   系统在计算上是单数的:倒数条件数= 1.6507e-131"

我已经仔细研究了这个论坛中的所有问题并阅读了奇点的概念,但是当IV是绝对的时候还没有找到解决这种情况的问题,我也不知道如何解决这个问题。

通过迭代过程,我相信我已经缩小了可能导致奇点问题的因变量,但1)我不确定,2)不确定我知道如何处理它。它是实际土地覆盖的分类变量,我不能将其排除在模型之外。

以下是仅包含DV(lu.srs.c4.age)和IV(lc1)的基本代码:

set.seed(12976)
library(mlogit)

ncdat<-dget("ncdat")
str(ncdat)
table(ncdat$lu_srs_c4_age, ncdat$lc1)

#REGRESSION OF lu_srs_c4_age VS lc1
lu4.log <- mlogit(lu_srs_c4_age ~ 1 | lc1,data=ncdat, shape="wide", 
                  reflevel="developed_no_age_given")
summary(lu4.log)

这是表格:

agri/pasture barren developed emerge wetland forest grass/herb shrub water
  ag_range_no_age_given              2612      2       197             12    277        100    74     8
  cover_unk_no_age_given                6      1         2              5    165          1     2     0
  developed_no_age_given              445     21      1136              7    410        134    76     7
  oth_for_closed                        0      0         0              5     29          0     0     0
  oth_for_neg_age                       0      0         0              0      6          0     0     0
  oth_for_new                           0      0         0              0      2          0     0     0
  oth_for_sapling                       0      0         0              1      4          0     0     0
  oth_for_seedling                      0      0         0              3      4          0     0     0
  timber_closed                       108      3       186             50   6598        124   262    10
  timber_neg_age                        3      0         1              1     28          6     6     0
  timber_new                            1      0         1              2     27          7     3     0
  timber_sapling                       55      5        16             15    612         45   158     3
  timber_seedling                      48      4        17             34    265        125   144     0
  water_no_age_given                    4     22        11             21     45          6     5  1431
  wetland_beach_no_age_given            1      5         2            108     18          3     0    14

我很抱歉没有正确排队 - 除了代码片段工具之外,我没有看到任何进入表格的方法。 lc1的土地覆盖类别是:1)农业/牧场,2)贫瘠,3)开发,4)出现湿地,5)森林,6)草/草,7)灌木,和8)水。

如果有人可以指导我如何正确加载表格以供查看,以及如何上传数据文件以配合此问题,以便有人可以运行该文件,我将非常感激。

我知道任何人都无法正确回答问题,直到他们能够获取数据并运行代码,但表中可能会有一些内容可能会提供一些解释,这可能会导致我成功使用这些数据。

0 个答案:

没有答案