R/cnmutHistoryMatrix.R
In SunPathLab/Butte: BoUnds of Time Till Expansion, inferring the SCNA arrival time in the somatic evolution toward the most recent common ancestor of tumor sample(s).
Defines functions
cnmutHistory
Documented in
cnmutHistory
#' generate history matrix in relating SCNA to SSNVs
#'
#' Given a SCNA configuration Nt (total copy) and Nb (minor copy)
#' this function produces the possible history matrices in relating
#' CN timing (time period for each stage) to the burden of SSNVs
#' at distinct allele states.
#'
#' @param nt total copy number
#' @param nb copy number of the minor allele
#' @return A list of possible matrices
#' @export
cnmutHistory <- function(nt, nb) {
cnmut <- list()
if (nt == 2 & nb == 0) {
f.con1 <- matrix(c(0, 2,
1, 0), nrow = 2, byrow = TRUE)
cnmut = list(f.con1)
} else if (nt == 3 & nb == 0){
f.con1 <- matrix(c(0, 1, 3,
0, 1, 0,
1, 0, 0), nrow = 3, byrow = TRUE)
cnmut = list(f.con1)
} else if (nt == 3 & nb == 1){
f.con1 <- matrix(c(1, 3,
1, 0), nrow = 2, byrow = TRUE)
cnmut = list(f.con1)
} else if (nt == 4 & nb == 0){
f.con1 <- matrix(c(0, 1, 2, 4,
0, 0, 1, 0,
0, 1, 0, 0,
1, 0, 0, 0), nrow = 4, byrow = TRUE)
f.con2 <- matrix(c(0, 0, 2, 4,
0, 2, 1, 0,
0, 0, 0, 0,
1, 0, 0, 0), nrow = 4, byrow = TRUE) #missing 3/4
cnmut = list(f.con1, f.con2)
} else if (nt == 4 & nb == 1){
f.con1 <- matrix(c(1, 2, 4,
0, 1, 0,
1, 0, 0), nrow = 3, byrow = TRUE)
cnmut = list(f.con1)
} else if (nt == 4 & nb == 2){
f.con1 <- matrix(c(0, 2, 4,
2, 1, 0), nrow = 2, byrow = TRUE)
#f.con2 <- matrix(c(0, 4,
# 2, 0), nrow = 2, byrow = TRUE) #syn gain
#cnmut = list(f.con1, f.con2)
cnmut = list(f.con1)
} else if (nt == 5 & nb == 0){
f.con1 <- matrix(c(0, 1, 2, 3, 5,
0, 0, 0, 1, 0,
0, 0, 1, 0, 0,
0, 1, 0, 0, 0,
1, 0, 0, 0, 0), nrow = 5, byrow = TRUE)
f.con2 <- matrix(c(0, 1, 1, 3, 5,
0, 0, 2, 1, 0,
0, 0, 0, 0, 0,
0, 1, 0, 0, 0,
1, 0, 0, 0, 0), nrow = 5, byrow = TRUE) #missing 3/5
f.con3 <- matrix(c(0, 0, 1, 3, 5,
0, 1, 2, 1, 0,
0, 1, 0, 0, 0,
0, 0, 0, 0, 0,
1, 0, 0, 0, 0), nrow = 5, byrow = TRUE) #missing 4/5
f.con4 <- matrix(c(0, 0, 2, 3, 5,
0, 1, 0, 1, 0,
0, 1, 1, 0, 0,
0, 0, 0, 0, 0,
1, 0, 0, 0, 0), nrow = 5, byrow = TRUE) #missing 4/5
cnmut = list(f.con1, f.con2, f.con3, f.con4)
} else if (nt == 5 & nb == 1){
f.con1 <- matrix(c(1, 2, 3, 5,
0, 0, 1, 0,
0, 1, 0, 0,
1, 0, 0, 0), nrow = 4, byrow = TRUE)
f.con2 <- matrix(c(1, 1, 3, 5,
0, 2, 1, 0,
0, 0, 0, 0,
1, 0, 0, 0), nrow = 4, byrow = TRUE) #missing 3/5
cnmut = list(f.con1, f.con2)
} else if (nt == 5 & nb == 2){
f.con1 <- matrix(c(0, 2, 3, 5,
1, 0, 1, 0,
1, 1, 0, 0), nrow = 3, byrow = TRUE)
f.con2 <- matrix(c(0, 1, 3, 5,
1, 2, 1, 0,
1, 0, 0, 0), nrow = 3, byrow = TRUE)
cnmut = list(f.con1, f.con2)
} else if (nt == 6 & nb == 0){
f.con1 <- matrix(c(0, 1, 2, 3, 4, 6,
0, 0, 0, 0, 1, 0,
0, 0, 0, 1, 0, 0,
0, 0, 1, 0, 0, 0,
0, 1, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0), nrow = 6, byrow = TRUE)
f.con2 <- matrix(c(0, 1, 2, 2, 4, 6,
0, 0, 0, 2, 1, 0,
0, 0, 0, 0, 0, 0,
0, 0, 1, 0, 0, 0,
0, 1, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0), nrow = 6, byrow = TRUE) #missing 3/6
f.con3 <- matrix(c(0, 1, 1, 3, 4, 6,
0, 0, 1, 0, 1, 0,
0, 0, 1, 1, 0, 0,
0, 0, 0, 0, 0, 0,
0, 1, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0), nrow = 6, byrow = TRUE) #missing 4/6
f.con4 <- matrix(c(0, 1, 1, 2, 4, 6,
0, 0, 1, 2, 1, 0,
0, 0, 1, 0, 0, 0,
0, 0, 0, 0, 0, 0,
0, 1, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0), nrow = 6, byrow = TRUE) #missing 4/6
f.con5 <- matrix(c(0, 0, 1, 3, 4, 6,
0, 0, 1, 0, 1, 0,
0, 2, 1, 1, 0, 0,
0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0), nrow = 6, byrow = TRUE) #missing 4/6 and 5/6
f.con6 <- matrix(c(0, 0, 1, 2, 4, 6,
0, 0, 1, 2, 1, 0,
0, 2, 1, 0, 0, 0,
0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0), nrow = 6, byrow = TRUE) #missing 4/6 and 5/6
f.con7 <- matrix(c(0, 0, 1, 3, 4, 6,
0, 1, 1, 0, 1, 0,
0, 0, 1, 1, 0, 0,
0, 1, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0), nrow = 6, byrow = TRUE) #missing 5/6
f.con8 <- matrix(c(0, 0, 2, 3, 4, 6,
0, 1, 0, 0, 1, 0,
0, 0, 0, 1, 0, 0,
0, 1, 1, 0, 0, 0,
0, 0, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0), nrow = 6, byrow = TRUE) #missing 5/6
f.con9 <- matrix(c(0, 0, 1, 2, 4, 6,
0, 1, 1, 2, 1, 0,
0, 0, 1, 0, 0, 0,
0, 1, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0), nrow = 6, byrow = TRUE) #missing 5/6
f.con10 <- matrix(c(0, 0, 2, 2, 4, 6,
0, 1, 0, 2, 1, 0,
0, 0, 0, 0, 0, 0,
0, 1, 1, 0, 0, 0,
0, 0, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0), nrow = 6, byrow = TRUE) #missing 3/6 and 5/6
f.con11 <- matrix(c(0, 0, 0, 2, 4, 6,
0, 1, 3, 2, 1, 0,
0, 0, 0, 0, 0, 0,
0, 1, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0), nrow = 6, byrow = TRUE) #missing 3/6 and 5/6
cnmut = list(f.con1, f.con2, f.con3, f.con4, f.con5, f.con6, f.con7, f.con8, f.con9, f.con10, f.con11)
} else if (nt == 6 & nb == 1){
f.con1 <- matrix(c(1, 2, 3, 4, 6,
0, 0, 0, 1, 0,
0, 0, 1, 0, 0,
0, 1, 0, 0, 0,
1, 0, 0, 0, 0), nrow = 5, byrow = TRUE)
f.con2 <- matrix(c(1, 2, 2, 4, 6,
0, 0, 2, 1, 0,
0, 0, 0, 0, 0,
0, 1, 0, 0, 0,
1, 0, 0, 0, 0), nrow = 5, byrow = TRUE) #missing 3/6
f.con3 <- matrix(c(1, 1, 2, 4, 6,
0, 1, 2, 1, 0,
0, 1, 0, 0, 0,
0, 0, 0, 0, 0,
1, 0, 0, 0, 0), nrow = 5, byrow = TRUE) #missing 4/6
f.con4 <- matrix(c(1, 1, 3, 4, 6,
0, 1, 0, 1, 0,
0, 1, 1, 0, 0,
0, 0, 0, 0, 0,
1, 0, 0, 0, 0), nrow = 5, byrow = TRUE) #missing 4/6
cnmut = list(f.con1, f.con2, f.con3, f.con4)
} else if (nt == 6 & nb == 2){
f.con1 <- matrix(c(0, 1, 3, 4, 6,
1, 1, 0, 1, 0,
0, 1, 1, 0, 0,
1, 0, 0, 0, 0), nrow = 4, byrow = TRUE)
f.con2 <- matrix(c(0, 2, 3, 4, 6,
1, 0, 0, 1, 0,
0, 0, 1, 0, 0,
1, 1, 0, 0, 0), nrow = 4, byrow = TRUE)
f.con3 <- matrix(c(0, 1, 2, 4, 6,
1, 1, 2, 1, 0,
0, 1, 0, 0, 0,
1, 0, 0, 0, 0), nrow = 4, byrow = TRUE)
f.con4 <- matrix(c(0, 2, 2, 4, 6,
1, 0, 2, 1, 0,
0, 0, 0, 0, 0,
1, 1, 0, 0, 0), nrow = 4, byrow = TRUE) #missing 3/6
f.con5 <- matrix(c(0, 0, 2, 4, 6,
1, 3, 2, 1, 0,
0, 0, 0, 0, 0,
1, 0, 0, 0, 0), nrow = 4, byrow = TRUE) #missing 3/6
cnmut = list(f.con1, f.con2, f.con3, f.con4, f.con5)
} else if (nt == 6 & nb == 3){
f.con1 <- matrix(c(0, 1, 3, 4, 6,
0, 1, 0, 1, 0,
2, 1, 1, 0, 0), nrow = 3, byrow = TRUE)
f.con2 <- matrix(c(0, 1, 2, 4, 6,
0, 1, 2, 1, 0,
2, 1, 0, 0, 0), nrow = 3, byrow = TRUE)
cnmut = list(f.con1, f.con2)
} else if (nt == 7 & nb == 1){
f.con1 <- matrix(c(1, 2, 3, 4, 5, 7,
0, 0, 0, 0, 1, 0,
0, 0, 0, 1, 0, 0,
0, 0, 1, 0, 0, 0,
0, 1, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0), nrow = 6, byrow = TRUE)
f.con2 <- matrix(c(1, 2, 3, 3, 5, 7,
0, 0, 0, 2, 1, 0,
0, 0, 0, 0, 0, 0,
0, 0, 1, 0, 0, 0,
0, 1, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0), nrow = 6, byrow = TRUE) #missing 3/7
f.con3 <- matrix(c(1, 2, 2, 4, 5, 7,
0, 0, 1, 0, 1, 0,
0, 0, 1, 1, 0, 0,
0, 0, 0, 0, 0, 0,
0, 1, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0), nrow = 6, byrow = TRUE) #missing 4/7
f.con4 <- matrix(c(1, 2, 2, 3, 5, 7,
0, 0, 1, 2, 1, 0,
0, 0, 1, 0, 0, 0,
0, 0, 0, 0, 0, 0,
0, 1, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0), nrow = 6, byrow = TRUE) #missing 4/7
f.con5 <- matrix(c(1, 1, 2, 4, 5, 7,
0, 0, 1, 0, 1, 0,
0, 2, 1, 1, 0, 0,
0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0), nrow = 6, byrow = TRUE) #missing 4/7 and 5/7
f.con6 <- matrix(c(1, 1, 2, 3, 5, 7,
0, 0, 1, 2, 1, 0,
0, 2, 1, 0, 0, 0,
0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0), nrow = 6, byrow = TRUE) #missing 4/7 and 5/7
f.con7 <- matrix(c(1, 1, 2, 4, 5, 7,
0, 1, 1, 0, 1, 0,
0, 0, 1, 1, 0, 0,
0, 1, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0), nrow = 6, byrow = TRUE) #missing 5/7
f.con8 <- matrix(c(1, 1, 3, 4, 5, 7,
0, 1, 0, 0, 1, 0,
0, 0, 0, 1, 0, 0,
0, 1, 1, 0, 0, 0,
0, 0, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0), nrow = 6, byrow = TRUE) #missing 5/7
f.con9 <- matrix(c(1, 1, 2, 3, 5, 7,
0, 1, 1, 2, 1, 0,
0, 0, 1, 0, 0, 0,
0, 1, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0), nrow = 6, byrow = TRUE) #missing 5/7
f.con10 <- matrix(c(1, 1, 3, 3, 5, 7,
0, 1, 0, 2, 1, 0,
0, 0, 0, 0, 0, 0,
0, 1, 1, 0, 0, 0,
0, 0, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0), nrow = 6, byrow = TRUE) #missing 3/7 and 5/7
f.con11 <- matrix(c(1, 1, 1, 3, 5, 7,
0, 1, 3, 2, 1, 0,
0, 0, 0, 0, 0, 0,
0, 1, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0), nrow = 6, byrow = TRUE) #missing 3/7 and 5/7
cnmut = list(f.con1, f.con2, f.con3, f.con4, f.con5, f.con6, f.con7, f.con8, f.con9, f.con10, f.con11)
} else if (nt == 7 & nb == 2) {
f.con1 <- matrix(c(0, 1, 2, 3, 5, 7,
1, 1, 1, 2, 1, 0,
0, 0, 1, 0, 0, 0,
0, 1, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0), nrow = 5, byrow = TRUE)
f.con2 <- matrix(c(0, 1, 2, 4, 5, 7,
1, 1, 1, 0, 1, 0,
0, 0, 1, 1, 0, 0,
0, 1, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0), nrow = 5, byrow = TRUE)
f.con3 <- matrix(c(0, 1, 3, 4, 5, 7,
1, 1, 0, 0, 1, 0,
0, 0, 0, 1, 0, 0,
0, 1, 1, 0, 0, 0,
1, 0, 0, 0, 0, 0), nrow = 5, byrow = TRUE)
f.con4 <- matrix(c(0, 2, 3, 4, 5, 7,
1, 0, 0, 0, 1, 0,
0, 0, 0, 1, 0, 0,
0, 0, 1, 0, 0, 0,
1, 1, 0, 0, 0, 0), nrow = 5, byrow = TRUE)
f.con5 <- matrix(c(0, 1, 1, 3, 5, 7,
1, 1, 3, 2, 1, 0,
0, 0, 0, 0, 0, 0,
0, 1, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0), nrow = 5, byrow = TRUE) #missing 3/7
f.con6 <- matrix(c(0, 1, 3, 3, 5, 7,
1, 1, 0, 2, 1, 0,
0, 0, 0, 0, 0, 0,
0, 1, 1, 0, 0, 0,
1, 0, 0, 0, 0, 0), nrow = 5, byrow = TRUE) #missing 3/7
f.con7 <- matrix(c(0, 2, 3, 3, 5, 7,
1, 0, 0, 2, 1, 0,
0, 0, 0, 0, 0, 0,
0, 0, 1, 0, 0, 0,
1, 1, 0, 0, 0, 0), nrow = 5, byrow = TRUE) #missing 3/7
f.con8 <- matrix(c(0, 0, 2, 4, 5, 7,
1, 2, 1, 0, 1, 0,
0, 1, 1, 1, 0, 0,
0, 0, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0), nrow = 5, byrow = TRUE) #missing 4/7
f.con9 <- matrix(c(0, 2, 2, 3, 5, 7,
1, 0, 1, 2, 1, 0,
0, 0, 1, 0, 0, 0,
0, 0, 0, 0, 0, 0,
1, 1, 0, 0, 0, 0), nrow = 5, byrow = TRUE) #missing 4/7
f.con10 <- matrix(c(0, 2, 2, 4, 5, 7,
1, 0, 1, 0, 1, 0,
0, 0, 1, 1, 0, 0,
0, 0, 0, 0, 0, 0,
1, 1, 0, 0, 0, 0), nrow = 5, byrow = TRUE) #missing 4/7
f.con11 <- matrix(c(0, 0, 2, 3, 5, 7,
1, 2, 1, 2, 1, 0,
0, 1, 1, 0, 0, 0,
0, 0, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0), nrow = 5, byrow = TRUE) #missing 4/7
f.con12 <- matrix(c(0, 0, 1, 3, 5, 7,
1, 2, 3, 2, 1, 0,
0, 1, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0), nrow = 5, byrow = TRUE) #missing 4/7
cnmut = list(f.con1, f.con2, f.con3, f.con4, f.con5, f.con6, f.con7, f.con8, f.con9, f.con10, f.con11, f.con12)
} else if (nt == 7 & nb == 3) {
f.con1 <- matrix(c(0, 1, 2, 4, 5, 7,
0, 0, 1, 0, 1, 0,
1, 2, 1, 1, 0, 0,
1, 0, 0, 0, 0, 0), nrow = 4, byrow = TRUE)
f.con2 <- matrix(c(0, 1, 2, 3, 5, 7,
0, 0, 1, 2, 1, 0,
1, 2, 1, 0, 0, 0,
1, 0, 0, 0, 0, 0), nrow = 4, byrow = TRUE)
f.con3 <- matrix(c(0, 1, 2, 4, 5, 7,
0, 1, 1, 0, 1, 0,
1, 0, 1, 1, 0, 0,
1, 1, 0, 0, 0, 0), nrow = 4, byrow = TRUE)
f.con4 <- matrix(c(0, 1, 3, 4, 5, 7,
0, 1, 0, 0, 1, 0,
1, 0, 0, 1, 0, 0,
1, 1, 1, 0, 0, 0), nrow = 4, byrow = TRUE)
f.con5 <- matrix(c(0, 1, 2, 3, 5, 7,
0, 1, 1, 2, 1, 0,
1, 0, 1, 0, 0, 0,
1, 1, 0, 0, 0, 0), nrow = 4, byrow = TRUE)
f.con6 <- matrix(c(0, 0, 2, 4, 5, 7,
0, 2, 1, 0, 1, 0,
1, 1, 1, 1, 0, 0,
1, 0, 0, 0, 0, 0), nrow = 4, byrow = TRUE)
f.con7 <- matrix(c(0, 1, 3, 3, 5, 7,
0, 1, 0, 2, 1, 0,
1, 0, 0, 0, 0, 0,
1, 1, 1, 0, 0, 0), nrow = 4, byrow = TRUE)
f.con8 <- matrix(c(0, 0, 2, 3, 5, 7,
0, 2, 1, 2, 1, 0,
1, 1, 1, 0, 0, 0,
1, 0, 0, 0, 0, 0), nrow = 4, byrow = TRUE)
f.con9 <- matrix(c(0, 0, 1, 3, 5, 7,
0, 2, 3, 2, 1, 0,
1, 1, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0), nrow = 4, byrow = TRUE)
f.con10 <- matrix(c(0, 1, 1, 3, 5, 7,
0, 1, 3, 2, 1, 0,
1, 0, 0, 0, 0, 0,
1, 1, 0, 0, 0, 0), nrow = 4, byrow = TRUE)
cnmut = list(f.con1, f.con2, f.con3, f.con4, f.con5, f.con6, f.con7, f.con8, f.con9, f.con10)
} else if (nt == 7 & nb == 0) {
f.con1 <- matrix(c(0, 1, 2, 3, 4, 5, 7,
0, 0, 0, 0, 0, 1, 0,
0, 0, 0, 0, 1, 0, 0,
0, 0, 0, 1, 0, 0, 0,
0, 0, 1, 0, 0, 0, 0,
0, 1, 0, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0, 0), nrow = 7, byrow = TRUE)
f.con2 <- matrix(c(0, 1, 2, 3, 3, 5, 7,
0, 0, 0, 0, 2, 1, 0,
0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 1, 0, 0, 0,
0, 0, 1, 0, 0, 0, 0,
0, 1, 0, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0, 0), nrow = 7, byrow = TRUE) #missing 3/7
f.con3 <- matrix(c(0, 1, 2, 2, 3, 5, 7,
0, 0, 0, 1, 2, 1, 0,
0, 0, 0, 1, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0,
0, 0, 1, 0, 0, 0, 0,
0, 1, 0, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0, 0), nrow = 7, byrow = TRUE) #missing 4/7
f.con4 <- matrix(c(0, 1, 1, 2, 3, 5, 7,
0, 0, 1, 1, 2, 1, 0,
0, 0, 0, 1, 0, 0, 0,
0, 0, 1, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0,
0, 1, 0, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0, 0), nrow = 7, byrow = TRUE) #missing 4/7
f.con5 <- matrix(c(0, 0, 1, 2, 3, 5, 7,
0, 1, 1, 1, 2, 1, 0,
0, 0, 0, 1, 0, 0, 0,
0, 0, 1, 0, 0, 0, 0,
0, 1, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0, 0), nrow = 7, byrow = TRUE) #missing 5/7
f.con6 <- matrix(c(0, 1, 2, 2, 3, 5, 7,
0, 0, 0, 1, 2, 1, 0,
0, 0, 0, 1, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0,
0, 0, 1, 0, 0, 0, 0,
0, 1, 0, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0, 0), nrow = 7, byrow = TRUE) #missing 4/7
f.con7 <- matrix(c(0, 1, 2, 2, 4, 5, 7,
0, 0, 0, 1, 0, 1, 0,
0, 0, 0, 1, 1, 0, 0,
0, 0, 0, 0, 0, 0, 0,
0, 0, 1, 0, 0, 0, 0,
0, 1, 0, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0, 0), nrow = 7, byrow = TRUE) #missing 4/7
f.con8 <- matrix(c(0, 1, 1, 1, 3, 5, 7,
0, 0, 1, 3, 2, 1, 0,
0, 0, 0, 0, 0, 0, 0,
0, 0, 1, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0,
0, 1, 0, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0, 0), nrow = 7, byrow = TRUE) #missing 3/7 and 5/7
f.con9 <- matrix(c(0, 1, 1, 2, 4, 5, 7,
0, 0, 1, 1, 0, 1, 0,
0, 0, 0, 1, 1, 0, 0,
0, 0, 1, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0,
0, 1, 0, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0, 0), nrow = 7, byrow = TRUE) #missing 5/7
f.con10 <- matrix(c(0, 0, 1, 2, 4, 5, 7,
0, 1, 1, 1, 0, 1, 0,
0, 0, 0, 1, 1, 0, 0,
0, 0, 1, 0, 0, 0, 0,
0, 1, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0, 0), nrow = 7, byrow = TRUE) #missing 6/7
f.con11 <- matrix(c(0, 1, 1, 2, 4, 5, 7,
0, 0, 0, 1, 0, 1, 0,
0, 0, 2, 1, 1, 0, 0,
0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0,
0, 1, 0, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0, 0), nrow = 7, byrow = TRUE) #missing 4/7 and 5/7
f.con12 <- matrix(c(0, 0, 0, 1, 3, 5, 7,
0, 1, 2, 3, 2, 1, 0,
0, 0, 1, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0,
0, 1, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0, 0), nrow = 7, byrow = TRUE) #missing 4/7 and 6/7
f.con13 <- matrix(c(0, 1, 1, 3, 3, 5, 7,
0, 0, 1, 0, 2, 1, 0,
0, 0, 0, 0, 0, 0, 0,
0, 0, 1, 1, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0,
0, 1, 0, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0, 0), nrow = 7, byrow = TRUE) #missing 3/7 and 5/7
f.con14 <- matrix(c(0, 1, 1, 3, 4, 5, 7,
0, 0, 1, 0, 0, 1, 0,
0, 0, 0, 0, 1, 0, 0,
0, 0, 1, 1, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0,
0, 1, 0, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0, 0), nrow = 7, byrow = TRUE) #missing 5/7
f.con15 <- matrix(c(0, 1, 1, 2, 3, 5, 7,
0, 0, 0, 1, 2, 1, 0,
0, 0, 2, 1, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0,
0, 1, 0, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0, 0), nrow = 7, byrow = TRUE) #missing 4/7 and 5/7
f.con16 <- matrix(c(0, 0, 1, 1, 3, 5, 7,
0, 1, 1, 3, 2, 1, 0,
0, 0, 0, 0, 0, 0, 0,
0, 0, 1, 0, 0, 0, 0,
0, 1, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0, 0), nrow = 7, byrow = TRUE) #missing 3/7 and 6/7
f.con17 <- matrix(c(0, 0, 0, 2, 4, 5, 7,
0, 0, 2, 1, 0, 1, 0,
0, 1, 1, 1, 1, 0, 0,
0, 1, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0, 0), nrow = 7, byrow = TRUE) #missing 5/7 and 6/7
f.con18 <- matrix(c(0, 0, 0, 2, 4, 5, 7,
0, 1, 2, 1, 0, 1, 0,
0, 0, 1, 1, 1, 0, 0,
0, 0, 0, 0, 0, 0, 0,
0, 1, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0, 0), nrow = 7, byrow = TRUE) #missing 4/7 and 6/7
f.con19 <- matrix(c(0, 0, 0, 2, 3, 5, 7,
0, 0, 2, 1, 2, 1, 0,
0, 1, 1, 1, 0, 0, 0,
0, 1, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0, 0), nrow = 7, byrow = TRUE) #missing 5/7 and 6/7
f.con20 <- matrix(c(0, 0, 1, 3, 4, 5, 7,
0, 1, 1, 0, 0, 1, 0,
0, 0, 0, 0, 1, 0, 0,
0, 0, 1, 1, 0, 0, 0,
0, 1, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0, 0), nrow = 7, byrow = TRUE) #missing 6/7
f.con21 <- matrix(c(0, 0, 1, 3, 3, 5, 7,
0, 1, 1, 0, 2, 1, 0,
0, 0, 0, 0, 0, 0, 0,
0, 0, 1, 1, 0, 0, 0,
0, 1, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0, 0), nrow = 7, byrow = TRUE) #missing 3/7 and 6/7
f.con22 <- matrix(c(0, 0, 1, 2, 4, 5, 7,
0, 0, 1, 1, 0, 1, 0,
0, 1, 1, 1, 1, 0, 0,
0, 1, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0, 0), nrow = 7, byrow = TRUE) #missing 5/7 and 6/7
f.con23 <- matrix(c(0, 0, 1, 2, 4, 5, 7,
0, 0, 0, 1, 0, 1, 0,
0, 1, 2, 1, 1, 0, 0,
0, 1, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0, 0), nrow = 7, byrow = TRUE) #missing 5/7 and 6/7
f.con24 <- matrix(c(0, 0, 1, 2, 3, 5, 7,
0, 0, 0, 1, 2, 1, 0,
0, 1, 2, 1, 0, 0, 0,
0, 1, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0, 0), nrow = 7, byrow = TRUE) #missing 5/7 and 6/7
f.con25 <- matrix(c(0, 0, 0, 2, 3, 5, 7,
0, 1, 2, 1, 2, 1, 0,
0, 0, 1, 1, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0,
0, 1, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0, 0), nrow = 7, byrow = TRUE) #missing 4/7 and 6/7
f.con26 <- matrix(c(0, 0, 0, 1, 3, 5, 7,
0, 0, 2, 3, 2, 1, 0,
0, 1, 1, 0, 0, 0, 0,
0, 1, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0, 0), nrow = 7, byrow = TRUE) #missing 5/7 and 6/7
f.con27 <- matrix(c(0, 0, 1, 1, 3, 5, 7,
0, 0, 1, 3, 2, 1, 0,
0, 1, 0, 0, 0, 0, 0,
0, 1, 1, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0, 0), nrow = 7, byrow = TRUE) #missing 5/7 and 6/7
f.con28 <- matrix(c(0, 0, 1, 2, 3, 5, 7,
0, 0, 1, 1, 2, 1, 0,
0, 1, 0, 1, 0, 0, 0,
0, 1, 1, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0, 0), nrow = 7, byrow = TRUE) #missing 5/7 and 6/7
f.con29 <- matrix(c(0, 0, 2, 3, 4, 5, 7,
0, 1, 0, 0, 0, 1, 0,
0, 0, 0, 0, 1, 0, 0,
0, 0, 0, 1, 0, 0, 0,
0, 1, 1, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0, 0), nrow = 7, byrow = TRUE) #missing 6/7
f.con30 <- matrix(c(0, 0, 2, 3, 3, 5, 7,
0, 1, 0, 0, 2, 1, 0,
0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 1, 0, 0, 0,
0, 1, 1, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0, 0), nrow = 7, byrow = TRUE) #missing 3/7 and 6/7
f.con31 <- matrix(c(0, 0, 2, 2, 3, 5, 7,
0, 1, 0, 1, 2, 1, 0,
0, 0, 0, 1, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0,
0, 1, 1, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0, 0), nrow = 7, byrow = TRUE) #missing 4/7 and 6/7
f.con32 <- matrix(c(0, 0, 2, 2, 4, 5, 7,
0, 1, 0, 1, 1, 1, 0,
0, 0, 0, 1, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0,
0, 1, 1, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0, 0), nrow = 7, byrow = TRUE) #missing 4/7 and 6/7
f.con33 <- matrix(c(0, 0, 1, 2, 4, 5, 7,
0, 0, 0, 1, 0, 1, 0,
0, 1, 2, 1, 1, 0, 0,
0, 1, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0, 0), nrow = 7, byrow = TRUE) #missing 5/7 and 6/7
cnmut = list(f.con1, f.con2, f.con3, f.con4, f.con5, f.con6, f.con7, f.con8, f.con9, f.con10, f.con11, f.con12, f.con13, f.con14, f.con15, f.con16, f.con17, f.con18, f.con19, f.con20, f.con21, f.con22, f.con23, f.con24, f.con25, f.con26, f.con27, f.con28, f.con29, f.con30, f.con31, f.con32, f.con33)
} else if (nt == 8 & nb == 1) {
f.con1 <- matrix(c(1, 2, 3, 4, 5, 6, 8,
0, 0, 0, 0, 0, 1, 0,
0, 0, 0, 0, 1, 0, 0,
0, 0, 0, 1, 0, 0, 0,
0, 0, 1, 0, 0, 0, 0,
0, 1, 0, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0, 0), nrow = 7, byrow = TRUE)
f.con2 <- matrix(c(1, 2, 3, 4, 4, 6, 8,
0, 0, 0, 0, 2, 1, 0,
0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 1, 0, 0, 0,
0, 0, 1, 0, 0, 0, 0,
0, 1, 0, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0, 0), nrow = 7, byrow = TRUE) #missing 3/7
f.con3 <- matrix(c(1, 2, 3, 4, 4, 6, 8,
0, 0, 0, 1, 2, 1, 0,
0, 0, 0, 1, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0,
0, 0, 1, 0, 0, 0, 0,
0, 1, 0, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0, 0), nrow = 7, byrow = TRUE) #missing 4/7
f.con4 <- matrix(c(1, 2, 2, 3, 4, 6, 8,
0, 0, 1, 1, 2, 1, 0,
0, 0, 0, 1, 0, 0, 0,
0, 0, 1, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0,
0, 1, 0, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0, 0), nrow = 7, byrow = TRUE) #missing 4/7
f.con5 <- matrix(c(1, 1, 2, 3, 4, 6, 8,
0, 1, 1, 1, 2, 1, 0,
0, 0, 0, 1, 0, 0, 0,
0, 0, 1, 0, 0, 0, 0,
0, 1, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0, 0), nrow = 7, byrow = TRUE) #missing 5/7
f.con6 <- matrix(c(1, 2, 3, 3, 4, 6, 8,
0, 0, 0, 1, 2, 1, 0,
0, 0, 0, 1, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0,
0, 0, 1, 0, 0, 0, 0,
0, 1, 0, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0, 0), nrow = 7, byrow = TRUE) #missing 4/7
f.con7 <- matrix(c(1, 2, 3, 3, 5, 6, 8,
0, 0, 0, 1, 0, 1, 0,
0, 0, 0, 1, 1, 0, 0,
0, 0, 0, 0, 0, 0, 0,
0, 0, 1, 0, 0, 0, 0,
0, 1, 0, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0, 0), nrow = 7, byrow = TRUE) #missing 4/7
f.con8 <- matrix(c(1, 2, 2, 2, 4, 6, 8,
0, 0, 1, 3, 2, 1, 0,
0, 0, 0, 0, 0, 0, 0,
0, 0, 1, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0,
0, 1, 0, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0, 0), nrow = 7, byrow = TRUE) #missing 3/7 and 5/7
f.con9 <- matrix(c(1, 2, 2, 3, 5, 6, 8,
0, 0, 1, 1, 0, 1, 0,
0, 0, 0, 1, 1, 0, 0,
0, 0, 1, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0,
0, 1, 0, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0, 0), nrow = 7, byrow = TRUE) #missing 5/7
f.con10 <- matrix(c(1, 1, 2, 3, 5, 6, 8,
0, 1, 1, 1, 0, 1, 0,
0, 0, 0, 1, 1, 0, 0,
0, 0, 1, 0, 0, 0, 0,
0, 1, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0, 0), nrow = 7, byrow = TRUE) #missing 6/7
f.con11 <- matrix(c(1, 2, 2, 3, 5, 6, 8,
0, 0, 0, 1, 0, 1, 0,
0, 0, 2, 1, 1, 0, 0,
0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0,
0, 1, 0, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0, 0), nrow = 7, byrow = TRUE) #missing 4/7 and 5/7
f.con12 <- matrix(c(1, 1, 1, 2, 4, 6, 8,
0, 1, 2, 3, 2, 1, 0,
0, 0, 1, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0,
0, 1, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0, 0), nrow = 7, byrow = TRUE) #missing 4/7 and 6/7
f.con13 <- matrix(c(1, 2, 2, 4, 4, 6, 8,
0, 0, 1, 0, 2, 1, 0,
0, 0, 0, 0, 0, 0, 0,
0, 0, 1, 1, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0,
0, 1, 0, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0, 0), nrow = 7, byrow = TRUE) #missing 3/7 and 5/7
f.con14 <- matrix(c(1, 2, 2, 4, 5, 6, 8,
0, 0, 1, 0, 0, 1, 0,
0, 0, 0, 0, 1, 0, 0,
0, 0, 1, 1, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0,
0, 1, 0, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0, 0), nrow = 7, byrow = TRUE) #missing 5/7
f.con15 <- matrix(c(1, 2, 2, 3, 4, 6, 8,
0, 0, 0, 1, 2, 1, 0,
0, 0, 2, 1, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0,
0, 1, 0, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0, 0), nrow = 7, byrow = TRUE) #missing 4/7 and 5/7
f.con16 <- matrix(c(1, 1, 2, 2, 4, 6, 8,
0, 1, 1, 3, 2, 1, 0,
0, 0, 0, 0, 0, 0, 0,
0, 0, 1, 0, 0, 0, 0,
0, 1, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0, 0), nrow = 7, byrow = TRUE) #missing 3/7 and 6/7
f.con17 <- matrix(c(1, 1, 1, 3, 5, 6, 8,
0, 0, 2, 1, 0, 1, 0,
0, 1, 1, 1, 1, 0, 0,
0, 1, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0, 0), nrow = 7, byrow = TRUE) #missing 5/7 and 6/7
f.con18 <- matrix(c(1, 1, 1, 3, 5, 6, 8,
0, 1, 2, 1, 0, 1, 0,
0, 0, 1, 1, 1, 0, 0,
0, 0, 0, 0, 0, 0, 0,
0, 1, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0, 0), nrow = 7, byrow = TRUE) #missing 4/7 and 6/7
f.con19 <- matrix(c(1, 1, 1, 3, 4, 6, 8,
0, 0, 2, 1, 2, 1, 0,
0, 1, 1, 1, 0, 0, 0,
0, 1, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0, 0), nrow = 7, byrow = TRUE) #missing 5/7 and 6/7
f.con20 <- matrix(c(1, 1, 2, 4, 5, 6, 8,
0, 1, 1, 0, 0, 1, 0,
0, 0, 0, 0, 1, 0, 0,
0, 0, 1, 1, 0, 0, 0,
0, 1, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0, 0), nrow = 7, byrow = TRUE) #missing 6/7
f.con21 <- matrix(c(1, 1, 2, 4, 4, 6, 8,
0, 1, 1, 0, 2, 1, 0,
0, 0, 0, 0, 0, 0, 0,
0, 0, 1, 1, 0, 0, 0,
0, 1, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0, 0), nrow = 7, byrow = TRUE) #missing 3/7 and 6/7
f.con22 <- matrix(c(1, 1, 2, 3, 5, 6, 8,
0, 0, 1, 1, 0, 1, 0,
0, 1, 1, 1, 1, 0, 0,
0, 1, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0, 0), nrow = 7, byrow = TRUE) #missing 5/7 and 6/7
f.con23 <- matrix(c(1, 1, 2, 3, 5, 6, 8,
0, 0, 0, 1, 0, 1, 0,
0, 1, 2, 1, 1, 0, 0,
0, 1, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0, 0), nrow = 7, byrow = TRUE) #missing 5/7 and 6/7
f.con24 <- matrix(c(1, 1, 2, 3, 4, 6, 8,
0, 0, 0, 1, 2, 1, 0,
0, 1, 2, 1, 0, 0, 0,
0, 1, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0, 0), nrow = 7, byrow = TRUE) #missing 5/7 and 6/7
f.con25 <- matrix(c(1, 1, 1, 3, 4, 6, 8,
0, 1, 2, 1, 2, 1, 0,
0, 0, 1, 1, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0,
0, 1, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0, 0), nrow = 7, byrow = TRUE) #missing 4/7 and 6/7
f.con26 <- matrix(c(1, 1, 1, 2, 4, 6, 8,
0, 0, 2, 3, 2, 1, 0,
0, 1, 1, 0, 0, 0, 0,
0, 1, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0, 0), nrow = 7, byrow = TRUE) #missing 5/7 and 6/7
f.con27 <- matrix(c(1, 1, 2, 2, 4, 6, 8,
0, 0, 1, 3, 2, 1, 0,
0, 1, 0, 0, 0, 0, 0,
0, 1, 1, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0, 0), nrow = 7, byrow = TRUE) #missing 5/7 and 6/7
f.con28 <- matrix(c(1, 1, 2, 3, 4, 6, 8,
0, 0, 1, 1, 2, 1, 0,
0, 1, 0, 1, 0, 0, 0,
0, 1, 1, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0, 0), nrow = 7, byrow = TRUE) #missing 5/7 and 6/7
f.con29 <- matrix(c(1, 1, 3, 4, 5, 6, 8,
0, 1, 0, 0, 0, 1, 0,
0, 0, 0, 0, 1, 0, 0,
0, 0, 0, 1, 0, 0, 0,
0, 1, 1, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0, 0), nrow = 7, byrow = TRUE) #missing 6/7
f.con30 <- matrix(c(1, 1, 3, 4, 4, 6, 8,
0, 1, 0, 0, 2, 1, 0,
0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 1, 0, 0, 0,
0, 1, 1, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0, 0), nrow = 7, byrow = TRUE) #missing 3/7 and 6/7
f.con31 <- matrix(c(1, 1, 3, 3, 4, 6, 8,
0, 1, 0, 1, 2, 1, 0,
0, 0, 0, 1, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0,
0, 1, 1, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0, 0), nrow = 7, byrow = TRUE) #missing 4/7 and 6/7
f.con32 <- matrix(c(1, 1, 3, 3, 5, 6, 8,
0, 1, 0, 1, 1, 1, 0,
0, 0, 0, 1, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0,
0, 1, 1, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0, 0), nrow = 7, byrow = TRUE) #missing 4/7 and 6/7
f.con33 <- matrix(c(1, 1, 2, 3, 5, 6, 8,
0, 0, 0, 1, 0, 1, 0,
0, 1, 2, 1, 1, 0, 0,
0, 1, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0, 0), nrow = 7, byrow = TRUE) #missing 5/7 and 6/7
cnmut = list(f.con1, f.con2, f.con3, f.con4, f.con5, f.con6, f.con7, f.con8, f.con9, f.con10, f.con11, f.con12, f.con13, f.con14, f.con15, f.con16, f.con17, f.con18, f.con19, f.con20, f.con21, f.con22, f.con23, f.con24, f.con25, f.con26, f.con27, f.con28, f.con29, f.con30, f.con31, f.con32, f.con33)
}
return(cnmut)
}
SunPathLab/Butte documentation built on Sept. 19, 2023, 9:42 a.m.
R Package Documentation
Browse R Packages
We want your feedback!
Note that we can't provide technical support on individual packages. You should contact the package authors for that.
Defines functions cnmutHistory
Documented in cnmutHistory
#' generate history matrix in relating SCNA to SSNVs
#'
#' Given a SCNA configuration Nt (total copy) and Nb (minor copy)
#' this function produces the possible history matrices in relating
#' CN timing (time period for each stage) to the burden of SSNVs
#' at distinct allele states.
#'
#' @param nt total copy number
#' @param nb copy number of the minor allele
#' @return A list of possible matrices
#' @export
cnmutHistory <- function(nt, nb) {
cnmut <- list()
if (nt == 2 & nb == 0) {
f.con1 <- matrix(c(0, 2,
1, 0), nrow = 2, byrow = TRUE)
cnmut = list(f.con1)
} else if (nt == 3 & nb == 0){
f.con1 <- matrix(c(0, 1, 3,
0, 1, 0,
1, 0, 0), nrow = 3, byrow = TRUE)
cnmut = list(f.con1)
} else if (nt == 3 & nb == 1){
f.con1 <- matrix(c(1, 3,
1, 0), nrow = 2, byrow = TRUE)
cnmut = list(f.con1)
} else if (nt == 4 & nb == 0){
f.con1 <- matrix(c(0, 1, 2, 4,
0, 0, 1, 0,
0, 1, 0, 0,
1, 0, 0, 0), nrow = 4, byrow = TRUE)
f.con2 <- matrix(c(0, 0, 2, 4,
0, 2, 1, 0,
0, 0, 0, 0,
1, 0, 0, 0), nrow = 4, byrow = TRUE) #missing 3/4
cnmut = list(f.con1, f.con2)
} else if (nt == 4 & nb == 1){
f.con1 <- matrix(c(1, 2, 4,
0, 1, 0,
1, 0, 0), nrow = 3, byrow = TRUE)
cnmut = list(f.con1)
} else if (nt == 4 & nb == 2){
f.con1 <- matrix(c(0, 2, 4,
2, 1, 0), nrow = 2, byrow = TRUE)
#f.con2 <- matrix(c(0, 4,
# 2, 0), nrow = 2, byrow = TRUE) #syn gain
#cnmut = list(f.con1, f.con2)
cnmut = list(f.con1)
} else if (nt == 5 & nb == 0){
f.con1 <- matrix(c(0, 1, 2, 3, 5,
0, 0, 0, 1, 0,
0, 0, 1, 0, 0,
0, 1, 0, 0, 0,
1, 0, 0, 0, 0), nrow = 5, byrow = TRUE)
f.con2 <- matrix(c(0, 1, 1, 3, 5,
0, 0, 2, 1, 0,
0, 0, 0, 0, 0,
0, 1, 0, 0, 0,
1, 0, 0, 0, 0), nrow = 5, byrow = TRUE) #missing 3/5
f.con3 <- matrix(c(0, 0, 1, 3, 5,
0, 1, 2, 1, 0,
0, 1, 0, 0, 0,
0, 0, 0, 0, 0,
1, 0, 0, 0, 0), nrow = 5, byrow = TRUE) #missing 4/5
f.con4 <- matrix(c(0, 0, 2, 3, 5,
0, 1, 0, 1, 0,
0, 1, 1, 0, 0,
0, 0, 0, 0, 0,
1, 0, 0, 0, 0), nrow = 5, byrow = TRUE) #missing 4/5
cnmut = list(f.con1, f.con2, f.con3, f.con4)
} else if (nt == 5 & nb == 1){
f.con1 <- matrix(c(1, 2, 3, 5,
0, 0, 1, 0,
0, 1, 0, 0,
1, 0, 0, 0), nrow = 4, byrow = TRUE)
f.con2 <- matrix(c(1, 1, 3, 5,
0, 2, 1, 0,
0, 0, 0, 0,
1, 0, 0, 0), nrow = 4, byrow = TRUE) #missing 3/5
cnmut = list(f.con1, f.con2)
} else if (nt == 5 & nb == 2){
f.con1 <- matrix(c(0, 2, 3, 5,
1, 0, 1, 0,
1, 1, 0, 0), nrow = 3, byrow = TRUE)
f.con2 <- matrix(c(0, 1, 3, 5,
1, 2, 1, 0,
1, 0, 0, 0), nrow = 3, byrow = TRUE)
cnmut = list(f.con1, f.con2)
} else if (nt == 6 & nb == 0){
f.con1 <- matrix(c(0, 1, 2, 3, 4, 6,
0, 0, 0, 0, 1, 0,
0, 0, 0, 1, 0, 0,
0, 0, 1, 0, 0, 0,
0, 1, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0), nrow = 6, byrow = TRUE)
f.con2 <- matrix(c(0, 1, 2, 2, 4, 6,
0, 0, 0, 2, 1, 0,
0, 0, 0, 0, 0, 0,
0, 0, 1, 0, 0, 0,
0, 1, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0), nrow = 6, byrow = TRUE) #missing 3/6
f.con3 <- matrix(c(0, 1, 1, 3, 4, 6,
0, 0, 1, 0, 1, 0,
0, 0, 1, 1, 0, 0,
0, 0, 0, 0, 0, 0,
0, 1, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0), nrow = 6, byrow = TRUE) #missing 4/6
f.con4 <- matrix(c(0, 1, 1, 2, 4, 6,
0, 0, 1, 2, 1, 0,
0, 0, 1, 0, 0, 0,
0, 0, 0, 0, 0, 0,
0, 1, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0), nrow = 6, byrow = TRUE) #missing 4/6
f.con5 <- matrix(c(0, 0, 1, 3, 4, 6,
0, 0, 1, 0, 1, 0,
0, 2, 1, 1, 0, 0,
0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0), nrow = 6, byrow = TRUE) #missing 4/6 and 5/6
f.con6 <- matrix(c(0, 0, 1, 2, 4, 6,
0, 0, 1, 2, 1, 0,
0, 2, 1, 0, 0, 0,
0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0), nrow = 6, byrow = TRUE) #missing 4/6 and 5/6
f.con7 <- matrix(c(0, 0, 1, 3, 4, 6,
0, 1, 1, 0, 1, 0,
0, 0, 1, 1, 0, 0,
0, 1, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0), nrow = 6, byrow = TRUE) #missing 5/6
f.con8 <- matrix(c(0, 0, 2, 3, 4, 6,
0, 1, 0, 0, 1, 0,
0, 0, 0, 1, 0, 0,
0, 1, 1, 0, 0, 0,
0, 0, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0), nrow = 6, byrow = TRUE) #missing 5/6
f.con9 <- matrix(c(0, 0, 1, 2, 4, 6,
0, 1, 1, 2, 1, 0,
0, 0, 1, 0, 0, 0,
0, 1, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0), nrow = 6, byrow = TRUE) #missing 5/6
f.con10 <- matrix(c(0, 0, 2, 2, 4, 6,
0, 1, 0, 2, 1, 0,
0, 0, 0, 0, 0, 0,
0, 1, 1, 0, 0, 0,
0, 0, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0), nrow = 6, byrow = TRUE) #missing 3/6 and 5/6
f.con11 <- matrix(c(0, 0, 0, 2, 4, 6,
0, 1, 3, 2, 1, 0,
0, 0, 0, 0, 0, 0,
0, 1, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0), nrow = 6, byrow = TRUE) #missing 3/6 and 5/6
cnmut = list(f.con1, f.con2, f.con3, f.con4, f.con5, f.con6, f.con7, f.con8, f.con9, f.con10, f.con11)
} else if (nt == 6 & nb == 1){
f.con1 <- matrix(c(1, 2, 3, 4, 6,
0, 0, 0, 1, 0,
0, 0, 1, 0, 0,
0, 1, 0, 0, 0,
1, 0, 0, 0, 0), nrow = 5, byrow = TRUE)
f.con2 <- matrix(c(1, 2, 2, 4, 6,
0, 0, 2, 1, 0,
0, 0, 0, 0, 0,
0, 1, 0, 0, 0,
1, 0, 0, 0, 0), nrow = 5, byrow = TRUE) #missing 3/6
f.con3 <- matrix(c(1, 1, 2, 4, 6,
0, 1, 2, 1, 0,
0, 1, 0, 0, 0,
0, 0, 0, 0, 0,
1, 0, 0, 0, 0), nrow = 5, byrow = TRUE) #missing 4/6
f.con4 <- matrix(c(1, 1, 3, 4, 6,
0, 1, 0, 1, 0,
0, 1, 1, 0, 0,
0, 0, 0, 0, 0,
1, 0, 0, 0, 0), nrow = 5, byrow = TRUE) #missing 4/6
cnmut = list(f.con1, f.con2, f.con3, f.con4)
} else if (nt == 6 & nb == 2){
f.con1 <- matrix(c(0, 1, 3, 4, 6,
1, 1, 0, 1, 0,
0, 1, 1, 0, 0,
1, 0, 0, 0, 0), nrow = 4, byrow = TRUE)
f.con2 <- matrix(c(0, 2, 3, 4, 6,
1, 0, 0, 1, 0,
0, 0, 1, 0, 0,
1, 1, 0, 0, 0), nrow = 4, byrow = TRUE)
f.con3 <- matrix(c(0, 1, 2, 4, 6,
1, 1, 2, 1, 0,
0, 1, 0, 0, 0,
1, 0, 0, 0, 0), nrow = 4, byrow = TRUE)
f.con4 <- matrix(c(0, 2, 2, 4, 6,
1, 0, 2, 1, 0,
0, 0, 0, 0, 0,
1, 1, 0, 0, 0), nrow = 4, byrow = TRUE) #missing 3/6
f.con5 <- matrix(c(0, 0, 2, 4, 6,
1, 3, 2, 1, 0,
0, 0, 0, 0, 0,
1, 0, 0, 0, 0), nrow = 4, byrow = TRUE) #missing 3/6
cnmut = list(f.con1, f.con2, f.con3, f.con4, f.con5)
} else if (nt == 6 & nb == 3){
f.con1 <- matrix(c(0, 1, 3, 4, 6,
0, 1, 0, 1, 0,
2, 1, 1, 0, 0), nrow = 3, byrow = TRUE)
f.con2 <- matrix(c(0, 1, 2, 4, 6,
0, 1, 2, 1, 0,
2, 1, 0, 0, 0), nrow = 3, byrow = TRUE)
cnmut = list(f.con1, f.con2)
} else if (nt == 7 & nb == 1){
f.con1 <- matrix(c(1, 2, 3, 4, 5, 7,
0, 0, 0, 0, 1, 0,
0, 0, 0, 1, 0, 0,
0, 0, 1, 0, 0, 0,
0, 1, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0), nrow = 6, byrow = TRUE)
f.con2 <- matrix(c(1, 2, 3, 3, 5, 7,
0, 0, 0, 2, 1, 0,
0, 0, 0, 0, 0, 0,
0, 0, 1, 0, 0, 0,
0, 1, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0), nrow = 6, byrow = TRUE) #missing 3/7
f.con3 <- matrix(c(1, 2, 2, 4, 5, 7,
0, 0, 1, 0, 1, 0,
0, 0, 1, 1, 0, 0,
0, 0, 0, 0, 0, 0,
0, 1, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0), nrow = 6, byrow = TRUE) #missing 4/7
f.con4 <- matrix(c(1, 2, 2, 3, 5, 7,
0, 0, 1, 2, 1, 0,
0, 0, 1, 0, 0, 0,
0, 0, 0, 0, 0, 0,
0, 1, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0), nrow = 6, byrow = TRUE) #missing 4/7
f.con5 <- matrix(c(1, 1, 2, 4, 5, 7,
0, 0, 1, 0, 1, 0,
0, 2, 1, 1, 0, 0,
0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0), nrow = 6, byrow = TRUE) #missing 4/7 and 5/7
f.con6 <- matrix(c(1, 1, 2, 3, 5, 7,
0, 0, 1, 2, 1, 0,
0, 2, 1, 0, 0, 0,
0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0), nrow = 6, byrow = TRUE) #missing 4/7 and 5/7
f.con7 <- matrix(c(1, 1, 2, 4, 5, 7,
0, 1, 1, 0, 1, 0,
0, 0, 1, 1, 0, 0,
0, 1, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0), nrow = 6, byrow = TRUE) #missing 5/7
f.con8 <- matrix(c(1, 1, 3, 4, 5, 7,
0, 1, 0, 0, 1, 0,
0, 0, 0, 1, 0, 0,
0, 1, 1, 0, 0, 0,
0, 0, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0), nrow = 6, byrow = TRUE) #missing 5/7
f.con9 <- matrix(c(1, 1, 2, 3, 5, 7,
0, 1, 1, 2, 1, 0,
0, 0, 1, 0, 0, 0,
0, 1, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0), nrow = 6, byrow = TRUE) #missing 5/7
f.con10 <- matrix(c(1, 1, 3, 3, 5, 7,
0, 1, 0, 2, 1, 0,
0, 0, 0, 0, 0, 0,
0, 1, 1, 0, 0, 0,
0, 0, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0), nrow = 6, byrow = TRUE) #missing 3/7 and 5/7
f.con11 <- matrix(c(1, 1, 1, 3, 5, 7,
0, 1, 3, 2, 1, 0,
0, 0, 0, 0, 0, 0,
0, 1, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0), nrow = 6, byrow = TRUE) #missing 3/7 and 5/7
cnmut = list(f.con1, f.con2, f.con3, f.con4, f.con5, f.con6, f.con7, f.con8, f.con9, f.con10, f.con11)
} else if (nt == 7 & nb == 2) {
f.con1 <- matrix(c(0, 1, 2, 3, 5, 7,
1, 1, 1, 2, 1, 0,
0, 0, 1, 0, 0, 0,
0, 1, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0), nrow = 5, byrow = TRUE)
f.con2 <- matrix(c(0, 1, 2, 4, 5, 7,
1, 1, 1, 0, 1, 0,
0, 0, 1, 1, 0, 0,
0, 1, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0), nrow = 5, byrow = TRUE)
f.con3 <- matrix(c(0, 1, 3, 4, 5, 7,
1, 1, 0, 0, 1, 0,
0, 0, 0, 1, 0, 0,
0, 1, 1, 0, 0, 0,
1, 0, 0, 0, 0, 0), nrow = 5, byrow = TRUE)
f.con4 <- matrix(c(0, 2, 3, 4, 5, 7,
1, 0, 0, 0, 1, 0,
0, 0, 0, 1, 0, 0,
0, 0, 1, 0, 0, 0,
1, 1, 0, 0, 0, 0), nrow = 5, byrow = TRUE)
f.con5 <- matrix(c(0, 1, 1, 3, 5, 7,
1, 1, 3, 2, 1, 0,
0, 0, 0, 0, 0, 0,
0, 1, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0), nrow = 5, byrow = TRUE) #missing 3/7
f.con6 <- matrix(c(0, 1, 3, 3, 5, 7,
1, 1, 0, 2, 1, 0,
0, 0, 0, 0, 0, 0,
0, 1, 1, 0, 0, 0,
1, 0, 0, 0, 0, 0), nrow = 5, byrow = TRUE) #missing 3/7
f.con7 <- matrix(c(0, 2, 3, 3, 5, 7,
1, 0, 0, 2, 1, 0,
0, 0, 0, 0, 0, 0,
0, 0, 1, 0, 0, 0,
1, 1, 0, 0, 0, 0), nrow = 5, byrow = TRUE) #missing 3/7
f.con8 <- matrix(c(0, 0, 2, 4, 5, 7,
1, 2, 1, 0, 1, 0,
0, 1, 1, 1, 0, 0,
0, 0, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0), nrow = 5, byrow = TRUE) #missing 4/7
f.con9 <- matrix(c(0, 2, 2, 3, 5, 7,
1, 0, 1, 2, 1, 0,
0, 0, 1, 0, 0, 0,
0, 0, 0, 0, 0, 0,
1, 1, 0, 0, 0, 0), nrow = 5, byrow = TRUE) #missing 4/7
f.con10 <- matrix(c(0, 2, 2, 4, 5, 7,
1, 0, 1, 0, 1, 0,
0, 0, 1, 1, 0, 0,
0, 0, 0, 0, 0, 0,
1, 1, 0, 0, 0, 0), nrow = 5, byrow = TRUE) #missing 4/7
f.con11 <- matrix(c(0, 0, 2, 3, 5, 7,
1, 2, 1, 2, 1, 0,
0, 1, 1, 0, 0, 0,
0, 0, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0), nrow = 5, byrow = TRUE) #missing 4/7
f.con12 <- matrix(c(0, 0, 1, 3, 5, 7,
1, 2, 3, 2, 1, 0,
0, 1, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0), nrow = 5, byrow = TRUE) #missing 4/7
cnmut = list(f.con1, f.con2, f.con3, f.con4, f.con5, f.con6, f.con7, f.con8, f.con9, f.con10, f.con11, f.con12)
} else if (nt == 7 & nb == 3) {
f.con1 <- matrix(c(0, 1, 2, 4, 5, 7,
0, 0, 1, 0, 1, 0,
1, 2, 1, 1, 0, 0,
1, 0, 0, 0, 0, 0), nrow = 4, byrow = TRUE)
f.con2 <- matrix(c(0, 1, 2, 3, 5, 7,
0, 0, 1, 2, 1, 0,
1, 2, 1, 0, 0, 0,
1, 0, 0, 0, 0, 0), nrow = 4, byrow = TRUE)
f.con3 <- matrix(c(0, 1, 2, 4, 5, 7,
0, 1, 1, 0, 1, 0,
1, 0, 1, 1, 0, 0,
1, 1, 0, 0, 0, 0), nrow = 4, byrow = TRUE)
f.con4 <- matrix(c(0, 1, 3, 4, 5, 7,
0, 1, 0, 0, 1, 0,
1, 0, 0, 1, 0, 0,
1, 1, 1, 0, 0, 0), nrow = 4, byrow = TRUE)
f.con5 <- matrix(c(0, 1, 2, 3, 5, 7,
0, 1, 1, 2, 1, 0,
1, 0, 1, 0, 0, 0,
1, 1, 0, 0, 0, 0), nrow = 4, byrow = TRUE)
f.con6 <- matrix(c(0, 0, 2, 4, 5, 7,
0, 2, 1, 0, 1, 0,
1, 1, 1, 1, 0, 0,
1, 0, 0, 0, 0, 0), nrow = 4, byrow = TRUE)
f.con7 <- matrix(c(0, 1, 3, 3, 5, 7,
0, 1, 0, 2, 1, 0,
1, 0, 0, 0, 0, 0,
1, 1, 1, 0, 0, 0), nrow = 4, byrow = TRUE)
f.con8 <- matrix(c(0, 0, 2, 3, 5, 7,
0, 2, 1, 2, 1, 0,
1, 1, 1, 0, 0, 0,
1, 0, 0, 0, 0, 0), nrow = 4, byrow = TRUE)
f.con9 <- matrix(c(0, 0, 1, 3, 5, 7,
0, 2, 3, 2, 1, 0,
1, 1, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0), nrow = 4, byrow = TRUE)
f.con10 <- matrix(c(0, 1, 1, 3, 5, 7,
0, 1, 3, 2, 1, 0,
1, 0, 0, 0, 0, 0,
1, 1, 0, 0, 0, 0), nrow = 4, byrow = TRUE)
cnmut = list(f.con1, f.con2, f.con3, f.con4, f.con5, f.con6, f.con7, f.con8, f.con9, f.con10)
} else if (nt == 7 & nb == 0) {
f.con1 <- matrix(c(0, 1, 2, 3, 4, 5, 7,
0, 0, 0, 0, 0, 1, 0,
0, 0, 0, 0, 1, 0, 0,
0, 0, 0, 1, 0, 0, 0,
0, 0, 1, 0, 0, 0, 0,
0, 1, 0, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0, 0), nrow = 7, byrow = TRUE)
f.con2 <- matrix(c(0, 1, 2, 3, 3, 5, 7,
0, 0, 0, 0, 2, 1, 0,
0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 1, 0, 0, 0,
0, 0, 1, 0, 0, 0, 0,
0, 1, 0, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0, 0), nrow = 7, byrow = TRUE) #missing 3/7
f.con3 <- matrix(c(0, 1, 2, 2, 3, 5, 7,
0, 0, 0, 1, 2, 1, 0,
0, 0, 0, 1, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0,
0, 0, 1, 0, 0, 0, 0,
0, 1, 0, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0, 0), nrow = 7, byrow = TRUE) #missing 4/7
f.con4 <- matrix(c(0, 1, 1, 2, 3, 5, 7,
0, 0, 1, 1, 2, 1, 0,
0, 0, 0, 1, 0, 0, 0,
0, 0, 1, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0,
0, 1, 0, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0, 0), nrow = 7, byrow = TRUE) #missing 4/7
f.con5 <- matrix(c(0, 0, 1, 2, 3, 5, 7,
0, 1, 1, 1, 2, 1, 0,
0, 0, 0, 1, 0, 0, 0,
0, 0, 1, 0, 0, 0, 0,
0, 1, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0, 0), nrow = 7, byrow = TRUE) #missing 5/7
f.con6 <- matrix(c(0, 1, 2, 2, 3, 5, 7,
0, 0, 0, 1, 2, 1, 0,
0, 0, 0, 1, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0,
0, 0, 1, 0, 0, 0, 0,
0, 1, 0, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0, 0), nrow = 7, byrow = TRUE) #missing 4/7
f.con7 <- matrix(c(0, 1, 2, 2, 4, 5, 7,
0, 0, 0, 1, 0, 1, 0,
0, 0, 0, 1, 1, 0, 0,
0, 0, 0, 0, 0, 0, 0,
0, 0, 1, 0, 0, 0, 0,
0, 1, 0, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0, 0), nrow = 7, byrow = TRUE) #missing 4/7
f.con8 <- matrix(c(0, 1, 1, 1, 3, 5, 7,
0, 0, 1, 3, 2, 1, 0,
0, 0, 0, 0, 0, 0, 0,
0, 0, 1, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0,
0, 1, 0, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0, 0), nrow = 7, byrow = TRUE) #missing 3/7 and 5/7
f.con9 <- matrix(c(0, 1, 1, 2, 4, 5, 7,
0, 0, 1, 1, 0, 1, 0,
0, 0, 0, 1, 1, 0, 0,
0, 0, 1, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0,
0, 1, 0, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0, 0), nrow = 7, byrow = TRUE) #missing 5/7
f.con10 <- matrix(c(0, 0, 1, 2, 4, 5, 7,
0, 1, 1, 1, 0, 1, 0,
0, 0, 0, 1, 1, 0, 0,
0, 0, 1, 0, 0, 0, 0,
0, 1, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0, 0), nrow = 7, byrow = TRUE) #missing 6/7
f.con11 <- matrix(c(0, 1, 1, 2, 4, 5, 7,
0, 0, 0, 1, 0, 1, 0,
0, 0, 2, 1, 1, 0, 0,
0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0,
0, 1, 0, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0, 0), nrow = 7, byrow = TRUE) #missing 4/7 and 5/7
f.con12 <- matrix(c(0, 0, 0, 1, 3, 5, 7,
0, 1, 2, 3, 2, 1, 0,
0, 0, 1, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0,
0, 1, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0, 0), nrow = 7, byrow = TRUE) #missing 4/7 and 6/7
f.con13 <- matrix(c(0, 1, 1, 3, 3, 5, 7,
0, 0, 1, 0, 2, 1, 0,
0, 0, 0, 0, 0, 0, 0,
0, 0, 1, 1, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0,
0, 1, 0, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0, 0), nrow = 7, byrow = TRUE) #missing 3/7 and 5/7
f.con14 <- matrix(c(0, 1, 1, 3, 4, 5, 7,
0, 0, 1, 0, 0, 1, 0,
0, 0, 0, 0, 1, 0, 0,
0, 0, 1, 1, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0,
0, 1, 0, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0, 0), nrow = 7, byrow = TRUE) #missing 5/7
f.con15 <- matrix(c(0, 1, 1, 2, 3, 5, 7,
0, 0, 0, 1, 2, 1, 0,
0, 0, 2, 1, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0,
0, 1, 0, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0, 0), nrow = 7, byrow = TRUE) #missing 4/7 and 5/7
f.con16 <- matrix(c(0, 0, 1, 1, 3, 5, 7,
0, 1, 1, 3, 2, 1, 0,
0, 0, 0, 0, 0, 0, 0,
0, 0, 1, 0, 0, 0, 0,
0, 1, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0, 0), nrow = 7, byrow = TRUE) #missing 3/7 and 6/7
f.con17 <- matrix(c(0, 0, 0, 2, 4, 5, 7,
0, 0, 2, 1, 0, 1, 0,
0, 1, 1, 1, 1, 0, 0,
0, 1, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0, 0), nrow = 7, byrow = TRUE) #missing 5/7 and 6/7
f.con18 <- matrix(c(0, 0, 0, 2, 4, 5, 7,
0, 1, 2, 1, 0, 1, 0,
0, 0, 1, 1, 1, 0, 0,
0, 0, 0, 0, 0, 0, 0,
0, 1, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0, 0), nrow = 7, byrow = TRUE) #missing 4/7 and 6/7
f.con19 <- matrix(c(0, 0, 0, 2, 3, 5, 7,
0, 0, 2, 1, 2, 1, 0,
0, 1, 1, 1, 0, 0, 0,
0, 1, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0, 0), nrow = 7, byrow = TRUE) #missing 5/7 and 6/7
f.con20 <- matrix(c(0, 0, 1, 3, 4, 5, 7,
0, 1, 1, 0, 0, 1, 0,
0, 0, 0, 0, 1, 0, 0,
0, 0, 1, 1, 0, 0, 0,
0, 1, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0, 0), nrow = 7, byrow = TRUE) #missing 6/7
f.con21 <- matrix(c(0, 0, 1, 3, 3, 5, 7,
0, 1, 1, 0, 2, 1, 0,
0, 0, 0, 0, 0, 0, 0,
0, 0, 1, 1, 0, 0, 0,
0, 1, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0, 0), nrow = 7, byrow = TRUE) #missing 3/7 and 6/7
f.con22 <- matrix(c(0, 0, 1, 2, 4, 5, 7,
0, 0, 1, 1, 0, 1, 0,
0, 1, 1, 1, 1, 0, 0,
0, 1, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0, 0), nrow = 7, byrow = TRUE) #missing 5/7 and 6/7
f.con23 <- matrix(c(0, 0, 1, 2, 4, 5, 7,
0, 0, 0, 1, 0, 1, 0,
0, 1, 2, 1, 1, 0, 0,
0, 1, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0, 0), nrow = 7, byrow = TRUE) #missing 5/7 and 6/7
f.con24 <- matrix(c(0, 0, 1, 2, 3, 5, 7,
0, 0, 0, 1, 2, 1, 0,
0, 1, 2, 1, 0, 0, 0,
0, 1, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0, 0), nrow = 7, byrow = TRUE) #missing 5/7 and 6/7
f.con25 <- matrix(c(0, 0, 0, 2, 3, 5, 7,
0, 1, 2, 1, 2, 1, 0,
0, 0, 1, 1, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0,
0, 1, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0, 0), nrow = 7, byrow = TRUE) #missing 4/7 and 6/7
f.con26 <- matrix(c(0, 0, 0, 1, 3, 5, 7,
0, 0, 2, 3, 2, 1, 0,
0, 1, 1, 0, 0, 0, 0,
0, 1, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0, 0), nrow = 7, byrow = TRUE) #missing 5/7 and 6/7
f.con27 <- matrix(c(0, 0, 1, 1, 3, 5, 7,
0, 0, 1, 3, 2, 1, 0,
0, 1, 0, 0, 0, 0, 0,
0, 1, 1, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0, 0), nrow = 7, byrow = TRUE) #missing 5/7 and 6/7
f.con28 <- matrix(c(0, 0, 1, 2, 3, 5, 7,
0, 0, 1, 1, 2, 1, 0,
0, 1, 0, 1, 0, 0, 0,
0, 1, 1, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0, 0), nrow = 7, byrow = TRUE) #missing 5/7 and 6/7
f.con29 <- matrix(c(0, 0, 2, 3, 4, 5, 7,
0, 1, 0, 0, 0, 1, 0,
0, 0, 0, 0, 1, 0, 0,
0, 0, 0, 1, 0, 0, 0,
0, 1, 1, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0, 0), nrow = 7, byrow = TRUE) #missing 6/7
f.con30 <- matrix(c(0, 0, 2, 3, 3, 5, 7,
0, 1, 0, 0, 2, 1, 0,
0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 1, 0, 0, 0,
0, 1, 1, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0, 0), nrow = 7, byrow = TRUE) #missing 3/7 and 6/7
f.con31 <- matrix(c(0, 0, 2, 2, 3, 5, 7,
0, 1, 0, 1, 2, 1, 0,
0, 0, 0, 1, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0,
0, 1, 1, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0, 0), nrow = 7, byrow = TRUE) #missing 4/7 and 6/7
f.con32 <- matrix(c(0, 0, 2, 2, 4, 5, 7,
0, 1, 0, 1, 1, 1, 0,
0, 0, 0, 1, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0,
0, 1, 1, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0, 0), nrow = 7, byrow = TRUE) #missing 4/7 and 6/7
f.con33 <- matrix(c(0, 0, 1, 2, 4, 5, 7,
0, 0, 0, 1, 0, 1, 0,
0, 1, 2, 1, 1, 0, 0,
0, 1, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0, 0), nrow = 7, byrow = TRUE) #missing 5/7 and 6/7
cnmut = list(f.con1, f.con2, f.con3, f.con4, f.con5, f.con6, f.con7, f.con8, f.con9, f.con10, f.con11, f.con12, f.con13, f.con14, f.con15, f.con16, f.con17, f.con18, f.con19, f.con20, f.con21, f.con22, f.con23, f.con24, f.con25, f.con26, f.con27, f.con28, f.con29, f.con30, f.con31, f.con32, f.con33)
} else if (nt == 8 & nb == 1) {
f.con1 <- matrix(c(1, 2, 3, 4, 5, 6, 8,
0, 0, 0, 0, 0, 1, 0,
0, 0, 0, 0, 1, 0, 0,
0, 0, 0, 1, 0, 0, 0,
0, 0, 1, 0, 0, 0, 0,
0, 1, 0, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0, 0), nrow = 7, byrow = TRUE)
f.con2 <- matrix(c(1, 2, 3, 4, 4, 6, 8,
0, 0, 0, 0, 2, 1, 0,
0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 1, 0, 0, 0,
0, 0, 1, 0, 0, 0, 0,
0, 1, 0, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0, 0), nrow = 7, byrow = TRUE) #missing 3/7
f.con3 <- matrix(c(1, 2, 3, 4, 4, 6, 8,
0, 0, 0, 1, 2, 1, 0,
0, 0, 0, 1, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0,
0, 0, 1, 0, 0, 0, 0,
0, 1, 0, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0, 0), nrow = 7, byrow = TRUE) #missing 4/7
f.con4 <- matrix(c(1, 2, 2, 3, 4, 6, 8,
0, 0, 1, 1, 2, 1, 0,
0, 0, 0, 1, 0, 0, 0,
0, 0, 1, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0,
0, 1, 0, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0, 0), nrow = 7, byrow = TRUE) #missing 4/7
f.con5 <- matrix(c(1, 1, 2, 3, 4, 6, 8,
0, 1, 1, 1, 2, 1, 0,
0, 0, 0, 1, 0, 0, 0,
0, 0, 1, 0, 0, 0, 0,
0, 1, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0, 0), nrow = 7, byrow = TRUE) #missing 5/7
f.con6 <- matrix(c(1, 2, 3, 3, 4, 6, 8,
0, 0, 0, 1, 2, 1, 0,
0, 0, 0, 1, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0,
0, 0, 1, 0, 0, 0, 0,
0, 1, 0, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0, 0), nrow = 7, byrow = TRUE) #missing 4/7
f.con7 <- matrix(c(1, 2, 3, 3, 5, 6, 8,
0, 0, 0, 1, 0, 1, 0,
0, 0, 0, 1, 1, 0, 0,
0, 0, 0, 0, 0, 0, 0,
0, 0, 1, 0, 0, 0, 0,
0, 1, 0, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0, 0), nrow = 7, byrow = TRUE) #missing 4/7
f.con8 <- matrix(c(1, 2, 2, 2, 4, 6, 8,
0, 0, 1, 3, 2, 1, 0,
0, 0, 0, 0, 0, 0, 0,
0, 0, 1, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0,
0, 1, 0, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0, 0), nrow = 7, byrow = TRUE) #missing 3/7 and 5/7
f.con9 <- matrix(c(1, 2, 2, 3, 5, 6, 8,
0, 0, 1, 1, 0, 1, 0,
0, 0, 0, 1, 1, 0, 0,
0, 0, 1, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0,
0, 1, 0, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0, 0), nrow = 7, byrow = TRUE) #missing 5/7
f.con10 <- matrix(c(1, 1, 2, 3, 5, 6, 8,
0, 1, 1, 1, 0, 1, 0,
0, 0, 0, 1, 1, 0, 0,
0, 0, 1, 0, 0, 0, 0,
0, 1, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0, 0), nrow = 7, byrow = TRUE) #missing 6/7
f.con11 <- matrix(c(1, 2, 2, 3, 5, 6, 8,
0, 0, 0, 1, 0, 1, 0,
0, 0, 2, 1, 1, 0, 0,
0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0,
0, 1, 0, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0, 0), nrow = 7, byrow = TRUE) #missing 4/7 and 5/7
f.con12 <- matrix(c(1, 1, 1, 2, 4, 6, 8,
0, 1, 2, 3, 2, 1, 0,
0, 0, 1, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0,
0, 1, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0, 0), nrow = 7, byrow = TRUE) #missing 4/7 and 6/7
f.con13 <- matrix(c(1, 2, 2, 4, 4, 6, 8,
0, 0, 1, 0, 2, 1, 0,
0, 0, 0, 0, 0, 0, 0,
0, 0, 1, 1, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0,
0, 1, 0, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0, 0), nrow = 7, byrow = TRUE) #missing 3/7 and 5/7
f.con14 <- matrix(c(1, 2, 2, 4, 5, 6, 8,
0, 0, 1, 0, 0, 1, 0,
0, 0, 0, 0, 1, 0, 0,
0, 0, 1, 1, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0,
0, 1, 0, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0, 0), nrow = 7, byrow = TRUE) #missing 5/7
f.con15 <- matrix(c(1, 2, 2, 3, 4, 6, 8,
0, 0, 0, 1, 2, 1, 0,
0, 0, 2, 1, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0,
0, 1, 0, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0, 0), nrow = 7, byrow = TRUE) #missing 4/7 and 5/7
f.con16 <- matrix(c(1, 1, 2, 2, 4, 6, 8,
0, 1, 1, 3, 2, 1, 0,
0, 0, 0, 0, 0, 0, 0,
0, 0, 1, 0, 0, 0, 0,
0, 1, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0, 0), nrow = 7, byrow = TRUE) #missing 3/7 and 6/7
f.con17 <- matrix(c(1, 1, 1, 3, 5, 6, 8,
0, 0, 2, 1, 0, 1, 0,
0, 1, 1, 1, 1, 0, 0,
0, 1, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0, 0), nrow = 7, byrow = TRUE) #missing 5/7 and 6/7
f.con18 <- matrix(c(1, 1, 1, 3, 5, 6, 8,
0, 1, 2, 1, 0, 1, 0,
0, 0, 1, 1, 1, 0, 0,
0, 0, 0, 0, 0, 0, 0,
0, 1, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0, 0), nrow = 7, byrow = TRUE) #missing 4/7 and 6/7
f.con19 <- matrix(c(1, 1, 1, 3, 4, 6, 8,
0, 0, 2, 1, 2, 1, 0,
0, 1, 1, 1, 0, 0, 0,
0, 1, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0, 0), nrow = 7, byrow = TRUE) #missing 5/7 and 6/7
f.con20 <- matrix(c(1, 1, 2, 4, 5, 6, 8,
0, 1, 1, 0, 0, 1, 0,
0, 0, 0, 0, 1, 0, 0,
0, 0, 1, 1, 0, 0, 0,
0, 1, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0, 0), nrow = 7, byrow = TRUE) #missing 6/7
f.con21 <- matrix(c(1, 1, 2, 4, 4, 6, 8,
0, 1, 1, 0, 2, 1, 0,
0, 0, 0, 0, 0, 0, 0,
0, 0, 1, 1, 0, 0, 0,
0, 1, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0, 0), nrow = 7, byrow = TRUE) #missing 3/7 and 6/7
f.con22 <- matrix(c(1, 1, 2, 3, 5, 6, 8,
0, 0, 1, 1, 0, 1, 0,
0, 1, 1, 1, 1, 0, 0,
0, 1, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0, 0), nrow = 7, byrow = TRUE) #missing 5/7 and 6/7
f.con23 <- matrix(c(1, 1, 2, 3, 5, 6, 8,
0, 0, 0, 1, 0, 1, 0,
0, 1, 2, 1, 1, 0, 0,
0, 1, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0, 0), nrow = 7, byrow = TRUE) #missing 5/7 and 6/7
f.con24 <- matrix(c(1, 1, 2, 3, 4, 6, 8,
0, 0, 0, 1, 2, 1, 0,
0, 1, 2, 1, 0, 0, 0,
0, 1, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0, 0), nrow = 7, byrow = TRUE) #missing 5/7 and 6/7
f.con25 <- matrix(c(1, 1, 1, 3, 4, 6, 8,
0, 1, 2, 1, 2, 1, 0,
0, 0, 1, 1, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0,
0, 1, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0, 0), nrow = 7, byrow = TRUE) #missing 4/7 and 6/7
f.con26 <- matrix(c(1, 1, 1, 2, 4, 6, 8,
0, 0, 2, 3, 2, 1, 0,
0, 1, 1, 0, 0, 0, 0,
0, 1, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0, 0), nrow = 7, byrow = TRUE) #missing 5/7 and 6/7
f.con27 <- matrix(c(1, 1, 2, 2, 4, 6, 8,
0, 0, 1, 3, 2, 1, 0,
0, 1, 0, 0, 0, 0, 0,
0, 1, 1, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0, 0), nrow = 7, byrow = TRUE) #missing 5/7 and 6/7
f.con28 <- matrix(c(1, 1, 2, 3, 4, 6, 8,
0, 0, 1, 1, 2, 1, 0,
0, 1, 0, 1, 0, 0, 0,
0, 1, 1, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0, 0), nrow = 7, byrow = TRUE) #missing 5/7 and 6/7
f.con29 <- matrix(c(1, 1, 3, 4, 5, 6, 8,
0, 1, 0, 0, 0, 1, 0,
0, 0, 0, 0, 1, 0, 0,
0, 0, 0, 1, 0, 0, 0,
0, 1, 1, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0, 0), nrow = 7, byrow = TRUE) #missing 6/7
f.con30 <- matrix(c(1, 1, 3, 4, 4, 6, 8,
0, 1, 0, 0, 2, 1, 0,
0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 1, 0, 0, 0,
0, 1, 1, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0, 0), nrow = 7, byrow = TRUE) #missing 3/7 and 6/7
f.con31 <- matrix(c(1, 1, 3, 3, 4, 6, 8,
0, 1, 0, 1, 2, 1, 0,
0, 0, 0, 1, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0,
0, 1, 1, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0, 0), nrow = 7, byrow = TRUE) #missing 4/7 and 6/7
f.con32 <- matrix(c(1, 1, 3, 3, 5, 6, 8,
0, 1, 0, 1, 1, 1, 0,
0, 0, 0, 1, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0,
0, 1, 1, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0, 0), nrow = 7, byrow = TRUE) #missing 4/7 and 6/7
f.con33 <- matrix(c(1, 1, 2, 3, 5, 6, 8,
0, 0, 0, 1, 0, 1, 0,
0, 1, 2, 1, 1, 0, 0,
0, 1, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0, 0), nrow = 7, byrow = TRUE) #missing 5/7 and 6/7
cnmut = list(f.con1, f.con2, f.con3, f.con4, f.con5, f.con6, f.con7, f.con8, f.con9, f.con10, f.con11, f.con12, f.con13, f.con14, f.con15, f.con16, f.con17, f.con18, f.con19, f.con20, f.con21, f.con22, f.con23, f.con24, f.con25, f.con26, f.con27, f.con28, f.con29, f.con30, f.con31, f.con32, f.con33)
}
return(cnmut)
}
R Package Documentation
Browse R Packages
We want your feedback!
Note that we can't provide technical support on individual packages. You should contact the package authors for that.
Embedding an R snippet on your website
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.