On Sun, 26 Apr 2009, Uwe Ligges wrote:

Works for me with

library(MASS)
plot(lda(Species~., data=iris))

hence you may want to profide the data to enable us to reproduce your problem...

He is trying to plot the results from a cross-validation. As the help page clearly states, that is a list (and has no assigned class). Plotting an arbitrary list makes little sense (and not does plotting the results of an LDA cross-validation).


Uwe Ligges



pgseye wrote:
Hi,

I've performed an lda and obtained a classification table for some of my
data:

efa.dfa<-lda(groups~.,efa.scores.8,CV=T)
str(efa.dfa)
List of 5
 $ class    : Factor w/ 2 levels "1","2": 1 2 1 2 1 1 2 2 1 2 ...
 $ posterior: num [1:160, 1:2] 0.99083 0.00852 0.93983 0.23186 0.85931 ...
  ..- attr(*, "dimnames")=List of 2
  .. ..$ : chr [1:160] "1" "2" "3" "4" ...
  .. ..$ : chr [1:2] "1" "2"
 $ terms    :Classes 'terms', 'formula' length 3 groups ~ Comp.1 + Comp.2 +
Comp.3 + Comp.4 + Comp.5 + Comp.6 +      Comp.7 + Comp.8 + Comp.9 + Comp.10
+ Comp.11 + Comp.12 +  ...
.. ..- attr(*, "variables")= language list(groups, Comp.1, Comp.2, Comp.3,
Comp.4, Comp.5, Comp.6,      Comp.7, Comp.8, Comp.9, Comp.10, Comp.11,
Comp.12, Comp.13,  ...
  .. ..- attr(*, "factors")= int [1:35, 1:34] 0 1 0 0 0 0 0 0 0 0 ...
  .. .. ..- attr(*, "dimnames")=List of 2
  .. .. .. ..$ : chr [1:35] "groups" "Comp.1" "Comp.2" "Comp.3" ...
  .. .. .. ..$ : chr [1:34] "Comp.1" "Comp.2" "Comp.3" "Comp.4" ...
  .. ..- attr(*, "term.labels")= chr [1:34] "Comp.1" "Comp.2" "Comp.3"
"Comp.4" ...
  .. ..- attr(*, "order")= int [1:34] 1 1 1 1 1 1 1 1 1 1 ...
  .. ..- attr(*, "intercept")= int 1
  .. ..- attr(*, "response")= int 1
.. ..- attr(*, ".Environment")=<environment: R_GlobalEnv> .. ..- attr(*, "predvars")= language list(groups, Comp.1, Comp.2, Comp.3,
Comp.4, Comp.5, Comp.6,      Comp.7, Comp.8, Comp.9, Comp.10, Comp.11,
Comp.12, Comp.13,  ...
  .. ..- attr(*, "dataClasses")= Named chr [1:35] "numeric" "numeric"
"numeric" "numeric" ...
.. .. ..- attr(*, "names")= chr [1:35] "groups" "Comp.1" "Comp.2" "Comp.3"
...
 $ call     : language lda(formula = groups ~ ., data = efa.scores.8, CV =
T)
 $ xlevels  : list()
table(groups, Classified=efa.dfa$class)
      Classified
groups  1  2
     1 59 21
     2 10 70

but when I try to plot the results I get:

plot(efa.dfa)
Error in plot.window(...) : need finite 'xlim' values
In addition: Warning messages:
1: In min(x) : no non-missing arguments to min; returning Inf
2: In max(x) : no non-missing arguments to max; returning -Inf
3: In min(x) : no non-missing arguments to min; returning Inf
4: In max(x) : no non-missing arguments to max; returning -Inf

anyone have any ideas?

Thanks a lot,

Paul

______________________________________________
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.


--
Brian D. Ripley,                  rip...@stats.ox.ac.uk
Professor of Applied Statistics,  http://www.stats.ox.ac.uk/~ripley/
University of Oxford,             Tel:  +44 1865 272861 (self)
1 South Parks Road,                     +44 1865 272866 (PA)
Oxford OX1 3TG, UK                Fax:  +44 1865 272595

______________________________________________
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.

Reply via email to