annotate rmarkdown_feature_counts.Rmd @ 13:eb51eca6cdf1 draft

update
author mingchen0919
date Fri, 29 Dec 2017 15:00:59 -0500
parents 238d57ce7d74
children bfa6440d01a3
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
0
8bcae2a95fe7 planemo upload
mingchen0919
parents:
diff changeset
1 ---
11
mingchen0919
parents: 10
diff changeset
2 title: 'Feature Counts'
0
8bcae2a95fe7 planemo upload
mingchen0919
parents:
diff changeset
3 output:
8bcae2a95fe7 planemo upload
mingchen0919
parents:
diff changeset
4 html_document:
8bcae2a95fe7 planemo upload
mingchen0919
parents:
diff changeset
5 number_sections: true
8bcae2a95fe7 planemo upload
mingchen0919
parents:
diff changeset
6 toc: true
8bcae2a95fe7 planemo upload
mingchen0919
parents:
diff changeset
7 theme: cosmo
8bcae2a95fe7 planemo upload
mingchen0919
parents:
diff changeset
8 highlight: tango
8bcae2a95fe7 planemo upload
mingchen0919
parents:
diff changeset
9 ---
8bcae2a95fe7 planemo upload
mingchen0919
parents:
diff changeset
10
8bcae2a95fe7 planemo upload
mingchen0919
parents:
diff changeset
11 ```{r setup, include=FALSE, warning=FALSE, message=FALSE}
8bcae2a95fe7 planemo upload
mingchen0919
parents:
diff changeset
12 knitr::opts_chunk$set(
7
mingchen0919
parents: 0
diff changeset
13 echo = opt$echo,
0
8bcae2a95fe7 planemo upload
mingchen0919
parents:
diff changeset
14 error = TRUE
8bcae2a95fe7 planemo upload
mingchen0919
parents:
diff changeset
15 )
8bcae2a95fe7 planemo upload
mingchen0919
parents:
diff changeset
16 ```
8bcae2a95fe7 planemo upload
mingchen0919
parents:
diff changeset
17
8bcae2a95fe7 planemo upload
mingchen0919
parents:
diff changeset
18
11
mingchen0919
parents: 10
diff changeset
19 # User input
0
8bcae2a95fe7 planemo upload
mingchen0919
parents:
diff changeset
20
11
mingchen0919
parents: 10
diff changeset
21 ```{r 'user input'}
7
mingchen0919
parents: 0
diff changeset
22 opt
0
8bcae2a95fe7 planemo upload
mingchen0919
parents:
diff changeset
23 ```
8bcae2a95fe7 planemo upload
mingchen0919
parents:
diff changeset
24
11
mingchen0919
parents: 10
diff changeset
25 # Calculate feature counts
mingchen0919
parents: 10
diff changeset
26
0
8bcae2a95fe7 planemo upload
mingchen0919
parents:
diff changeset
27 ```{r 'ste[ 2'}
10
1cff3a3e2b94 working version
mingchen0919
parents: 7
diff changeset
28 res = featureCounts(
11
mingchen0919
parents: 10
diff changeset
29 files = strsplit(opt$input_bam_paths, ',')[[1]],
10
1cff3a3e2b94 working version
mingchen0919
parents: 7
diff changeset
30 # annotation
1cff3a3e2b94 working version
mingchen0919
parents: 7
diff changeset
31 annot.inbuilt=opt$annot_inbuilt,
1cff3a3e2b94 working version
mingchen0919
parents: 7
diff changeset
32 annot.ext=opt$annot_ext,
1cff3a3e2b94 working version
mingchen0919
parents: 7
diff changeset
33 isGTFAnnotationFile=opt$isGTFAnnotationFile,
1cff3a3e2b94 working version
mingchen0919
parents: 7
diff changeset
34 GTF.featureType=opt$gtf_feature_type,
1cff3a3e2b94 working version
mingchen0919
parents: 7
diff changeset
35 GTF.attrType=opt$gtf_attr_type,
1cff3a3e2b94 working version
mingchen0919
parents: 7
diff changeset
36 chrAliases=opt$chr_aliases,
1cff3a3e2b94 working version
mingchen0919
parents: 7
diff changeset
37
1cff3a3e2b94 working version
mingchen0919
parents: 7
diff changeset
38 # level of summarization
1cff3a3e2b94 working version
mingchen0919
parents: 7
diff changeset
39 useMetaFeatures=opt$use_meta_features,
1cff3a3e2b94 working version
mingchen0919
parents: 7
diff changeset
40
1cff3a3e2b94 working version
mingchen0919
parents: 7
diff changeset
41 # overlap between reads and features
1cff3a3e2b94 working version
mingchen0919
parents: 7
diff changeset
42 allowMultiOverlap=opt$allow_multi_overlap,
1cff3a3e2b94 working version
mingchen0919
parents: 7
diff changeset
43 minOverlap=opt$min_overlap,
1cff3a3e2b94 working version
mingchen0919
parents: 7
diff changeset
44 largestOverlap=opt$largest_overlap,
1cff3a3e2b94 working version
mingchen0919
parents: 7
diff changeset
45 readExtension5=opt$read_extension_5,
1cff3a3e2b94 working version
mingchen0919
parents: 7
diff changeset
46 readExtension3=opt$read_extension_3,
1cff3a3e2b94 working version
mingchen0919
parents: 7
diff changeset
47 read2pos=opt$read_2_pos,
1cff3a3e2b94 working version
mingchen0919
parents: 7
diff changeset
48
1cff3a3e2b94 working version
mingchen0919
parents: 7
diff changeset
49 # multi-mapping reads
1cff3a3e2b94 working version
mingchen0919
parents: 7
diff changeset
50 countMultiMappingReads=opt$count_multi_mapping_reads,
1cff3a3e2b94 working version
mingchen0919
parents: 7
diff changeset
51 fraction=opt$fraction,
0
8bcae2a95fe7 planemo upload
mingchen0919
parents:
diff changeset
52
10
1cff3a3e2b94 working version
mingchen0919
parents: 7
diff changeset
53 # read filtering
1cff3a3e2b94 working version
mingchen0919
parents: 7
diff changeset
54 minMQS=opt$min_mqs,
1cff3a3e2b94 working version
mingchen0919
parents: 7
diff changeset
55 splitOnly=opt$split_only,
1cff3a3e2b94 working version
mingchen0919
parents: 7
diff changeset
56 nonSplitOnly=opt$non_split_only,
1cff3a3e2b94 working version
mingchen0919
parents: 7
diff changeset
57 primaryOnly=opt$primary_only,
1cff3a3e2b94 working version
mingchen0919
parents: 7
diff changeset
58 ignoreDup=opt$ignore_dup,
1cff3a3e2b94 working version
mingchen0919
parents: 7
diff changeset
59
1cff3a3e2b94 working version
mingchen0919
parents: 7
diff changeset
60 # strandness
1cff3a3e2b94 working version
mingchen0919
parents: 7
diff changeset
61 strandSpecific=opt$strand_specific,
1cff3a3e2b94 working version
mingchen0919
parents: 7
diff changeset
62
1cff3a3e2b94 working version
mingchen0919
parents: 7
diff changeset
63 # exon-exon junctions
1cff3a3e2b94 working version
mingchen0919
parents: 7
diff changeset
64 juncCounts=opt$junc_counts,
1cff3a3e2b94 working version
mingchen0919
parents: 7
diff changeset
65 genome=opt$genome,
1cff3a3e2b94 working version
mingchen0919
parents: 7
diff changeset
66
1cff3a3e2b94 working version
mingchen0919
parents: 7
diff changeset
67 # parameters specific to paired end reads
1cff3a3e2b94 working version
mingchen0919
parents: 7
diff changeset
68 isPairedEnd=opt$is_paired_end,
1cff3a3e2b94 working version
mingchen0919
parents: 7
diff changeset
69 requireBothEndsMapped=opt$require_both_ends_mapped,
1cff3a3e2b94 working version
mingchen0919
parents: 7
diff changeset
70 checkFragLength=opt$check_frag_length,
1cff3a3e2b94 working version
mingchen0919
parents: 7
diff changeset
71 minFragLength=opt$min_frag_length,
1cff3a3e2b94 working version
mingchen0919
parents: 7
diff changeset
72 maxFragLength=opt$max_frag_length,
1cff3a3e2b94 working version
mingchen0919
parents: 7
diff changeset
73 countChimericFragments=opt$count_chimeric_fragments,
1cff3a3e2b94 working version
mingchen0919
parents: 7
diff changeset
74 autosort=opt$auto_sort,
1cff3a3e2b94 working version
mingchen0919
parents: 7
diff changeset
75
1cff3a3e2b94 working version
mingchen0919
parents: 7
diff changeset
76 # miscellaneous
1cff3a3e2b94 working version
mingchen0919
parents: 7
diff changeset
77 nthreads=opt$n_threads,
1cff3a3e2b94 working version
mingchen0919
parents: 7
diff changeset
78 maxMOp=opt$max_mop,
1cff3a3e2b94 working version
mingchen0919
parents: 7
diff changeset
79 reportReads=opt$report_reads
1cff3a3e2b94 working version
mingchen0919
parents: 7
diff changeset
80 )
11
mingchen0919
parents: 10
diff changeset
81 ```
mingchen0919
parents: 10
diff changeset
82
mingchen0919
parents: 10
diff changeset
83 # Write counts into CSV file
mingchen0919
parents: 10
diff changeset
84
mingchen0919
parents: 10
diff changeset
85 ```{r}
mingchen0919
parents: 10
diff changeset
86 colnames(res$counts) = strsplit(opt$input_bam_names, ',')[[1]]
mingchen0919
parents: 10
diff changeset
87 # write count into csv file
13
mingchen0919
parents: 11
diff changeset
88 write.table(res$counts, file = 'feature_counts.txt')
mingchen0919
parents: 11
diff changeset
89 ```
mingchen0919
parents: 11
diff changeset
90
mingchen0919
parents: 11
diff changeset
91 Display the first 100 rows.
mingchen0919
parents: 11
diff changeset
92
mingchen0919
parents: 11
diff changeset
93 ```{r}
mingchen0919
parents: 11
diff changeset
94 datatable(head(res$counts, 100))
11
mingchen0919
parents: 10
diff changeset
95 ```
mingchen0919
parents: 10
diff changeset
96
mingchen0919
parents: 10
diff changeset
97 # Save results into RData file
mingchen0919
parents: 10
diff changeset
98
mingchen0919
parents: 10
diff changeset
99 ```{r}
mingchen0919
parents: 10
diff changeset
100 save(res, file = 'feature_counts.RData')
10
1cff3a3e2b94 working version
mingchen0919
parents: 7
diff changeset
101 str(res)
0
8bcae2a95fe7 planemo upload
mingchen0919
parents:
diff changeset
102 ```
8bcae2a95fe7 planemo upload
mingchen0919
parents:
diff changeset
103
8bcae2a95fe7 planemo upload
mingchen0919
parents:
diff changeset
104