[R] empty density plot for point pattern in spatstat

[Adrian Baddeley]

Dear Pablo,

> Yes, effectively utilizando useRaster = FALSE, the plot is printed. 

Muy bien!
 
> Only the ribbon appears without color.

Try setting ribargs=list(useRaster=TRUE).

> I'm using R version 3.0.2, Windows 7 system, and spatstat 1.33-0 package version.

Things may improve if you upgrade to the current versions, R 3.1.0 and spatstat 1.37-0.

A

Prof Adrian Baddeley FAA
University of Western Australia
________________________________________
From: Pablo Ramón [paramon at utpl.edu.ec]
Sent: Thursday, 19 June 2014 9:12 PM
To: Adrian Baddeley
Cc: Rolf Turner; Christopher W Ryan; R-help
Subject: Re: [R] empty density plot for point pattern in spatstat

Dear Adrian,
Yes, effectively utilizanddo useRaster = FALSE, the plot is printed. Only the ribbon appears without color.
I'm using R version 3.0.2, Windows 7 system, and spatstat 1.33-0 package version.

Best wishes
Pablo



2014-06-18 19:50 GMT-05:00 Adrian Baddeley <adrian.baddeley at uwa.edu.au<mailto:adrian.baddeley at uwa.edu.au>>:
Dear Pablo

Please provide details of the system you are using (Windows - version?  R version? spatstat version?)

Also, when you issue the plot command, what graphics device are you using?
[To find out, type dev.cur() after performing the plot, and read the printout]

Try changing the colour map by running 'example(plot.im<http://plot.im>)' or
setting col=grey(seq(0,1,length=32)) in the call to plot.im<http://plot.im>

Try setting 'useRaster=FALSE' in the call to plot.im<http://plot.im>.

Let me know if any of these work. If not, let's talk on skype
to figure it out.

regards
Adrian


Prof Adrian Baddeley FAA
University of Western Australia
________________________________________
From: Pablo Ramón [paramon at utpl.edu.ec<mailto:paramon at utpl.edu.ec>]
Sent: Wednesday, 18 June 2014 9:36 PM
To: Rolf Turner
Cc: Christopher W Ryan; R-help; Adrian Baddeley
Subject: Re: [R] empty density plot for point pattern in spatstat

Dear Rolf,

The problem persists with the function plot.im<http://plot.im><http://plot.im>(), or plot (density ()), similar to with(bei.extra, plot (elev)). This problem occurs only with my personal machine, to another machine this problem does not occur. However, when I use another function, such kde2d () in the MASS package and perform the plot, the graph is obtained without problem.

Best wishes
Pablo





2014-06-17 16:23 GMT-05:00 Rolf Turner <r.turner at auckland.ac.nz<mailto:r.turner at auckland.ac.nz><mailto:r.turner at auckland.ac.nz<mailto:r.turner at auckland.ac.nz>>>:
On 18/06/14 03:49, Christopher W Ryan wrote:
I upgraded to version 3.1.0 of R, still on Win XP, but that did not
solve the problem:

library(spatstat)
data(redwood)
str(redwood)
plot(redwood)  # works OK
plot(Kest(redwood))  # works OK
plot(density(redwood))  # produces an "empty" plot
                                        # with just the figure border
                                        # and a thin vertical box for
the color ramp
                                        # but no color ramp in it

Any other ideas? I'd be much obliged.

I have had an inquiry of this nature previously and sad to say was completely unable to find a solution.  I made numerous suggestions to the inquirer, e.g.:

*looking for a plot.im<http://plot.im><http://plot.im>() function masking the "real"
plot.im<http://plot.im><http://plot.im>() function from spatstat

* cleaning out the global environment
(.RData) in case there were some "artefacts" hanging around that had an impact on plot.im<http://plot.im><http://plot.im>()

* completely *removing* .RData and starting again,...

Nothing worked.

It is very difficult to diagnose a problem like this remotely.  If you could find an R guru who was able to come and sit with you at your computer it is ***just possible*** that they could track down what is going wrong.

However it seems to me that the problem is due to some sort of recalcitrance of the Windoze system and is thus probably insoluble.

I am cc-ing this to the person who made the previous inquiry.  Pablo, have there been any further developments since we last communicated about this problem?

cheers,

Rolf Turner

P. S. BTW Chris:  Does plot.im<http://plot.im><http://plot.im>() screw up with other images?  E.g. what happens if you do

with(bei.extra,plot(elev))

???

R. T.




--
Mat. Pablo Ramón
Sección de Ecología
Departamento de Ciencias Naturales
Universidad Técnica Particular de Loja (http://www.utpl.edu.ec/)
Telf. 593 7 2570275 Ext. 2505, 2258



--
Mat. Pablo Ramón
Sección de Ecología
Departamento de Ciencias Naturales
Universidad Técnica Particular de Loja (http://www.utpl.edu.ec/)
Telf. 593 7 2570275 Ext. 2505, 2258