R/hct_beta.R

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 2, 2019, 6:41 a.m.