R/binomialissima.R

In Fioroni1863179/LabMec: Prova Individuale LabMec

#' Stima del valore della probabilità di ogni successo
#'
#' La funzione fa inferenza sulla probabilità p di ogni processo di Bernoulli a partire dal numero di successi ottenuti su \emph{n} misure
#'
#' @param x
#' Il numero di successi ottenuti
#' @param n
#' Il numero di tentativi effettuati

binomialissima <- function(x, n){
  if(missing(x)) stop("manca il numero di successi")
  if(missing(n)) stop("manca il numero di prove")

  p <- (x+1)/(n+2)
  sigmap <- sqrt(p*(1-p)/(n+3))
  cat("Probabiliotà di successo:", p, "+-", sigmap, "\n")
  cat("Moda:", x/n)
}