R/FindBestTransform.R

In PoiClaClu: Classification and Clustering of Sequencing Data Based on a Poisson Model

Defines functions FindBestTransform

Documented in FindBestTransform

FindBestTransform <-
function(x){
  alphas <- seq(.01, 1, len=50)
  gof <- rep(NA, length(alphas))
  for(alpha in alphas){
    gof[alphas==alpha] <- GoodnessOfFit(x^alpha,type="mle")
  }
  return(alphas[which.min(abs(gof-(nrow(x)-1)*(ncol(x)-1)))])
}

Try the PoiClaClu package in your browser