Displaying 20 results from an estimated 7000 matches similar to: "Periodic B-spline surface"
2006 Feb 27
3
how to use the basis matrix of "ns" in R? really confused by multi-dim spline filtering?
Hi all,
Could anybody recommend some easy-to-understand and example based
notes/tutorials on how to use cubic splines to do filtering on
multi-dimension data?
I am confused by the 1-dimensional case, and more confused by
multi-dimensional case.
I found all the books suddenly become very abstract when it comes to this
subject.
They don't provide examples in R or Splus at all.
Specifically,
2011 Sep 20
2
Multivariate spline regression and predicted values
Hello,
I am trying to estimate a multivariate regression of Y on X with
regression splines. Y is (nx1), and X is (nxd), with d>1. I assume the
data is generated by some unknown regression function f(X), as in Y =
f(X) + u, where u is some well-behaved regression error. I want to
estimate f(X) via regression splines (tensor product splines). Then, I
want to get the predicted values for some new
2003 Apr 07
3
spline with multiple predictor vars?
Hi, is there a way in R to generate a polynomial spline with multiple predictor
variables? I have one response and two predictors and I'm trying to fit a
spline model for this...
Please cc me on the reply..
Thanks,
nirmal
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
2010 May 14
1
Cubic B-spline, how to numerically integrate?
(corrected version of previous posting)
I fit a GAM to turtle growth data following methods in Limpus & Chaloupka
1997 (http://www.int-res.com/articles/meps/149/m149p023.pdf).
I want to obtain figures similar to Fig 3 c & f in Limpus & Chaloupka
(1997), which according to the figure legend are "expected size-at-age
functions" obtained by numerically integrating
2009 Oct 13
2
How to choose a proper smoothing spline in GAM of mgcv package?
Hi, there,
I have 5 datasets. I would like to choose a basis spline with same knots in
GAM function in order to obtain same basis function for 5 datasets.
Moreover, the basis spline is used to for an interaction of two covarites.
I used "cr" in one covariate, but it can only smooth w.r.t 1 covariate. Can
anyone give me some suggestion about how to choose a proper smoothing spline
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 Oct 10
1
how to evaluate a cubic Bezier curve (B-spline?) given the four control points
I'm trying to use R to determine the quality of a cubic Bezier curve
approximation of an elliptical arc. I know the four control points
and I want to compute (x,y) coordinates of many points on the curve.
I can't find anything in either the base distribution or CRAN that
does this; all the spline-related packages seem to be about *fitting*
piecewise Bezier curves to a data set.
2009 Aug 24
2
Formulas in gam function of mgcv package
Dear R-experts,
I have a question on the formulas used in the gam function of the mgcv
package.
I am trying to understand the relationships between:
y~s(x1)+s(x2)+s(x3)+s(x4)
and
y~s(x1,x2,x3,x4)
Does the latter contain the former? what about the smoothers of all
interaction terms?
I have (tried to) read the manual pages of gam, formula.gam, smooth.terms,
linear.functional.terms but
2009 Aug 24
2
Formulas in gam function of mgcv package
Dear R-experts,
I have a question on the formulas used in the gam function of the mgcv
package.
I am trying to understand the relationships between:
y~s(x1)+s(x2)+s(x3)+s(x4)
and
y~s(x1,x2,x3,x4)
Does the latter contain the former? what about the smoothers of all
interaction terms?
I have (tried to) read the manual pages of gam, formula.gam, smooth.terms,
linear.functional.terms but
2005 Dec 12
2
Bivariate Splines in R
Hi..,
is there a function in R to fit bivariate splines
?
I came across 'polymars' (POLSPLINE) and 'mars' (mda)
packages. Are these the one to use or are there other
specific commands?
Thanks.
Harsh
2006 Nov 07
1
gamm(): nested tensor product smooths
I'd like to compare tests based on the mixed model representation of additive models, testing among others
y=f(x1)+f(x2) vs y=f(x1)+f(x2)+f(x1,x2)
(testing for additivity)
In mixed model representation, where X represents the unpenalized part of the spline functions and Z the "wiggly" parts, this would be:
y=X%*%beta+ Z_1%*%b_1+ Z_2%*%b_2
vs
y=X%*%beta+ Z_1%*%b_1+ Z_2%*%b_2 + Z_12
2009 Jun 23
1
Model fitting with GAM and "by" term
Hello R Users,
I have a question regarding fitting a model with GAM{mgcv}. I have data
from several predictor (X) variables I wish to use to develop a model to
predict one Y variable. I am working with ecological data, so have data
collected many times (about 20) over the course of two years. Plotting
data independently for each date there appears to be relationships
between Y (fish density)
2012 Jul 30
1
te( ) interactions and AIC model selection with GAM
Hello R users,
I'm working with a time-series of several years and to analyze it, I?m using
GAM smoothers from the package mgcv. I?m constructing models where
zooplankton biomass (bm) is the dependent variable and the continuous
explanatory variables are:
-time in Julian days (t), to creat a long-term linear trend
-Julian days of the year (t_year) to create an annual cycle
- Mean temperature
2013 Mar 11
1
Use pcls in "mgcv" package to achieve constrained cubic spline
Hello everyone,
Dr. wood told me that I can adapting his example to force cubic spline to pass through certain point.
I still have no idea how to achieve this. Suppose we want to force the cubic spline to pass (1,1), how can
I achieve this by adapting the following code?
# Penalized example: monotonic penalized regression spline .....
# Generate data from a monotonic truth.
2009 Jan 16
2
Smooth periodic splines
Hello group!
Is there a package that allows to fit smooth *periodic* splines to
data? I'm interested in a function which combines the functionality of
smooth.spline and splines::periodicSpline.
Thanks,
Andrey
2007 May 21
2
comparing fit of cubic spline
I want to compare the fit of a quadratic model to continuous data, with that
of a cubic spline fit. Is there a way of computing AIC from for e.g. a GAM
with a smoothing spine, and comparing this to AIC from a quadratic model?
Cheers
******************************************
Tom Reed
PhD Student
Institute of Evolutionary Biology
102 Ashworth Laboratories
Kings Buildings
University of
2010 Jun 04
1
package mgcv inconsistency in help files? cyclic P-spline "cs" not cyclic?
Dear all,
I'm a bit stunned by the behaviour of a gam model using cyclic
P-spline smoothers. I cannot provide the data, as I have about 61.000
observations from a time series.
I use the following model :
testgam <- gam(NO~s(x)+s(y,bs="cs")+s(DD,bs="cs")+s(TT),data=Final)
The problem lies with the cyclic smoother I use for seasonal trends.
The variable Final$y is a
2010 Aug 20
1
Surface generated as product of two curves
I have the following data table:
[,1] [,2] [,3] [,4] [,5]
[1,] 2575 1927 1754 581 354
[2,] 1156 810 730 541 237
[3,] 417 297 199 125 110
[4,] 281 132 132 47 58
[5,] 152 84 87 54 19
Each column represents a time interval expressed in nano-seconds
Each row represents an energy interval expressed m=n MeV
Each table value represents the average gamma rays intensity
2006 Aug 31
0
periodic spline in glm
I have been trying to find a command similiar to bs or ns to add a
periodic/or cyclical splie term to a regression model in glm.nb
I have been using ns to fit a spline to day of the year (1:365):
glm.nb(CNT ~ WKDY + ns(DY,df=6) + HOLIDAY_FLAG + Trend)
but I think a periodic spline might be more appropriate. Any suggestions
would be appreciated.
best,
Spencer Jones, M.Stat.
NLM Fellow
Dept.