[R] border of a polygon in contour.kriging - geoR

Paulo Justiniano Ribeiro Jr paulojus at est.ufpr.br
Mon Feb 23 14:51:05 CET 2004


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




More information about the R-help mailing list