critVal.target {geneRxCluster} | R Documentation |
critical region cutpoints
critVal.target(k, p0, target, posdiff = NULL, ns)
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 |
This version uses TFD and will find alpha implicitly
list of cutoffs and attributes
Charles Berry
gRxCluster
for how and why this function is
used
# 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)