Ich muss ein LMM mit einem Interaktions-Random-Effekt, aber ohne den Random-Random-Effekt, unter Verwendung des lme
-Befehls anpassen. Das heißt, ich möchte das Modell in oats.lmer
(siehe unten) passen, aber die Funktion lme
aus dem nlme
Paket verwenden.Verschachtelte zufällige Effekte in `lme {nlme}`
Der Code ist
require("nlme")
require("lme4")
oats.lmer <- lmer(yield~nitro + (1|Block:Variety), data = Oats)
summary(oats.lmer)
#Linear mixed model fit by REML ['lmerMod']
#Formula: yield ~ nitro + (1 | Block:Variety)
# Data: Oats
#
#REML criterion at convergence: 598.1
#
#Scaled residuals:
# Min 1Q Median 3Q Max
#-1.66482 -0.72807 -0.00079 0.56416 1.85467
#
#Random effects:
# Groups Name Variance Std.Dev.
# Block:Variety (Intercept) 306.8 17.51
# Residual 165.6 12.87
#Number of obs: 72, groups: Block:Variety, 18
#
#Fixed effects:
# Estimate Std. Error t value
#(Intercept) 81.872 4.846 16.90
#nitro 73.667 6.781 10.86
#
#Correlation of Fixed Effects:
# (Intr)
#nitro -0.420
Ich begann mit diesem
oats.lme <- lme(yield~nitro, data = Oats, random = (~1|Block/Variety))
summary(oats.lme)
#Linear mixed-effects model fit by REML
# Data: Oats
# AIC BIC logLik
# 603.0418 614.2842 -296.5209
#
#Random effects:
# Formula: ~1 | Block
# (Intercept)
#StdDev: 14.50596
#
# Formula: ~1 | Variety %in% Block
# (Intercept) Residual
#StdDev: 11.00468 12.86696
#
#Fixed effects: yield ~ nitro
# Value Std.Error DF t-value p-value
#(Intercept) 81.87222 6.945273 53 11.78819 0
#nitro 73.66667 6.781483 53 10.86291 0
# Correlation:
# (Intr)
#nitro -0.293
#
#Standardized Within-Group Residuals:
# Min Q1 Med Q3 Max
#-1.74380770 -0.66475227 0.01710423 0.54298809 1.80298890
#
#Number of Observations: 72
#Number of Groups:
# Block Variety %in% Block
# 6 18
spielen, aber das Problem ist, dass es auch eine marginale Zufallseffekt für Variety
setzt die ich weglassen soll.
Die Frage ist: Wie spezifiziert man die zufälligen Effekte in oats.lme
, so dass oats.lme
ist identisch (zumindest strukturell) zu oats.lmer
?
Das ist richtig. Ich habe nicht darüber nachgedacht, aber es funktioniert. Vielen Dank. – utobi
meinst du "eine marginale zufällige Wirkung für' Block' "? –