[GRASS-stats] Superimpose density plots in R

[Jasiewicz Jarosław]

Nikos Alexandris pisze:
> Hi all!
>
> I started playing around with R and GRASS-maps.
>
> I've a couple or raster maps (MODIS bands) in R and what I am trying to
> achieve is to get density plots superimposed in one graph (for the sake
> of easy visual comparison).
>
> The only function that I have found and is able to perform this is
> "sm.density.compare" ( belongs to the "sm" package). Although this
> function works with *other* small (in terms of number of rows) examples,
> for my MODIS bands (~87.000 rows each) I get an empty graph.
>
> Any pointers for this "comparison" task?
>
> Thank you, Nikos
>
>
>
> _______________________________________________
> grass-stats mailing list
> grass-stats at lists.osgeo.org
> http://lists.osgeo.org/mailman/listinfo/grass-stats
>   
Hi
simply use density function with the same bandwith=, from=, to= for each 
raster and next use lines to add curves for vuisual comparing

for example:
plot(density(bronze$a_slope,bw=0.1,from=0, to=12), col=2)
lines(density(neolitic$a_slope,bw=0.1,from=0, to=12), col=1)

Jarek