Skip to content

Instantly share code, notes, and snippets.

@knausb
Last active May 19, 2017 22:50
Show Gist options
  • Save knausb/01a3d8efb21e0a4dbac0735270d147af to your computer and use it in GitHub Desktop.
Save knausb/01a3d8efb21e0a4dbac0735270d147af to your computer and use it in GitHub Desktop.
Test for Population Genomics Workshop dependencies
myFile <- "apstest.txt"
writeLines(capture.output(sessionInfo()), myFile)
writeLines(capture.output(sessionInfo()))
cat("\n", file = myFile, append = TRUE)
cat("\n", file = myFile, append = TRUE)
cat("\n")
cat("\n")
myString <- "Testing to see if you are ready for the Population Genomics in R Workshop!"
cat(myString)
cat("\n")
cat("\n")
cat(myString, file = myFile, append = TRUE)
cat("\n", file = myFile, append = TRUE)
cat("\n", file = myFile, append = TRUE)
pkgTest <- function(x){
if( suppressPackageStartupMessages( require(x, quietly = TRUE, character.only = TRUE)) ){
myString <- paste(x, "is installed.")
# print(myString)
# write(myString)
cat(myString)
cat("\n")
cat(myString, file = myFile, append = TRUE)
cat("\n", file = myFile, append = TRUE)
} else {
myString <- paste(x, "is not installed.")
cat(myString)
cat("\n")
cat(myString, file = myFile, append = TRUE)
cat("\n", file = myFile, append = TRUE)
}
flush.console()
}
myPackages <- c("adegenet", "ape", "cowplot", "dplyr", "ggplot2", "hierfstat",
"igraph", "knitr", "lattice", "magrittr", "mmod", "pegas", "pinfsc50",
"poppr", "RColorBrewer", "reshape2", "treemap", "vcfR", "vegan",
"viridisLite")
for(i in 1:length(myPackages)){
pkgTest(myPackages[i])
}
cat("\n", file = myFile, append = TRUE)
myString <- "Testing to see if you can load your data."
cat("\n")
#cat("\n")
cat(myString)
cat("\n")
cat(myString, file = myFile, append = TRUE)
cat("\n", file = myFile, append = TRUE)
#cat("\n", file = myFile, append = TRUE)
dataTest <- function(x){
vcf <- read.vcfR(x, verbose = FALSE)
myShow <- capture.output(vcf)
cat("\n")
cat("\n", file = myFile, append = TRUE)
for(i in 1:length(myShow)){
cat(myShow[i])
cat("\n")
cat(myShow[i], file = myFile, append = TRUE)
cat("\n", file = myFile, append = TRUE)
}
}
myData <- c("pinfsc50_filtered.vcf.gz", "prubi_gbs.vcf.gz")
for(i in 1:length(myData)){
dataTest(myData[i])
}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment