![]() ![]() One of the few exceptions in the pile of useless emails I have archived. Thanks to Etienne Laliberté, who sent me an email in early 2010 pointing out that I had to use doMC. Gratuitous picture: Trees at 8 pm illuminated by bonfire and full moon (Photo: Luis). # Run the bivariate example in parallel for two pairs of sitesįoreach(trial1=rep("FR38_1",2), trial2=c("FR38_2", "FR38_3")) %dopar% # Returns only log-Likelihood for this example # heterogeneous variances and diagonal matrix for heterogeneousĭen.bi <- asreml(cbind(den1,den2) ~ trait, # Fits overall mean, random correlated additive effects with T2$den2 <- ifelse(t2$Trial_id = trial2, t2$DEN, NA) ![]() T2$den1 <- ifelse(t2$Trial_id = trial1, t2$DEN, NA) # (This one uses the same model for all sites # Wrap call to a generic bivariate analysis # Inverse of the numerator relationship matrix Ped <- read.table('pedigree.txt', header = TRUE) RegisterDoMC() # to start a "parallel backend"ĭat <- read.table('density09.txt', header = TRUE) Library(asreml) # Multivariate mixed models ![]()
0 Comments
Leave a Reply. |
AuthorWrite something about yourself. No need to be fancy, just an overview. ArchivesCategories |