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_main_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_QxEC(Qeff, EC = EC, env_id = c('CIAM', 'TUM', 'INRA', 'KWS'),
QTL = 2, EC_id = 'cum rain', trait_id = 'DMY')
}
Run the code above in your browser using DataLab