get.uri.2d: Compute the correspondence profile (Psi(t, v)) and the...

In idr: Irreproducible Discovery Rate

Compute the correspondence profile (Psi(t, v)) and the derivative of the correspondence profile (Psi'(t, v))

Description

An internal function to compute the correspondence profile (Psi(t, v)) and the derivative of the correspondence profile (Psi'(t, v)) for t and v between 0 and 1

Usage

get.uri.2d(x1, x2, tt, vv, spline.df = NULL)

Arguments

x1
x2
tt	A vector of right-tail percentages for x1. It is a numeric vector between 0 and 1, e.g. seq(0.01, 0.99, by=0.01)
vv	A vector of right-tail percentages for x2. It is a numeric vector between 0 and 1, e.g. seq(0.01, 0.99, by=0.01)
spline.df	The degree of freedom to control the smoothness of the spline to be superimposed on the points

Author(s)

Qunhua Li

References

Q. Li, J. B. Brown, H. Huang and P. J. Bickel. (2011) Measuring reproducibility of high-throughput experiments. Annals of Applied Statistics, Vol. 5, No. 3, 1752-1779.

See Also

comp.uri, get.correspondence

Examples

# salmon data from Kallenberg and Ledwina, 1999
data(salmon)

# get.correspondence() needs large ranks have high correlation and
# small ranks have low correlation. In this dataset, small values
# have high correlation and large values have low correlation.
# Ranking negative values makes the data follow the structure required
# by get.correspondence().
# There are 28 observations in this data set.

rank.x <- rank(-salmon$spawners)
rank.y <- rank(-salmon$recruits)
t <- seq(0.01, 0.99, by=1/28)
psi.dpsi <- get.uri.2d(rank.x, rank.y, t, t, spline.df=6.4)

idr documentation built on June 21, 2022, 9:05 a.m.