我很难进行累积链接混合模型。
有一个名为clm function package ordinal “Hessian is numerically singular: parameters are not uniquely determined的SO问题,但尚未得到解答。
我的因变量是序数因子(1到4,不同级别具有特定含义),我的两个自变量(活动变量和速度变量)也是。我还有两个自变量,分别是时间和治疗因素,以及一个随机的数字变量(但也可以是一个因子)。
下面是我的data.frame的摘要,看起来像Christensen在CRAN ordinal package pdf中概述的data.frames。
'data.frame': 158 obs. of 6 variables:
$ fish_id: num 1 2 3 5 7 8 13 14 15 16 ...
$ react : Ord.factor w/ 4 levels "1"<"2"<"3"<"4": 3 3 3 3 3 2 4 3 2 3 ...
$ treat : Factor w/ 3 levels "C","S","T": 2 2 2 2 2 2 2 2 2 2 ...
$ time : Factor w/ 2 levels "1","2": 1 1 1 1 1 1 1 1 1 1 ...
$ active : Ord.factor w/ 2 levels "1"<"2": 2 1 1 1 2 1 1 2 1 1 ...
$ speed : Ord.factor w/ 3 levels "1"<"2"<"3": 2 2 2 3 3 2 2 3 2 2 ...
Ordinal factor values may be used as independent variables,并且我已经检查了每个级别的等级,以确保它们合适:
> levels(reaction_123$treat)
[1] "C" "S" "T"
> levels(reaction_123$react)
[1] "1" "2" "3" "4"
> levels(reaction_123$active)
[1] "1" "2"
> levels(reaction_123$speed)
[1] "1" "2" "3"
运行模型时:
r<-clm(react ~ treat + time + speed + active + (1|fish_id),
data=reaction_123)
我收到警告:
(1)粗麻布在数字上是单数的:参数不是唯一确定的。另外:满足绝对收敛准则,但不满足相对准则
尽管这是一个警告,但我的摘要(r)指出
系数(由于奇异而未定义1)
我的变量都不包含单个值(所有变量都至少包含两个),所以如果有人可以帮助我确定为什么会发生这种情况,我将不胜感激。