Takes the output of HPDplot(...,plot=FALSE), adds two columns of adjusted p-values (z and mcmc) Disregards the entries corresponding to control genes
padj.qpcr(data, controls = NULL, method = "BH")
Output of HPDplot(...,plot=FALSE); may be several several such outputs concatenated with rbind
A vector of control gene names
p-value correction method (see function p.adjust), default is Benjamini-Hochberg
The dataframe derived from the input, with added columns "padj.z" and "padj.mcmc". See tutorial for examples.
Matz MV, Wright RM, Scott JG (2013) No Control Genes Required: Bayesian Analysis of qRT-PCR Data. PLoS ONE 8(8): e71448. doi:10.1371/journal.pone.0071448