given a plating matrix, unwrap it, then perform association testing
p1 <- read.csv('platesample.csv',header=F)
p1_unwrap <- as.data.frame(unlist(p1),stringsAsFactors=F)
p1_unwrap$row <- rep(paste(c(1:8),sep=""),12)
p1_unwrap$col <- sort(rep(c(1:12),8))
p1_unwrap$case <- substr(p1_unwrap[,1],1,4)
p1_unwrap$visit <- substr(p1_unwrap[,1],5,9)
p1_unwrap$sample <- substr(p1_unwrap[,1],10,14)
p1_unwrap2 <- p1_unwrap[-which(p1_unwrap$case == "CEPH"),]
p1_unwrap2[which(p1_unwrap2$case == "Crtl"),]$case <- "Ctrl"
summary(glm(as.factor(case) ~ as.factor(row) ,data=p1_unwrap2),family='binomial')
table(p1_unwrap2$case,p1_unwrap2$row)
table(p1_unwrap2$case,p1_unwrap2$col)
table(p1_unwrap2$visit,p1_unwrap2$row)
table(p1_unwrap2$visit,p1_unwrap2$col)
Case Prep 65 CEPHDNA Crtl Post 31 Case Prep 29 Ctrl Prep 42 Case Post 64 Case Post 83 Ctrl Prep 75 Case Prep 85 Ctrl Post 81 Ctrl Post 116 Case Prep 107
Case Post 65 Ctrl Prep 179 Ctrl Prep 31 Case Post 29 Ctrl Post 42 Case Prep 64 Case Prep 83 Ctrl Post 75 Case Post 85 Ctrl Prep 81 Ctrl Prep 116 Case Post 107
Ctrl Post 28 Case Prep 72 Ctrl Prep 37 Case Post 93 Case Post 48 Ctrl Prep 56 Case Prep 49 Ctrl Post 66 Ctrl Post 68 Case Prep 95 Case Prep 104 Ctrl Post 136
Ctrl Prep 28 Case Post 72 Ctrl Post 37 Case Prep 93 Case Prep 48 Ctrl Post 56 Case Post 49 Ctrl Prep 66 Ctrl Prep 68 Case Post 95 Case Post 104 Ctrl Prep 136
Ctrl Prep 185 Case Post 90 Case Post 74 Ctrl Prep 39 Case Prep 77 Ctrl Post 60 Ctrl Post 58 Case Prep 87 Ctrl Prep 1 Case Post 97 Case Post 102 Ctrl Prep 73
Ctrl Post 185 Case Prep 90 Case Prep 74 Ctrl Post 39 Case Post 77 Ctrl Prep 60 Ctrl Prep 58 Case Post 87 Ctrl Post 1 Case Prep 97 Case Prep 102 Ctrl Post 73
Case Post 69 Ctrl Prep 35 Case Prep 215 Ctrl Post 41 Ctrl Post 54 Case Prep 78 Ctrl Prep 53 Case Post 59 Case Post 113 Ctrl Prep 44 Ctrl Prep 70 Case Post 119
Case Prep 69 Ctrl Post 35 Case Post 215 Ctrl Prep 41 Ctrl Prep 54 Case Post 78 Ctrl Post 53 Case Prep 59 Case Prep 113 Ctrl Post 44 Ctrl Post 70 Case Prep 119