Marino Taussig De Bodonia, Agnese
2010-Aug-03 16:04 UTC
[R] "glmulti": defining which intractions between variables are to be included
Hello,
I'm using the "glmulti" package to run models of all the possible
combinations of my variables. However, I am only interested in a few
interactions between my variables.
I have tried the equivalent of:
mod1<-lm(y~a+b+c+a:b)
glmulti(mod1, level=1)
mod2<-lm(y~a+b+c+a:b)
glmulti(mod2, level=2)
and
glmulti("y", c("a", "b", "c"),
exclude=c("a:c", "b:c" ), data=mydata)
where "mydata" is an object reading a text file.
However, the first code only gives me combinations of main effects (without
interactions), the second code doesn't allow me to exclude the a:c and b:c
interactions, and the third code doesn't allow me to define my variables as
factors becuase it uses character strings (and therefore runs a lot more models
than necessary).
Is there a way to use the third code, but tell R that some of the character
strings are factors?
Otherwise, since there are a lot more incteractions that I do not wish to test
than interactions that I do wish to test, is there an alternative by which
instead of having to exclude the interactions that I do not wish to test, I can
simply run models with the interaction terms that I include in my lm fromula?
Thank you in advance for your time,
Agnese
Marino Taussig De Bodonia, Agnese
2010-Aug-03 16:28 UTC
[R] FW: "glmulti": defining which intractions between variables are to be included
Hello,
I'm using the "glmulti" package to run models of all the possible
combinations of my variables. However, I am only interested in a few
interactions between my variables.
I have tried the equivalent of:
mod1<-lm(y~a+b+c+a:b)
glmulti(mod1, level=1)
mod2<-lm(y~a+b+c+a:b)
glmulti(mod2, level=2)
and
glmulti("y", c("a", "b", "c"),
exclude=c("a:c", "b:c" ), data=mydata)
where "mydata" is an object reading a text file.
However, the first code only gives me combinations of main effects (without
interactions), the second code doesn't allow me to exclude the a:c and b:c
interactions, and the third code doesn't allow me to define my variables as
factors becuase it uses character strings (and therefore runs a lot more models
than necessary).
Is there a way to use the third code, but tell R that some of the character
strings are factors?
Otherwise, since there are a lot more incteractions that I do not wish to test
than interactions that I do wish to test, is there an alternative by which
instead of having to exclude the interactions that I do not wish to test, I can
simply run models with only the interaction terms that I include in my lm
fromula?
Thank you in advance for your time,
Agnese