library(dplyr)
library(tidyr)
d = readr::read_delim("data-raw/CorrectedGenomeResearchData.txt",
delim = '\t',
col_names = TRUE)
Paschold2012 = d %>%
gather(tmp,count,-GeneID) %>%
mutate(tmp_ = gsub("\\.","_",tmp)) %>%
separate(tmp_, c("replicate","genotype","line")) %>%
mutate(genotype = plyr::revalue(genotype, c("B"="B73","M"="Mo17","BM"="B73xMo17","MB"="Mo17xB73")),
line = plyr::revalue(line, c("b"="B73","m"="Mo17")),
replicate = substr(replicate,4,4),
flow_cell = ceiling(as.numeric(replicate)/2)) %>%
select(GeneID, genotype, replicate, flow_cell, line, count) %>%
arrange(GeneID, genotype, replicate) %>%
spread(line,count)
devtools::use_data(Paschold2012, overwrite=TRUE, compress='xz')
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.