CLM(累积链接混合模型)中的数值奇异性问题

时间:2018-06-20 18:38:25

标签: mixed-models

我很难进行累积链接混合模型。

有一个名为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)

我的变量都不包含单个值(所有变量都至少包含两个),所以如果有人可以帮助我确定为什么会发生这种情况,我将不胜感激。

0 个答案:

没有答案