Re: [R] Finding Lower Envelope of Points on a Plot?

From: Ted Harding <Ted.Harding_at_manchester.ac.uk>
Date: Wed, 30 Jun 2010 22:26:57 +0100 (BST)


I think this has reached the point where it is absolutely necessary for Asha Sharma to state *exactly* what is meant by a "linear lower envelope" of a set of points. There are too many possibilities for the problem as stated so far. For instance (possibilities which have not been considered so far):

[1] A straight line which is such that any point which falls below it is "censored" and so would not be present (but this would seem to be inconsistent with Asha's proviso of "a specified percentage of points allowed to lie outside it");

[2] A line parallel to the regression line, but displaced downwards in such a way as to exclude a ceretain percantage of ths points (i.e. such that that percentage lie below the line); but this does not match my own usual interpretation of "envelope" (which is such that it is a boundary for the positions of the points, as in [1]).

And so on ...
Ted.

On 30-Jun-10 20:58:45, David Winsemius wrote:

> On Jun 30, 2010, at 4:02 PM, Asha Sharma wrote:

>> Hi,
>> Thanks to both of you for taking the time to answer my question. I
>> was maybe not very clear in the way I framed my question. By plot, I
>> meant an x-y plot with a cloud of points which should have a linear
>> lower envelope. Is there a way to both plot as well as get the
>> parameters of the lower envelope (intercept, slope, etc.) and to
>> also set the percentage of points outside it?
> 
> I suspect you are requesting a scatterplot with a superimposed linear  
> fit and 5% and 95% prediction intervals. Not much point in pursuing  
> until you "confirm or deny".
> 
> -- 
> David

>>
>> Thanks!
>>
>> Asha
>>
>>
>> Stephan Kolassa wrote:
>>> Hi,
>>>
>>> one possibility would be to calculate the convex hull using
>>> chull(). I believe that the hull points are returned by chull() in
>>> a clockwise order (?), so the points between the rightmost and the
>>> leftmost point in the chull() result are the lower half of the
>>> convex hull. Remove these points from the original dataset (a
>>> variant of "peeling convex hulls") and iterate until you have
>>> removed your prespecified percentage of points - these all lie
>>> outside of the final lower hull (though the percentage will, of
>>> course, only be approximated, but you should be able to modify this
>>> to taste).
>>>
>>> HTH
>>> Stephan
>>>
>>>
>>> David Winsemius schrieb:
>>>>
>>>> On Jun 30, 2010, at 2:05 PM, Asha Sharma wrote:
>>>>
>>>>> Hi,
>>>>>
>>>>> I am looking for a way to find the lower envelope of points on a
>>>>> plot, preferably specifying what percentage of points should be
>>>>> allowed to lie outside the envelope. There must be a
>>>>> straightforward way to do this, but I do not seem to be able to
>>>>> find it. I would greatly appreciate any help.
>>>>
>>>> You probably want something like the "lower half" of the convex
>>>> hull. You should find quite a bit of code with your favorite r
>>>> search engine on the topic of "convex hull".
>>>
>>
> 
> David Winsemius, MD
> West Hartford, CT

--------------------------------------------------------------------
E-Mail: (Ted Harding) <Ted.Harding_at_manchester.ac.uk> Fax-to-email: +44 (0)870 094 0861
Date: 30-Jun-10                                       Time: 22:26:54
------------------------------ XFMail ------------------------------

______________________________________________
R-help_at_r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. Received on Wed 30 Jun 2010 - 21:31:14 GMT

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.2.0, at Wed 30 Jun 2010 - 21:40:44 GMT.

Mailing list information is available at https://stat.ethz.ch/mailman/listinfo/r-help. Please read the posting guide before posting to the list.

list of date sections of archive