Skip to content

Instantly share code, notes, and snippets.

@tgirke
tgirke / Jekyll_Doc_Theme.md
Last active December 8, 2022 19:28
How to set up and maintain Jekyll Documentation Theme
@tgirke
tgirke / PeakSeq.sh
Last active May 29, 2018 00:50
How to run PeakSeq
##################################################
## 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
@tgirke
tgirke / reactiveUrl.R
Last active July 26, 2019 13:37
Shiny app example with reactive URL to query database with IDs
###############################################
## 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)
@tgirke
tgirke / BiocGitHub.md
Last active June 11, 2023 23:09
New Bioc-GitHub Source Control

Bioconductor-GitHub Package Sync

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).


Table of Content

@tgirke
tgirke / _Nvim-R-Tmux.md
Last active October 13, 2023 05:29
Nvim-R-Tmux: An Integrated Working Environment for R

Nvim-R-Tmux: Neovim-based IDE for R

!!! This Gist tutorial is deprecated. Its new version is available here !!!

This is some red text.

  • Author: Thomas Girke
  • Last update: 18-Nov-2020

###################################################################
## 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
@tgirke
tgirke / kegg_cmp_download.R
Last active March 8, 2018 15:35
KEGG CMP Download into R
####################################################
## 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 ##
################################
@tgirke
tgirke / crouton.md
Last active September 10, 2022 23:29
Crouton
@tgirke
tgirke / hmmscan.R
Last active May 28, 2018 21:57
Map Pfam domain with HMMER3
##############################################
## 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'
@tgirke
tgirke / motifEnrich.R
Last active April 11, 2021 23:58
Motif enrichment for ChIP-Seq
###############################
## 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'))