извлечение стандартизированных бета-коэффициентов из GAM со случайными эффектами в R - PullRequest
Новая
       118

извлечение стандартизированных бета-коэффициентов из GAM со случайными эффектами в R

1 голос
/ 22 марта 2020

Я ищу способ расчета стандартизированных бета-коэффициентов по моей модели гаммы, построенной с использованием mgcv. Вопрос в этой ссылке полезен, но только для лм. извлечение стандартизированных коэффициентов из lm в R Я не знаю, как будут обрабатываться случайные эффекты (вписываются в сплайны с "re"). Я бы предпочел простую функцию, такую ​​как подход lm.beta. 

library(codyn)
library(tidyr) ## for spread()
names(pplots)

pplots.spread <- as.data.frame(spread(data = pplots, 
                        key = species, 
                        value = relative_cover))
names(pplots.spread)<-make.names(names(pplots.spread),unique = TRUE)

mod1 <- mgcv::gam(sorghastrum.nutans~
                    s(solidago.missouriensis,m=3,k=5,bs="bs",fx=T) +
                    s(year,bs="cr",k=3) + 
                    s(treatment,bs="re"),
                  family=quasibinomial(link="logit"),
                  data=pplots.spread, select=T, method="REML") 
summary(mod1)

Можно ли извлечь стандартизированные бета-коэффициенты без переоснащения модели?

...