Thanks Bill,
Yes in my data too I get some different posteriors also.  I guess what I was
getting at was that when I'm dealing with multi-elemental data with 16
variables a probability of 1 is a bit more "certain" an assignment than I
was expecting.  I'm a student, so perhaps I don't have quite the handle on
the theory that I should.  I'm attempting to explain some of the
classifications in my write-up so was just looking for some discussion
really.

Thanks for your help.


Gareth Campbell
PhD Candidate
The University of Auckland

P    +649 815 3670
M    +6421 256 3511
E   gareth.campb...@esr.cri.nz
    gcam...@gmail.com


2009/4/12 <bill.venab...@csiro.au>

> If you try
>
> library(MASS)
> example(lda)
> predict(z)$posterior
>
> You will see there are very many posterior probabilities there, too, which
> differe from 1 by round-off error at most.  So your data are not unique.
>
> I don't get why you don't get it.
>
>
> Bill Venables
> http://www.cmis.csiro.au/bill.venables/
>
>
> -----Original Message-----
> From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org]
> On Behalf Of gcam032
> Sent: Sunday, 12 April 2009 11:16 AM
> To: r-help@r-project.org
> Subject: [R] predict.lda posteriors equal to 1
>
>
> Hi all,
>
> I was wondering if anyone was able to explain to me why when I run
> predict.lda() on some new data I get lots of posterior probabilities = 1.
> This seems markedly unrealistic.  What do the posteriors mean in the
> context
> of the predict.lda() function?
>
> Thanks
>
> Gareth
> --
> View this message in context:
> http://www.nabble.com/predict.lda-posteriors-equal-to-1-tp23007165p23007165.html
> Sent from the R help mailing list archive at Nabble.com.
>
> ______________________________________________
> 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.
>

        [[alternative HTML version deleted]]

______________________________________________
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