search for: bspline

...eed to compute tensor product of B-spline defined over equi-spaced break-points. I wrote my own program (it works in a 2-dimensional setting) library(splines) # set the break-points Knots = seq(-1,1,length=10) # number of splines M = (length(Knots)-4)^2 # short cut to splineDesign function bspline = function(x) splineDesign(Knots,x,outer.ok = T) # bivariate tensor product of bspline btens = function(x) t(bspline(x[1]))%*%bspline(x[2]) # numebr of points to plot ng = 51 # create vectors for plotting xgr = seq(-1,1,length=ng) xgr2 = expand.grid(xgr,xgr) # generate random coef. of...

...) ##- Error in predict.npolySpline(isP, xo, deriv = der) : ##- subscript out of bounds ## Same for B-spline (instead of Polynom.): for(der in 0:3) # deriv=3 fails! print(formatC(try(predict(isB, xo, deriv = der)$y), wid=7,format="f"), quote = FALSE) ##- Error in predict.nbSpline(isB, xo, deriv = der) : ##- subscript out of bounds ## Here, also for deriv= 1 & 2, ## the value left to first proper knot is NA (and should not) <<< ! >>> ------ Martin Maechler <maechler@stat.math.ethz.ch> http://stat.ethz.ch/~maechler/ Seminar fuer Statistik, ETH-Z...

Hello, I am using library fda and I can not run a lot of functions because I receive the error: Error in bsplineS(evalarg, breaks, norder, nderiv) : couldn't find function "spline.des" do you know how I can fix that? Thnaks. Liliana

On Sat, Aug 17, 2013 at 7:08 AM, Tom Browder <tom.browder at gmail.com> wrote: > On Sat, Aug 17, 2013 at 6:15 AM, Tom Browder <tom.browder at gmail.com>wrote: > >> On Fri, Aug 16, 2013 at 2:58 PM, Tom Browder <tom.browder at gmail.com>wrote: >> >>> On Fri, Aug 16, 2013 at 2:55 PM, Tom Browder <tom.browder at gmail.com>wrote: >>>

...erroneous attempts are below. Can you suggest corrections? library(fda) domain<-matrix(c(1,1.1,1.2,2,2.1,2.1,3,3.1,3.2,4,4.1,4.2,5,5.1,5.2), nrow = 5, ncol=3, byrow=TRUE) range<-sin(domain)+matrix(rnorm(15,sd=0.1),nrow=5, ncol=3) myfd<-Data2fd(argvals=domain, y=range) #Error in create.bspline.basis(argvals) : # rangeval is not a vector; class(rangeval) = matrix mybasis<-create.bspline.basis(rangeval=c(0.5,5.2), nbasis=6) myfd<-Data2fd(argvals=domain, y=range,basisobj=mybasis) #Error in smooth.basis(argvals, y, fdP, wtvec = w, fdnames = fdnames) : # 'y' is not the sam...

....0 OS: Mac OS X Submission from: (NULL) (128.135.239.11) I have recently installed the version 2.8.0 of R along with package fda (v 2.0.2) and its dependencies (including package splines v. 2.8.0). Here are my problems: 1) The package splines should feature functions such a predict.bs, predict.bSpline and such and it does not! I can make calls to bs, ns, and interpSpline but not to any predicting function. Example: > library(splines) > predict.bs Erreur : objet "predict.bs" non trouv? # in english: object "predict.bs" not found Also, I cannot track the package splines...

...-spline/smooth.basis derivative > matrices > > > Hello again,(my message was truncated the first time!) > > Thank you, vito for these functions: I've been able to look > at the smoothing spline matrix and its derivatives that I > constucted with functions "create.bspline.basis" and > "smooth.basis" > > Understanding from Hastie and Tishbirani that the > S(lambda)-- the penalty matrix > > decreases in magnitude compared to the matrix [y[i,j] - > f(x)] would it not be of fundamental importance to > > have the penalty matri...

...looks perfectly reasonable. The prediction from the inverted spline matches the curve over part of the range but goes crazy elsewhere. I would have expected it to be reasonably close to this well-behaved curve over the whole range. I have looked at the docs for interpSpline, backSpline, predict.bSpline, ... and nothing jumps out at me. Please be gentle, if possible, in explaining to me what I'm missing ... thanks, Ben Bolker

Dear All, I am researching financial market microstructure and have approx 4 x 10^7 multivariate 2D data samples which I have counted into a 250 x 390 bin matrix (frequency counted 2D histogram) in order to more efficiently manage the volume of data. I now wish to construct a smooth kernel density estimate (Gaussian kernel function) using this binned data. Does anyone in the R community know

Hello, We would like to apply the smooth monotone function to our data which correspond to a non-linear function. We follow the example posted on the web, but in our case it did not apply. We always get a straight line in response. Which parameters we should change. ind.basis = create.bspline.basis(c(min(time),max(time)),nbasis=38,norder=4) Wfdob = data2fd(rep(0,length(time)),time,ind.basis) fd = data2fd(alphaest,time,ind.basis) solution= smooth.monotone(time,alphaest,Wfdobj=Wfdob,lambda=0.1, conv=0.01) Alpha = solution$beta[1]+solution$beta[2]*eval.monfd(time,solution$Wfdobj)...

Hi. I have been using the splines package for my work, in particular, the bs() function and associated predict() method. I now find myself in need of the derivatives of this beast. In the man page for predict.bSpline I found `predict(object, x, nseg=50, deriv=0, ...)' but alas this is not implemented (deriv= is ignored). I contacted the package maintainers who were most helpful. Bill Venables kindly replied "Curiously, the S-PLUS code did have that feature, I recall, but as I no longer have S-PLUS I...

Hi Can any one inform me an r package used to fit Cubic M-splines or r functions used to have design and penalty matrices in fitting cubic  M-splines please? Thanks, Tigist [[alternative HTML version deleted]]

Hello, I've noticed that SPLUS seems to have a function for evaluating derivative matrices of splines. I've found the R function that evaluates matrices from 'smooth.spline'; maybe someone has written something to do the same with smooth.basis? regards, s

Dear all, I'm trying to build some fda objects using smooth.pos() from fda package. Here is what I try: >X <-c(0.4541139, 0.4802537, 0.5091228, 0.3894931, 2.1512258, >7.9301281, 62.3876109, 54.1168247, 67.8769904, 91.0670615, 20.0682440) >basis <- create.bspline.basis(c(-100, 100), 7, 5) >MDA.fdPar <- fdPar(basis, 2, 0.01) >MDA.fda.joint <- smooth.pos(argvals=seq(from=-100, to=100, length=11), >y=X, WfdParobj=MDA.fdPar) Iter. PENSSE Grad Length 0 1780.78 -1780.78 1 1753.525 14.318 2 1749.296 15.0611 3...

...39 boys from the Berkeley growth study, each one measured 31 times at the ages of 1.00 1.25 1.50 1.75 2.00 3.00 4.00 5.00 6.00 7.00 8.00 8.50 9.00 9.50 10.00 10.50 11.00 11.50 12.00 12.50 13.00 13.50 14.00 14.50 15.00 15.50 16.00 16.50 17.00 17.50 18.00 from experience we know that a bspline basis with at least 12 to 15 knots is needed to model a complete growth curve. I take the spline basis matrix as both fixed and random effect (not forgetting the -1, since the intercept is included in the splines). When I use 3 knots (5 parameters), computation time is about 5 Minutes on a 750MHz...

...t; * The way the DC classes are designed are great but: > ??? - FXDCPrint is just way too buggy and incomplete.? It probably > should be removed from the release (or a word of warning in the docs > should be present) > ??? - FXDC is in serious need of supporting: > ??? ??? * bezier/bspline drawing > ??? ??? * hit testing > ??? ??? * FXDCWindow''s with double buffering (sans opengl) > ????? Three things that are relatively common in other toolkits (Tk, > FLTK, etc). These are again shortcomings of FOX itself. You can report them to the foxgui-users list if you l...

...atex.Rd had unbalanced braces in section "Supported Options". o merge() with multiple and differently-named match columns failed. (PR#1510) o NAs in right-justified unquoted character matrices were being mis-aligned. This mainly affected printing data frames. o predict.*bSpline() bugs extrapolating for deriv >= 1 (PR#1473), and predict.[bn]s bug e.g. for bs(x), reported by Ch.Sangiorgio. o qr.X was failing if n < p. Fixed, but only in the case when pivoting does not occur. (PR#1519) o xx[, 1:3] was returning a list if xx had only one row, even though x...

...atex.Rd had unbalanced braces in section "Supported Options". o merge() with multiple and differently-named match columns failed. (PR#1510) o NAs in right-justified unquoted character matrices were being mis-aligned. This mainly affected printing data frames. o predict.*bSpline() bugs extrapolating for deriv >= 1 (PR#1473), and predict.[bn]s bug e.g. for bs(x), reported by Ch.Sangiorgio. o qr.X was failing if n < p. Fixed, but only in the case when pivoting does not occur. (PR#1519) o xx[, 1:3] was returning a list if xx had only one row, even though x...

R-1.7.0 built on NetBSD 1.6, but the validation test suite failed: Machinetype: Intel Pentium III (600 MHz); NetBSD 1.6 (GENERIC) Remote gcc version: gcc (GCC) 3.2.2 Remote g++ version: g++ (GCC) 3.2.2 Configure environment: CC=gcc CXX=g++ LDFLAGS=-Wl,-rpath,/usr/local/lib make[5]: Entering directory `/local/build/R-1.7.0/src/library' >>> Building/Updating

...ion: create D: file: /usr/local/lib/R/library/splines/R-ex/ns.R action: create D: file: /usr/local/lib/R/library/splines/R-ex/periodicSpline.R action: create D: file: /usr/local/lib/R/library/splines/R-ex/polySpline.R action: create D: file: /usr/local/lib/R/library/splines/R-ex/predict.bSpline.R action: create D: file: /usr/local/lib/R/library/splines/R-ex/predict.bs.R action: create D: file: /usr/local/lib/R/library/splines/R-ex/splineDesign.R action: create D: file: /usr/local/lib/R/library/splines/R-ex/splineKnots.R action: create D: file: /usr/local/lib/R/library/splines/...