gm = ~al:an:st
vcc = list(~me+st, ~ve+an, ~al)
ecc = list(~me:ve+me:al, ~ve:al+al:st)
data(math)
m1 <- rcox(gm=gm, vcc=vcc, ecc=ecc, data=math, fit=FALSE)
f1 <- matching(m1)
## Use f1$K as starting value
scoring(m1, K0=f1$K)
ipm(m1, K0=f1$K)
Run the code above in your browser using DataLab