Displaying 1 result from an estimated 1 matches for "fpi1".
Did you mean:
fp1
2008 Apr 06
1
lme cant get parameter estimated correctly
...t.seed(76589437887)
fph = 0.4
Sigh = sqrt(0.0002)
Sigi = sqrt(0.04)
ci = 1
fpi = matrix(,7200,3)
for (i in 1:90) {
fph = rnorm(1, fph, Sigh)
for (k in 1:80) {
fpi[ci,1:3] = matrix(c(i, k, rnorm(1, fph, Sigi)),1)
ci = ci+1
}
}
colnames(fpi) = c("hospid", "empid", "fpi1")
dta = as.data.frame(fpi)
lme = lme(fpi1 ~ 1, dta, ~1|hospid)
summary(lme)
lme = lme(fpi1 ~ 1, dta, ~1|hospid)
summary(lme)
Linear mixed-effects model fit by REML
Data: dta
AIC BIC logLik
-2555.416 -2534.771 1280.708
Random effects:
Formula: ~1 | hospid...