我在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)水。
如果有人可以指导我如何正确加载表格以供查看,以及如何上传数据文件以配合此问题,以便有人可以运行该文件,我将非常感激。
我知道任何人都无法正确回答问题,直到他们能够获取数据并运行代码,但表中可能会有一些内容可能会提供一些解释,这可能会导致我成功使用这些数据。