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)
m1
comparecc(m1, type="vcc")
comparecc(m1, type="ecc")
Run the code above in your browser using DataLab