R/globalFstat.R

In WRI: Wasserstein Regression and Inference

#' An internal function used in bootstrap global F test
#' @name globalFstat
#' @keywords internal
globalFstat <- function(xfit, Qobs, t_vec) {

        Qmat_lm <- stats::lm(Qobs ~ xfit)
        Qfitted = fitted(Qmat_lm)
        n = nrow(Qobs)
        Qmean = colMeans(Qobs)
        ssr = sum(sapply(1:n, function(i) fdapace::trapzRcpp(t_vec, (Qmean - Qfitted[i, ])^2)))
        return(ssr)
}