R/estimation.R

In mo-mahdavi/RandomGenerator:

get.vector <- function(address){
  a = read.table(address)
  return(as.numeric(a[1,]))
}

duest <- function(address){
  f = get.vector(address)
  return(c(min(f), max(f)))
}

brest <- function(address){
  f = get.vector(address)
  return(mean(f))
}

biest <- function(address){
  f = get.vector(address)
  miu = mean(f)
  variance = mean(f ^ 2) - (miu ^ 2)
  p = 1 - (variance / miu)
  n = ceiling(miu / p)
  return(c(n, p))
}

geest <- function(address){
  f = get.vector(address)
  return(1/mean(f))
}

expest <- function(address){
  return(geest(address))
}

gaest <- function(address){
  f = get.vector(address)
  miu = mean(f)
  variance = mean(f ^ 2) - (miu ^ 2)
  beta = miu / variance
  alpha = miu * beta
  #s = log(miu) - mean(log(f))
  #alpha = (3 - s + sqrt(((s - 3) ^ 2) + 24 * s)) / (12 * s)
  #beta = alpha / miu
  return(c(alpha, beta))
}

poest <- function(address){
  return(brest(address))
}

noest <- function(address){
  f = get.vector(address)
  return(c(mean(f), mean(f ^ 2) - (mean(f) ^ 2)))
}