Формальная цитата гладких случайных членов в смешанной модели mgcv :: gam - PullRequest
Новая
       24

Формальная цитата гладких случайных членов в смешанной модели mgcv :: gam

0 голосов
/ 14 декабря 2018

У меня есть смешанная модель mgcv :: gam вида: 

m1 <- gam(Y ~ A + s(B, bs = "re"), data = dataframe, family = gaussian,
 method = "REML")

Случайный член s(B, bs = "re") указан в summary(m1) как, например, 

Approximate significance of smooth terms:
           # edf Ref.df      F p-value   
s(B)  4.486      5 97.195 6.7e-08 ***

Мой вопрос: как бы я процитировал этот результат (статистику и значение P) в официальном документе, например, в техническом отчете или в документе?

Например, одна возможность -

F [4.486,5] = 97,195, P = 6,7e-08

Однако, опровергая эту идею, «обратное проектирование» результата с использованием 

pf(q= 97.195, df1= 4.486, df2= 5, lower.tail=FALSE)

даетневерное значение p: 

[1] 5.931567e-05

Буду очень признателен за ваш совет.Большое спасибо за вашу помощь!

1 Ответ

0 голосов
/ 14 декабря 2018

Статистика F, о которой идет речь, на самом деле не следует за F с указанными вами степенями свободы.Значение Ref df относится к тесту, но вам нужно прочитать и понять, как работает Wood (2013), чтобы полностью понять, как получены степени свободы для теста.

Я бы просто процитировал статистикуи значение p, а затем процитируйте статью Саймона, если кто-то хочет знать, как они были вычислены.Я не думаю, что вы легко можете достичь тех степеней свободы, которые действительно используются.(ну, не без отладки кода summary.gam() и просмотра того, как они вычисляются.)

Ссылки

Wood, SN 2013. Простой тест на случайные эффекты в регрессионных моделях.Биометрика 100: 1005–1010.DOI: 10,1093 / Biomet / ast038

...