R/generateResults.R

In dschlauch/stego: Similarity Test for Estimation of Genetic Outliers

generateToyResults <- function(){
    data(toyGenotypes)
    toyResults$toyAllTogether <-  run_stego(toyGenotypes, simFun=cor)
    toyResults$toyUnphasedAllTogether <-   run_stego(toyGenotypes[,rep(c(T,F),100),,with=F]+toyGenotypes[,rep(c(F,T),100),with=F], groups="all.together", phased=F, simFun=cov)
    
    labels <- paste("Group",c(LETTERS[rep(1:5,20)]))
    toyResults$toyEachSeparately <- run_stego(toyGenotypes, groups="each.separately", labels=labels, simFun=cor)
    
    labels <- paste("Group",c(LETTERS[rep(1:5,10)],LETTERS[rep(6:10,10)]))
    super <- c(rep("Super A",50), rep("Super B",50))
    toyResults$toyPairwiseWithinSuperpop <- run_stego(toyGenotypes, groups="pairwise.within.superpop", labels=labels, super=super, simFun=cor)

    save(toyResults, file="./data/toyResults.rda")
}