tests/testthat/helper-00-freqs.R

In DNAtools: Tools for Analysing Forensic Genetic DNA Data

# From DNAtools01022, i.e. version 0.1-22 from before refactoring
if (FALSE) {
  if (!requireNamespace('DNAtools01022', quietly = TRUE)) {
    install.packages("DNAtools01022_0.1-22.tar.gz", repos = NULL)
  }
  
  library(DNAtools01022)
  
  set.seed(1)
  freqs <-  structure(replicate(10, {g = rgamma(n=10,scale=4,shape=3); g/sum(g)}, simplify=FALSE), .Names=paste("locus",1:10,sep="."))
  dput(freqs)
}


freqs <- structure(list(locus.1 = c(0.0474946168330123, 0.149039228357251, 
0.145256669862042, 0.094498400655289, 0.198907997603563, 0.112374968690082, 
0.103204475390989, 0.0602827479229705, 0.0412507951226619, 0.0476900995621393
), locus.2 = c(0.074146387385448, 0.0879018061731534, 0.0466247431670957, 
0.143534968752236, 0.127620629358631, 0.150666185736466, 0.140732142776328, 
0.0947701351097818, 0.0124442728814128, 0.121558728659448), locus.3 = c(0.177479220550563, 
0.0600659281787312, 0.0708557662895187, 0.0842339258845552, 0.0920591454047341, 
0.0620993146775552, 0.0313437854899729, 0.169743616701942, 0.0734468733928604, 
0.178672423429567), locus.4 = c(0.156425873241834, 0.0912264835619842, 
0.0745255553645018, 0.140400119644217, 0.0621098494751476, 0.0259389163884027, 
0.15685312211086, 0.0944298043729314, 0.119320678049069, 0.0787695977910527
), locus.5 = c(0.112274428730273, 0.0378042227780907, 0.10651892081973, 
0.130739060067879, 0.0817605799974613, 0.127404479579307, 0.0838283473463391, 
0.0908370172922858, 0.135730814319751, 0.0931021290688829), locus.6 = c(0.051866159519443, 
0.0267485368211831, 0.189830514072871, 0.0974354379929345, 0.25226766401127, 
0.117296575986563, 0.0533051308520478, 0.0820815868972172, 0.0322894566855689, 
0.0968789371609022), locus.7 = c(0.075755387188565, 0.0856165759021539, 
0.0540968703966548, 0.159405678012481, 0.126513276887885, 0.0707758329472124, 
0.115312830325823, 0.0998448640716341, 0.145916945816844, 0.0667617384507473
), locus.8 = c(0.102163683236432, 0.140173152152202, 0.0368963169328276, 
0.1315805688531, 0.105055927510634, 0.158821828794473, 0.0974876945093594, 
0.0250394065539911, 0.0809873607876016, 0.121794060669379), locus.9 = c(0.120772854816577, 
0.0441798872016384, 0.0393508153581742, 0.0320528442665698, 0.18866118822539, 
0.116815807036734, 0.128758914786706, 0.0976464038788485, 0.18220278429176, 
0.049558500137602), locus.10 = c(0.0968160027785965, 0.0969431498835093, 
0.0209293175086413, 0.0730796172945451, 0.116577979448978, 0.127429617413566, 
0.103780999758667, 0.0882562623297492, 0.193398734210837, 0.0827883193729112
)), .Names = c("locus.1", "locus.2", "locus.3", "locus.4", "locus.5", 
"locus.6", "locus.7", "locus.8", "locus.9", "locus.10"))