R/mcpv.function.R

mcpv.function <-
function(obs.stat,ms.stat,direction)
{
	D=get(direction)
	pv=apply(t(t(obs.stat)),1,function(x) (sum(D(ms.stat,x))+1)/(length(ms.stat)+1))
	return(pv=pv)
}

Try the SNscan package in your browser

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

SNscan documentation built on May 2, 2019, 8:26 a.m.