search for: norder

Dear all, I'm trying to use the Pspline add-on package to fit a quintic spline (norder =3), but I keep running into a Singularity error. > traj.spl <- smooth.Pspline(time, x, norder=3 ) Error in smooth.Pspline(time, x, norder = 3) : Singularity error in solving equations > Playing around with the other parameters produces an "unused arguments" error: > tra...

...39;t know whether this is really a bug with the Pspline package or only a problem with my installation. Things work fine in Linux but not in Mac 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...

Hi All, I want to fit a piecewise spline of degree 1, i.e. a spline consisting of a straight line over each piece. I downloaded the R package pspline, then I have following questions: 1. in the program, the degree of the spline is specified by 2*norder-1. Why do they adopt such scheme that we can only fit a spline with odd degree? 2. norder cannot be set to 1. Is there any specific reason for doing so? Maybe I need to code this by myself. So I'd like to hear your advices. Thanks in advance, Gang _________________________________________...

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

...e n=4476, but not for anything larger (say, n=4477). For example: THIS WORKS: ----------------------------- random = array(0,c(4476,2)) random[,1] = runif(4476,0,1) random[,2] = runif(4476,0,1) random = random[order(random[,1]),] plot(random[,1],random[,2]) fit2 = sm.spline(random[,1],random[,2], norder=2, cv=FALSE) lines(fit2$x,fit2$y) THIS FAILS: ----------------------------- random = array(0,c(4477,2)) random[,1] = runif(4477,0,1) random[,2] = runif(4477,0,1) random = random[order(random[,1]),] plot(random[,1],random[,2]) fit2 = sm.spline(random[,1],random[,2], norder=2, cv=FALSE) lines(fit2$x...

...e "Pspline" 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 a...

I'm a Beginner on R-language and I' have tried to use some packages but the following message have appeared: "Error in .Fortran("pspline", as.integer(n), as.integer(nvar), as.integer(norder), : C/Fortran function name not in load table" I'm using a Windows version and I'd like to know if it is a installation problem or problems with paths? Thanks, Marlon!!! -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Rea...

...its derivative" (https://stat.ethz.ch/pipermail/r-help/2004-July/052376.html) the following code could be use to estimate (and plot) a smooth baseline hazard: > t<-seq(min(data$start),max(data$stop),length=100) > lines(t, predict(sm.spline(x=basehaz(coxfit)[,2],y=basehaz(coxfit)[,1],norder=2), t,1)) #there is a problem with this code. One should add the contraint that the baseline hazard cannot be negative. The following computes the parametric part of the cox model. > risk <- predict(coxfit, type='risk') # gives exp(X'b) something like > baseline.hazard*risk w...

...ave more useful defaults and they return an object of class 'fdSmooth', with methods for 'plot', 'lines', and 'plotfit'. In addition, there have also been improvements to the 'plotfit' function and the 'CanadianWeather' data. Also, a function 'norder' has been added, which returns the order of a B-spline. Spencer Graves _______________________________________________ R-packages mailing list R-packages at r-project.org https://stat.ethz.ch/mailman/listinfo/r-packages

...oth 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) plot(Alpha) Thanks a lot. Best, Rafa...

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 spline was LINEAR outside the knots. I did this by plotting the predicted values from the spline fitting in the following way: pl...

...ave more useful defaults and they return an object of class 'fdSmooth', with methods for 'plot', 'lines', and 'plotfit'. In addition, there have also been improvements to the 'plotfit' function and the 'CanadianWeather' data. Also, a function 'norder' has been added, which returns the order of a B-spline. Spencer Graves _______________________________________________ R-packages mailing list R-packages at r-project.org https://stat.ethz.ch/mailman/listinfo/r-packages

...= 0 munmap(0x7f4310e06000, 4096) = 0 open("/etc/host.conf", O_RDONLY|O_CLOEXEC) = 3 fstat(3, {st_mode=S_IFREG|0644, st_size=63, ...}) = 0 mmap(NULL, 4096, PROT_READ|PROT_WRITE, MAP_PRIVATE|MAP_ANONYMOUS, -1, 0) = 0x7f4310e06000 read(3, "#\n# /etc/host.conf\n#\n\norder host"..., 4096) = 63 read(3, "", 4096) = 0 close(3) = 0 munmap(0x7f4310e06000, 4096) = 0 getpid() = 30534 open("/etc/resolv.conf", O_RDONLY|O_CLOEXEC) = 3 fstat(3, {st_mode=S_IFRE...

Hi, I'm trying to get the so-called "single instance store" (I think cyrus has got the name for the first time) with dovecot --version = 2.0.19 binary package installed from ubuntu 12.04 lts official repo. I have checked that "maildir_copy_with_hardlinks" is enabled ("dovecot -a|grep hard" shows "yes") then I have installed and enabled the lmtp

I got a little way into this. The two attached patches are preliminary work. My proposed split is: libguestfs.git common -> git submodule libguestfs-common.git generator/ lib/ all language bindings C based tools (eg. virt-df, virt-edit, guestfish) guestfs-tools.git common -> git submodule libguestfs-common.git