R/fun.beta1.R

"fun.beta1" <-
function(a, b){
index<-as.logical((a>0)*(b>0))

result<-gamma(a[index])*gamma(b[index])/gamma(a[index] + b[index])
result[!index]<-NA
return(result)
}

Try the GLDEX package in your browser

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

GLDEX documentation built on Aug. 21, 2023, 9:08 a.m.