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]]

