BOL: Related items

plumber:An R package that converts your existing R code to a web API

BioJoker — Wed, 13 Mar 2019 19:20:10 -0500

plumber allows you to create a REST API by merely decorating your existing R source code with special comments. Take a look at an example.

# plumber.R

#* Echo back the input
#* @param msg The message to echo
#* @get /echo
function(msg=""){
  list(msg = paste0("The message is: '", msg, "'"))
}

#* Plot a histogram
#* @png
#* @get /plot
function(){
  rand <- rnorm(100)
  hist(rand)
}

#* Return the sum of two numbers
#* @param a The first number to add
#* @param b The second number to add
#* @post /sum
function(a, b){
  as.numeric(a) + as.numeric(b)
}

Address of the bookmark: https://www.rplumber.io/

Github replacement !

Rahul Nayak — Thu, 26 Sep 2019 03:42:39 -0500

For a number of reasons researchers have been trying out GitLab as a replacement for for both GitHub and various continuous integration systems, and have been exploring configurations useful for model-fitting pipelines. Researchers turned one of these into an example repository that shows how to use GitLab together with the Rocker Docker images and the drake build system to reproducibly run a project pipeline, using the cacheing functionality across all three tools to make things reasonably speedy and enable both local and remote builds.

Check it out at https://gitlab.com/ecohealthalliance/drake-gitlab-docker-example.

Scientist Bioinformatics Positions

Thu, 30 Jan 2020 06:53:40 -0600

Bioinformatics-Multi_Omics_Integration

https://www.researchgate.net/job/939073_Senior_Scientist_Bioinformatics-Multi_Omics_Integration

Senior_Scientist_Bioinformatics-Transcriptomics_Analysis

https://www.researchgate.net/job/939075_Senior_Scientist_Bioinformatics-Transcriptomics_Analysis-Belgium_France_Switzerland_The_Netherlands

Senior Scientist Bioinformatics - Network Analytics

https://www.researchgate.net/job/939070_Senior_Scientist_Bioinformatics-Network_Analytics_Belgium_France_Switzerland_the_Netherlands

Team Leader Bioinformatics Data Sciences - Mechelen, Belgium

https://www.researchgate.net/job/938787_Team_Leader_Bioinformatics_Data_Sciences-Mechelen_Belgium

Instructions for Creating Your Own R Package

Jit — Wed, 19 Feb 2020 01:22:48 -0600

The following is a step-by-step guide to creating your own R package. Even beyond this course, youmay find this useful for storing functions you create for your own research or for editing existingR packages to suit your needs.

This guide contains three different sets of instructions. If you use RStudio, you can follow the “Ba-sic Instructions” in Section 2 which involve using RStudio’s interface.

Address of the bookmark: http://web.mit.edu/insong/www/pdf/rpackage_instructions.pdf

PheWAS: R package is designed to provide an accessible interface to the phenome wide association study

BioStar — Thu, 30 Jul 2020 22:06:37 -0500

The PheWAS R package is designed to provide an accessible interface to the phenome wide association study. For a description of the methods available and some simple examples, please see the package vignette or the R documentation. For installation help, see below. ##Installing the PheWAS Package The PheWAS package can be installed using the devtools package. The following code when executed in R will get you started:

install.packages("devtools")
#It may be necessary to install required as not all package dependencies are installed by devtools:
install.packages(c("dplyr","tidyr","ggplot2","MASS","meta","ggrepel","DT"))
devtools::install_github("PheWAS/PheWAS")
library(PheWAS)

Address of the bookmark: https://github.com/PheWAS/PheWAS

R Shiny in Life Sciences – Top 7 Dashboard Examples

Neel — Fri, 01 Apr 2022 19:05:03 -0500

R Shiny is one of the easiest ways for developers to make production-ready dashboards when speed and functionality are crucial. Shiny is approachable with a lot of documentation available, and because of this, a lot of developers/researchers with non-coding backgrounds are able to produce some impressive results. The whole ecosystem is easy to get your head around and pretty much limitless with regard to what you can do.

Address of the bookmark: https://www.r-bloggers.com/2022/03/r-shiny-in-life-sciences-top-7-dashboard-examples/

GraphMap - A highly sensitive and accurate mapper for long, error-prone reads

Jit — Wed, 07 Jun 2017 04:18:16 -0500

GraphMap - A highly sensitive and accurate mapper for long, error-prone reads http://www.nature.com/ncomms/2016/160415/ncomms11307/full/ncomms11307.html

Features

    Mapping position agnostic to alignment parameters.
    Consistently very high sensitivity and precision across different error profiles, rates and sequencing technologies even with default parameters.
    Circular genome handling to resolve coverage drops near ends of the genome.
    E-value.
    Meaningful mapping quality.
    Various alignment strategies (semiglobal bit-vector and Gotoh, anchored).
    Overlapping of reads for de novo assembly.
    Transcriptome mapping through internal construction of a transcriptome from a given genomic reference and a GTF file.
    ...and much more.

GraphMap is also used as an overlapper in a new de novo genome assembly project called Ra (https://github.com/mariokostelac/ra-integrate).
Ra attempts to create de novo assemblies from raw nanopore and PacBio reads without requiring error correction, for which a highly sensitive overlapper is required.

Currently, development of a new spliced-alignment mode for mapping RNA-seq reads is under way.
Description of the current effort as well as how to reach the experimental implementation can be found here: doc/rnaseq.md.

Address of the bookmark: https://github.com/isovic/graphmap

MashMap: a fast and approximate software for mapping long reads (PacBio/ONT) or assembly to reference genome(s)

Jit — Tue, 12 Dec 2017 17:23:31 -0600

MashMap is a fast and approximate software for mapping long reads (PacBio/ONT) or assembly to reference genome(s). It maps a query sequence against a reference region if and only if its estimated alignment identity is above a specified threshold. It does not compute the alignments explicitly, but rather estimates a k-mer based Jaccard similarity using a combination of Winnowing and MinHash. This is then converted to an estimate of sequence identity using the Mash distance. An appropriate k-mer sampling rate is automatically determined given minimum local alignment length and identity thresholds. The efficiency of the algorithm improves as both of these thresholds are increased.

Address of the bookmark: https://github.com/marbl/MashMap

ALPACA: A hybrid strategy for assembly of genomic DNA shotgun sequencing reads.

Seema Singh — Mon, 30 Apr 2018 04:38:40 -0500

ALPACA requires Celera Assembler 8.3 or later. It is recommended to build Celera Assembler from source. (Why? The pre-built binaries CA_8.3rc1 and CA8.3rc2 will work for any large data set.

Detail paper at https://bmcgenomics.biomedcentral.com/articles/10.1186/s12864-017-3927-8

Address of the bookmark: https://github.com/VicugnaPacos/ALPACA

BlasR Mapping single molecule sequencing reads using Basic Local Alignment with Successive Refinement (BLASR): Theory and Application,

Jit — Wed, 23 May 2018 06:54:32 -0500

BLASR (Basic Local Alignment with Successive Refinement) for mapping Single Molecule Sequencing (SMS) reads that are thousands to tens of thousands of bases long with divergence between the read and genome dominated by insertion and deletion error.

Here is how I use the blasr to align PacBio reads to the contigs (target.fasta). The “target.fasta.sa” is the suffix array from “target.fasta” generated by sawriter.

blasr query.fa ./target.fasta -sa ./target.fasta.sa -bestn 40 -maxScore -500 -m 4 -nproc 24 -out target.m4 -maxLCPLength 15

the output format option “-m 4″ generate the alignment coordinate. Not fully documented, but I can explain that to you.

I use a 24 cores / 48G ram server for the alignment. It took about 2 to 3 hours aligning 3G PacBio Reads to 10^6 sequences of short read contigs with a mean 3.5kbp length.

Address of the bookmark: http://bix.ucsd.edu/projects/blasr/