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)
getecc(m1)
getSlot(m1,"type")
fitInfo(m1)
fitInfo(m1,"K")
Run the code above in your browser using DataLab