Coëfficiënten extraheren
We willen vaak weten welke waarden een model als coëfficiënten schat. Hoewel de functie summary() de modelresultaten laat zien, willen we soms ook direct toegang tot deze outputs.
De schattingen voor fixed effects kun je direct opvragen met de functie fixef(). De schattingen voor random effects kun je direct opvragen met de functie ranef(). We kunnen ook betrouwbaarheidsintervallen voor de fixed effects extraheren met de functie confint().
Het pakket broom.mixed bevat ook tidy-methoden om modelresultaten uit lmer()-modellen te halen, namelijk de functie tidy(). Deze resultaten zijn echter complexer en minder tidy dan veel andere tidy-outputs, vanwege de complexiteit van mixed-effectmodellen.
Er is al een lmer-model voor je gefit en opgeslagen als out.
Deze oefening maakt deel uit van de cursus
Hiërarchische en Mixed-Effects-modellen in R
Oefeninstructies
- Extraheer de coëfficiënten van de fixed effects met
fixef()op het opgeslagen modelout. - Extraheer de coëfficiënten van de random effects met
ranef()op het opgeslagen modelout. - Schat de 95%-betrouwbaarheidsintervallen met de functie
confint()op het opgeslagen modelout. - Gebruik
tidy()metoutenconf.int = TRUEom je vorige drie aanroepen te herhalen met één tidy-opdracht.
Praktische interactieve oefening
Probeer deze oefening eens door deze voorbeeldcode in te vullen.
# Extract the fixed-effect coefficients
___
# Extract the random-effect coefficients
___
# Estimate the confidence intervals
___
# Use the tidy function
___