nullp: Probability Weighting Function

Description Usage Arguments Details Value Author(s) References See Also Examples

View source: R/nullp.R


Calculates a Probability Weighting Function for a set of genes based on a given set of biased data (usually gene length) and each genes status as differentially expressed or not.


nullp(DEgenes, genome, id,,



A named binary vector where 1 represents DE, 0 not DE and the names are gene IDs.


A string identifying the genome that genes refer to. For a list of supported organisms run supportedGenomes.


A string identifying the gene identifier used by genes. For a list of supported gene IDs run supportedGeneIDs.

A numeric vector containing the data on which the DE may depend. Usually this is the median transcript length of each gene in bp. If set to NULL nullp will attempt to fetch length using getlength.

Plot the PWF or not? Calls plotPWF with default values if TRUE.


It is essential that the entire analysis pipeline, from summarizing raw reads through to using goseq be done in just one gene identifier format. If your data is in a different format you will need to obtain the gene lengths and supply them to the nullp function using the arguement. Converting to a supported format from another format should be avoided whenever possible as this will almost always result in data loss.

NAs are allowed in the vector if you do not have information about a certain gene. Setting a gene to NA is preferable to removing it from the analysis.

If is left as NULL, nullp attempts to use getlength to fetch GO catgeory to gene identifier mappings.

It is recommended you review the fit produced by the nullp function before proceeding by leaving as TRUE.


A data frame with 3 columns, named "DEgenes", "" and "pwf" with the rownames set to the gene names. Each row corresponds to a gene with the DEgenes column specifying if the gene is DE (1 for DE, 0 for not DE), the column giving the numeric value of the DE bias being accounted for (usually the gene length or number of counts) and the pwf column giving the genes value on the probability weighting function. This object is usually passed to goseq to calculate enriched categories or plotPWF for further plotting.


Matthew D. Young [email protected]


Young, M. D., Wakefield, M. J., Smyth, G. K., Oshlack, A. (2010) Gene ontology analysis for RNA-seq: accounting for selection bias Genome Biology Date: Feb 2010 Vol: 11 Issue: 2 Pages: R14

See Also

supportedGenomes, supportedGeneIDs, goseq, getlength


pwf <- nullp(genes, 'hg19', 'ensGene')

goseq documentation built on Nov. 17, 2017, 12:24 p.m.