geneRxCluster (version 1.8.0)

critVal.power: critical regions

Description

critical region cutpoints

Usage

critVal.power(k, p0, target, pwr = 0.8, odds = 7)

Arguments

k

- window width(s)

p0

- length 2 probabilities

target

- false discoveries wanted

pwr

- desired power

odds

- alternative odds ratio

Value

list of cutoffs and attributes

Details

This version uses power and TFD and will limit windows screened

See Also

gRxCluster for how and why this function is used

Examples

# symmetric odds:
crit <-
    critVal.power(5:25,c(1,1),5,pwr=0.8,odds=7)
crit[[1]]
sapply(crit,c)
# 5:1 odds
asymmetric.crit <-
    critVal.power(5:25,c(1,5),5,pwr=0.8,odds=7)
# show the critical regions
par(mfrow=c(1,2))
gRxPlot(crit,method="critical")
gRxPlot(asymmetric.crit,method="critical")
rm(crit,asymmetric.crit)