R/Pvero_ratio.R

Pvero_ratio <-
function (x,n,rho) {
ene=sum(n);a1=sum(x)
b=ene-(n[2]-x[2])+(ene-(n[1]-x[1]))*rho
raizp=b^2-4*ene*a1*rho
pvero=(b-sqrt(raizp))/(2*ene*rho)
pvero[pvero>1]=1;pvero[pvero<0]=0
pvero}

Try the prop.comb.RR package in your browser

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

prop.comb.RR documentation built on May 2, 2019, 12:41 p.m.