# NOT RUN {
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","D"))
# make a simple fit object - OLS-b only
sheep.fit1 <- dmm(sheep.mdf, Ymat ~ 1 + Year + Sex)
# look at model plus fixed effects and components for all traits
print(sheep.fit1)
# }
# NOT RUN {
# can do the same thing without saving fit object
dmm(sheep.mdf, Ymat ~ 1 + Year + Tb + Sex)
# so this is the default print method for an object of class 'dmm'
# }
# NOT RUN {
#cleanup
rm(sheep.fit1)
rm(sheep.mdf)
rm(sheep.df)
# }
Run the code above in your browser using DataLab