geneRxCluster (version 1.8.0)

critVal.alpha: critical regions

Description

critical region cutpoints

Usage

critVal.alpha(k, p0, alpha, posdiff)

Arguments

k

- window width(s)

p0

- length 2 probabilities

alpha

- two tailed

posdiff

- position difference matrix

Value

list of cutoffs and attributes

Details

This version uses alpha and will find TFD

See Also

gRxCluster for how and why this function is used

Examples

# symmetric odds:
crit <- critVal.alpha(5:25,c(1,1)/2,alpha=0.05,
                      matrix(1,nr=50,nc=21))
crit[[1]]
sapply(crit,c)
# 5:1 odds
asymmetric.crit <- critVal.alpha(5:25,c(1,5)/6,alpha=0.05,
                      matrix(1,nr=50,nc=21))
# show the critical regions
par(mfrow=c(1,2))
gRxPlot(crit,method="critical")
gRxPlot(asymmetric.crit,method="critical")
rm(crit,asymmetric.crit)