Re: [GRASS-user] r.quantile: what it exaclty do?

2009-01-02 Thread Dylan Beaudette
On Fri, Jan 2, 2009 at 4:55 AM, Jarek Jasiewicz jar...@amu.edu.pl wrote:
 Hi
 I tirend used r.quantile on SRTM data for Poland (min value about to 0, max
 value about to 2500 m a.s.l)

 r.quantile input=pol...@permanent quantiles=4
 percentiles=0.001,0.01,0,1,0.25,0.50,0.75,0.90,0.99,0.999 bins=100

 and I recived:

 0:0.00:0.00
 1:0.001000:0.00
 2:0.01:0.00
 3:0.25:0.00
 4:0.50:1.00
 5:0.75:2.00
 6:0.90:3.00
 7:0.99:3.00
 8:0.999000:3.00
 9:1.00:3.00

 It looks like the results are divided by 1000 and rounded to the nearest
 integer.

 Did I something wrong?

 Jarek


I thought that the user supplies one of [quantiles=] | [percentiles=]
... could that be related to the odd output. I have verified (with R)
that r.quantile can compute correct quantiles... however, on my system
I need to set bins=100 or so, as the higher i set 'bins' odd things
would happen.

Glynn should know for sure.

Dylan
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


Re: [GRASS-user] r.quantile: what it exaclty do?

2009-01-02 Thread Jarek Jasiewicz



Dylan Beaudette pisze:

On Fri, Jan 2, 2009 at 4:55 AM, Jarek Jasiewicz jar...@amu.edu.pl wrote:
  

Hi
I tirend used r.quantile on SRTM data for Poland (min value about to 0, max
value about to 2500 m a.s.l)

r.quantile input=pol...@permanent quantiles=4
percentiles=0.001,0.01,0,1,0.25,0.50,0.75,0.90,0.99,0.999 bins=100

and I recived:

0:0.00:0.00
1:0.001000:0.00
2:0.01:0.00
3:0.25:0.00
4:0.50:1.00
5:0.75:2.00
6:0.90:3.00
7:0.99:3.00
8:0.999000:3.00
9:1.00:3.00

It looks like the results are divided by 1000 and rounded to the nearest
integer.

Did I something wrong?

Jarek




I thought that the user supplies one of [quantiles=] | [percentiles=]
... could that be related to the odd output. I have verified (with R)
that r.quantile can compute correct quantiles... however, on my system
I need to set bins=100 or so, as the higher i set 'bins' odd things
would happen.

Glynn should know for sure.

Dylan
  

Thanks
it really seems that the problem was in default number of bins
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


Re: [GRASS-user] r.quantile: what it exaclty do?

2009-01-02 Thread Markus Neteler
On Fri, Jan 2, 2009 at 8:40 PM, Dylan Beaudette
dylan.beaude...@gmail.com wrote:
 On Fri, Jan 2, 2009 at 4:55 AM, Jarek Jasiewicz jar...@amu.edu.pl wrote:
 Hi
 I tirend used r.quantile on SRTM data for Poland (min value about to 0, max
 value about to 2500 m a.s.l)

 r.quantile input=pol...@permanent quantiles=4
 percentiles=0.001,0.01,0,1,0.25,0.50,0.75,0.90,0.99,0.999 bins=100

 and I recived:

 0:0.00:0.00
 1:0.001000:0.00
 2:0.01:0.00
 3:0.25:0.00
 4:0.50:1.00
 5:0.75:2.00
 6:0.90:3.00
 7:0.99:3.00
 8:0.999000:3.00
 9:1.00:3.00

 It looks like the results are divided by 1000 and rounded to the nearest
 integer.

 Did I something wrong?

 Jarek


 I thought that the user supplies one of [quantiles=] | [percentiles=]
 ... could that be related to the odd output. I have verified (with R)
 that r.quantile can compute correct quantiles... however, on my system
 I need to set bins=100 or so, as the higher i set 'bins' odd things
 would happen.

 Glynn should know for sure.

Can you please file a bug report on this? Preferable with a
Spearfish/NC reproducible example?

Thanks
Markus
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user