|
18 | 18 | #' @returns |
19 | 19 | #' A dataframe with all possible variant-environmental pairs and their estimated interaction effect |
20 | 20 | #' @examples |
21 | | -#' g_vec <- matrix(0, nrow = 100000, ncol = 3) |
| 21 | +#' N_run <- 50000 |
| 22 | +#' g_vec <- matrix(0, nrow = N_run, ncol = 3) |
22 | 23 | #' freqs <- runif(ncol(g_vec), min = 0, max = 1) |
23 | | -#' env_vec <- matrix(0, nrow = 100000, ncol = 3) |
| 24 | +#' env_vec <- matrix(0, nrow = N_run, ncol = 3) |
24 | 25 | #' for(i in 1:ncol(g_vec)){ |
25 | | -#' g_vec[, i] <- rbinom(100000, 2, freqs[i]) |
| 26 | +#' g_vec[, i] <- rbinom(N_run, 2, freqs[i]) |
26 | 27 | #' } |
27 | 28 | #' for( i in 1:ncol(env_vec)){ |
28 | | -#' env_vec[, i] <- round(runif(100000,min=0,max=6)) |
| 29 | +#' env_vec[, i] <- round(runif(N_run,min=0,max=6)) |
29 | 30 | #' } |
30 | | -#' cc_vec <- rbinom(100000,1,0.1 * (1.05 ^ g_vec[, 1]) * |
| 31 | +#' cc_vec <- rbinom(N_run,1,0.1 * (1.05 ^ g_vec[, 1]) * |
31 | 32 | #' (1.06 ^ env_vec[,1]) * (0.95 ^ g_vec[, 2]) * |
32 | 33 | #' (1.1^(g_vec[, 1] * env_vec[, 1]))) |
33 | 34 | #' res <- pairwise_env_int_CC.calc(cc_vec, g_vec, env_vec) |
|
0 commit comments