fitNull fits the null distribution to a set of negative controls

fitNull(logRr, seLogRr)

## Arguments

logRr A numeric vector of effect estimates on the log scale The standard error of the log of the effect estimates. Hint: often the standard error = (log() - log())/qnorm(0.025)

## Value

An object containing the parameters of the null distribution.

## Details

This function fits a Gaussian function to the negative control estimates as described in Schuemie et al (2014).

## References

Schuemie MJ, Ryan PB, Dumouchel W, Suchard MA, Madigan D. Interpreting observational studies: why empirical calibration is needed to correct p-values. Statistics in Medicine 33(2):209-18,2014

## Examples

data(sccs)
negatives <- sccs[sccs$groundTruth == 0, ] null <- fitNull(negatives$logRr, negatives\$seLogRr)
null#> Estimated null distribution
#>
#>      Estimate
#> Mean   0.7921
#> SD     0.2834