Mercurial > repos > greg > multilocus_genotype

comparison multilocus_genotype.R @ 14:96ee9122823e draft

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Uploaded
author greg
date Sun, 09 Dec 2018 09:44:09 -0500
parents ea2914ddea50
children 62ad61eac1ff
comparison
equal deleted inserted replaced
13:6b8061a95c04 14:96ee9122823e
299 labels <- paste(c("missing data", "mixed", "reference", "alternative"), " (", round(spy, 1), "%)", sep=""); 299 labels <- paste(c("missing data", "mixed", "reference", "alternative"), " (", round(spy, 1), "%)", sep="");
300 col <- c("GREY", "#006DDB", "#24FF24", "#920000"); 300 col <- c("GREY", "#006DDB", "#24FF24", "#920000");
301 main <- "Average breakdown of SNP assignments across all samples"; 301 main <- "Average breakdown of SNP assignments across all samples";
302 pie(spy, labels=labels, radius=0.60, col=col, main=main, cex.main=.75); 302 pie(spy, labels=labels, radius=0.60, col=col, main=main, cex.main=.75);
303 par(mfrow=c(3, 2)); 303 par(mfrow=c(3, 2));
304 col <- c("GREY", "#006DDB", "#24FF24", "#920000");
304 for (i in 1:96) { 305 for (i in 1:96) {
305 labels <- paste(labels, " (", round(tdt1_matrix[,i], 1), "%)", sep=""); 306 tmp_labels <- paste(labels, " (", round(tdt1_matrix[,i], 1), "%)", sep="");
306 col <- c("GREY", "#006DDB", "#24FF24", "#920000");
307 main <- paste("Breakdown of SNP assignments for", tdt1[1, i]); 307 main <- paste("Breakdown of SNP assignments for", tdt1[1, i]);
308 pie(tdt1_matrix[,i], labels=labels, radius=0.90, col=col, main=main, cex.main=.85, cex=0.75); 308 pie(tdt1_matrix[,i], labels=tmp_labels, radius=0.90, col=col, main=main, cex.main=.85, cex=0.75);
309 } 309 }
310 dev.off() 310 dev.off()
311 311