Displaying 20 results from an estimated 2000 matches similar to: "smooth.spline"
2009 Jan 17
2
DierckxSpline segfault
I've just encountered a segfault when using DierckxSpline::percur
function. Below is the minimal example which triggers the error:
---
library(DierckxSpline)
x <- 1:10
y <- rep(0, 10)
pspline <- percur(x, y)
---
*** caught segfault ***
address (nil), cause 'memory not mapped'
Traceback:
1: .Fortran("percur", iopt = as.integer(iopt), m = as.integer(m),
x =
2009 Apr 03
1
DierckxSpline fitting with different sets of y-values in one time
Dear "R" users,
I have a question about the Package DierckxSpline. I have tried to find the answer by myself but it didn't worked out.
I wondered if Dierckxspline can use different sets of y values in one time to fit a line with knot. I have different sets of Y values representing the same thing for different voxels (in an fmri image). I have already fitted the data in different
2012 Feb 24
1
B-spline/smooth.basis derivative matrices
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
2012 Mar 12
1
Fwd: Re[2]: B-spline/smooth.basis derivative matrices
--- 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
>
>
>
> -------- ???????????? ?????????
> --------
> ?? ????:
2017 Jun 01
1
Function with multi return path?
Hi folks,
I'm new to LLVM and non-expert in hardware architecture, so my question
might be dumb.
I'm looking for the possibility to define/call a function with multi return
path, it might look like this:
```
; Define a function with 1 alternate return path
define i32 @f(i1 %i) fork 1 {
br i1 %i, label %noraml, label %alternate
noraml:
...
setpath ret i32 42 ; take normal return
2005 Oct 03
3
spline.des
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
2012 Feb 24
1
evaluating derivative matrices of spline functions
Hello,
I've noticed that SPLUS has a function for evaluating the derivative matrices of splines.
I've also noticed that there's a function in R for evaluating matrices from smooth.spline;
maybe someone knows if something has been written to evaluate matrices from 'smooth.basis'?
regards,
s
2006 Nov 07
1
multivariate splines
Hi,
I am looking for an R package that would calculate multivarite (mostly
2d and 3d, tensor) cubic interpolating splines, so that I could
evaluate these splines (and their derivatives) at many points (unkown
at the time of calculating the spline polynomials) repeatedly.
To make things concrete, I have an array V with
dim(V) = k
and gridpoint vectors grid=list(...), length(grid[[i]])==k[i],
2008 Jul 29
1
tensor product of equi-spaced B-splines in the unit square
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
2006 Dec 13
2
caching frequently used values
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
2012 Aug 02
2
Rd] Numerics behind splineDesign
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
2016 Mar 04
2
R 3.2.4 rc issue
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)
*
2006 Apr 28
1
gauss.quad.prob
I've written a series of functions that evaluates an integral from -inf to a or b to +inf using equally spaced quadrature points along a normal distribution from -10 to +10 moving in increments of .01. These functions are working and give very good approximations, but I think they are computationally wasteful as I am evaluating the function at *many* points.
Instead, I would prefer to use
2008 May 01
1
Optimal knot locations for splines
Suppose I have two variables, x and y. For a fixed number of knots, I want
to create a spline transformation of x such that a loss function is
minimized. Presumably, this loss function would be least squares, i.e. sum
(f(x)-y)^2. The spline transformations would be linear, quadratic or
cubic. I know I can solve this problem using some optimization function in
R, but I was wondering if anyone
2006 Nov 15
1
splineDesign and not-a-knot conditions
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
2008 May 06
2
gfortran: Command not found
Hi, All:
How do I get information about the "R CMD check / build / install"
process, especially regarding what software to install and how to
configure it under Windows XP?
I ask, because "R CMD check" complained, "gfortran: Command not
found". I found 'gfortran-sjlj.exe' in 'C:\Program
Files\R\Rtools\MinGW\bin', which is in the
2003 May 23
2
predict.smooth.spline
I'm using R 1.7.0 on linux. With this version of R the package modreg is
automatically loaded at start of session. However attempting to use
predict.smooth.spline() produces Error: couldn't find function
predict.smooth.spline.
The function smooth.spline() is OK. What am I missing?
======================================
I.White
ICAPB, University of Edinburgh
Ashworth Laboratories, West
2011 Aug 25
3
Application of results from smooth.spline outside R
Hi,
I want to use the result from smooth.spline outside R.
I take my data ,which is 180 point stored in x and y
s <- smooth(x,y)
I can know use to e.g. find the interpolated value at e.g. x=500
predict (s,500)
My problem is, that i don't know how to implement the predict function. I
have looked at literature, but i cannot connect the output of the
smooth.spline() to an actual spline
2006 Mar 17
1
smooth.spline
I have noticed a slightly puzzling behaviour exhibited by
smooth.spline(). If I do
sss <- smooth.spline(x,y)
for a certain pair of data vectors x and y, and then do
length(sss$x)
I get the result ``18''. However if I do
length(unique(x))
I get ``27''. Trying to force smooth.spline() to use more knots I
tried
sss <- smooth.spline(x,y,all.knots=TRUE)
but again
2006 Jan 20
1
Hardwiring a Tellabs echo canceller - help req
Hi All, Greg has been a huge help getting me going with this tellabs echo
can, but I'm still having some problems getting it to work... I suspect I
wired it up incorrectly, so I thought I'd see if anyone can point me in the
right direction. Digium tech support pointed me to this doc for a standard
T1 cable:
http://www.arcelect.com/RJ48C_and_RJ48S_8_position_jack_.htm
which looks like it