similar to: documentation mistake

Thanks for the info, Dirk. The tarball builds don't run make check (because of a policy decision that it is better to have the sources available on all platforms for testing than to have none if it breaks on a single platform). However the build as of tonight has no problem with make check on the build machine. Did you by any chance forget that Matrix is a recommended package and expected to

I generally run 'make; make check' (with more settings) when building the Debian package. Running 3.2.4 rc from last night, I see a lot of package loading issues during 'make check'. Here is splines as one examples: checking package 'splines' * using log directory '/build/r-base-3.2.3.20160303/tests/splines.Rcheck' * using R version 3.2.4 RC (2016-03-02 r70270) *

On 4 March 2016 at 09:11, peter dalgaard wrote: | Thanks for the info, Dirk. | | The tarball builds don't run make check (because of a policy decision that it is better to have the sources available on all platforms for testing than to have none if it breaks on a single platform). However the build as of tonight has no problem with make check on the build machine. Did you by any chance forget

Hello I'm using function gam from package mgcv to fit splines. ?When I try to make a prediction slightly beyond the original 'x' range, I get this error: > A = runif(50,1,149) > B = sqrt(A) + rnorm(50) > range(A) [1] 3.289136 145.342961 > > > fit1 = gam(B ~ s(A, bs="ps"), outer.ok=TRUE) > predict(fit1, newdata=data.frame(A=149.9), outer.ok=TRUE) Error

Dear all, I need 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

Hi, I am trying to find an elegant way to compute and store some frequently used matrices "on demand". The Matrix package already uses something like this for storing decompositions, but I don't know how to do it. The actual context is the following: A list has information about a basis of a B-spline space (nodes, order) and gridpoints at which the basis functions would be

--- On Mon, 3/12/12, aleksandr shfets <a_shfets at mail.ru> wrote: > From: aleksandr shfets <a_shfets at mail.ru> > Subject: Fwd: Re[2]: [R] B-spline/smooth.basis derivative matrices > To: "Vassily Shvets" <shv736 at yahoo.com> > Received: Monday, March 12, 2012, 5:15 PM > > > > -------- ???????????? ????????? > -------- > ?? ????:

Laurent 14 juin 12:02 afficher les options De : "Laurent" <eddy_l... at hotmail.com> - Rechercher les messages de cet auteur Date : Tue, 14 Jun 2005 03:02:37 -0700 Local : Mar 14 juin 2005 12:02 Objet : bs() function of the splines package R??pondre | R??pondre ?? l'auteur | Transf??rer | Imprimer | Message individuel | Afficher l'original | Retirer | Signaler un

On 08/02/2012 05:00 AM, r-devel-request at r-project.org wrote: > Now I just have to grovel over the R code in ns() and bs() to figure > out how exactly they pick knots and handle boundary conditions, plus > there is some code that I don't understand in ns() that uses qr() to > postprocess the output from spline.des. I assume this is involved > somehow in imposing the boundary

Hello, I'm implementing a function using non uniform B-Splines. Looking at the code of the bs() function, I realized that if the intercept was set to FALSE, the behavior of the function was the following (df is the number of degrees of freedom that I believe can be interpreted as the number of control points): - Compute df- ord + 1 _internal_ knots (ord is the order of the spline) - Add ord

Hi, I would like to fit an (interpolating) spline to data where the derivatives at the endpoints of the interval are nonzero, thus the natural spline endpoint-specification does not make sense. Books (de Boor, etc) suggest that in this case I use not-a-knot splines. I know what not-a-knot splines are (so if I were solving for the coefficients directly I knew how to do this), but I don't

Goodmorning, This is a documentation related question about the B-spline function in R. In the help file it is stated that: "df degrees of freedom; one can specify df rather than knots; bs() then chooses df-degree-1 knots at suitable quantiles of x (which will ignore missing values)." So if one were to specify a spline with 6 degrees of freedom (and no intercept) then a basis

I've got a problem with the R function spline.des. I use the following arguments: x_seq(0,1,length=200) knot_quantile(x,(1:8)/9) ord_3 k_sort(c(rep(range(x),ord,knot))) derivs_rep(2,200) When I do spline.des(k,x,ord,derivs)$design on Splus and on R, I don't have the same results. However, if I take an order ord=4 (ie a spline of degree 3), then, I have the same

Hello, I am using R and the smooth.Pspline function in the pspline package to smooth some data by using natural cubic splines. After fitting a sufficiently smooth spline using the following call: (ps=smooth.Pspline(x,y,norder=2,spar=0.8,method=1) [the values of x are age in years from 1 to 100] I tried to check that R in fact had fitted a natural cubic spline by checking that the resulting

Hi, I'm using natural cubic splines from splines::ns() in survival regression (regressing inter-arrival times of patients to a queue on queue size). The queue size fluctuates between 3600 and 3900. I would like to be able to run predict.survreg() for sizes <3600 and >3900 by assuming that the rate for <3600 is the same as for 3600 and that for >4000 it's the same as for

Hi all, I have a vector of proportions (post_op_prw) such that >summary(amb$post_op_prw) Min. 1st Qu. Median Mean 3rd Qu. Max. NA's 0.0000 0.0000 0.0000 0.3985 0.9134 0.9962 1.0000 > summary(cut2(amb$post_op_prw,0.0001)) [0.0000,0.0001) [0.0001,0.9962] NA's 1904 1672 1

Hi r-helpers. Please forgive my ignorance, but I would like to plot a smoothing spline (smooth.spline) from package "stats", and show the knots in the plot, and I can't seem to figure out where smooth.spline has located the knots (when I use nknots). Unfortunately, I don't know a lot about splines, but I know that they provide me an easy way to estimate the location of local

--AZaLVt6Pw+ Content-Type: text/plain; charset=us-ascii Content-Description: message body text Content-Transfer-Encoding: 7bit Dear all, I was looking at "splinefun" and the underlying C code and believe that there is a memory access error in the C routine "spline_eval". Specifically, on line 368 and following the following code appears: if(ul < x[i] || x[i+1] < ul)

Dear R developers, The trace of the hat matrix H~(n,n) is computed as follows: tr(H) = tr(BS^-1B') = tr(S^-1B'B) := tr(X) = sum(diag(X)) with B~(n,p), S~(p,p). Since p is of the order 10^3 but S is sparse I would like to employ Taucs linear solver ( http://www.tau.ac.il/~stoledo/taucs/ ) on SX = B'B. (Further improvement by implying a looping over i=1,...,p, calling

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