zscore {limma} | R Documentation

## Z-score Equivalents

### Description

Compute z-score equivalents of non-normal random deviates.

### Usage

    zscore(q, distribution, ...)
    zscoreGamma(q, shape, rate = 1, scale = 1/rate)
    zscoreT(x, df, approx=FALSE)
    tZscore(x, df)
    zscoreHyper(q, m, n, k)

### Arguments

`q, x` | numeric vector or matrix giving deviates of a random variable
`distribution` | character name of probabability distribution for which a cumulative distribution function exists
`...` | other arguments specify distributional parameters and are passed to the cumulative distribution function
`shape` | gamma shape parameter (>0)
`rate` | gamma rate parameter (>0)
`scale` | gamma scale parameter (>0)
`df` | degrees of freedom (>0 for `zscoreT` or >=1 for `tZscore`)
`approx` | logical, if `TRUE` then a fast approximation is used to convert t-statistics into z-scores. If `FALSE`, z-scores will be exact.
`m` | as for `qhyper`
`n` | as for `qhyper`
`k` | as for `qhyper`

### Details

These functions compute the standard normal deviates which have the same quantiles as the given values in the specified distribution. For example, if `z <- zscoreT(x,df=df)` then `pnorm(z)` equals `pt(x,df=df)`.

`zscore` works for any distribution for which a cumulative distribution function (like `pnorm`) exists in R. The argument `distribution` is the name of the cumulative distribution function with the `"p"` removed.

`zscoreGamma`, `zscoreT` and `zscoreHyper` are specific functions for the gamma, t and hypergeometric distributions respectively.

`tZscore` is the inverse of `zscoreT`, and computes t-distribution equivalents for standard normal deviates.

The transformation to z-scores is done by converting to log tail probabilities, and then using `qnorm`. For numerical accuracy, the left or right tail is used, depending on which is likely to be smaller.

If `approx=TRUE`, then the approximation from Hill (1970) is used to convert t-statistics to z-scores directly without computing tail probabilities. Brophy (1987) showed this to be most accurate of a variety of possible closed-form transformations.

### Value

Numeric vector giving equivalent deviates from the standard normal distribution. The exception is `tZscore` which gives deviates from the specified t-distribution.

### Author(s)

Gordon Smyth

### References

Hill, GW (1970). Algorithm 395: Student's t-distribution. Communications of the ACM 13, 617-620.

Brophy, AL (1987). Efficient estimation of probabilities in the t distribution. Behavior Research Methods 19, 462–466.

### See Also

`qnorm`, `pgamma`, `pt` in the stats package.

### Examples

    # First three are equivalent
    zscore(c(1,2.5), dist="gamma", shape=0.5, scale=2)
    zscore(c(1,2.5), dist="chisq", df=1)
    zscoreGamma(c(1,2.5), shape=0.5, scale=2)

    zscoreT(2, df=3)
    tZscore(2, df=3)

[Package limma version 3.36.5 Index]