Blog posts

2025

BGLR with Reps

less than 1 minute read

Published:

To fit a model that has reps or any unbalanced design in BGLR, you need to first design the Z matrix, you can use model.matrix function for that. Then, in the ETA object you need to provide a ZGZ’ or ZM instead of standard X or K matrices. More on it here. Example below:

2020

Common cases that the GRM doesn’t invert

6 minute read

Published:

Here, I show common genetic reasons where can lead to genomic relationship matrix (GRM) with problems in the inversion. They are commonly based on population structure or repeated (or highly similar) information. The cases are not always true, but they can represent why you have a strong linear dependence in your G matrix (Van Raden 2008) and a not unique inverse for the matrix. After identifying the reason, you can take action in order to manage the data before building the G matrix.

2018

R Parallel with Genomic Prediction Application

8 minute read

Published:

Material developed for a workshop on R Parallel with Genomic Selection Application for the Genomic Selection Studies group at University of Florida, Gainesville, FL.

2017

2016

Introduction to R

26 minute read

Published:

For beginners with focus on Bioinformatics and Genetics

Better Graphics on R

7 minute read

Published:

The idea of this course is to introduce several cool graphics over different packages and how to build it step-by-step. Hopefully you will get some tools to run away from the basic R graphics.