# Load simulated data
data(test);
# Fit multivariate ridge regression
fit = pegs(Y,X);
# Heritability
print(fit$h2);
# Genetic correlations
print(fit$GC);
# Accuracy
print(mean(diag(cor(fit$hat,simu$tbv))));
Run the code above in your browser using DataLab