R/tbsb.R

Defines functions Tbsb

Tbsb <- function(c, p) {
    ksiint <- function(s) {
        (2^s) * gamma(s + p/2) * pgamma(c^2/2, s + p/2)/gamma(p/2)
    }
    
    y1 <- ksiint(1) * 3/c - ksiint(2) * 3/(c^3) + ksiint(3)/(c^5)
    y2 <- c * (1 - pchisq(c^2, p))

    y1 + y2
}

Try the riv package in your browser

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

riv documentation built on May 24, 2018, 9:04 a.m.