Re: [R] Quelplot

2011-09-22 Thread Den Alpin
Hi Hadley,
have a look at:

http://www.riani.it/pub/zrc-csda98.pdf
and some Gauss code:
http://www.riani.it/Gauss/procedures/BOXPLOTB.G

Best regards,
Daniele


2011/9/21 Hadley Wickham had...@rice.edu:
 Hi all,

 Does anyone have an R implementation of the queplot (K. M. Goldberg
 and B. Iglewicz. Bivariate extensions of the boxplot. Technometrics,
 34(3):pp. 307–320, 1992)?  I'm struggling with the estimation of the
 asymmetry parameters.

 Hadley

 --
 Assistant Professor / Dobelman Family Junior Chair
 Department of Statistics / Rice University
 http://had.co.nz/

 __
 R-help@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.


__
R-help@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.


Re: [R] Quelplot

2011-09-22 Thread Richard M. Heiberger
Hadley,

I forwarded your request to Boris and Ken.  Here is Boris' response.

Rich


Dear Richard:

   I looked a bit more and found the second edition of book referenced, Rand
R Wilcox, Introduction to Robust Estimation and Hypothesis Testing, Second
Edition, Academic Press, 2005. It has over 600 quotes on google scholar. If
my memory is correct, Wilcox developed code for relplots and possibly
quelplots.

Best Wishes,
Boris

On Wed, Sep 21, 2011 at 3:11 PM, Hadley Wickham had...@rice.edu wrote:

 Hi all,

 Does anyone have an R implementation of the queplot (K. M. Goldberg
 and B. Iglewicz. Bivariate extensions of the boxplot. Technometrics,
 34(3):pp. 307–320, 1992)?  I'm struggling with the estimation of the
 asymmetry parameters.

 Hadley

 --
 Assistant Professor / Dobelman Family Junior Chair
 Department of Statistics / Rice University
 http://had.co.nz/

 __
 R-help@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.


[[alternative HTML version deleted]]

__
R-help@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.


Re: [R] Quelplot

2011-09-22 Thread Hadley Wickham
Thanks Rich - unfortunately, as far as I can, tell Rand Wilcox only
developed relplot code.
Hadley

On Thu, Sep 22, 2011 at 12:52 PM, Richard M. Heiberger r...@temple.edu wrote:
 Hadley,

 I forwarded your request to Boris and Ken.  Here is Boris' response.

 Rich


 Dear Richard:

    I looked a bit more and found the second edition of book referenced, Rand
 R Wilcox, Introduction to Robust Estimation and Hypothesis Testing, Second
 Edition, Academic Press, 2005. It has over 600 quotes on google scholar. If
 my memory is correct, Wilcox developed code for relplots and possibly
 quelplots.

 Best Wishes,
 Boris

 On Wed, Sep 21, 2011 at 3:11 PM, Hadley Wickham had...@rice.edu wrote:

 Hi all,

 Does anyone have an R implementation of the queplot (K. M. Goldberg
 and B. Iglewicz. Bivariate extensions of the boxplot. Technometrics,
 34(3):pp. 307–320, 1992)?  I'm struggling with the estimation of the
 asymmetry parameters.

 Hadley

 --
 Assistant Professor / Dobelman Family Junior Chair
 Department of Statistics / Rice University
 http://had.co.nz/

 __
 R-help@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.





-- 
Assistant Professor / Dobelman Family Junior Chair
Department of Statistics / Rice University
http://had.co.nz/

__
R-help@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.


Re: [R] Quelplot

2011-09-22 Thread Michael Friendly

On 9/21/2011 3:11 PM, Hadley Wickham wrote:

Hi all,

Does anyone have an R implementation of the queplot (K. M. Goldberg
and B. Iglewicz. Bivariate extensions of the boxplot. Technometrics,
34(3):pp. 307–320, 1992)?  I'm struggling with the estimation of the
asymmetry parameters.


Do you have some code you have tried or don't you know how to calculate 
P_1 and P_2?  The description is not particularly explicit, but the

major/minor axes of the ellipse are defined in terms of standard scores
for X and Y as
Z1 = Y + X
Z2 = Y - X

If you can figure out what they mean by proportion of the total 
standard deviation due to residuals..., then something like


P_1 = sum(Z_10) / n
P_2 = sum(Z_20) / n

Maybe Boris can help.

Thanks for bringing this to our attention.  If the bagplot is a totally
non-parametric analog of a bivariate boxplot, the quelplot might
be a 2df parametric version.

-Michael



--
Michael Friendly Email: friendly AT yorku DOT ca
Professor, Psychology Dept.
York University  Voice: 416 736-5115 x66249 Fax: 416 736-5814
4700 Keele StreetWeb:   http://www.datavis.ca
Toronto, ONT  M3J 1P3 CANADA

__
R-help@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.


Re: [R] Quelplot

2011-09-22 Thread Goldberg, Ken [JRDUS]
Hi Hadley,

Thank-you for your interest in the quelplot. I agree that Rand Wilcox only 
published relplot code. I do not know of any other R or S implementations. I 
wrote all of my code in HTBasic (Rocky Mountain BASIC) originally. I later 
copied some of it to SAS proc iml, but I think that was just the multivariate 
analog of the relplot (including the bivariate relplot). 

I wonder if there is a typo in the Technometrics article, because the quelplot 
always worked fine for me but I recall somebody else having trouble with 
programming it. I think I sent that person my HTBasic code back in the early 
1990's. Perhaps that was Rand. I think I did send it to somebody who asked for 
it in Australia.

I will have to search for my code. I know I have it on my PC in HTBasic format, 
but it has been a long time since I used it and may have to work a bit even to 
print it or to get it out into a DOS text file. I may have a DOS text file of 
it at home on a 3.5 diskette, but even if I could find it, I doubt if I could 
find a disk drive to read it! When I get it in a useable format, I can send you 
that code, or at least the parts of it relevant to the quelplot, if you would 
like it. Perhaps just browsing it would show the problem you are having with my 
Technometrics article, whether it was a typo or just something that I wrote 
poorly. There may just be a few subroutines that would need translating from 
HTBasic into R, and that when added to Rand Wilcox's library would give his 
software the quelplot capacity.

Hadley, I took your R/Splus Fundamentals and Programming Techniques course 5 
years ago in Washington DC, but have used R just a couple of times since. I've 
rarely programmed in R, but this would be a good project for me to start with. 
It may be a few weeks though until I can get started on any of this though, but 
I will keep you informed.

Regards, 
 Ken 
 
Kenneth M. Goldberg, Ph.D. 
Associate Director, Nonclinical Statistics 
Johnson  Johnson Pharmaceutical Research  Development, LLC
965 Chesterbrook Boulevard, Mailstop C-4-1
Wayne, PA 19087
Tel  610-651-6924   Fax 610-651-6717
kgold...@its.jnj.com

Confidentiality Notice:  This e-mail transmission may contain confidential or 
legally privileged information that is intended only for the individual or 
entity named in the e-mail address.  If you are not the intended recipient, you 
are hereby notified that any disclosure, copying, distribution, or reliance 
upon the contents of this e-mail is strictly prohibited. If you have received 
this e-mail transmission in error, please reply to the sender, so that Centocor 
can arrange for proper delivery, and then please delete the message from your 
inbox.  Thank you.






Does anyone have an R implementation of the queplot (K. M. Goldberg
 and B. Iglewicz. Bivariate extensions of the boxplot. Technometrics,
 34(3):pp. 307-320, 1992)?  I'm struggling with the estimation of the
 asymmetry parameters.

-Original Message-
From: h.wick...@gmail.com [mailto:h.wick...@gmail.com] On Behalf Of Hadley 
Wickham
Sent: Thursday, September 22, 2011 1:56 PM
To: Richard M. Heiberger
Cc: R-help; BORIS IGLEWICZ; Goldberg, Ken [JRDUS]
Subject: Re: [R] Quelplot

Thanks Rich - unfortunately, as far as I can, tell Rand Wilcox only
developed relplot code.
Hadley

On Thu, Sep 22, 2011 at 12:52 PM, Richard M. Heiberger r...@temple.edu wrote:
 Hadley,

 I forwarded your request to Boris and Ken.  Here is Boris' response.

 Rich


 Dear Richard:

    I looked a bit more and found the second edition of book referenced, Rand
 R Wilcox, Introduction to Robust Estimation and Hypothesis Testing, Second
 Edition, Academic Press, 2005. It has over 600 quotes on google scholar. If
 my memory is correct, Wilcox developed code for relplots and possibly
 quelplots.

 Best Wishes,
 Boris

 On Wed, Sep 21, 2011 at 3:11 PM, Hadley Wickham had...@rice.edu wrote:

 Hi all,

 Does anyone have an R implementation of the queplot (K. M. Goldberg
 and B. Iglewicz. Bivariate extensions of the boxplot. Technometrics,
 34(3):pp. 307-320, 1992)?  I'm struggling with the estimation of the
 asymmetry parameters.

 Hadley

 --
 Assistant Professor / Dobelman Family Junior Chair
 Department of Statistics / Rice University
 http://had.co.nz/

 __
 R-help@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.





-- 
Assistant Professor / Dobelman Family Junior Chair
Department of Statistics / Rice University
http://had.co.nz/

__
R-help@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

[R] Quelplot

2011-09-21 Thread Hadley Wickham
Hi all,

Does anyone have an R implementation of the queplot (K. M. Goldberg
and B. Iglewicz. Bivariate extensions of the boxplot. Technometrics,
34(3):pp. 307–320, 1992)?  I'm struggling with the estimation of the
asymmetry parameters.

Hadley

-- 
Assistant Professor / Dobelman Family Junior Chair
Department of Statistics / Rice University
http://had.co.nz/

__
R-help@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.