lda output missing
That's odd. You don't provide a reproducible example, but using
a built-in dataset (from the help for lda) I get the Proportion of Trace
given by the print.lda method.
library(MASS)
Iris <- data.frame(rbind(iris3[,,1], iris3[,,2], iris3[,,3]), Sp =
rep(c("s","c","v"), rep(50,3)))
train <- sample(1:150, 75)
z <- lda(Sp ~ ., Iris, prior = c(1,1,1)/3, subset = train)
print(z)
Call:
lda(Sp ~ ., data = Iris, prior = c(1, 1, 1)/3, subset = train)
Prior probabilities of groups:
c s v
0.3333333 0.3333333 0.3333333
Group means:
Sepal.L. Sepal.W. Petal.L. Petal.W.
c 5.947826 2.786957 4.326087 1.3652174
s 5.050000 3.380769 1.465385 0.2346154
v 6.811538 3.023077 5.700000 2.0653846
Coefficients of linear discriminants:
LD1 LD2
Sepal.L. -1.6665885 1.104460
Sepal.W. -0.8681909 1.872188
Petal.L. 2.4670375 -1.499070
Petal.W. 4.2156275 2.948759
Proportion of trace:
LD1 LD2
0.9916 0.0084
Here's my sessionInfo(). Are you certain everything is up-to-date?
What OS are you running? Do you have other packages loaded that could
be interfering? What happens when you run the example I used?
sessionInfo()
R version 2.14.0 (2011-10-31) Platform: x86_64-redhat-linux-gnu (64-bit) locale: [1] LC_CTYPE=en_US.UTF-8 LC_NUMERIC=C [3] LC_TIME=en_US.UTF-8 LC_COLLATE=en_US.UTF-8 [5] LC_MONETARY=en_US.UTF-8 LC_MESSAGES=en_US.UTF-8 [7] LC_PAPER=C LC_NAME=C [9] LC_ADDRESS=C LC_TELEPHONE=C [11] LC_MEASUREMENT=en_US.UTF-8 LC_IDENTIFICATION=C attached base packages: [1] stats graphics grDevices utils datasets methods base other attached packages: [1] MASS_7.3-16 loaded via a namespace (and not attached): [1] tools_2.14.0
On Thu, Dec 8, 2011 at 10:29 AM, David Lutz <lutzda at wfu.edu> wrote:
Hello everyone, I am working on a linear discriminant analysis and am having issues finding the full output of my lda. Specifically, there is no reporting of the Proportion of Trace that is a normal output of the procedure. I'm using a csv file and everything is reading in correctly. I've looked and looked and can't figure out why my output is not complete. Is it something simple that I am overlooking? I am working in R 2.14.0 with an updated MASS package. Here is my syntax: lda(CC ~ sumstem + maxdbh + maxh + cdlai + sumlai + sumbas + ? ?totalbio + ratio, data = practice) and the output is: Prior probabilities of groups: normal ? warm ? 0.5 ? ?0.5 Group means: ? ? ? ?sumstem ? maxdbh ? ? maxh ? ? ? cdlai ? sumlai ? sumbas totalbio normal 901.3722 16.66994 15.70927 0.010393258 1.672247 5.471812 27.67875 warm ? 972.4916 22.27247 19.64740 0.002429775 3.181994 9.892683 47.68511 ? ? ? ? ? ?ratio normal 0.01162921 warm ? 0.85721910 Coefficients of linear discriminants: ? ? ? ? ? ? ? ? ?LD1 sumstem ? 0.001111176 maxdbh ? ?0.034299258 maxh ? ? -0.287343783 cdlai ? ?-2.021350057 sumlai ? ?1.407215702 sumbas ? -0.419422181 totalbio ?0.110403369 ratio ? ? 0.001657649
Sarah Goslee http://www.functionaldiversity.org