pfbar <- function(x, df1, df2, wt.bar) {
if (x <= 0) {
return(0)
}
zed <- df1 == 0
cdf <- ifelse(any(zed), wt.bar[zed], 0)
cdf <- cdf + sum(pf(x/df1[!zed], df1[!zed], df2) * wt.bar[!zed])
return(cdf)
}
pchibar <- function(x, df, wt.bar) {
if (x <= 0) {
return(0)
}
zed <- df == 0
cdf <- ifelse(any(zed), wt.bar[zed], 0)
cdf <- cdf + sum(pchisq(x, df[!zed]) * wt.bar[!zed])
return(cdf)
}
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.