Can you krige with a nominal response variable
I might consider using gam() in R to smooth the 0/1 indicator for each class individually based on a Bernoulli model - doing each class individually. A downside of course is that the resulting smoothed maps of the proportions for each class would not add to one at an individual location, but they would probably be reasonably close and good enough for some basic analysis. If you're talking about truly modeling a joint distribution for the 10 classes spatially, it's starting to sound like a more complicated modeling problem based on a hierarchical statistical model. Not aware of any standard recipe or software for this. chris ------------------------------------------------------------------------------- Chris Paciorek / Asst. Professor Email: paciorek at hsph.harvard.edu Department of Biostatistics Voice: 617-432-4912 Harvard School of Public Health Fax: 617-432-5619 655 Huntington Av., Bldg. 2-407 WWW: www.biostat.harvard.edu/~paciorek Boston, MA 02115 USA Permanent forward: paciorek at alumni.cmu.edu
"Pierce, Ken" <ken.pierce at oregonstate.edu> 12/06/06 2:31 PM >>>
Is there a method for kriging a response with multiple states, for instance vegetation class in 10 classes? I know indicator kriging would work for a boolean variable. Is there an extension for multiple responses? Kenneth B. Pierce Jr. Research Ecologist Landscape Ecology, Modeling, Mapping and Analysis Team PNW Research Station - USDA- FS 3200 SW Jefferson Way, Corvallis, OR 97331 ken.pierce at oregonstate.edu 541 750- 7393 http://www.fsl.orst.edu/lemma/gnnfire _______________________________________________ R- sig- Geo mailing list R- sig- Geo at stat.math.ethz.ch https://stat.ethz.ch/mailman/listinfo/r- sig- geo