critVal.target: critical regions

In geneRxCluster: gRx Differential Clustering

Description

critical region cutpoints

Usage

  critVal.target(k, p0, target, posdiff = NULL, ns)

Arguments

k	window width(s)
p0	length 2 probabilities
target	- two tailed
posdiff	- position difference matrix
ns	the number of windows passing filter at each k

Details

This version uses TFD and will find alpha implicitly

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.target(5:25,c(1,1),1,ns=rep(10,21))
crit[[1]]
sapply(crit,c)
# 5:1 odds
asymmetric.crit <- critVal.target(5:25,c(1,5),1,ns=rep(10,21))
# 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.