Mercurial > repos > bgruening > methylkit

annotate test-data/test.r @ 0:872907d4a311 draft default tip

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/methylkit commit 6054d0d3c7038ffc0303169a70bac4ff37040b98
author bgruening
date Tue, 20 Dec 2016 11:26:47 -0500
parents
children
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
0
872907d4a311 planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/methylkit commit 6054d0d3c7038ffc0303169a70bac4ff37040b98
bgruening
parents:
diff changeset
1 library("methylKit")
872907d4a311 planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/methylkit commit 6054d0d3c7038ffc0303169a70bac4ff37040b98
bgruening
parents:
diff changeset
2
872907d4a311 planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/methylkit commit 6054d0d3c7038ffc0303169a70bac4ff37040b98
bgruening
parents:
diff changeset
3 file.list = list("input_test1.myCpG.txt", "input_test2.myCpG.txt", "input_control1.myCpG.txt", "input_control2.myCpG.txt")
872907d4a311 planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/methylkit commit 6054d0d3c7038ffc0303169a70bac4ff37040b98
bgruening
parents:
diff changeset
4
872907d4a311 planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/methylkit commit 6054d0d3c7038ffc0303169a70bac4ff37040b98
bgruening
parents:
diff changeset
5 myobj=methRead( file.list,
872907d4a311 planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/methylkit commit 6054d0d3c7038ffc0303169a70bac4ff37040b98
bgruening
parents:
diff changeset
6 sample.id=list("test 1","test 2","control 1","control 2"),assembly="hg18",pipeline="amp",treatment=c(1,1,0,0))
872907d4a311 planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/methylkit commit 6054d0d3c7038ffc0303169a70bac4ff37040b98
bgruening
parents:
diff changeset
7
872907d4a311 planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/methylkit commit 6054d0d3c7038ffc0303169a70bac4ff37040b98
bgruening
parents:
diff changeset
8 pdf('output_statistics.pdf')
872907d4a311 planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/methylkit commit 6054d0d3c7038ffc0303169a70bac4ff37040b98
bgruening
parents:
diff changeset
9 for (obj in myobj){
872907d4a311 planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/methylkit commit 6054d0d3c7038ffc0303169a70bac4ff37040b98
bgruening
parents:
diff changeset
10 getMethylationStats(obj,plot=TRUE,both.strands=FALSE)
872907d4a311 planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/methylkit commit 6054d0d3c7038ffc0303169a70bac4ff37040b98
bgruening
parents:
diff changeset
11 getCoverageStats(obj,plot=TRUE,both.strands=FALSE)
872907d4a311 planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/methylkit commit 6054d0d3c7038ffc0303169a70bac4ff37040b98
bgruening
parents:
diff changeset
12 }
872907d4a311 planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/methylkit commit 6054d0d3c7038ffc0303169a70bac4ff37040b98
bgruening
parents:
diff changeset
13 devname = dev.off()
872907d4a311 planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/methylkit commit 6054d0d3c7038ffc0303169a70bac4ff37040b98
bgruening
parents:
diff changeset
14
872907d4a311 planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/methylkit commit 6054d0d3c7038ffc0303169a70bac4ff37040b98
bgruening
parents:
diff changeset
15 # unite function
872907d4a311 planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/methylkit commit 6054d0d3c7038ffc0303169a70bac4ff37040b98
bgruening
parents:
diff changeset
16 methidh = unite(myobj)
872907d4a311 planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/methylkit commit 6054d0d3c7038ffc0303169a70bac4ff37040b98
bgruening
parents:
diff changeset
17
872907d4a311 planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/methylkit commit 6054d0d3c7038ffc0303169a70bac4ff37040b98
bgruening
parents:
diff changeset
18 pdf("output_correlation.pdf")
872907d4a311 planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/methylkit commit 6054d0d3c7038ffc0303169a70bac4ff37040b98
bgruening
parents:
diff changeset
19 getCorrelation(object = methidh, plot=TRUE, method = "pearson")
872907d4a311 planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/methylkit commit 6054d0d3c7038ffc0303169a70bac4ff37040b98
bgruening
parents:
diff changeset
20 devname = dev.off()
872907d4a311 planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/methylkit commit 6054d0d3c7038ffc0303169a70bac4ff37040b98
bgruening
parents:
diff changeset
21
872907d4a311 planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/methylkit commit 6054d0d3c7038ffc0303169a70bac4ff37040b98
bgruening
parents:
diff changeset
22 # differential methylation
872907d4a311 planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/methylkit commit 6054d0d3c7038ffc0303169a70bac4ff37040b98
bgruening
parents:
diff changeset
23 myDiff = calculateDiffMeth(methidh, overdispersion="none",
872907d4a311 planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/methylkit commit 6054d0d3c7038ffc0303169a70bac4ff37040b98
bgruening
parents:
diff changeset
24 adjust="SLIM", effect="wmean", test="Chisq",
872907d4a311 planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/methylkit commit 6054d0d3c7038ffc0303169a70bac4ff37040b98
bgruening
parents:
diff changeset
25 slim=FALSE, weighted.mean=FALSE)
872907d4a311 planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/methylkit commit 6054d0d3c7038ffc0303169a70bac4ff37040b98
bgruening
parents:
diff changeset
26
872907d4a311 planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/methylkit commit 6054d0d3c7038ffc0303169a70bac4ff37040b98
bgruening
parents:
diff changeset
27 bedgraph(myDiff, file.name="output_myDiff.bedgraph", col.name="meth.diff",
872907d4a311 planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/methylkit commit 6054d0d3c7038ffc0303169a70bac4ff37040b98
bgruening
parents:
diff changeset
28 unmeth=FALSE, log.transform=FALSE, negative=FALSE, add.on="")
872907d4a311 planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/methylkit commit 6054d0d3c7038ffc0303169a70bac4ff37040b98
bgruening
parents:
diff changeset
29
872907d4a311 planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/methylkit commit 6054d0d3c7038ffc0303169a70bac4ff37040b98
bgruening
parents:
diff changeset
30 MethPerChr = diffMethPerChr(myDiff, plot=FALSE,
872907d4a311 planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/methylkit commit 6054d0d3c7038ffc0303169a70bac4ff37040b98
bgruening
parents:
diff changeset
31 qvalue.cutoff=0.01,
872907d4a311 planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/methylkit commit 6054d0d3c7038ffc0303169a70bac4ff37040b98
bgruening
parents:
diff changeset
32 meth.cutoff=25)
872907d4a311 planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/methylkit commit 6054d0d3c7038ffc0303169a70bac4ff37040b98
bgruening
parents:
diff changeset
33 write.table(MethPerChr, sep="\t", row.names=FALSE, quote=FALSE, file="output_MethPerChr.tsv")
872907d4a311 planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/methylkit commit 6054d0d3c7038ffc0303169a70bac4ff37040b98
bgruening
parents:
diff changeset
34
872907d4a311 planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/methylkit commit 6054d0d3c7038ffc0303169a70bac4ff37040b98
bgruening
parents:
diff changeset
35 MethylDiff = getMethylDiff(myDiff, difference=25,
872907d4a311 planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/methylkit commit 6054d0d3c7038ffc0303169a70bac4ff37040b98
bgruening
parents:
diff changeset
36 qvalue=0.01, type="all")
872907d4a311 planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/methylkit commit 6054d0d3c7038ffc0303169a70bac4ff37040b98
bgruening
parents:
diff changeset
37 bedgraph(MethylDiff, file.name="output_MethylDiff.bedgraph", col.name="meth.diff",
872907d4a311 planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/methylkit commit 6054d0d3c7038ffc0303169a70bac4ff37040b98
bgruening
parents:
diff changeset
38 unmeth=FALSE,log.transform=FALSE,negative=FALSE,add.on="")
872907d4a311 planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/methylkit commit 6054d0d3c7038ffc0303169a70bac4ff37040b98
bgruening
parents:
diff changeset
39
872907d4a311 planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/methylkit commit 6054d0d3c7038ffc0303169a70bac4ff37040b98
bgruening
parents:
diff changeset
40 pdf( "output_clustering.pdf" )
872907d4a311 planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/methylkit commit 6054d0d3c7038ffc0303169a70bac4ff37040b98
bgruening
parents:
diff changeset
41 methClust = clusterSamples(methidh, dist="correlation", method="ward")
872907d4a311 planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/methylkit commit 6054d0d3c7038ffc0303169a70bac4ff37040b98
bgruening
parents:
diff changeset
42 devname = dev.off()
872907d4a311 planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/methylkit commit 6054d0d3c7038ffc0303169a70bac4ff37040b98
bgruening
parents:
diff changeset
43
872907d4a311 planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/methylkit commit 6054d0d3c7038ffc0303169a70bac4ff37040b98
bgruening
parents:
diff changeset
44 pdf( "output_PCA.pdf" )
872907d4a311 planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/methylkit commit 6054d0d3c7038ffc0303169a70bac4ff37040b98
bgruening
parents:
diff changeset
45 PCASamples(methidh)
872907d4a311 planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/methylkit commit 6054d0d3c7038ffc0303169a70bac4ff37040b98
bgruening
parents:
diff changeset
46 devname = dev.off()
872907d4a311 planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/methylkit commit 6054d0d3c7038ffc0303169a70bac4ff37040b98
bgruening
parents:
diff changeset
47
872907d4a311 planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/methylkit commit 6054d0d3c7038ffc0303169a70bac4ff37040b98
bgruening
parents:
diff changeset
48 ## methSeg works for methylRaw or methylDiff with resolution region,
872907d4a311 planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/methylkit commit 6054d0d3c7038ffc0303169a70bac4ff37040b98
bgruening
parents:
diff changeset
49 ## so methylBase has to be tiled before
872907d4a311 planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/methylkit commit 6054d0d3c7038ffc0303169a70bac4ff37040b98
bgruening
parents:
diff changeset
50 tileRaw = tileMethylCounts(myobj[[1]])
872907d4a311 planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/methylkit commit 6054d0d3c7038ffc0303169a70bac4ff37040b98
bgruening
parents:
diff changeset
51 tileBase = tileMethylCounts(methidh)
872907d4a311 planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/methylkit commit 6054d0d3c7038ffc0303169a70bac4ff37040b98
bgruening
parents:
diff changeset
52 tileDiff = calculateDiffMeth(tileBase)
872907d4a311 planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/methylkit commit 6054d0d3c7038ffc0303169a70bac4ff37040b98
bgruening
parents:
diff changeset
53
872907d4a311 planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/methylkit commit 6054d0d3c7038ffc0303169a70bac4ff37040b98
bgruening
parents:
diff changeset
54 ## methseg generates Granges
872907d4a311 planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/methylkit commit 6054d0d3c7038ffc0303169a70bac4ff37040b98
bgruening
parents:
diff changeset
55 segRaw = methSeg(tileRaw, diagnostic.plot = FALSE)
872907d4a311 planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/methylkit commit 6054d0d3c7038ffc0303169a70bac4ff37040b98
bgruening
parents:
diff changeset
56 segDiff = methSeg(tileDiff, diagnostic.plot = FALSE)
872907d4a311 planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/methylkit commit 6054d0d3c7038ffc0303169a70bac4ff37040b98
bgruening
parents:
diff changeset
57
872907d4a311 planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/methylkit commit 6054d0d3c7038ffc0303169a70bac4ff37040b98
bgruening
parents:
diff changeset
58 ## and can be exported as BED
872907d4a311 planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/methylkit commit 6054d0d3c7038ffc0303169a70bac4ff37040b98
bgruening
parents:
diff changeset
59 methSeg2bed(segments = segRaw, filename = "output_seg_raw.bed")
872907d4a311 planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/methylkit commit 6054d0d3c7038ffc0303169a70bac4ff37040b98
bgruening
parents:
diff changeset
60 methSeg2bed(segments = segDiff, filename = "output_seg_diff.bed")