Thank you akk.
I know it is not statistically sounded to check the distribution of
response before glm.
I will check the distribution of xmodel$residuals later on.
About the program problem.
It can print summary(xmodel) but not confint(xmodel) by amending my
code as suggested by Bill Venables.
R
csiro.au> writes:
>
> Finally, I'm a bit puzzled why you use glm() when the simpler lm() would
> have done the job. You are fitting a linear model and do not need the
> extra paraphernaila that generalized linear models require.
>
> Bill Venables.
>
Perhaps the original poster is confuse
ECTED] On Behalf Of Chung-hong Chan
Sent: Sunday, 6 May 2007 6:47 PM
To: r-help@stat.math.ethz.ch
Subject: [R] 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.
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")