library(dmm)
data(sheep.df)
sheep.mdf <- mdf(sheep.df,pedcols=c(1:3),factorcols=c(4:6),ycols=c(7:9),
sexcode=c("M","F"),relmat=c("E","A"))
# make a simple fit object - OLS-b only
sheep.fit1 <- dmm(sheep.mdf, Ymat ~ 1 + Year + Sex,
components=c("VarE(I)","VarG(Ia)"))
# compute a 'summary.dmm' object, use all the defaults
sheep.sum <- summary(sheep.fit1)
# print the summary of genetic parameters
print(sheep.sum)
if (FALSE) {
# can do the same thing without saving response object
summary(sheep.fit1)
# so this is the default print method for an object of class 'summary.dmm'
}
#cleanup
rm(sheep.fit1)
rm(sheep.sum)
rm(sheep.mdf)
rm(sheep.df)
Run the code above in your browser using DataLab