Skip to content
Prev 299408 / 398503 Next

Declaring a density function with for loop

David L Carlson 
dnewdis <- function(x, a, b, n) {
    i <- 0:exp(10)
    term <-sum((((-1)**i)*(choose(b-1,i))*(beta(x+a+a*i,n-x+1))))
    a*b*choose(n,x)*term
}

----------------------------------------------
David L Carlson
Associate Professor of Anthropology
Texas A&M University
College Station, TX 77843-4352

Thread (5 messages)

chamilka Declaring a density function with for loop Jul 7 Berend Hasselman Declaring a density function with for loop Jul 7 chamilka Declaring a density function with for loop Jul 7 David L Carlson Declaring a density function with for loop Jul 7 chamilka Declaring a density function with for loop Jul 8