project_rnsq: Create directory for RNA-seq pipeline

Description Usage Arguments Examples

View source: R/project_rnsq.R

Description

Create Project directories for RNA-seq pipeline using hisat2-stringtie. Then fastq files move to './project/fastq' directory, followning execution of 'rskoseq::rep_hisat2' or 'rskoseq::rep_rsem', for read mapping.

Usage

1
project_rnsq(prjd, alnd, algnr)

Arguments

prjd

character: path of a project directory name

alnd

character: path of a alignment directory name, the default is "alignment1"

algnr

character: select an alignment program from "hisat2"(default), "rsem", "others"

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
# # create new project
# prj <- "~/pub/sampledata/rnaseq/project1"
# project_rnsq(prjd = prj, alnd = "test.h2", algnr = "hisat2")
# project_rnsq(prjd = prj, alnd = "test.rsm", algnr = "rsem")
#
# # now there is a project directory, create another alignment directory.
# system(paste("tree", prj))
# # ./project1
# # ├── test.h2
# # │   ├── res_hisat2
# # │   │   └── failalign
# # │   └── res_stringtie
# # │       ├── ballgown
# # │       ├── gff
# # │       ├── mgff
# # │       └── tab
# # └── test.rsm
# #     ├── resdir
# #     ├── sortbam
# #     └── stats

shkonishi/rskoseq documentation built on Dec. 6, 2018, 6:18 p.m.