R/hct_beta.R

Defines functions hct_beta

Documented in hct_beta

hct_beta <-
function(pvalue, p, n)
{

	j=1:length(pvalue)
	nj1=p+1-j
	beta_pvalue=apply(cbind(pvalue, j, nj1), 1, beta_prob)
	k=which.max(beta_pvalue)
	pvalue[k]


}

Try the HDDesign package in your browser

Any scripts or data that you put into this service are public.

HDDesign documentation built on May 29, 2017, 8:13 p.m.