It's not clear exactly what you mean by 'automate' but you can simplify
a bit by fitting a multivariate linear model to all the responses
together, and using . on the RHS of the formula to represent all
other variables in the data set as independent variables,
m.all <- glm(cbind(O3, temp) ~ ., data=ozone)
(assuming that only humidity, ibh and ibt remain; otherwise, use
data=subset(ozone, ...))
-Michael
On 11/4/2013 2:55 AM, Kumar Raj wrote:
I want to estimate the effect of several independent variables on several
dependent
variables. In the example below I wanted to estimate the
effect of three independent variables on ozone and temperature. My aim is
to create a list of dependent and independent variables and automate the
process rather than writing every dependent and independent variable in
each model as I have done below.
Example data is provided by the following library:
library(faraway)
data(ozone)
mo3 <- glm(O3 ~ humidity + ibh + ibt, data=ozone)
mtemp<- glm(temp ~ humidity + ibh + ibt, data=ozone)
Thanks
[[alternative HTML version deleted]]
--
Michael Friendly Email: friendly AT yorku DOT ca
Professor, Psychology Dept. & Chair, Quantitative Methods
York University Voice: 416 736-2100 x66249 Fax: 416 736-5814
4700 Keele Street Web: http://www.datavis.ca
Toronto, ONT M3J 1P3 CANADA
______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
