data(math)
gm = ~al:an:st
vcc = list(~me + st, ~ve + an, ~al)
ecc = list(~me:ve + me:al, ~ve:al + al:st)
m1 <- rcox(gm=gm, vcc=vcc, ecc=ecc, data=math, fit=FALSE)
fit(m1, method="matching")
fit(m1, method="scoring")
fit(m1, method="ipm")
Run the code above in your browser using DataLab