inst/example.r

#example
data(jtrain1)
(fm.fem <- fem(log(scrap)~d88+d89+grant+grant.1,id=fcode,data=jtrain1)) #from wood p430
(fm.bem <- bem(log(scrap)~d88+d89+grant+grant.1,id=fcode,data=jtrain1))

library(nlme)
summary(gls(log(scrap)~d88+d89+grant+grant.1,cor=corCompSymm(form=~1|fcode),data=jtrain1,na.action="na.exclude"))

Try the bpm package in your browser

Any scripts or data that you put into this service are public.

bpm documentation built on May 2, 2019, 5:46 p.m.