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]]