similar to: ssanova help

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

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.fit) part.fit=ssanova(y~x, partial=~cbind(as.numeric(x1),x2))

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:

Hi all, I'm fairly new to R, so I'm still trying to feel out what is available to me. I would like to be able to plot joint density in a two dimensional plot where density is indicated by color or darkness gradients, like a 2d color coded topographic map. Ideally, the output would be something I could then plot other points or lines on. Currently, I'm calculating joint density with

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

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

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

Hi, Have anyone using speex experience the following effect: When I say "hello", it will be become "hellololololololo....". This does not always happened. The chance of happen is quite random. This effect does not happen in my LAN. It only happen when ping time is about 500 milisecond or more. Is it a speex problem? Any suggestion how to resolve this issue? Thanks in

Hi, I am running a negative binomial model using Gamlss and when I try to include random effect, I get the following message: Warning messages: 1: In vcov.gamlss(object, "all") : addive terms exists in the mu formula standard errors for the linear terms maybe are not appropriate 2: In vcov.gamlss(object, "all") : addive terms exists in the sigma formula standard

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

** High Priority ** Hello R users My issue is both theorical and technical. I would like to run a penalised discriminant analysis with the fda() function, but I don''t know all the details of splines theory. I try on the example of the phonems from the article "Penalised Discriminant Analysis" of Hastie, Buja and Tibshirani 1994 : 5 groups and 256 variables. The 256

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

Hi http://mailman.ds9a.nl/pipermail/lartc/2002q3/004977.html I have the same problem to discuss as in the above link. I want to allocate say X MBit per individual connection regardless of the number of connection . KIndly could anyonen suggest me how to proceed. I have tried with SFq but is doesnot yeild my requirement.. Thanks, Namitha. _______________________________________________ LARTC

Hellos, I have date/times in 64-bit Julian format, that is the number of milliseconds since 1 jan 1970 (or something like that). E.g., "1159884877406" For the life of me I can't figure out how to render these in a more readable format, e.g., "2006-10-04 23:12:93.191" (and I have tried to do my best searching the archives and help etc ...) Thank you most

On Thu, Dec 7, 2017 at 12:18 PM, Timothy B. Terriberry <tterribe at xiph.org> wrote: > [...] > > Sorry, I can't reply to the original patches because I didn't actually get > that e-mail due to local trouble with my mail server. I could pull the > patches from the list archive, however. Thanks for the reports. > Thanks for recovering them and having a look. I

James Zern wrote: > On Tue, Nov 28, 2017 at 3:22 PM, James Zern <jzern at google.com> wrote: >> On Mon, Nov 20, 2017 at 1:07 PM, James Zern <jzern at google.com> wrote: >>> Just an attempt to avoid overflows with an explicit check, I don't know if >>> there's a better way to identify corrupt input here. >>> >>> James Zern (2):