aneuploid.info <- matrix(4, nrow=tetra.solcap$n.ind, ncol = 12)
set.seed(8080)
aneuploid.info[sample(1:length(aneuploid.info), round((4*length(aneuploid.info))/100),0)] <- 3
aneuploid.info[sample(1:length(aneuploid.info), round((4*length(aneuploid.info))/100),0)] <- 5
colnames(aneuploid.info) <- paste0(1:12)
aneuploid.info <- cbind(inds = tetra.solcap$ind.names, aneuploid.info)
filt.dat <- filter_aneuploid(input.data = tetra.solcap,
aneuploid.info = aneuploid.info, ploidy = 4)
Run the code above in your browser using DataLab