Commit 516a97e3 authored by Guillaume Devailly's avatar Guillaume Devailly
Browse files

Export of CpG density in bigwig format

parent 6ac8c19a
......@@ -18,6 +18,12 @@ cpgDens <- cpgDensityCalc(windows, organism=Sscrofa, w.function="linear", window
# windows_genome <- readRDS("/work2/genphyse/genepi/guillaume/rosepigs/windows_full_genome.RDS")
# load("/work2/genphyse/genepi/guillaume/rosepigs/cpgDens.RData")
cpgdens <- windows_genome
cpgdens$score <- cpgDens
seqlengths(cpgdens) <- seqlengths(Sscrofa)
export(cpgdens, "/work2/genphyse/genepi/guillaume/rosepigs/cpgDens.Sscrofa.11.ensembl.bigWig", format = "bigWig")
methyl1 <- readRDS("/work2/genphyse/genepi/guillaume/rosepigs/annotationBlockCountRDS/Tem-methyl_R1.rds")
methyl2 <- readRDS("/work2/genphyse/genepi/guillaume/rosepigs/annotationBlockCountRDS/Tem-methyl_R2.rds")
methyl3 <- readRDS("/work2/genphyse/genepi/guillaume/rosepigs/annotationBlockCountRDS/Tem-methyl_R3.rds")
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment