Dredge в MuMin (R) сохраняет модели с членами более высокого порядка без их соответствующих членов более низкого порядка - PullRequest
Новая
       22

Dredge в MuMin (R) сохраняет модели с членами более высокого порядка без их соответствующих членов более низкого порядка

0 голосов
/ 20 марта 2019

Я использую функцию Dredge из пакета MuMin для выбора модели.В моей модели есть многочлен.У меня сложилось впечатление (https://www.rdocumentation.org/packages/MuMIn/versions/1.42.1/topics/dredge, см. «Взаимодействия»), что MuMin как стандартные модели канав, которые включают член более высокого порядка без соответствующего члена более низкого порядка;однако, когда я использую это, верхняя модель включает только член более высокого порядка.

Вот код, который я использую для создания глобальной модели, которая включает hour: как в первом, так и во втором порядке 

sc.mod.env.hour<-clmm(seedcat~hour+I(hour^2)+MeanSpeed+RH+TKE+AirTemp+MeanUpdraft+(1|jdayfac), data=sc.data, na.action = na.fail, Hess =TRUE)

Однако, когда я запускаю Dredge: 

dro.env.hour<-dredge(sc.mod.env.hour)

Выходные данные включают модели, которые включают I(hour^2) без hour

Model selection table 
    (Intrc)   ArTmp    hour hour^2 MnSpd    MnUpd     RH      TKE df   logLik   AICc  delta weight
46        + -1.2170         -3.340 1.894          -4.219          13 -571.772 1170.5   0.00  0.154
62        + -1.3410         -3.478 1.964  0.30990 -4.269          14 -570.882 1170.9   0.37  0.128
45        +                 -2.882 2.052          -3.214          12 -573.245 1171.3   0.81  0.103
47        +         -0.3595 -2.670 1.942          -3.601          13 -572.374 1171.7   1.20  0.084
63        +         -0.4495 -2.712 2.010  0.33870 -3.645          14 -571.350 1171.8   1.30  0.080

Я также пытался использовать poly(hour,2) для определения глобальной модели, но это приводит только к включению одного термина длячас.

Я использую clmm, но попробовал более простую модель с использованием lm и получил тот же результат.

Любое руководство приветствуется, спасибо.Я не был уверен, что это должно быть на перекрестной проверке, но это проблема программирования, а не статистика, поэтому подумал, что это должно быть здесь.

РЕДАКТИРОВАТЬ: Я решил проблему с помощью подмножества: 

dro.env.hour<-dredge(sc.mod.env.hour, subset=(dc(hour, I(hour^2)))

Хотя неясно, почему это проблема в первую очередь.

...