inst/snippet/prob-dist02.R

In fastR: Foundations and Applications of Statistics Using R

require(MASS)                       # for fractions()
kernel <- function(x) { x*(x-3) * as.numeric(x >=0 & x <= 3) }
k <- 1 / integrate(kernel,0,3)$value; k
fractions(k)
g <- function(x) { k * kernel(x) }
integrate(g,0,3)                    # check that we have pdf
integrate(g,0,1)
fractions(integrate(g,0,1)$value)
integrate(g,0,2)
fractions(integrate(g,0,2)$value)
integrate(g,1,2)
fractions(integrate(g,1,2)$value)