similar to: Crash report: projection pursuit & predict

Dear R users, I am looking for the source code of the implementation of ppr (Projection Pursuit Regression) in R. It will be great if citations of the source papers on which the implementation is based, are also provided. Thank you, Arvind Iyer, Grad student, Deptt. of Biomedical Engineering Viterbi School of Engineering University of Southern California, Los Angeles [[alternative HTML

Hi, Recently, I installed the 1.1.0 version of R (for Windows), since it includes an implementation of Projection Pursuit (I failed to write my own version of PP as a standalone C++ program). As far as I know, R offers two interfaces/sintax for the ppr() function. The first one requieres a regression formula and a data frame. The other requieres X, a matrix with the explanatory variables, and Y,

Sorry about the lack of detail. I am running R v.1.2.2. I can recast my question (which I think I have partially answered) more succinctly as follows: 1. This seems to work (note that group takes values 1,2,3,4, or 5): my.newfun <- function(x) myfile <- lm(award ~ ilogemp + ilogage, x) test.by <- by(wintemp, as.factor(wintemp$group), my.newfun) 2. This does not work (leaving aside

I used the supersmoother function in the modreg package as follows: super <- supsmu(ilogemp,award) Then I decided that I might want additional explanatory variables (other than ilogemp) in my model. The ppr function in modreg seemed a logical extension of supsmu from univariate to multidimensional explanatory variables. As a "check" I ran the following: pprest <-

Hello, Just a quick question about ppr in library modreg. I have looked at Ripley and Venables 2002 and it says that projection pursuit works "by projecting X in M carefully chosen directions" I want to know how it choses the directions? I presume it moves around the high-dimensional space of unit vectors finding ones that separate the response variables, but how. I looked at the

Dear R list: I want to make separate estimates for each level of the variable "group." After consulting many sources I am stumped as to why the following does not work: > wintemp <- subset(alltemp, winner==1) > my.ppr <- function(x) + { + if(nrow(x) >= 50) { + pprfile <- ppr(award~ilogemp, data=x,nterms=5,max.terms=10,optlevel=3) + summary(pprfile) +

Dear R listers -- The program below does the following tasks: 1. It creates a file (wintemp4) that is a subset of alldata4 consisting of "winner" records in 50 industry groups (about 5400 obs); 2. It defines a function (myppr1) that runs the ppr function in modreg once to generate goodness of fit (sum of squared errors) measures by number of terms included in model and then reruns

Dear R listers -- I have several questions about using the ppr command in the modreg module. I discovered -- quite by accident -- that if I re-order the data, I obtain different results. The output below shows what I mean. I have two datasets (dataset1 and dataset2) that are identical (tested using proc compare in SAS) except for the fact that the records are in different order. Below I have

Hi, I am solving a projection pursuit regression problem, of the form y = \sum_i f_i (a_i^T x), where a_i are unknown directions, while f_i are unknown univariate link functions. The following is known about each f_i: 1. f_i (0) = 0  (that is, each f_i passes through the origin) 2. f_i is monotonic. Is there a way to ensure that the function ppr() in R produces solutions that respect the

Buenas. El Paquete pdp es muy fácil de usar, pero cuando se lo aplico a mis datos me da: Error in eval(stats::getCall(object)$data) : object 'x.data' not found. Os copio abajo un ejemplo de aplicación a un RF. El mio es de un boosted regression trees (paquete gbm). No sé si esa puede ser la razón del error. En el paquete pdp no especifica que sea solo para RF, aunque en los

Hi all, I am doing a projection pursuit regression using the ppr() function from modreg. I would also like to use predict.ppr(). However, I cannot find any information about it in the help files. There is a link to predict.ppr in the index for modreg, but that link is to the help for ppr(). Has predict.ppr() not been implemented? If not, does anyone have a suggestion as to how to implement

The pursuit projection packages ppr is an excellent contribution to R. It is great for one-to-three ridge fits, often somewhat intuitive, and for multi-ridge fits, where it at least describes a lot of variance. Like many folk, I need to report the fits obtained from ppr to the greater, outside, non-R world. It is fairly obvious how to use the terms alpha and beta to report on

Dear R-developers, I am trying to experiment with projection pursuit (PP), and different indexes for the same, especially using the tourr package. However, I've noticed that a PP index in the said package is only a function of the projected data. Could I modify the function so that the index sees the projection bases instead (or in addition to the data) ? Alternatively, if some other package

Hello, I am using projection pursuit regression parameters selection. Does anyone has experience on the range to test for df parameter (spline kernel) and gcvpen (gcvspline kernel)? I don't find any information about this. Thanks in advance. Joao Moreira [[alternative HTML version deleted]]

Dear helpers Does R have a package that performs projection pursuit density estimation? Or anyone knows code in Matlab or C for example to do this? Thank you all Luis

Hello everybody, I'm a new user of R and I'm working right now with the ACE function from the acepack library. I Have a question: Is there a way to predict new responses using ACE? What I mean is doing something similar to the following code that uses PPR (Projection Pursuit Regression): library(MASS) x <- runif(20, 0, 1) xnew <- runif(2000, 0, 1) y <- sin(x) a <- ppr(x, y,

Dear all, I have been doing tests using SVM, random forests and PPR. The data is from a data stream (that is, the data for training and for test is always increasing / changing). With SVM and random forests everything is ok, but with ppr there are situations where it crashes. For the examples I have used I noticed that if one of the variables has just one value (it can happen), it crashes for

Full_Name: Yves Rosseel Version: 2.2.1 OS: i686-pc-linux-gnu Submission from: (NULL) (157.193.116.152) Dear developers, For the single-model case, the anova.mlm() function does not seem to handle multi-parameter predictors (eg factors) correctly. A toy example illustrates the problem: Y <- cbind(rnorm(100),rnorm(100),rnorm(100)) A <- factor(rep(c(1,2,3,4), each=25)) fit <- lm(Y ~ A)

I have been running into difficulties with dispatching on an S4 class defined in the SparseM package, when the method calls are inside a function passed as the f= argument to optimize() in functions in the spdep package. The S4 methods are typically defined as: setMethod("det","matrix.csr", function(x, ...) det(chol(x))^2) that is within setMethod() rather than by name before

Full_Name: Hugh Miller Version: 2.8.1 OS: XP Submission from: (NULL) (128.250.24.101) Hi there, I've been looking at approaches that use the projection pursuit regression function fairly (ppr) heavily, and have discovered that it crashes my R system on occasion. It only happens with the inputs are pathological in some way I don't understand. I have pasted such an example below. Any