I would like to run an anova on a set of raster variables imported from
GRASS into R but am confused about how R handles the raster data. So if
I have a model with a treatment with 2 levels, each raster is a binary
file with cells coded for each level of the treatment. I have 1 response
raster layer that corresponds spatially with the treatment layers. So
let's say I have 5 rasters representing a different rep of the treatment
how would I code this in R? I started out trying this:
treatment<-readRAST6(c("rep1","rep2","rep3","rep4","rep5"),cat=TRUE)
response<-readRAST6("response",cat=FALSE)
test.aov<-aov(response~treatment)
but I don't think this is right. Is it possible to do what I am trying?
My data is actually a bit more complicated than this but if I can clear
this up than I think I can do the rest on my own.
Thanks for your help,
Sam
--
****************************************************
Sam Veloz
Postdoctoral Researcher
Department of Environmental Science and Policy
University of California, Davis
[EMAIL PROTECTED]
_______________________________________________
grass-stats mailing list
[email protected]
http://lists.osgeo.org/mailman/listinfo/grass-stats
