R/dneymanA.R

dneymanA <- Vectorize(function(x, lambda1, lambda2, K, log=FALSE) {
  j <- 0:K
  dna <- exp(-lambda1)*lambda2^x/factorial(x)*sum(lambda1*exp(-lambda2)^j*j^x/(factorial(j)))
  if(log) return(log(dna)) else return(dna)
}, "x")

Try the jointNmix package in your browser

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

jointNmix documentation built on May 2, 2019, 8:18 a.m.