In asreml-r version 2, it was possible to do multiple predict statements at the same time by making 'classify' to be a list:
Code:
p2 <- predict(m2, classify=list("hyb", "hyb:enc"),
present=list(hyb=c("hyb","enc"),
"hyb:enc"=c("loc","enc")),
ignore=list(hyb=c("hyb:year","hyb:year:enc"),
"hyb:enc"=c("hyb:year","hyb:year:enc")))
The predictions were then extracted using the name from 'classify':
Code:
p2$predictions$hyb
p2$predictions$"hyb:enc"
Has this ability been deprecated for version 3? Only one predict is allowed at a time? If so, it seems like a good move. The original functionality made things harder more often than it made things easier.
The example above was a proof-of-concept for me to figure out how to do multiple predictions and not something I ever actually used in practice. The big change is that I often wrote things like
Code:
p2$predictions$"hyb:enc"$pvals
which now will simplify to become
Code:
p2$predictions$pvals
Any other major changes to watch out for?
Kevin Wright
-------------------- m2f --------------------
Sent using Mail2Forum (http://www.mail2forum.com).
Read this topic online here:
http://www.vsni.co.uk/forum/viewtopic.php?p=645#645
-------------------- m2f --------------------
Received on Fri Jul 09 2009 - 20:11:30 EST
This webpage is part of the ASReml-l discussion list archives 2004-2010. More information on ASReml can be found at the VSN website. This discussion list is now deprecated - please use the VSN forum for discussion on ASReml. (These online archives were generated using the hypermail package.)