This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
################################################## | |
## Create mappability file for reference genome ## | |
################################################## | |
## Thomas Girke | |
## Date: May 16, 2018 | |
## PeakSeq docs: | |
## https://github.com/gersteinlab/PeakSeq | |
## http://gensoft.pasteur.fr/docs/PeakSeq/1.1/PeakSeq.readme | |
## Note: mappability step should be skipped in newer version of PeakSeq |
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
############################################### | |
## Shiny App Example with Reactive URL Input ## | |
############################################### | |
## Author: Yuzhu Duan & Thomas Girke | |
## Last update: 06-Jul-17 | |
## Usage: save this file under name app.R to your shiny app directory, e.g. reactiveUrl | |
## Run in parent directory the following from within R: | |
## > library(shiny) | |
## > runApp(reactiveUrl) |
This page outlines how to maintain and sync R packages to both GitHub and Bioconductor's new git source control system.
Note: the master
branch has been renamed to devel
branch following the instructions here (Mar-2023).
!!! This Gist tutorial is deprecated. Its new version is available here !!!
This is some red text.
- Author: Thomas Girke
- Last update: 18-Nov-2020
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
################################################################### | |
## Alignment Stats with Support for FASTQ Files with >500M Reads ## | |
################################################################### | |
alignStats <- function(args) { | |
fqpaths <- infile1(args) | |
bampaths <- outpaths(args) | |
bamexists <- file.exists(bampaths) | |
fqpaths <- fqpaths[bamexists] | |
bampaths <- bampaths[bamexists] | |
## Obtain total read number from FASTQ files |
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
#################################################### | |
## Import of KEGG Compounds into SDFset Container ## | |
#################################################### | |
## Date: March 07, 2018 | |
## Motivation: response to request on Bioc support list: https://support.bioconductor.org/p/106712/ | |
################################ | |
## Download KEGG CMPs via URL ## | |
################################ |
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
############################################## | |
## Map Pfam Domains to Proteins with HMMER3 ## | |
############################################## | |
## Author: Thomas Girke | |
## Date: May 11, 2018 | |
## Utility: mapping of Pfam domains to protein sequences. | |
## The module load and Pfam database paths given below are specific to the HPCC/biocluster system. | |
## For details consult the man page for hmmscan from the command-line with 'hmmscan -h' |
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
############################### | |
## PWM matching in sequences ## | |
############################### | |
## Date: May 22, 2020 | |
## Starter code for ChIP-Seq1 project | |
library(MotifDb) | |
as.list(query(MotifDb, 'athaliana')) | |
as.list(query(MotifDb, 'hsapiens')) |