Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- library(mgcv)
- library(arm)
- example<-gamm(mag ~ s(depth), random=list(stations=~1), data=quakes)
- summary(example$gam)
- #Family: gaussian
- #Link function: identity
- #Formula:
- # mag ~ s(depth)
- #Parametric coefficients:
- # Estimate Std. Error t value Pr(>|t|)
- #(Intercept) 5.02300 0.04608 109 <2e-16 ***
- # ---
- # Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
- #Approximate significance of smooth terms:
- # edf Ref.df F p-value
- #s(depth) 3.691 3.691 43.12 <2e-16 ***
- # ---
- # Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
- #
- #R-sq.(adj) = 0.0725
- #Scale est. = 0.036163 n = 1000
- ranef(example$lme)$stations[[1]] # extract random intercepts
- #se.ranef(example$lme) # extract se of random intercepts - Problem line - doesn't work?
Add Comment
Please, Sign In to add comment