Displaying 1 result from an estimated 1 matches for "sngl".
Did you mean:
sgl
2007 Jun 10
1
{nlme} Multilevel estimation heteroscedasticity
...normexam ~ standLRT,
data = Exam,
weights = varIdent(form = ~ 1 | type),
random = ~ 1 | school)
This gives me nice and clear output, part of which is shown below:
Variance function:
Structure: Different standard deviations per stratum
Formula: ~normexam | type
Parameter estimates:
Mxd Sngl
1.000000 1.034607
Number of Observations: 4059
Number of Groups: 65
Though, should I hypothesize that the variance on the normexam-
variable is larger on schools that have a higher average score on
intake-exams (schavg), I run into troubles. I'd use weights = varIdent
(form = ~ 1 | schavg...