Learn R Programming

MAST

Model-based Analysis of Single-cell Transcriptomics

Installation Instructions

If you have previously installed the package SingleCellAssay you will want to remove it as MAST superscedes SingleCellAssay. (If both MAST and SingleCellAssay are attached, opaque S4 dispatch errors will result.) Remove it with:

 remove.packages('SingleCellAssay')

Then you may install or update MAST with: install.packages('devtools') library(devtools) install_github('RGLab/MAST') # or if you don't have a working latex setup install_github(RGLab/'MAST', build_vignettes=FALSE) vignette('MAST-intro')

New Features

  • gseaAfterBoot for competitive geneset analysis under variance inflation
  • Support tests of arbitrary contrasts using LRT/zlm.SingleCellAssay

Copy Link

Version

Version

0.931

License

GPL(>= 2)

Maintainer

Last Published

December 22nd, 2017

Functions in MAST (0.931)

Drop

Drop specified dimension from an array
FluidigmAssay

Constructor for a FluidigmAssay
LMlike-class

Linear Model-like Class
LRT

Likelihood Ratio Tests for SingleCellAssays
Hypothesis

Describe a linear model hypothesis to be tested
LMERlike-class

Wrapper for lmer/glmer
FromMatrix

Construct a SingleCellAssay from a matrix or array of expression
GLMlike-class

Wrapper for regular glm/lm
BayesGLMlike-class

Wrapper for bayesian GLM
DataLayer-class

DataLayer class
SingleCellAssay

SingleCellAssay: A constructor for an object of type SingleCellAssay.
ThresholdedNanoString-class

Holds output and diagnostics from thresholdNanoString Not intended to be called by the user.
applyFlat

Apply a vectorized binary operation recycling over last dimension
bootVcov1

Bootstrap a zlmfit
condSd

Report standard deviation of et, for positive et for each gene
condmean

Report the mean et value for each gene
RNASeqAssay-class

RNASeqAssay class. Doesn't require ncells
outlierDetection

Filter an RNASeqAssay
combine,SingleCellAssay,ANY-method

Combine a SingleCellAssay and a vector, data.frame or AnnotatedDataFrame
SingleCellAssay-class

SingleCellAssay class
SingleCellAssay-package

"Tools for Single Cell Assay Analysis"
coef

Coefficients of zero-inflated
SCASet-class

SCASet is a set of SingleCellAssay objects or objects of its subclasses (i.e. FluidigmAssay) The constructor SCASet should be used to make objects of this class.
calcZ

Get Z or T statistics and P values after running gseaAfterBoot
cellData

Accessor for cellData AnnotatedDataFrame
combine,SCASet,missing-method

Combines a SCASet into a unified SingleCellAssay
defaultPrior

Initialize a prior to be used a prior for BayeGLMlike/BayesGLMlike2
deviance_residuals_hook

Standardized deviance residuals hook
ebayes

Estimate hyperparameters for hierarchical variance model for continuous component
expavg

Exponential average
getss

getss
SCASet

Constructs a SCASet
conform

Do two objects conform in dimension and type?
continuous_residuals_hook

Hook to get the discrete residuals
filterLowExpressedGenes

Filter low-expressing genes
fit

fit a zero-inflated regression
freq

Report the proportion of expression for each gene
getConcordance

getConcordance
ZlmFit-class

An S4 class to hold the output of a call to zlm
addlayer

Layer management methods
burdenOfFiltering

Burden of filtering
combine,SingleCellAssay,SingleCellAssay-method

Combine two SingleCellAssay or derived classes
combined_residuals_hook

Hook to get the combined residuals
computeEtFromCt

computeEtFromCt
exprs,DataLayer-method

Get or set a matrix of measurement values in a SingleCellAssay
fData

fData
hushWarning

Selectively muffle warnings based on output
fit,LMERlike,missing-method

##'
getwss

getwss
gseaAfterBoot

Gene set analysis for hurdle model
getrc

getrc
logFC

Calculate log-fold changes from hurdle model components
logmean

Log mean
cData

Accessor for cellData data.frame
discrete_residuals_hook

hHook to get the continuous residuals
dof

Degrees of freedom of Zero inflated model
getwellKey

Accessor for wellKey
layer

Which layer is active?
layername<-

Set name of active layer
influence.bayesglm

influence for bayesglm objects.
layername

Get name of active layer
length,SCASet-method

Get the length of an SCASet
plotSCAConcordance

Concordance plots of filtered single vs n-cell assays
plotlrt

Plot a likelihood ratio test object
split,SingleCellAssay-method

Split into SCASet
model.matrix

Model matrix accessor
myBiplot

Makes a nice BiPlot
sampleNames,SCASet-method

get the sample names in an SCASet
score_2_hook

Hook to compute the combined residuals with treatment effects added back.
nlayer

How many layers does a DataLayer have?
numexp

Report number of expressing cells per gene
read.fluidigm

read.fluidigm
lrTest,ZlmFit,character-method

Likelihood ratio test
lrTest

Run a likelihood-ratio test
pbootVcov1

Bootstrap a zlmfit
plot.thresholdSCRNACountMatrix

Plot cutpoints and densities for thresholding
show,Mapping-method

Show methods
split,DataLayer-method

Split a DataLayer correctly
summary,ZlmFit-method

Summarize model features from a ZlmFit object
summary

Return a human-readable summary
waldTest

Run a Wald test
xform

Make matrix of continuous expression values, orthogonal to discrete
featureData

Accessor for featureData AnnotatedDataFrame
filter

Filter a SingleCellAssay
lapply,SCASet-method

Apply a function to each member of a SCASet
rstandard.bayesglm

rstandard for bayesglm objects.
[[<-,DataLayer-method

Replace a section of a DataLayer
subset

Subset a SingleCellAssay by data in cellData
[,DataLayer,ANY,ANY-method

Subset a DataLayer
[[,DataLayer-method

Extract from a DataLayer
[[,SCASet-method

Index into an SCASet object
layer<-

Set active layer
melt.SingleCellAssay

Melt a rectangular array
model.matrix<-

Replace model matrix
summary.thresholdSCRNACountMatrix

Summarize the effect of thresholding
thresholdNanoString

thresholds for positive expression
zlm

Convenience function for running a zero-inflated regression
zlm.SingleCellAssay

Zero-inflated regression for SingleCellAssay
[,SCASet,ANY,ANY-method

Subset an SCASet to a smaller SCASet
[,SingleCellAssay,ANY,ANY-method

Subset a SingleCellAssay
vcov

Variance-covariance matrix for zero inflated
vbeta

Vbeta Data Set
vbetaFA

Vbeta Data Set, FluidigmAssay
waldTest,ZlmFit,matrix-method

Wald test
primerAverage

Average within duplicated genes/primers
print.summaryZlmFit

Print summary of a ZlmFit
score_hook

Used to compute module "scores"
shalek_weights

Hook to return p_hat from the model
summarize

Return programmatically useful summary of a fit
summary,SingleCellAssay-method

Summarize expression parameters
thresholdSCRNACountMatrix

Threshold a count matrix using an adaptive threshold.
todt

Convert a pair of character names to a list that can be placed into a data.table