critVal.power: critical regions

In geneRxCluster: gRx Differential Clustering

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

Details

This version uses power and TFD and will limit windows screened

Value

list of cutoffs and attributes

Author(s)

Charles Berry

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)

geneRxCluster documentation built on Nov. 8, 2020, 4:56 p.m.