Skip to content

Fitting Hastie's principal surfaces in R

4 messages · Neverstop ., Bert Gunter, David Winsemius +1 more

#
Hello,

I need to summarize a three-dimensional dataset through a principal surface that passes through the middle of the data. Principal surfaces are non-linear generalization of the plane created by the first two principal components and provide a non-linear summary of p-dimensional dataset. Principal surfaces are described in this 1989 article by Hastie and Stuetzle: https://web.stanford.edu/~hastie/Papers/Principal_Curves.pdf . They were introduced by Trevor Hastie in his Ph.D dissertation: http://www.slac.stanford.edu/cgi-wrap/getdoc/slac-r-276.pdf

I'm looking for a package to fit principal surfaces with R.
I've come across the package princurve created by TrevorHastie, but it allows to fit principal curves only. How can I fit two-dimensional principal surfaces in R?

Thank you.
#
I couldn't find anything, but you might try searching on "thin plate
splines" on rseek.org. I realize these are different than principal
surfaces, but they might nevertheless be useful to you.   Or not.

Cheers,
Bert


Bert Gunter

"The trouble with having an open mind is that people keep coming along
and sticking things into it."
-- Opus (aka Berkeley Breathed in his "Bloom County" comic strip )
On Sun, Jan 1, 2017 at 2:56 AM, Neverstop . <neverstop at hotmail.it> wrote:
#
What are the operational requirements for satisfaction with this assignment? I wonder if a regression "surface" of some other sort would be sufficient?  Assuming the three dimensions are X,Y, Z, then a loess (or locfit) "surface" or a regression estimate that used a spline surface for  Z ~ X+Y would not be a principal surface, but should satisfy some of your other requirements. This question would become clearer if you would offer data for analysis (the lack of which is one reason that duplicates of this question were closed as off-topic at both CrossValidated.com and StackOverflow.)
R-help is a plain-text mailing list.
David Winsemius
Alameda, CA, USA
1 day later
#
Did you search for the princurve package?  Sounds like it may be what you
want.
See https://cran.r-project.org/web/packages/princurve/index.html

Best, MEH

Mark E. Hall, PhD
Assistant Field Manager
Black Rock Field Office
Winnemucca District Office
775-623-1529.
On Sun, Jan 1, 2017 at 2:56 AM, Neverstop . <neverstop at hotmail.it> wrote: