Skip to content

Instantly share code, notes, and snippets.

@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 / processFileInChunks.R
Last active December 25, 2018 22:44
Streaming through large tabular files in R
################################################
## Streaming through large tabular files in R ##
################################################
## Author: Thomas Girke
## Last update: 21-Dec-2018
## Utility: the function 'processFileInLineChunks' streams through a file in
## batches of lines and applies to each imported line batch a function assigned
## to the 'myFct' argument. The number of lines processed in each iteration can
## be defined under the 'n_rows' argument. As importer functions, both 'fread'
@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 / Bioconductor_Git_Mirror.md
Last active April 15, 2020 00:36
How to maintain Bioconductor package on GitHub and keep it in sync with SVN repository on Bioconductor
@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'))
## Arranging plots with patchwork. For details see this manual page: https://ggplot2-book.org/arranging-plots.html
## Load libraries and create some sample data
library(ggplot2); library(reshape2); library(patchwork)
iris_mean <- aggregate(iris[,1:4], by=list(Species=iris$Species), FUN=mean)
df_mean <- melt(iris_mean, id.vars=c("Species"), variable.name = "Samples", value.name="Values")
## Create lots of plots
p1 <- ggplot(cbind(iris, Samples="Sepal.Length"), aes(Species, Sepal.Length)) + geom_boxplot(aes(fill = Species)) + ylim(4, 8) + theme(axis.title.x=element_blank(), axis.title.y=element_blank(), legend.position = "none") + facet_grid(rows = ~Samples) + ggtitle("(A)") + theme(plot.title = element_text(hjust=-0.3))
p2 <- ggplot(df_mean[df_mean$Samples=="Sepal.Width",], aes(Species, Values)) + geom_bar(position="dodge", aes(fill = Species), stat="identity") + yli
@tgirke
tgirke / crouton.md
Last active September 10, 2022 23:29
Crouton
@tgirke
tgirke / Jekyll_Doc_Theme.md
Last active December 8, 2022 19:28
How to set up and maintain Jekyll Documentation Theme
@tgirke
tgirke / gist-markdown-template.md
Last active April 4, 2023 00:51
Gist Markdown Template

Project Name

TODO: Write a project description

Installation

TODO: Describe the installation process

Usage

@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