[R] Re: projection pursuit
Liaw, Andy
andy_liaw at merck.com
Thu Mar 18 17:13:25 CET 2004
> From: Jim_Garrett at bd.com
>
> Luis,
>
> See the fastICA package, in particular the final example in
> the function
> fastICA's help page. This doesn't leave you with density
> estimates, but
> with projection-pursuit directions; you still have to figure
> out how to fit
> a density estimate to the rotated data. Actually, as I understand it,
> XGobi also finds directions but does not fit a density. However, any
> multivariate density estimator ought to be applicable. I'm
> not aware at
> the moment of the tools R offers for multivariate density
> estimation, but
> I'm sure there are multiple possibilities.
I believe locfit() can handle up to 5 dimensions.
> _The Elements of Statistical Learning_ by Hastie, Tibshirani,
> and Friedman
> mention a "trick" to use classification tools (which model class
> probabilities) to estimate density. Fundamentally, generate
> data from a
> reference distribution, and use the classification tool to estimate
> probability of observed data (as opposed to generated data)
> as a function
> of the inputs. These probabilities, normalized to integrate
> to 1, form a
> density estimate. Since there are so very many classification tools
> available, this trick offers a lot of flexibility.
Leo Breiman had talked about this `trick' in his early presentations on
random forest. However, from my limited experiment, this is an extremely
poor way of estimating density. Has anyone else tried it?
Cheers,
Andy
> Good luck.
>
> Jim Garrett
> Baltimore, Maryland, USA
>
>
> **********************************************************************
> This message is intended only for the designated
> recipient(s...{{dropped}}
>
> ______________________________________________
> R-help at stat.math.ethz.ch mailing list
> https://www.stat.math.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide!
> http://www.R-project.org/posting-guide.html
>
>
------------------------------------------------------------------------------
Notice: This e-mail message, together with any attachments,...{{dropped}}
More information about the R-help
mailing list