hypsometric integral from ecdf curve

[Pierluigi De Rosa]

Some javascript/style in this post has been disabled (why?)

Dear All,

I need to evaluate an Hypsometric integral but i should "rotate" the graphical output of

plot(ecdf(elev$elev))

I want in x the fraction of observations (pixel value of elevation) that are above the elevation in y

the  command: plot(ecdf(elev$elev)) show in x the elevation and in y the fraction of observation that are below the value of elevation in x

The value of 1 in plot(ecdf(elev$elev)) correspond to the maximum value of elevation. I want the it correspond to the value of minimum.

thanks

Pierluigi
_______________________________________________
grass-stats mailing list
[hidden email]
http://lists.osgeo.org/mailman/listinfo/grass-stats

[Roger Bivand]

On Wed, 28 Jan 2009, pierluigi de rosa wrote:

> Dear All,
>
> I need to evaluate an Hypsometric integral but i should "rotate" the
> graphical output of
>
> plot(ecdf(elev$elev))
>
> I want in x the fraction of observations (pixel value of elevation) that
> are above the elevation in y
>
> the  command: plot(ecdf(elev$elev)) show in x the elevation and in y the
> fraction of observation that are below the value of elevation in x
>
> The value of 1 in plot(ecdf(elev$elev)) correspond to the maximum value
> of elevation. I want the it correspond to the value of minimum.

It looks as though:

res <- plot.stepfun(ecdf(elev$elev))
plot(res$y[-1], res$t[-c(1, length(res$t))], xlim=c(1, 0), type="l",
ylab="elev", xlab="F(x)")

gets you there. I can't see an argument in the first step to turn plotting
off.

Best wishes,

Roger

>
> thanks
>
> Pierluigi
>

--
Roger Bivand
Economic Geography Section, Department of Economics, Norwegian School of
Economics and Business Administration, Helleveien 30, N-5045 Bergen,
Norway. voice: [phone]; fax [phone]
e-mail: [hidden email]

_______________________________________________
grass-stats mailing list
[hidden email]
http://lists.osgeo.org/mailman/listinfo/grass-stats