R/poissonissima.R

In Fioroni1863179/LabMec: Prova Individuale LabMec

#' Stima del valore di lambda e del rate
#'
#' La funzione fa inferenza su lambda e su r a partire da due vettori contenenti rispettivamente i conteggi effettuati e i tempi di esposizione
#'
#' @param x
#' Il vettore contenente i conteggi
#' @param t
#' Il vettore contenente i tempi di esposizione per ogni misurazione
#'
poissonissima <- function(x, t){
  if(missing(x)) stop("manca il vettore dei conteggi")
  if(missing(t)) stop("manca il vettore dei tempi")

  count <- sum(x)+1
  ttot <- sum(t)
  cat("Lambda:", count, "+-", sqrt(count), "\n")
  cat("Moda di lambda:", count-1, "\n")
  cat("rate:", count/ttot, "+-", sqrt(count)/ttot, "\n")
}