similar to: a question about spline A(ns or bs)

"A. Tsiolakidis" <tsioltas@cs.keele.ac.uk> writes: > OK, thanks, it looks like I can use the splines library instead (I had no > idea how, but now I do!). Should we call it a bug then? > > Question No. 2: > > >From the same S-plus program, I was forced to change a line that was giving > me Segmentation Error problems (a crash on Solaris 2.5): > >

I am fittling a spline to a variable in a regression model, I am then using the predict.glm funtion to make some predictions. When I use bs to fit the spline I don't have any problems using the predict.glm function however when I use ns I get the following error: Error in model.frame(formula, rownames, variables, varnames, extras, extranames, : variable lengths differ (found for

BUG IN SPLINE()? Version R-1.0.1, system i486,linux If the spline(x,y,method="natural") function is given values outside the range of the data, it does not give a warning. Moreover, the extrapolated value reported is not the ordinate of the natural spline defined by (x,y). Example. Let x <- c(2,5,8,10) and y <- c(1.2266,-1.7606,-0.5051,1.0390). Then interpolate/extrapolate with

Hi, I am a total beginner with this whole thing so please have patience! I am trying to run an S-plus program with a certain line: spline(1:nrow(y), y[,1],n=100); This crashes with: Error: NAs in foreign function call (arg 8) Apparently, this is caused by the last command of spline: u <- seq(xmin, xmax, length.out = n) .C("spline_eval", z$method, length(u), x = u, y =

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

Dear Sir/Madam: this is shuangge Ma, graduate student in UW-Madison statistics department. I have a computational question. I have a function f(x,y). I want to find the y(x) that maximize f(x,y) under the constraint y(x) is a non-decreasing step function. Is there any R package or algorithm I can use for this purpose? thanks a lot for your time and help, Sincerely, Shuangge Ma

--- 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 > > > > -------- ???????????? ????????? > -------- > ?? ????:

Hi, I am looking at the change in N concentration in plant roots over 4 time points and I have fit a spline to the data using ns and lme: fit10 <- lme( N~ns(day, 3), data = rcn10G) I may want to adjust the model a little bit, but for now, let's assume it's good. I get output for the fixed effects: Fixed: N ~ ns(day, 3) (Intercept) ns(day, 3)1 ns(day, 3)2 ns(day, 3)3 1.15676524

Hi, I am looking at the change in N concentration in plant roots over 4 time points and I have fit a spline to the data using ns and lme: fit10 <- lme( N~ns(day, 3), data = rcn10G) I may want to adjust the model a little bit, but for now, let's assume it's good. I get output for the fixed effects: Fixed: N ~ ns(day, 3) (Intercept) ns(day, 3)1 ns(day, 3)2 ns(day, 3)3

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

mtext(..., outer = TRUE) does not seem to center the text in the outer margins as it did previously. Reproduce with example(plot.profile.nls, package = "nls") On versions of the development sources before November 17 this produced main and subtitles that were centered horizontally on the page. After that the main title (from mtext("..", side = 3, outer = TRUE)) is centered

In the following case, library(splines) tt <- c(55, 251, 380, 289, 210, 385, 361, 669) nn <- rep(0:7,tt) ns(nn,4) ## knots are located at (0.25,0.5,0.75); quantiles = (2,5,7) we get the error Error in qr.default(t(const)) : NA/NaN/Inf in foreign function call (arg 1) because the 75th quantile (the location of the last "interior" knot) ends up on the boundary. As a

Dear R'ers, I've scanned available documentation and the web, but I can't seem to figure out where I've gone wrong in adding numbers to scale the colorkey in levelplot (Lattice package). For example, levelplot(z, contour=T, labels=T, cuts=10, region=T,... colorkey=list(T, space="bottom", at=seq(0,2000,length=20),

R folks, I can't seem to find the instructions in the help files for the lattice package that explain how to add lines, such as with lines() or ?, to a levelplot. I'd be grateful if someone could point me in the proper direction. Cheers, Jeff ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ Jeff Jorgensen Center for Limnology jcjorgensen at wisc.edu

Does anyone know what could be causing the following error message on startup of R: Error in .addMethodFrom(def,argName[1],class[1],fromClass) : object "*tmp*" not found. I'm using R 1.7.0 and Windows 2000. This happened shortly after I installed 1.7.0. Once it occurs, R freezes. >From then on R will always freeze with this error message. I reinstalled 1.7.0 and it worked

As described below, external network connections to the campus of the University of Wisconsin - Madison will be either inaccessible or slow from 20:00 CST 2000/01/03 until 01:00 CST 2000/01/04. CST is our local time zone: Central Standard Time in North America. These local times correspond to 04:00 - 07:00 GMT Tuesday 2000/01/04. This planned outage will affect the r-project.org, omegahat.org,

Sorry to bring up this old topic again but ... The 2.0 release of Debian Linux is about to happen and the release team requires copyright statements on any code included in the release. Presently there are 4 packages that I have created from R source code and libraries and contributed to the release. I think this will help with the visibility of R when there is a well-regarded release of an

# I need to generate some data. I'm modeling some time series that follow a # negative exponential decay (mostly). I have 20 samples that can easily be fit with cubic splines. # What I want to do is generate many thousands of similar samples using the parameters from the splines # For instance one data sample looks not unlike this: foo.curve <- 1 * exp(-0.01 * 1:500) + 0.5

I want to fit a series of lme() regression models that differ only in the degrees of freedom of a ns() spline. I want to use a wrapper function to do this. The models will be of the form y ~ ns(x, df=splineDF) where splineDF is passed as an argument to a wrapper function. This works fine if the regression function is lm(). But with lme(), I get an error. fitfunction() below demonstrates this.

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