if (FALSE) {
data(mppData_GE)
Qpos <- c("PZE.105068880", "PZE.106098900")
EC <- matrix(c(180, 310, 240, 280), 4, 1)
rownames(EC) <- c('CIAM', 'TUM', 'INRA', 'KWS')
colnames(EC) <- 'cum_rain'
Qeff <- QTL_effect_QxEC(mppData = mppData_GE,
trait = c('DMY_CIAM', 'DMY_TUM', 'DMY_INRA_P', 'DMY_KWS'),
env_id = c('CIAM', 'TUM', 'INRA', 'KWS'),
QTL = Qpos, EC = EC)
pl <- plot_QTLxEC(Qeff, Q_id = 2, RP = "UH007", EC_id = 'cum rain',
trait_id = 'DMY')
}
Run the code above in your browser using DataLab