Mercurial > repos > rnateam > methylkit

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