## Not run: ------------------------------------
# library(synbreedData)
# data(maize)
#
# # realized kinship matrix
# maizeC <- codeGeno(maize)
# U <- kin(maizeC,ret="realized")/2
#
# # solution with gpMod
# m <- gpMod(maizeC,kin=U,model="BLUP")
#
# # solution with MME
# diag(U) <- diag(U) + 0.000001 # to avoid singularities
# # determine shrinkage parameter
# lambda <- m$fit$sigma[2]/ m$fit$sigma[1]
# # multiply G with shrinkage parameter
# GI <- solve(U)*lambda
# y <- maizeC$pheno[,1,]
# n <- length(y)
# X <- matrix(1,ncol=1,nrow=n)
# mme <- MME(y=y,Z=diag(n),GI=GI,X=X,RI=diag(n))
#
# # comparison
# head(m$fit$predicted[,1]-m$fit$beta)
# head(mme$u)
## ---------------------------------------------
Run the code above in your browser using DataCamp Workspace