Displaying 1 result from an estimated 1 matches for "bmiz".
Did you mean:
biz
2007 May 07
1
Simple question about function with glm
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, li...