ich arbeite seit kurzem mit gam, gamm und gamm4 und bitte euch um Unterstützung beim Verständnis dieser Funktionen:
Ich habe ein multilineares Modell mit einer GAM komponente an einen Datensatz anpassen lassen:
Code: Alles auswählen
rawmod = gamm4(data = dfscaled, formula = bm ~ mnh_d_ar + sdi + bon_100 + maxh + api + s(DG, bs='cr', k = -1), random = ~(1|Tnr), REML=REML, weights = weig_lm_r)
Ich habe mir summary(modl) ausgeben lassen
Code: Alles auswählen
Random effects:
Groups Name Variance Std.Dev.
Tnr (Intercept) 10376213.2738 3221.2130
Xr s(DG) 10092728.7380 3176.9055
Residual 0.7849 0.8859
Number of obs: 7784, groups: Tnr, 4858; Xr, 8
Dafür habe ich mir die Werte des Effektes erzeugt wie folgt:
Code: Alles auswählen
p_l_ti0 = predict(modl, re.form = ~0) # no random effects
p_l_tnr = predict(modl, re.form = ~(1|Tnr)) # prediction on level of group Tnr
Code: Alles auswählen
sd(p_l_tnr - p_l_ti0)
Vielen Dank
Werner