| Demand {SASmixed} | R Documentation |
The Demand data frame has 77 rows and 8 columns of data on
per-capita demand deposits by state and year.
This data frame contains the following columns:
WA < FL < CA < TX < IL < DC < NY
1949 < ...{} < 1959
lme.
library(SASmixed)
data(Demand)
names( Demand )
formula( Demand )
# the "grp" factor is a dummy factor with only one level.
unique( Demand$grp )
# Crossed random-effects factors have to be created by pdIdent
# applied to the indicator variables and joined by pdBlocked.
fm1Demand <-
lme( log(d) ~ log(y) + log(rd) + log(rt) + log(rs), data = Demand,
random = list(grp = pdBlocked(list(pdIdent(~ State - 1),
pdIdent(~ Year - 1)))))
summary( fm1Demand ) # compare to output 3.13, p. 132
VarCorr( fm1Demand )