critVal.alpha: critical regions
In geneRxCluster: gRx Differential Clustering

Description Usage Arguments Details Value Author(s) See Also Examples

View source: R/Rutils.R

critical region cutpoints

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

`k`	- window width(s)
`p0`	- length 2 probabilities
`alpha`	- two tailed
`posdiff`	- position difference matrix

This version uses alpha and will find TFD

list of cutoffs and attributes

Charles Berry

gRxCluster for how and why this function is used

# 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)