Working with EMMA

It has been a while, but I need to make sure that EMMA works for Dr. Kleeberger’s project

Understand the population genetics

A hands on tutorial here.

Understand the Mixed model

plink provides many good tools.

Implementation with R — I encountered some hurdles with R programming and would like to document them here before I forget.

If fctr.cols are the names of your factor columns, to convert
them to characters

X[, fctr.cols] <- sapply(X[, fctr.cols], as.character)

 

distances <- matrix(1:25, nrow=5, ncol=5)
Now, 
apply(distances, c(1, 2), function(x) 0)
distances[] <- 0L
distances*0

Viewing/reporting the results

I found a very useful and handy help doc on Manhattan Plot