BatchProc.R
In permRand: Permutation Randomization

## ----include = FALSE----------------------------------------------------------
knitr::opts_chunk$set(
  collapse = TRUE,
  comment = "#>"
)

## ----libraries, message=FALSE-------------------------------------------------

library(dplyr)
library(magrittr)
library(permRand)

## ----testData-----------------------------------------------------------------
data("testR",package="permRand")

serumIDs <- testR[[1]]
serumLoc <- testR[[2]]
emptyQC <- testR[[3]]
motherQC <- testR[[4]]

## ----qualData-----------------------------------------------------------------
tests1 <- uniqueID(serumIDs,"serumID")
tests2 <- uniqueID(serumLoc,"serumID")
tests3 <- uniqueID(emptyQC,"serumID")
tests4 <- uniqueID(motherQC,"serumID")

## ----format-------------------------------------------------------------------
QCMaster <- randTest(dataMom=motherQC,dataChild=emptyQC,maxAliq=4, nEvent=c(28,27,28,30))

serumMaster <- formatRand(QCdata=QCMaster,serumIDR=serumIDs,serumPack=serumLoc)

## ----rand---------------------------------------------------------------------

serumMaster3 <- serumMaster %>% filter(event == 3)
serumRand <- allRand(dataR=serumMaster3,batchTot=c(40,44), numQC=2,
                     withinN=2,numMatch=2,chkRep=1)

## ----test1ar------------------------------------------------------------------
test <- testCCAcross(dataS=serumRand)

## ----test2ar------------------------------------------------------------------
test <- testQCmatch(dataS=serumRand,numQCs=4,numMatch=2)

## ----test3ar------------------------------------------------------------------
test <- uniqueID(serumRand,"serumID")

## ----test4ar------------------------------------------------------------------
test <- testPair(dataS=serumRand)

## ----test5ar------------------------------------------------------------------
test <- orderCases(dataI=serumRand,betW=4)

## ----test6ar------------------------------------------------------------------
test <- batchCount(dataS=serumRand,batchSizeT=84)

## ----test7ar------------------------------------------------------------------
test <- countQC(dataS=serumRand,QCN=4)

## ----packingar----------------------------------------------------------------
unBlind <- outputLab(dataOut=serumRand,blind=0,origP=1,maxRows=9,maxCols=9,newPack=1)
blind <- outputLab(dataOut=serumRand,blind=1,origP=0,maxRows=9,maxCols=9,newPack=1)

## ----switchR------------------------------------------------------------------
serumSwitch <- switchR(dataIn=serumRand,numqc=2,numqcM=2,batchS=44)

## ----checksr------------------------------------------------------------------
test1 <- testCCAcross(serumSwitch)
test2 <- testQCmatch(serumSwitch,2,2)
test3 <- uniqueID(serumSwitch,"serumID")
test4 <- testPair(serumSwitch)
test5 <- orderCases(serumSwitch,4)
test6 <- batchCount(serumSwitch,44)
test7 <- countQC(serumSwitch,2)

## ----packingarS---------------------------------------------------------------
serumSwitchP <- unBlind %>% 
                select(serumID,rack,row,col) %>%
                merge(.,serumSwitch,by='serumID')
unBlindSw <- outputLab(serumSwitchP,blind=0,origP=.,maxRows=.,maxCols=.,newPack=0)
blindSw <- outputLab(serumSwitchP,blind=1,origP=.,maxRows=.,maxCols=.,newPack=0)