Hello,

I have a file with two dependent variables (three and five) and one
independent variable. I do  i.mod <- lm(cbind(three, five) ~ species,
data=i.txt) and get the following output:


Coefficients:
             three   five
(Intercept)   9.949   9.586
species      -1.166  -1.156

I do a" i.can<-candisc(i.mod,data=i):

and get the following output:

Canonical Discriminant Analysis for species:

    CanRsq Eigenvalue Difference Percent Cumulative
1 0.096506    0.10681                100        100

Test of H0: The canonical correlations in the
current row and all that follow are zero
LR test stat approx F num Df den Df   Pr(> F)
1        0.903   63.875      1    598 6.859e-15 ***
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

this is different than the output I get with SAS:

         Eigenvalue Difference Proportion Cumulative      Ratio     F Value
Num DF Den DF Pr > F

       1     0.1068                1.0000     1.0000 0.90349416
31.88      2    597 <.0001


I am also wondering how to plot the can1*can1 like it is done in SAS.

proc plot;
    plot can1*can1=species;
    format species spechar.;
    title2 'Plot of Constits_vs_cassettes';
 run;

Thanks

	[[alternative HTML version deleted]]

