R help - Sep 2010 - Stats Question

I have site data with variables that vary across sites (like wind speed,
moisture content) and within sites I have experimental units (logs) with
associated variables (like decay class, suspension). With normal response y one
can use R to get the between and within site information using
aov(y~u+v+Error(site)) where u and v are independent variables. Problem is we
have a binary response (burn or not), so need to use glm. I need to know if it
is possible to use glm with multistrata and if so how.


	[[alternative HTML version deleted]]

On Sep 13, 2010, at 8:27 PM, <Jennifer.Hollis at csiro.au> wrote:
> I have site data with variables that vary across sites (like wind  
> speed, moisture content) and within sites I have experimental units  
> (logs) with associated variables (like decay class, suspension).  
> With normal response y one can use R to get the between and within  
> site information using aov(y~u+v+Error(site)) where u and v are  
> independent variables. Problem is we have a binary response (burn or  
> not), so need to use glm. I need to know if it is possible to use  
> glm with multistrata and if so how.
Logistic regression is done in R with glm using family="binomial". See
the help pages.


-- 
David Winsemius, MD
West Hartford, CT