[R] varclus function

Peter Malewski malewski.peter at mh-hannover.de
Fri Aug 2 15:19:15 CEST 2002


Hmisc. 
Peter


On Fri, Aug 02, 2002 at 01:47:08PM +0100, Huan Huang wrote:
> Hi, dear list
> 
> I am fitting a linear model for 50 variables. I know there are collinear
> predictors. So I want to use function varclus. Could any one tell me in
> which library I can find it? Thanks a lot.
> 
> Huan
> ----- Original Message -----
> From: "Shravan Vasishth" <vasishth at coli.uni-sb.de>
> To: <r-help at stat.math.ethz.ch>
> Sent: Tuesday, July 30, 2002 11:15 AM
> Subject: RE: [R] Barplot coloring question
> 
> 
> > I forgot to say in the summary below how mybarplot2 is called:
> >
> > mybarplot2(data2,        # or data3
> >         beside=FALSE,
> >         col= colormatrix,
> >         width = c(1,.55),
> >         space = c(2.75,.25),
> >         axes = F)
> >
> > --
> > Dr. Shravan Vasishth                           Phone: +49 (681) 302 4504
> > Computational Linguistics, Universit舩 des Saarlandes, Postfach 15 11 50
> > D-66041 Saarbr・瘢雹ken, Germany         http://www.coli.uni-sb.de/~vasishth
> >
> >
> > ---------- Forwarded message ----------
> > Date: Tue, 30 Jul 2002 11:33:08 +0200 (MEST)
> > From: Shravan Vasishth <vasishth at gnome.at.coli.uni-sb.de>
> > To: Marc Schwartz <mschwartz at medanalytics.com>
> > Cc: r-help at stat.math.ethz.ch
> > Subject: RE: [R] Barplot coloring question
> >
> > Hi all,
> >
> > I summarize below the eventual solution I used following the discussion in
> > this thread, in case others are interested. This also answers Marc's
> > question below (why the inner loop?).
> >
> > The problem:
> >
> > I need to build two kinds of plots automatically; I need to keep things as
> > automated as possible, since a lot of plots will be generated in the
> > future and we just want to see the output quickly without having to fiddle
> > with the plots individually.
> >
> > The first kind's initial data looks like this:
> >
> > > data2
> >   A  A  B  B  C  C
> > 0 5 10 15 20 70 90
> >
> > The second kind's data looks like this:
> >
> > > data3
> >    A  A  B  B  C  C
> > 1  5 10 15 20 70 90
> > 2 20 15 24 56 34 70
> > 3 40 54 23 67 12 78
> >
> > For uniformity, we want to keep the format of the data identical in both
> > cases.
> >
> > In both cases, I want the distance between the bars A-A, B-B, C-C, to be
> > closer to each other than the A-B, B-C distance, and I want to alternate
> > colors, so that the first A has color x and the second A has a different
> > color y. The x, y colors must alternate after that: the first B has x,
> > then the second B has y, etc.
> >
> > However, in data3, I also want the three rows for a given column to appear
> > *within* a single bar, in three colors, say, c1, c2, c3, and the column
> > next to this column should have colors c4, c5, c6. So, now we want the
> > triplets <c1, c2, c3>, <c4,c5,c6> alternating for each successive bar. The
> > reason for this is that we want to display two related but different
> > percentages in the same graph, and we want to do this because then we can
> > quickly eyeball the graph and get maximum information from a single graph,
> > rather than having to pore over dozens of them and to compare bars
> > pairwise across different graphs (a sure way to lose your mind).
> >
> > Solution:
> >
> > (for the reasons why this should be so, see Marc's previous responses)
> >
> > For data2, barplot is rewritten as mybarplot1 (thanks to Marc) as follows:
> >
> >  ...
> >
> >  if (beside)
> >             xyrect(0, w.l, c(height), w.r, horizontal = horiz,
> >                 col = col)
> >         else {
> >             for (i in 1:NC) {
> >                 xyrect(height[1:NR, i], w.l[i], height[-1, i],
> >                   w.r[i], horizontal = horiz, col = col[i])
> >             }
> >
> >  ...
> >
> >
> > For data3, barplot is rewritten as mybarplot2 (again, thanks to Marc) as
> > follows:
> >
> >    ...
> >
> >    if (beside)
> >             xyrect(0, w.l, c(height), w.r, horizontal = horiz,
> >                 col = col)
> >         else {
> >             for (i in 1:NC) {
> >                      for (j in 1:NR) {
> >                        mycol = col[,i]
> >                        xyrect(height[1:NR, i], w.l[i], height[-1, i],
> >                             w.r[i], horizontal = horiz, col = mycol)
> >                      }
> >             }
> >
> >    ...
> >
> >
> > Now, one nice thing is that we can just use the above modification, i.e.,
> > mybarplot2, for both data2 and data3, we don't need mybarplot1 any more.
> > Not a real big deal, but now we can use the same R script for all our
> > data, irrespective of whether it comes in in data2 or data3 format. Just
> > might make things easier for other people who will one day have to
> > read and/or modify this and related code at our end here.
> >
> > On Mon, 29 Jul 2002, Marc Schwartz wrote:
> >
> > > The question this whole process now raises is why the inner loop?
> >
> > Does the above explanation answer that question?
> >
> > > I am still a little confused as to how you want the colors structured.
> > > That is do you want each bar to contain the three colors, one for each
> > > segment, or do you want each of the three pairs of bars to be colored
> > > differently?
> >
> > I want them to alternate, for data3, as follows:
> >
> > c3 c6    c3 c6
> > c2 c5    c2 c5
> > c1 c4    c1 c4   ...
> > A  A     B  B
> >
> > Marc, I want to thank you again for the enormous amount of time you put
> > into this, and for your extraordinary patience as you educated me on
> > these revisions to the code. It would have taken much, much longer to
> > figure all this out without your help.
> >
> > --
> > Dr. Shravan Vasishth                           Phone: +49 (681) 302 4504
> > Computational Linguistics, Universit舩 des Saarlandes, Postfach 15 11 50
> > D-66041 Saarbr・瘢雹ken, Germany         http://www.coli.uni-sb.de/~vasishth
> >
> >
> >
> > -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.
> -.-.-
> > r-help mailing list -- Read
> http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html
> > Send "info", "help", or "[un]subscribe"
> > (in the "body", not the subject !)  To: r-help-request at stat.math.ethz.ch
> >
> _._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._.
> _._
> >
> 
> -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
> r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html
> Send "info", "help", or "[un]subscribe"
> (in the "body", not the subject !)  To: r-help-request at stat.math.ethz.ch
> _._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._

-- 
P.Malewski, Limmerstr.47, 30451 Hannover, 0511-2135008
At work: http://www.MH-Hannover.de 0511 532 3194 / Fax: 0511 532 3190, 
P.Malewski at tu-bs.de, peter.malewski at gmx.de, malewski.peter at mh-hannover.de.
-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html
Send "info", "help", or "[un]subscribe"
(in the "body", not the subject !)  To: r-help-request at stat.math.ethz.ch
_._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._



More information about the R-help mailing list