Subject: Re: spikes in contour and persp (PR#327)
From: Peter Dalgaard BSA (firstname.lastname@example.org)
Date: Tue 16 Nov 1999 - 23:15:59 EST
> > From: email@example.com
> > Date: Tue, 16 Nov 1999 08:57:00 +0100 (MET)
> > The following matrix of normed likelihoods should give a smooth
> > surface but instead gives a series of spikes in both persp and contour
> > (the dim labels are the axes values). I know an algorithm cannot be
> > infallible but it would be nice to have some parameters to control the
> > smoothing. (It takes close to an hour to produce this matrix on a
> > Pentium II 300mh. It was after that that it crashed with the Inf in my
> > previous bug message...)
> I am sorry, but I don't understand. contour does not do any
> smoothing, but interpolates assuming smoothness. The data
> appears to me to be a series of spikes just off the diagonal, and not
> what contour thinks of as a smooth surface. It is common to
> interpolate the data from such a coarse grid before contouring.
> persp just plots the data with visual linear interpolation.
To elaborate a little: The contouring is of course algorithm
Let's look at a section of the data:
1 17 15 1 0 0
0 3 34 21 1 0
0 0 6 57 26 1
0 0 0 12 79 29
One contouring method works by linear interpolation on triangles, so when
creating the line for, say z=25, it will draw four lines *around* 34,
then one connecting the 21-57 line, the 6-57, the 12-57 and the 12-79
and on the other end 21-26 (diag) 26-1 (vertical) 26-1 (horizontal) and 1-29.
This gives two disjoint contours and an apparent peak at 34. The
algorithm in R is pretty similar to this as far as I recall.
Alternatively one could have bilinear interpolation which would fit a
curve with a saddlepoint to the square
and draw one connected contour. However, bilinear interpolation is a
bit of a pain, because the contours are not lines, but pieces of a
hyperbola, which is not nearly as easy to draw. And if the contour
level is close enough to 34, you'll get the peak effect anyway.
In any case, this is not a bug, it's an algorithm with shortcomings.
-- O__ ---- Peter Dalgaard Blegdamsvej 3 c/ /'_ --- Dept. of Biostatistics 2200 Cph. N (*) \(*) -- University of Copenhagen Denmark Ph: (+45) 35327918 ~~~~~~~~~~ - (firstname.lastname@example.org) FAX: (+45) 35327907 -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-devel mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info", "help", or "[un]subscribe" (in the "body", not the subject !) To: email@example.com _._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._
This archive was generated by hypermail 2b25 : Tue 04 Jan 2000 - 14:16:10 EST