R/gaussianissima.R

In Fioroni1863179/LabMec: Prova Individuale LabMec

#' Stima del valore vero
#'
#' La funzione da la migliore stima con incertezza del valore vero di un misurando, a partire da un vettore passato di misure distribuite secondo una gaussiana. Disegna inoltre l'istogramma delle misure sovrapponendo l'andamento della gaussiana di parametri stimati
#'
#' @param x
#' Il vettore di misure

gaussianissima <- function(x) {
  if(missing(x)) stop("manca il vettore di dati")

  len <- length(x)
  devstd <- sd(x)
  cat("Miglior stima:", mean(x), "+-", devstd/sqrt(len), "\n")
  cat("Incertezza campionaria (sulla misura singola):", devstd, "\n")
  width <- max(x) - min(x)
  nbins <- round(width*2/devstd)
  hist(x, breaks = nbins, col = "grey", freq = FALSE)
  curve(dnorm(x, mean=mean(x),sd=devstd), min(x), max(x), add=T, col="blue")
}