nievergeltlab
10/31/2018 - 4:11 PM

Check plating

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