critVal.alpha {geneRxCluster} | R Documentation |
critical region cutpoints
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)