border of a polygon in contour.kriging - geoR
Dear Patrick border=NULL should work but there is a bug. I'm fixing this right now and a new version will be available later today P.J.
On Mon, 23 Feb 2004, Patrick Giraudoux wrote:
Dear all, When a conventional kriging and then a contour plot is limited with a polygon (as possible with krige.conv and contour.kriging), the polygon border is displayed in black by default.
kc<-krige.conv(CZdata,loc=pred.grid,borders=czpoly,krige=krige.control(obj.m=ls)) plot(CZcoord,xlab="x",ylab="y",type="n",asp=1) contour(kc,loc=pred.grid,add=T)
I would like to be capable to get the contour lines only (limited by the polygon, but without the polygon border displayed). I tried to pass the argument border=NULL, unsuccessfully:
contour(kc,loc=pred.grid,add=T,border=NULL)
Error in contour(as.double(x), as.double(y), z, as.double(levels), labels, :
dimension mismatch
In addition: Warning message:
Replacement length not a multiple of the elements to replace in matrix(...)
Any hint? Thanks in advance, Patrick Giraudoux
______________________________________________ R-help at stat.math.ethz.ch mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
Paulo Justiniano Ribeiro Jr Departamento de Estat?stica Universidade Federal do Paran? Caixa Postal 19.081 CEP 81.531-990 Curitiba, PR - Brasil Tel: (+55) 41 361 3471 Fax: (+55) 41 361 3141 e-mail: pj at est.ufpr.br http://www.est.ufpr.br/~paulojus