search for: ssanova

Dear all, I am trying to use ssanova of the gss package but met some error that I cannot figure out the answer for. Here is the code I am using to explain the problem. library(gss) set.seed(5732) x=(1:100)/100 y=1+3*sin(2*pi*x)+2*(x>0.7)+rnorm(x) x1=rnorm(100) x2=rnorm(100) part.fit=ssanova(y~x, partial=~cbind(x1,x2)) summary(part...

Hi all, I'm using the ssanova function from the gss package to fit smoothing spline anovas, and am running into some difficulty. For my data, I have measurements at 2 milisecond intervals for every observation. Every observation does not have the same duration, so I have scaled the times for each observation to a scale between...

Hi I'm using the gss package to fit thin plate splines. I've fitted the tps without problems and got a ssanova object. Then I wanted to do some prediction using a new set of data (latlon data) and I got an error message: > predict.ssanova(recs.spliney,recn.grid,se.fit=FALSE) Error in array(x, c(length(x), 1), if (!is.null(names(x))) list(names(x), : attempt to set an attribute on NULL Can some...

Hi everyone, I use this code to add the fitted line to a plot, however, I think I cannot access the fitted equation. How my I find the equation in the object ksi? thanks alot. library (gss) ####generate (simple linear) x and y x= 2* (1:40) y= x-1 plot (x,y) ## fit and draw ksi = ssanova0 (x~y, method= "m") lines( ksi$qwk, col = "red") -- View this message in context: http://r.789695.n4.nabble.com/ssanova-ssanova0-and-adding-the-fitted-line-to-a-plot-tp4264679p4264679.html Sent from the R help mailing list archive at Nabble.com.

Hello: I sent this question the other day with the wrong subject heading and couple typos, with no response. So, here I go again, having made those corrections. I would like to estimate, for lack of a better description, a partially additive non-parametric model with the following structure: z~ f(x,y):w1 + g(x,y):w2 + e In other words, I'd like to estimate the marginals with respect to

A new version of gss, version 0.8-2, is on CRAN now. Numerous new functionalities have been added since my last r-announce post. An ssanova1 suite has been added since version 0.7-4. It implements low-dimensional approximations of the smoothing spline ANOVA models of the ssanova suite. ssanova1 scales much better than ssanova with large sample sizes. A gssanova1 suite is added for non Gaussian regression. Similar to ssanova1, it pr...

Hi I would like to know if the "penalty associated with the fit" that is printed by the summary.ssanova function is the smoothing parameter (\lambda) of the criterion function that's minimized to find f. Thanks EJ -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info", &q...

A new version of gss, version 0.8-2, is on CRAN now. Numerous new functionalities have been added since my last r-announce post. An ssanova1 suite has been added since version 0.7-4. It implements low-dimensional approximations of the smoothing spline ANOVA models of the ssanova suite. ssanova1 scales much better than ssanova with large sample sizes. A gssanova1 suite is added for non Gaussian regression. Similar to ssanova1, it pr...

I use smooth.spline(x, y) in package modreg and I would like to get value of penalized log likelihood and preferable also its two parts. To make clear what I am asking for (and make sure that I am asking for the right thing) I clarify my problem trying to use the same notation as in help(smooth.spline): I want to find the natural cubic spline f(x) such that L(f) = \sum_{k=1}{n} w[k](y[k] -

...of the square second derivative of beta upon the frequencies and J(beta)=(beta)^T*OMEGA*beta No function for producing one dimensional penalty object is providing with gen.ridge function.. What is the best (and more efficient) function to use to fit OMEGA ? I have ever tried to use the function ssanova in the gss package but the solution of my program is not symetrical : I am definitly wrong in my approach. Thank you in advance for the users of mda package that would answer me. Damien Berlemont student at : DESS Methodes Scientifiques de Gestion - Paris-X Nanterre Damien.Berlemont at ccf.fr...

Hello, I'm a student at the University of Klagenfurt / Austria and I need some help ! I have to predict 24 daily load-values. Therefor I got a dataset with following colums: 24 past daily load-values 6 past daily temperature-values My goal is to find a model (GAM with thin plate splines) in R. I found the function "gam" in the R-library "mgcv", but it just fits

Hi, I was wondering if there is a function in R which performs bivariate non parametric smoothing which allows for the possibility of including some weights in the smoothing (for each data points in my grid I have some predefined weights that I would like to include in the smoothing). Thanks, Ben _________________________________________________________________ [[alternative HTML

Dear all, There has been some time since I asked about interpolation in higher (>2) dimensions, and I must admit I failed to write a function to do this myself the last time, but eventually ended up doing it in MATLAB. I tried to translate the MATLAB code, but MATLAB code is so much more opaque than R (S) code, so I failed that too, mainly because I could only get one MATLAB session, I would

I just heard a talk about a semi-parametric model. I was quite excited by the idea. This model is fitted y= xB + g(z) + e where x is a data matrix, B a column vector, z is another data matrix, and g is a smooth model fitted by a Kernel Smoothing regression (I got the idea any smoother would do as well). The speaker said that when z is considered as a "control" variable, and there is