similar to: Source code for ppr (Projection Pursuit Regression)

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,

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

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

Folks, The projection pursuit regression function in the base R seems to crash when the optimization level is set to zero, i.e. the initial ridge terms are accepted without refitting. I encountered this problem in an out-of-sample prediction exercise using predict. But further investigation suggests the issue is with the ppr fit and predict just sppeds up the crash. The other optlevels seem to be

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

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]]

Is there a way to create a 'bubble plot' in R? For example, if we define the following data frame containing the level of y observed for 5 patients at three time points: time<-c(rep('time 1',5),rep('time 2',5),rep('time 3',5))

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

G'day all, I had a look at the GLM code of R (1.4.1) and I believe that there are problems with the function "glm.fit" that may bite in rare circumstances. Note, I have no data set with which I ran into trouble. This report is solely based on having a look at the code. Below I append a listing of the glm.fit function as produced by my system. I have added line numbers so that I

I'm not understanding something. I'm trying to add xlab & ylab to a balloon plot of a table object. From docs I thought following should work: require(gplots) # From balloonplot example: # Create an example using table xnames <- sample( letters[1:3], 50, replace=2) ynames <- sample( 1:5, 50, replace=2) tab <- table(xnames, ynames) balloonplot(tab)

I received a floating point exception from "make check" on a Linux alpha (Red Hat 6.0) using egcs-1.1.2 compiler, and compiling from R-base-0.64.1-5.src.rpm packaged by Red Hat Contrib|Net. make check fails with return code 136 from modreg-Ex.Rout. The final output is: > plot(update(rock.ppr, sm.method="gcv", gcvpen=2), + main = "update(...,

I think this is another one of the same kind of bugs in ts: Version 0.64.1 (May 8, 1999) ... > z <- ts(matrix(1:20,10,2), start=c(1969,1), frequency=12) > (z > 5) | (z < 2) Error: invalid time series parameters specified > Paul -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-devel mailing list -- Read

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) +

Hi all, I have a data.frame with the following colnames pattern: x1 y11 x2 y21 y22 y23 x3 y31 y32 ... I.e. I have an x followed by a few y's. What I would like to do is turn this wide format into a tall format with two columns: "x", "y". The structure is that xi needs to be associated with yij (e.g. x1 should next to y11 and y12, x2 should be next to y21, y22, and

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

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)

Dear R-team As I didn't get any answer to my bug-report last week I have taken the effort and extracted a minimal data set from my data (see below) where the following bug occurs: > glm(SKR.ein.aus ~ ., family = binomial, data = bugdata, na.action = na.omit) Error in names<-.default(*tmp*, value = ynames) : names attribute must be the same length as the vector In addition: Warning

Hello, I've been looking through ?phantom and ?expression and this forum for examples of how I might be able to manipulate some of the names that appear on the y-axis of the barplot below. For example, the "gw" in "ECgw" would appear as a subscript...or "qr" would be the theta symbol followed by subscript "r". My attempts haven't even come close

Hi all, I'm wanting to be able to find and store the z-score of coxph below: - modz=coxph(Surv(TSURV,STATUS)~RAGE+DAGE+REG_WTIME_M+CLD_ISCH+POLY_VS, data=kidneyT,method="breslow") I know summary(modz) will give me this, but how do i extract the standard error or z-score values in a similar way to obtaining the coefficients by coef(modz) ? I think it must be something to do with

Hi, I want to test for a contrast from a regression where I am regressing the columns of a matrix. In short, the following. X <- matrix(rnorm(50),10,5) Y <- matrix(rnorm(50),10,5) lm(Y~X) Call: lm(formula = Y ~ X) Coefficients: [,1] [,2] [,3] [,4] [,5] (Intercept) 0.3350 -0.1989 -0.1932 0.7528 0.0727 X1 0.2007 -0.8505 0.0520