myY <- read.table("~/temp_file/GAPIT3-master/inst/extdata/mdp_traits.txt", head =TRUE) myGD=read.table(file="~/temp_file/GAPIT3-master/inst/extdata/mdp_numeric.txt",head=T) myGM=read.table(file="~/temp_file/GAPIT3-master/inst/extdata/mdp_SNP_information.txt",head=T)
myGAPIT_MLM <- GAPIT( Y=myY[,c(1,3)],#fist column is individual ID, the third columns is days to pollination GD=myGD, GM=myGM, PCA.total=3, model=c("GLM","MLM","CMLM","FarmCPU","Blink"), Multiple_analysis=TRUE)