library(EstimateBreed)
data("maize")
#Extract heterosis and heterobeltiosis
general <- with(maize,het(GEN,GM,GP,PR,REP,param="all"))
#Only extract heterosis
het <- with(maize,het(GEN,GM,GP,PR,REP,param = "het"))
#Extract only heterobeltiosis
hetb <- with(maize,het(GEN,GM,GP,PR,REP,param = "hetb"))
Run the code above in your browser using DataLab