Skip to content

Instantly share code, notes, and snippets.

@sgibb
sgibb / dispatch.R
Created July 17, 2019 14:45
s4 dispatching
bin <- function(x)x
setGeneric("bin", function(x) standardGeneric("bin"))
setMethod("bin", "matrix", function(x)x)
i <- 1:10
m <- matrix(1:10, nrow=2)
showMethods("bin")
# Function: bin (package .GlobalEnv)
@sgibb
sgibb / ondiskvec.R
Created March 3, 2019 21:27
simple on-disk vector
setClass("ondiskvec",
slots=list(path="character", n="numeric"),
prototype=list(path=character(), n=numeric())
)
ondiskvec <- function(x, path=tempfile()) {
writeBin(as.double(x), path, size=8L)
new("ondiskvec", path=path, n=length(x))
}
@sgibb
sgibb / class-draft.R
Created January 24, 2019 13:30
MSnbase3 - multiple MSnExp backends
library("mzR")
library("msdata")
library("rhdf5")
library("digest")
.sha1 <- function(x)vapply(x, digest::sha1, NA_character_, USE.NAMES=FALSE)
setClass(
"Backend",
contains="VIRTUAL"
)
@sgibb
sgibb / MassSpectrumOnDisk.R
Created November 6, 2018 21:22
Example of an on-disk MassSpectrum class for MALDIquant
library("MALDIquant")
setClass("MassSpectrumOnDisk",
slots=list(path="character", n="numeric"),
prototype=list(mass=numeric(), intensity=numeric(), metaData=list(),
path=character(), n=numeric()),
contains="AbstractMassObject")
## overwrite default validation method (checks for non-empty intensity etc)
## needs to be rewritten obviously
@sgibb
sgibb / figure1.mmd
Created October 17, 2017 19:23
topdownr paper figure1
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.

Keybase proof

I hereby claim:

  • I am sgibb on github.
  • I am sgibb (https://keybase.io/sgibb) on keybase.
  • I have a public key whose fingerprint is 22C0 E8AC A81A 3A13 C788 BCC4 ECAA 2C0A 2843 245B

To claim this, I am signing this object:

@sgibb
sgibb / memstat.sh
Created January 22, 2017 16:55
monitor synapter memory usage
#!/bin/sh
PROCESSPID=${1}
pidstat -urIh -p ${PROCESSPID} 1 1 | grep -v "^Linux\|^ *$" > memstat.txt && pidstat -urIh -p ${PROCESSPID} 60 | grep -v "^#\|^ *$\|^Linux" | tee -a memstat.txt
@sgibb
sgibb / acqu
Last active December 26, 2016 16:11
correctly read d10/1ref/acqu
##TITLE= XMASS Parameter file
##JCAMPDX= 5.0
##DATATYPE= CONTINUOUS MASS SPECTRUM
##ORIGIN= XMASS, Bruker-Daltonics (USA) and Bruker-Daltonik (GER)
##OWNER= TOF User
##SPECTROMETER/DATASYSTEM= Bruker Flex Series
##.SPECTROMETER TYPE= TOF
##.INLET= DIRECT
##.IONIZATION MODE= LD+
##$ACQMID= 0
@sgibb
sgibb / plotProteome.R
Created September 19, 2016 08:24
plot proteomes in a similar manner as Kuharev et al 2015 10.1002/pmic.201400396
library("synapter")
library("MSnbase")
load_all()
combMSnSet <- readRDS("refCombMSNSet.RDS")
plotProteomes <- function(x,
xlim=log2(range(exprs(x), na.rm=TRUE)),
ylim=c(-4, 4),
xlab=expression(log[2](B)),
@sgibb
sgibb / ionmobility.Rmd
Created April 29, 2014 14:57
synapter ion mobility example
# Ion Mobility Example
```{r setup, include=FALSE}
library(knitr)
opts_chunk$set(fig.height=10, fig.width=14, fig.path="figure/ionmobility/", cache=TRUE, autodep=TRUE, tidy=FALSE)
```
## Default Synapter Preprocessing
(same as in `synergise`)