search for: logtrans

...800 48990 65400 79000 285000 2 > library(MASS) > boxcox(price ~ doors+CC+KW+KG+LENGTH, lambda=seq(-0.1,0.5,length=30)) > boxcox(price ~ doors+CC+KW+KG+LENGTH, lambda=seq(-0.5,0.1,length=30)) This all works fine, I get a fairly sharp peak near lambda=-0.25, but then: > logtrans(price ~ doors+CC+KW+KG+LENGTH, alpha=seq(-10000,0,length=30)) Error in eval(expr, envir, enclos) : Object "price" not found Comments welcome, Murray -- Dr Murray Jorgensen http://www.stats.waikato.ac.nz/Staff/maj.html Department of Statistics, University of Waikato, Hamilton, New...

Dear all, I coded a function called u.glm u.glm <- function (x,ahi,age,bmiz,gender) { library(nortest) lil.rslt <- lillie.test(x) if (lil.rslt$p.value >0.05) { cat("Logtrans=0, lillie=",lil.rslt$p.value,"\n") xmodel<-glm(x~ahi+age+bmiz+as.factor(gender)) summary(xmodel) confint(xmodel) } else { cat("Logtrans=1, lillie=",lil.rslt$p.value,"\n") xmodel<-glm(x~ahi+age+bmiz+as.factor(gender)) summary(xmodel)...

Hi all, Does anyone know of a function that return the best lambda in box-cox transformation? many thanks for any help. -- Cristiane S. Rocha Laboratorio Genoma Funcional - Bioinform?tica Centro de Biologia Molecular e Engenharia Genetica Universidade Estadual de Campinas Campinas - SP - Brasil Tel:(19)3788-1119