R/residui.R

In Fioroni1863179/LabMec: Prova Individuale LabMec

#' Calcolo dei residui e della sigma
#'
#' La funzione accetta come parametro l'oggetto restituito dal fit lineare e calcola i residui e la migliore stima della sigma su ogni misura
#'
#' @param fitobj
#' L'oggetto restituito dalla funzione linearfit
#' @param verbose
#' DEFAULT = FALSE
#'
#' Se il valore è TRUE vengono stampati i valori di tutti i singoli residui

residui <- function(fitobj, verbose = FALSE){
  d = fitobj$y-fitobj$m*fitobj$x-fitobj$c
  if(verbose){
    cat("\n", "Singoli residui:", "\n", d, "\n")
  }
  cat("Sigma stimata dai residui:", sqrt(sum(d^2)/(length(fitobj$x)-2)))
}