plot.sparncp: plot an object of class sparncpt, i.e., semiparametric...

In gitlongor/pi0: Estimating the Proportion of True Null Hypotheses for FDR

Description

Plot the histogram of observed t-statistics together with its fitted density estimate; also plotted is the estimated density of noncentrality parameters.

Usage

## S3 method for class 'sparncpt'
plot(x,...)

Arguments

x	an object of class sparncpt
...	currently not used

Details

Left panel shows the density estimate of observed t-statistics, overlapped with a histogram; right panel shows the estimated density of noncentrality parameters. Solid line is the actual mean of the estimate; dashed line is located at zero.

Value

the invisible x itself

Author(s)

Long Qu

References

Qu L, Nettleton D, Dekkers JCM. (2012) Improved Estimation of the Noncentrality Parameter Distribution from a Large Number of $t$-statistics, with Applications to False Discovery Rate Estimation in Microarray Data Analysis. Biometrics, 68, 1178–1187.

See Also

parncpt, nparncpt, sparncpt

Examples

## Not run: 
data(simulatedTstat)
(npfit=nparncpt(tstat=simulatedTstat, df=8)); 
(pfit=parncpt(tstat=simulatedTstat, df=8, zeromean=FALSE)); plot(pfit)
(pfit0=parncpt(tstat=simulatedTstat, df=8, zeromean=TRUE)); plot(pfit0)
(spfit=sparncpt(npfit,pfit)); plot(spfit)

## End(Not run)

gitlongor/pi0 documentation built on May 17, 2019, 5:29 a.m.