R/find.beta.R

"find.beta" <-
function(m,alpha)
  {  
    quant.ecd <- -log(0.5*log(1/(1-alpha))) 
    c.n <- 2*log(log(m))+0.5*log(log(log(m)))-0.5*log(4*pi)
    b.n <- sqrt(2*log(log(m)))
    beta <- 1/sqrt(m)*(quant.ecd+c.n)/b.n 
  }

Try the howmany package in your browser

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

howmany documentation built on May 2, 2019, 4:18 a.m.