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
 wrote:
> On Fri, Jan 2, 2009 at 4:55 AM, Jarek Jasiewicz  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


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  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 Dylan Beaudette
On Fri, Jan 2, 2009 at 4:55 AM, Jarek Jasiewicz  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


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

2009-01-02 Thread Jarek Jasiewicz

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


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