In stranda/unpakathon: Code Used and Data Needed to work with unpak data

Three levels of experimental organization

library(VCA)
library(unpakathon)
library(dplyr)

pl = phenolong%>%filter(meta.experiment=="1", variable=="fruitnum") %>%
    mutate(institution=gsub("-GC.","",facility)) %>%
    mutate(knockout=grepl("SALK",accession)) %>%
    mutate(knockout=ifelse(grepl("60000",accession),F,knockout))

fit1 = anovaMM(value~(institution)+(institution/facility), Data=pl)

print(fit1)

stranda/unpakathon documentation built on Nov. 9, 2021, 7:48 a.m.