library(ctl)
data(multitrait) # Arabidopsis Thaliana (R/qtl cross object)
multitrait$pheno <- multitrait$pheno[,1:4]
slopes <- scanSlopes.cross(multitrait)
image(1:nrow(slopes), 1:ncol(slopes), -log10(slopes))
Run the code above in your browser using DataLab