| 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)