search for: spar

...users NA\connelmp and NA\guminssa both started getting messages that they were not part of the required group. Here is the log for you all to see... >From /var/log/apache2/error_log [Thu Feb 14 13:22:24 2008] [error] [client 192.63.212.40] CHKAUTH: is na\\huynhsv a member of NA\\USTR-LINUX-1-SPAR? [Thu Feb 14 13:22:24 2008] [error] [client 192.63.212.40] CHKAUTH: YES, na\\huynhsv is listed amongst the NA\\USTR-LINUX-1-SPAR group members [Thu Feb 14 13:22:24 2008] [error] [client 192.63.212.40] CHKAUTH: is na\\huynhsv a member of NA\\USTR-LINUX-1-SPAR? [Thu Feb 14 13:22:24 2008] [error] [cli...

I've had a play with this and, due to my own short-comings, remain none the wiser. In particular, I'm not sure what value of 'spar' is consistent with the magic lambda=1/1600 for quarterly data. I initially interpreted spar as lambda and tried setting spar=1/1600. This results in almost no smoothing while spar=1600 causes an error. The smooth.spline function seems to want something in the region of 0 to 1. The closer...

...mail) given _df_ . Unfortunately the documentation is not clear to me. Maybee someone can help to answer in my view the basic question: If the penalized log likelihood is L = (y - f)' W (y - f) + lambda c' Sigma c how the _lambda_ in the above equation is computed if _df_ is given and _spar_ not? And, Is there a way to define lambda directly? Yours. /E Hi, I am usign the smooth.spline function. I am not sure how the _df_ (degrees of freedom) parameter, if set, influences _lambda_ in eq: L = (y - f)' W (y - f) + lambda c' Sigma c Is _df_ substituting tr(Sigma), if de...

...bpage that is secured by a Domain Local group. Sometimes they get in, others they don't. Here are some examples from the logs. /var/log/apache2/error_log [Tue Feb 12 18:54:52 2008] [error] [client 172.xx.xxx.xxx] GROUP: NA\\selltc not in required group(s)., referer: https://ustr-linux-1/scm/spar/trac/browser/trunk/common/include/channels [Tue Feb 12 18:55:00 2008] [error] [client 172.xx.xxx.xxx] GROUP: NA\\selltc not in required group(s)., referer: https://ustr-linux-1/scm/spar/trac/browser/trunk/common/include/channels [Tue Feb 12 18:56:12 2008] [error] [client 172.xx.xxx.xxx] GROUP: NA\\...

...f there was an easier way. Thank you for any advice. Here is an example: # use an example panel function using smooth.spline; however, the issue relates to all panel functions # a panel function to fit smoothed lines through data panel.smooth.spline <- function(x, y, w=NULL, df, spar = NULL, cv = FALSE, lwd=plot.line$lwd, lty=plot.line$lty,col, col.line=plot.line$col, type, horizontal=FALSE,... ) { # Deepayan Sarkar code from: http://www.mail-archive.com/r-help at r-project.org/msg39654.html x <- as.numeric(x) y <- as.numeric(y...

...(this is part of a bigger program) the following fails as a function, but runs OK if we comment out the fnfn_ function() line. Any hint would be appreciated. -Yudi- R : Copyright 1998, The R Development Core Team Version 0.63.0 Beta (Nov 13, 1998) -- on WIndows3.11 fnfn _ function (m=10,n=10,spar=2) { fn _ function(u,v){ uc_ u-floor(m/2)-1 vc_ v-floor(n/2)-1 exp(-(uc^2 + vc^2)/2/spar^2) } x1_ 1:m ; x2_ 1:n x12 _ outer(x1,x2,FUN='fn') filt _ x12/sum(x12) image(filt) } > source('fnfn.r') > m_10; n_10; spar_2 > fnfn(m,n,spar) Error in ge...

...oth.spline in R with parameters that match my problem above. So I tried this sequence of calls in R: > x <- c(1.,5.,10., 15., 20., 25., 30., 35.) > y <- c(-999.98099, -1001.61875, -1007.9, -1019.36875, -1036.4, > -1059.21875, -1087.9, -1122.36875) > smooth.spline( x,y, w=NULL, spar=0.5, cv=TRUE, all.knots= TRUE ) I am unsure about spar being the smoothness parameter, about where to put the standard errors of the points, and about the return of the smooth.spline function: Smoothing Parameter spar= 0.5 lambda= 0.006833112 Equivalent Degrees of Freedom (Df): 3.221101 Pe...

...ng to compile R-2.2.1 on Solaris 2.9 with a 64-bit build. Following the instructions in "R Installation and Adminstration", I changed the following settings in "config.site": CC="gcc -m64" F77="g77 -64" CXX="g++ -m64" LDFLAGS="-L/usr/local/lib/sparcv9 -L/usr/local/lib" But I got the following error messages: configure:3987: gcc -m64 -I/usr/local/include -L/usr/local/lib/sparcv9 -L/usr/local/lib conftest.c >&5 /usr/ccs/bin/ld: skipping incompatible /work/net-local-b/sparc- sun-solaris2.9/bin/../lib/gcc/sparc-sun-solaris2.9/3.4.2...

...ument. Here is the custom panel, a wrapper for smooth spline. I copied panel.loess and replaced the loess arguments with smooth.spline(). [Note: I would like to use the cross-validation fitting properties of smooth.spline.] library(lattice) panel.smooth.spline<-function(x,y,w=NULL, df, spar = NULL, cv = FALSE, lwd=plot.line$lwd, lty=plot.line$lty,col, col.line=plot.line$col, type, horizontal=FALSE,... ){ x <- as.numeric(x) y <- as.numeric(y) ok <- is.finite(x) & is.finite(y) if (sum(ok) < 1) return() if (!missing(col)) {...

Can someone please show me how to smooth time series data that I have in the form of a zoo object? I have a monthly economies series and all I really need is to see a less jagged line when I plot it. If I do something like s <- smooth.spline(d.zoo$Y, spar = 0.2) plot(predict(s,index(d.zoo)), xlab = "Year") # not defined for Date objects and if I do something like plot(predict(s,as.numeric(index(d.zoo))), xlab = "Year") # one straight line, no matter the value of spar What am I doing wrong? (The unsmoothed series plots just fin...

Dear list, if I do smooth.spline(tmpSec, tmpT, all.knots=T) with the attached data, I get this error-message: Error in smooth.spline(tmpSec, tmpT, all.knots = T) : smoothing parameter value too small If I do smooth.spline(tmpSec[-single arbitrary number], tmpT[-single arbitrary number], all.knots=T) it works! I just don't see it. It works for hundrets other datasets, but not for

...OS X (Darwin). Both system run the latest public versions of R and Pspline. predict.smooth.Pspline produces only NaN instead of predicted values when norder>2: > library (Pspline) > tt <- seq (0,1,length=20) > xt <- tt^3 > > fit <- smooth.Pspline (tt, xt, norder=3,spar=0.0001, method=1) > predict.smooth.Pspline (fit, tt, nderiv=0)[,1] [1] NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN [20] NaN > smooth.Pspline seems to work fine since >plot (tt,xt) >lines (fit) produces the expected result (curve fits data closely). U...

...586179e-01, -1.189e+00, -2.712e-03, -2.126218e+00, 7.242e-01, -1.149e+00, 9.377e-01, -8.457e-03, -3.035807e-02, -1.697e-02, 2.279305e-02, -4.7e-07, -1.797e-07, -1.183e-07, -4.3494e-05) abunddst=c(Cnictdat$SegDistKm)*1000 esw=0.5 fit4=vgam(Cnict~factor(Hab)+s(Clr, spar=2)+s(Rd, spar=2)+s(YL, spar=2)+ s(Pop, spar=2)+offset(log(2*esw*abunddst)), zanegbinomial, data=Cnictdat, coefstart=coefnict, control=vgam.control(se.fit=T, exp=1e-4, maxit=500, trace=T)) Thanks, John

...smooth.spline in R with parameters that match my problem above. So I tried this sequence of calls in R: > x <- c(1.,5.,10., 15., 20., 25., 30., 35.) > y <- c(-999.98099, -1001.61875, -1007.9, -1019.36875, -1036.4, -1059.21875, -1087.9, -1122.36875) > smooth.spline( x,y, w=NULL, spar=0.5, cv=TRUE, all.knots= TRUE ) I am unsure about spar being the smoothness parameter, about where to put the standard errors of the points, and about the return of the smooth.spline function: Smoothing Parameter spar= 0.5 lambda= 0.006833112 Equivalent Degrees of Freedom (Df): 3.221101 Penali...

...<- sm.spline (time,x, norder=2, cv=FALSE) > summary (traj.spl) Length Class Mode norder 1 -none- numeric x 36 -none- numeric ysmth 36 -none- numeric lev 36 -none- numeric gcv 1 -none- numeric cv 1 -none- numeric df 1 -none- numeric spar 1 -none- numeric call 6 -none- call > Has anybody been able to fit quintic splines with Pspline or has some advise how to avoid the Singularity error? Thank you, Martin -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read...

...at the code for smooth.spline, but can't readily follow it. Here's a toy example: > x<-seq(from=0,to=4*pi,length=1002) > y<-sin(x) > ss<-smooth.spline(x,y=y,all.knots=F,nknots=25) > ss Call: smooth.spline(x = x, y = y, all.knots = F, nknots = 25) Smoothing Parameter spar= -0.4573636 lambda= 1.006117e-09 (14 iterations) Equivalent Degrees of Freedom (Df): 26.99935 Penalized Criterion: 3.027077e-06 GCV: 3.190666e-09 > str(ss) List of 15 $ x : num [1:1002] 0 0.0126 0.0251 0.0377 0.0502 ... $ y : num [1:1002] 2.88e-05 1.26e-02 2.51e-02 3.77e-02 5.02e-...

...spline" package are successfully installed while others are not. The code with which I'm working is more complicated, but the following highlights my problem. If I run the following code > tt <- seq (0,1,length=20) > xt <- tt^3 > fit <- smooth.Pspline (tt, xt, norder=3,spar=0.0001, method=1) > pred<- predict.smooth.Pspline (fit, tt, nderiv=0) I get this error: "Error: could not find function "predict.smooth.Pspline" " So it seems like "smooth.Pspline" is loading but "predict.smooth.Pspline" isn't. Any ideas about wh...

...e, whence it has a calculable derivative, but is there a (simple) way of getting at it? One workaround that I have thought of is to take a fine grid of points, evaluate the function returned by splinefun() at these points, put an interpolating spline through these points using smooth.spline() with spar=0, and then calculate the derivative with predict(...,deriv=1). This seems pretty kludgy, and fraught with a certain amount of peril. Is there a better way? cheers, Rolf Turner rolf at math.unb.ca

Dear All, I know how to export graphics as pdf files and then how to include them in LaTeX documents. However, I do not know how to do in order to have the text of the graphics written with the font selected for the LaTeX document. Is that possible? Thanks in advance, Paul

...lt;-as.array(data11, dimnames=dimnames(data11)) The analysis that I would like to do with R(CollocInfer) runs as follows(I print for the sake of introduction here though to run this line of code requires extensive setup) res11a = inneropt(coefs, times=times, data=data11a, lik=lik, proc=proc, pars=spars, in.meth='nlminb', control.in=control.out) The response I get is Error in 'colnames<-'('*tmp*', value = c("V","R")) : length of 'dimnames' [2] not equal to array extent It seems to me that there is a fundamental match of 'dimname...