R/get.alpha.R

In SMPracticals: Practicals for Use with Davison (2003) Statistical Models

Defines functions get.alpha

Documented in get.alpha

get.alpha <- function(d)
{  # estimate alpha from data
  rho <- d$y/d$x
  n <- length(d$y)
  mean(rho)^2/( (n-1)*var(rho)/n - mean(rho)*mean(1/d$x) )
}

Try the SMPracticals package in your browser