Nothing
## ----setup, include = FALSE---------------------------------------------------
knitr::opts_chunk$set(
collapse = TRUE,
comment = "#>"
)
## ---- fig.cap='PRAM workflow', out.width='237px', out.height='233px', echo=F----
knitr::include_graphics("workflow_noScreen.jpg")
## ----installFromGitHub, eval=FALSE--------------------------------------------
# devtools::install_github('pliu55/pram')
## ----installFromBioconductor, eval=FALSE--------------------------------------
# if (!requireNamespace("BiocManager", quietly = TRUE))
# install.packages("BiocManager")
# BiocManager::install("pram")
## ----runQuickPRAM, eval=FALSE-------------------------------------------------
# ##
# ## assuming the stringtie binary is in folder /usr/local/stringtie-1.3.3/
# ##
# pram::runPRAM( in_gtf, in_bamv, out_gtf, method='plst',
# stringtie='/usr/loca/stringtie-1.3.3/stringtie')
## ----accessExample------------------------------------------------------------
system.file('extdata/demo/in.gtf', package='pram')
## ----runPRAM, eval=FALSE------------------------------------------------------
# in_gtf = system.file('extdata/demo/in.gtf', package='pram')
#
# in_bamv = c(system.file('extdata/demo/SZP.bam', package='pram'),
# system.file('extdata/demo/TLC.bam', package='pram') )
#
# pred_out_gtf = tempfile(fileext='.gtf')
#
# ##
# ## assuming the stringtie binary is in folder /usr/local/stringtie-1.3.3/
# ##
# pram::runPRAM( in_gtf, in_bamv, pred_out_gtf, method='plst',
# stringtie='/usr/loca/stringtie-1.3.3/stringtie')
## ----echo=F-------------------------------------------------------------------
url_hg19_chromsize=
'http://hgdownload.cse.ucsc.edu/goldenpath/hg19/database/chromInfo.txt.gz'
## ----defIgRanges, warning=FALSE-----------------------------------------------
pram::defIgRanges(system.file('extdata/gtf/defIgRanges_in.gtf', package='pram'),
genome = 'hg38')
## ----prepIgBam----------------------------------------------------------------
finbam =system.file('extdata/bam/CMPRep2.sortedByCoord.raw.bam',
package='pram')
iggrs = GenomicRanges::GRanges('chr10:77236000-77247000:+')
foutbam = tempfile(fileext='.bam')
pram::prepIgBam(finbam, iggrs, foutbam)
## ----echo=F-------------------------------------------------------------------
url_cf_web = paste0('[Cufflinks, Cuffmerge]',
'(http://cole-trapnell-lab.github.io/cufflinks/)')
url_cf_lnx = paste0('[v2.2.1]',
'(http://cole-trapnell-lab.github.io/cufflinks/assets/',
'downloads/cufflinks-2.2.1.Linux_x86_64.tar.gz)')
url_cf_mac = paste0('[v2.1.1]',
'(http://cole-trapnell-lab.github.io/cufflinks/assets/',
'downloads/cufflinks-2.1.1.OSX_x86_64.tar.gz)')
cf_methods = '__plcf__, __cfmg__, __cftc__, and __cf__'
url_st_web = paste0('[StringTie, StringTie-merge]',
'(https://ccb.jhu.edu/software/stringtie/)')
url_st_lnx = paste0('[v1.3.3b]',
'(http://ccb.jhu.edu/software/stringtie/dl/',
'stringtie-1.3.3b.Linux_x86_64.tar.gz)')
url_st_mac = paste0('[v1.3.3b]',
'(http://ccb.jhu.edu/software/stringtie/dl/',
'stringtie-1.3.3b.OSX_x86_64.tar.gz)')
st_methods = '__plst__, __stmg__, and __st__'
url_tc_web = '[TACO](https://tacorna.github.io)'
url_tc_lnx = paste0('[v0.7.0]',
'(https://github.com/tacorna/taco/releases/download/',
'v0.7.0/taco-v0.7.0.Linux_x86_64.tar.gz)')
url_tc_mac = paste0('[v0.7.0]',
'(https://github.com/tacorna/taco/releases/download/',
'v0.7.0/taco-v0.7.0.OSX_x86_64.tar.gz)')
## ----buildModel, eval=FALSE---------------------------------------------------
# fbams = c( system.file('extdata/bam/CMPRep1.sortedByCoord.clean.bam',
# package='pram'),
# system.file('extdata/bam/CMPRep2.sortedByCoord.clean.bam',
# package='pram') )
#
# foutgtf = tempfile(fileext='.gtf')
#
# ##
# ## assuming the stringtie binary is in folder /usr/local/stringtie-1.3.3/
# ##
# pram::buildModel(fbams, foutgtf, method='plst',
# stringtie='/usr/loca/stringtie-1.3.3/stringtie')
## ----selModel-----------------------------------------------------------------
fin_gtf = system.file('extdata/gtf/selModel_in.gtf', package='pram')
fout_gtf = tempfile(fileext='.gtf')
pram::selModel(fin_gtf, fout_gtf)
## ----evalModel----------------------------------------------------------------
fmdl = system.file('extdata/benchmark/plcf.tsv', package='pram')
ftgt = system.file('extdata/benchmark/tgt.tsv', package='pram')
mdldt = data.table::fread(fmdl, header=TRUE, sep="\t")
tgtdt = data.table::fread(ftgt, header=TRUE, sep="\t")
pram::evalModel(mdldt, tgtdt)
## ----sessionInfo, echo=FALSE--------------------------------------------------
sessionInfo()
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.