 
Dear list, 

I'm using the coal ash percentage data to perform a median polish kriging.

Why the the residuals are different from data-overall effect-row effect-column effect?



library(gstat)
data(coalash)
coal.mat <- tapply(coalash$coalash, 
       list(factor(coalash$x),factor # percentages according to their
       (coalash$y)),function(x)x) # locations.
coal.mp <- medpolish(coal.mat, na.rm=T) 

Thanks in advance,
Gema




Gema Fernández-Avilés 

Department of Statistics
Faculty of Law and Social Sciences
University of Castilla-La Mancha at Toledo

http://www.uclm.es/profesorado/gemafaviles/


	[[alternative HTML version deleted]]

