R/con_pbar.R

In LeonardV/restriktor: Restricted Statistical Estimation and Inference for Linear Models

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)
}