R/covariance.R

In Fioroni1863179/LabMec: Prova Individuale LabMec

#' Calcolo della covarianza tra due vettori di misure
#'
#' Calcola la covarianza ed il coefficiente di correlazione tra i due vettori passati
#'
#' @param x,y
#' Vettori di misure

covariance <- function(x, y){
  if(missing(x)) stop("manca un vettore di misure")
  if(missing(y)) stop("manca un vettore di misure")

  meanx = mean(x)
  meany = mean(y)
  cov <- sum( (x-meanx)*(y-meany) )/(length(x)-1)
  cat("Covarianza: ", cov, "\n")
  ccor <- cov/(sd(x)*sd(y))
  cat("Coefficiente di correlazione ", ccor, "\n")
}