w****g
发帖数: 44 | 1 Hello:
any one know how to obtain the derivative of the cubic spline interpolation
results using a program ?
I am thinking about using a program to find the cubic spline interpolation
and the
derivative of the cubic spline interpolation results based on the input. the
example input file is listed below: a list of potential energies of a pair
particles as a function of the distance
r1 E1
r2 E2
r3 E3
The goal is to find the interpolated potential and force (negative
derivative of the cubic spline |
|
c*******h
发帖数: 1096 | 2 【 以下文字转载自 Computation 讨论区 】
发信人: cockroach (冬冬), 信区: Computation
标 题: spline interpolation
发信站: BBS 未名空间站 (Mon Sep 7 16:35:13 2009, 美东)
有没有哪本书有通用的算spline的系数的公式?
需要大于三次的spline。我见过的数值分析的书上都只有三次spline的公式。高于三次
的推导挺麻烦的 |
|
G***G
发帖数: 16778 | 3 in previous posts, it seems that most people use 0.124 screen spline.
just curious about whether there are anyone else here who used
bigger screen spline such as 0.140 or 0.160?
any experience of how difficult it is to install the bigger screen spline? |
|
c*******h
发帖数: 1096 | 4 有没有哪本书有通用的算spline的系数的公式?
需要大于三次的spline。我见过的数值分析的书上都只有三次spline的公式。高于三次
的推导挺麻烦的 |
|
w****g
发帖数: 44 | 5 Hello:
any one know how to obtain the derivative of the cubic spline interpolation
results using a program ?
I am thinking about using a program to find the cubic spline interpolation
and the
derivative of the cubic spline interpolation results based on the input. the
example input file is listed below: a list of potential energies of a pair
particles as a function of the distance
r1 E1
r2 E2
r3 E3
r4 E4
r5 E5
The goal is to find the interpolated potential and force (negative
derivative of the c |
|
f***a
发帖数: 329 | 6 你的问题看不懂，referent 你这里指什么？
我乱说一些，不对的地方大家指正。
cubic spline 大概就是选好knots，得到一组basis funtions，然后在保证Knots处1st
and 2nd derivatives continuous的情况下用data去fit，解出参数。 再算estimated
response 的standard error，得出C.I.
一般在boundary上用高阶曲线去fit结果variance太大，所以可以限定在boundary周围
是linear的，这样就是natural cubic spline了。
再如果懒得去选knots，干脆就每一个data point都当作knot来算，然后用penalized
residual sum of squares来限制Model complexity，这样就是smoothing spline了。 |
|
t****n
发帖数: 15 | 7 Sorry, should be basis since B-spline is represented as sum of polynomial
basis functions. Other spline may not be represented in this way if I
remember correctly. Google for details. |
|
jl
发帖数: 398 | 8 【 以下文字转载自 Statistics 讨论区 】
发信人: jl (打火石), 信区: Statistics
标 题: 那位能给出 spline的 一般定义
发信站: BBS 未名空间站 (Thu Jul 14 00:26:56 2005), 转信
什么是 spline,
能不能随便解释一下,
为什么要这么的定义. |
|
d********d
发帖数: 497 | 9 需要用Fortran实现多元/多变量Spline插值。
参考了Matlab里面的函数，用得好象是Tensor Product.但是Matlab用了Cell的格式写
的很简洁，但是比较难看懂。
不知道，有没有简单的使用Tensor Product实现Spline插值的程序或者文章。
Thanks! |
|
s******t
发帖数: 150 | 10 发信人: slarlett (木木), 信区: Computation
标 题: 请问B-spline中order 3和basic function的关系
发信站: BBS 未名空间站 (Sat Apr 4 21:51:07 2009)
已知B(x)=1/12(|x-2|^3-4|x-1|^3+6|x|^3-4|x+1|^3+|x+2|^3);
现在要求B-spline function of order 3,怎么求啊？
我是这样做的：就是求B(x)*B(x)*B(x)， 其中*表示convolution
直接用matlab函数conv(B,B)*h_step.
本来出来的结果应该是关于y-axis对称的，可是我出来的呢？根本就不对称。
请问是什么原因啊？是否我理解有问题啊？ |
|
c*****z
发帖数: 182 | 11 How to calculate the total arc length of a b spline curve, is there some
analytic solution. Another dual question is that, suppose I fix a point on a
b spline and I want to find the point that the arc length between these two
points are, say 10, is there analytic solution to that?
Thanks |
|
m****s
发帖数: 402 | 12 能否具体谈谈这个base的意义，
跟非base的spline有什么差别？ |
|
g*****n
发帖数: 56 | 13 假设要FIT 一个正旋波SIN(3X)， 如果用另外一个SIN(BX)去FIT，可以看到B=3处的
RESIDUAL 最小。但是用SMOOTHING SPLINE 去FIT为什么就达不到这个效果呢？请大家指
教。 |
|
s*****l
发帖数: 167 | 14 subroutine seval(n, u, x, y, b, c, d, f, df, d2f)
implicit none
integer n
double precision u, x(n), y(n), b(n), c(n), d(n)
double precision f, df, d2f
!
! this subroutine evaluates the cubic spline function
!
! seval = y(i) + b(i)*(u-x(i)) + c(i)*(u-x(i))**2 + d(i)*(u-x(i))**3
!
! where x(i) .lt. u .lt. x(i+1), using horner's rule
!
! if u .lt. x(1) then i = 1 is used.
! if u .ge. x(n) then i = n is used.
!
! input.. |
|
s******t
发帖数: 150 | 15 已知B(x)=1/12(|x-2|^3-4|x-1|^3+6|x|^3-4|x+1|^3+|x+2|^3);
现在要求B-spline function of order 3,怎么求啊？
我是这样做的：就是求B(x)*B(x)*B(x)， 其中*表示convolution
直接用matlab函数conv(B,B)*h_step.
本来出来的结果应该是关于y-axis对称的，可是我出来的呢？根本就不对称。
请问是什么原因啊？是否我理解有问题啊？ |
|
c*******h
发帖数: 1096 | 16 does anyone know if there's any library for cubic spline interpolation
using mpi? otherwise i'll just implement my own code |
|
y***z
发帖数: 4 | 17 用eviews做cubic spline怎么做呢？
我eviews太烂了，谢谢大家。 |
|
I****M
发帖数: 143 | 18 你用Matlab吗？那个做cubic spline很方便的。 |
|
I****M
发帖数: 143 | 19 那就在help里找spline，怎么用都写得很清楚的，good luck |
|
y***z
发帖数: 4 | 20 是这样的， 我一作业要求estimate money deamnd function for Japan
time range: 1970Q1 2008Q1
regression: log(m/P) = c i^-1 i^-2 i^-3 i i^2 i^3 log(GDP/P) lags
m=M2+CD or MB
i=short term rate
目的在于i=0， money demand 是不是存在。
所有变量都assume I(1), 以上是老师要求的。extension 自己想。
作出来的结果存在以下问题：
1）用mb做，1993之前的interest rate elasticiy 都大于零。
2）用m2+cd做，完全跟 i 没关系。
我觉得可能是因为money demand的方程变了， 于是乎想到了spline。
不过就我观察， 好像现在用这个做money demand的很少。
莫非这个方法不可行？？？
那么有什么其他的方法呢？ |
|
B****g
发帖数: 43 | 21 not difficult to calculate by yourself, the key is to calcualte the B-
spline coefficients, then you can compute gradient
interpolation
the |
|
|
l*******l
发帖数: 204 | 23 anyone familiar with spline regression. I know nothing about it. I try
google and did not find introductory material. Many thanks if you can
provide some introductory material. The question I try to answer is to find
the non-linear relationship between x and y and how this relationship
differs between two groups. |
|
n*****x
发帖数: 9 | 24 Is anyone here familiar with noparametric regression, spline regression
modelling? Especially SAS macro codes?
Thanks very much |
|
n*****x
发帖数: 9 | 25 Does anyone know how to calculate odds ratios of every knots in spline
regression? I know there is a SAS macro could do this. However, I am not
able to find on the website.
Thanks |
|
M****e
发帖数: 178 | 26 对B-spline basis functions 不太理解，希望能找本写得深入浅出的书，了解一下来
龙去脉，谢谢啦。 |
|
W**********E
发帖数: 242 | 27 比方说, binary logistic regression model, 如果只有1个变量X1，很好做出一个y
轴为log_odds的restricted cubic spline（ＲＣＳ）图来表示x1与log odds的
functional form.
但是通常加入其它的变量或confounders后，x1的几个rcs变量系数也会变，但画出这种
调整其它变量的rcs图好像很麻烦，因为不但只有x1的值还有其它变量的值如何确定。
有谁用图形表示过矫正其它变量后的某变量的RCS图？
谢谢， |
|
d******3
发帖数: 93 | 28 看到几篇文献，做subic spline regression并且plot了predicted value和95%
confidence interval (CI)， 但是整条CI曲线上都没有referent
这是怎么回事？谁有好的讲义或tutorial么？谢了 |
|
|
t*********e
发帖数: 143 | 30 Suppose we have two sequence of data (x, y), there are a lot of data points.
How to fit a continuous monotonic piecewise linear spline using SAS? |
|
s***h
发帖数: 487 | 31 比如对于实函数，我总可以把一个常函数，跟一个二次函数，在零点 spline，而由于
一维分配率，成功与商学院系花偶遇。
但复数上不行嘎，复数的两维分配率麻烦了，同时满足分配率和连续性，这 spline 偶
遇似乎很难。得有缘有份咋都不能错过。。。
: 我觉得对于 intuitive 图景，首先我们能自然想象实数函数能够 spline，包括
多变量
: spline 平面。Spline 在多次求导后总是会出现不可导的点。
: 而复函数必须不可能 Spline？
: 样，以
|
|
a*********a
发帖数: 3656 | 32 then the formula I gave should be it. if you don't require crossing the
points, then it is a different matter.
starting from 0th, you can pick an arbitrary tangent. you have now y0, y0'.
then going to the 1st point, the only thing you can do is a quadratic form.
3 equations for 3 unknowns like I laid out before. y0=a x0^2+b x0+c,y1=a x1^
2+b x1+c,y0'=2 a1 x0+b1)
This will fix the tangent at point 1. then you solve anther equation set of
3. (y1=a x1^2+b x1+c,y2=a x2^2+b x2+c,y1'=2 a1 x1+b1). lath... 阅读全帖 |
|
w*******y
发帖数: 60932 | 33 I've been looking for a Bosch Demo Hammer for awhile and trying to find the
right combo of size/power/price that I want. I have a 64lb Brute, but need
something in the handheld range for more delicate breaking. Additionally, I
'd like the capability of drilling some larger holes (than a 1" Makita
Rotary I have can do), but am stuck between going with a strictly demo
hammer (such as the 11316evs) or going with a combo (such as a 11264evs).
This morning I happened to be looking at CPO Bosch's si... 阅读全帖 |
|
p******s
发帖数: 2552 | 34 换纱窗一般要以下几样
screen 纱窗布
screen frame 外面的铝合金框，可以买来用hack saw 裁
screen corner，4个角
spline 橡胶条,自己买注意尺寸
spline roller,把spline挤进窗框缝的工具
如果用现成的窗框和spline，那么只要买screen和spline roller自己换上就好 |
|
q*********n
发帖数: 103 | 35 借楼主的帖子问下三次样条拟合。
Matlab里用的spline函数是not a knot的边界条件，不知道实际效果和自然边界(两头二阶
导数为零)有多大的区别？
还有一种 Kochanek-Bartels Hermite Spline interpolation，能通过调节三个参数
得到不同的拟合效果。这种spline是怎么算的？常用么？
(1) tension: 调节峰的尖锐程度
(2) continuity: 节点两边的平滑度
(3) bias: 好像是关于节点两侧对称程度
http://www.shaneaherne.com/research/splines.html
请问怎么为不同类型的曲线选spline的算法？ 比如，是不是光滑曲线要用自然三次样
条， 有很多峰值的曲线则应该用其他的方法？
多谢多谢~ |
|
R*******c
发帖数: 249 | 36 搞了三四天了，还是没成功，实在是崩溃了。。。只能上来请大牛们帮忙看看问题出在哪？
从网上下载了一个R的package（一个zip文件）,这个包是在mac上build的，在mac里运行没有
问题，我现在想在windows下用这个包，直接通过R里面的install package from local
zip file会有错误如下：
Error in gzfile(file, "r") : cannot open the connection
In addition: Warning messages:
1: In unzip(zipname, exdir = dest) : error 1 in extracting from zip
file
2: In gzfile(file, "r") :
cannot open compressed file 'curves2_1.0.1.tar.gz/DESCRIPTION',
probable reason 'No such file or directory'
所以决定先将zip文件解压缩成文件夹，然后在windows下把这个包buil... 阅读全帖 |
|
s***h
发帖数: 487 | 37 我觉得对于 intuitive 图景，首先我们能自然想象实数函数能够 spline，包括多变量
spline 平面。Spline 在多次求导后总是会出现不可导的点。
而复函数必须不可能 Spline？
: 而复数可导就得在分配率上有缘有份？不管从哪条路径走过去，复斜率要总是一
样，以
: 至于此生注定要跟美腿系花偶遇在那个外企写字楼？
: 不过还是没有解释怎么得出那就一定就高阶可导的图景？
: 现我和
|
|
d*j
发帖数: 756 | 38 刚刚看了这个链接关于 splinefit的例子：
http://www.mathworks.com/matlabcentral/fileexchange/13812-splin
请教下，里面的break 和 order 是啥意思啊？
就是这个例子里面的：
% x = linspace(0,2*pi,100);
% y = sin(x) + 0.1*randn(size(x));
% % Breaks
% breaks = [0:5,2*pi];
%
% % Fit a spline of order 5
% pp = splinefit(x,y,breaks,5);
%
% % Fit a spline of order 3 with periodic boundary conditions
% pp = splinefit(x,y,breaks,3,'p');
%
% % Constraints: y(0) = 0, y'(0) = 1 and y(3) + y"(3) = 0... 阅读全帖 |
|
d*j
发帖数: 756 | 39 刚刚看了这个链接关于 splinefit的例子：
http://www.mathworks.com/matlabcentral/fileexchange/13812-splin
请教下，里面的break 和 order 是啥意思啊？
就是这个例子里面的：
% x = linspace(0,2*pi,100);
% y = sin(x) + 0.1*randn(size(x));
% % Breaks
% breaks = [0:5,2*pi];
%
% % Fit a spline of order 5
% pp = splinefit(x,y,breaks,5);
%
% % Fit a spline of order 3 with periodic boundary conditions
% pp = splinefit(x,y,breaks,3,'p');
%
% % Constraints: y(0) = 0, y'(0) = 1 and y(3) + y"(3) = 0... 阅读全帖 |
|
h*********n
发帖数: 278 | 40 c-value 又是啥？查了一下，说是equivalent to roc, 这个前面就说了呀，85%, 被认
为太高了呀。
spline我不太懂, 原以为比较高深的东西，同事说有别的组就是这么做的(也是GLM
model)，我们就照搬过来了。基本上就是这样code:
if x=a, x1=x;
if x>a, x2=a, if x<=a, x2=x;
能不能讲讲怎么个无厘头法？当时我跟他们说这东西可能没那么简单，貌似被嘲笑了。
不过后来我咨询另外组的一个同胞，好象真是这么做的。话说回来，即使没有spline，
model一样fit得很好，spline就是把本来有一小段fit得不算非常好的也变成非常好了
，我说了，当时被push追求这个。所以即使没有spline，还是too good to be true.
x没有跟y一起的derived.
。。 |
|
e**y
发帖数: 51 | 41 陈不像是民科。
关于spline model的缺陷的讨论是对的，陈的这一点没错。大部分民科是不会站在一个
已知的理论缺陷上来提出自己的观点的。后面的修正部分，我很赞赏用概率推断的办法
去找临界点的思路。这一点，就我所知，在spline model里面还没有人讨论过。但是在
工程应用上，概率推断的办法已经广泛应用了。从这一点来说，陈的思路也并非非常新
颖。除非，陈能像whaba那样，建立一个reproducing kernel hilbert space的理论来
justify spline model一样来justify他的模型，那今天的这些所谓争论就可以归结为
粪土了。
陈虽然不是统计学家，但是就他自己所说，“我对统计学仅仅懂得一丁点的皮毛，尽管
我在这个领域干了20多年了，而且我的经历从数据输入员到现场调查员到调查设计者、
组织者和管理者再到数据库的建立与维护者最后到统计分析人员等几乎贯穿了一个应用
统计研究的全过程。”我完全相信，这个经历是板上绝大部分人所没有的。绝大部分人
甚至都不知道数据是怎么一个收集来的。而绝大部分人甚至还在相信，统计学的最有价
值部分在后面的分析，而不是前... 阅读全帖 |
|
e**y
发帖数: 51 | 42 如果你有过审稿经历，你就会知道大部分审稿人只看introduction，根本用不着看后面
的细节。我就只看懂了他的introduction，我就觉得够了。
另外，我不知道你是否真懂spline或者piecewise linear。陈的东西几乎就是讲spline
上的knot。请你好好翻翻jasa，看看讨论spline上knot估计的文章有多少。我不觉得陈
会去看jasa，但是他提出的问题确实实实在在的正确。 |
|
g********r
发帖数: 8017 | 43 OMG。你审稿不用负责推导是否正确，结论是否成立么？
我确实对spline理解肤浅。上课学过不做研究。但是你如果看了陈的文章，你会认为他
是从spline的思路出发的么？
spline |
|
r*****d
发帖数: 346 | 44 受教了。
一直以为decision tree就是decision tree, 原来CART是一种，CHAID也是一种。
还有就是spline regression,
from wiki, "I-splines can be used as basis splines for regression analysis
and data transformation when monotonicity is desired (constraining the
regression coefficients to be non-negative for a non-decreasing fit, and non
-positive for a non-increasing fit)." |
|
m*y
发帖数: 1778 | 45 大概搜了一下，Sequoia 和LC platform 和四驱好像还是有区别
LC platform 是单独系列，和其它款都不通用，frame 比唐抓和sequoia 要厚。
四驱的话，后轴一样，LC的前轴，前桥都要大很多，真正前轴能当后轴使的。30
spline vs 27 spline， CVs 和 front ring gear 都要大。
总之sequoia 是向full size truck 看齐， LC 比Full size 要更heavy duty
厂家对可靠性的要求也不一样 |
|
m********8
发帖数: 938 | 46 嗯，用的0.125的spline，我一般买套件，spline和FRAME都是配套的。 |
|
G***G
发帖数: 16778 | 47 which kind of screen spline do you use?
I assume you use 0.140.
0.125 spline is easy to push in.
thanks. |
|
p****o
发帖数: 1340 | 48 the problem lies in her spline, C6-7 and C1-2. she did MRI and
X-ray: they showed scoliosis -- minor according to the doc.
the spline also affects her muscles: her neck movement
is constrained, dizziness is now a common place, her left
arm is extremely weak, and the muscles around her chest
is tight too.
the massage does help a lot; temporary relief, i agree.
she is now under regular accupucture treatment, and we
are looking for a good accupucturist...
any recommendation? |
|
T*******4
发帖数: 13143 | 49 丫丫大厨要把food版的名ID们都吸引到花版不肯走了
在 SPLINE (spline) 的大作中提到: 】 |
|
z******o
发帖数: 1224 | 50 赞三多，看到音量线了，但图上的音量线怎么是step function，不是linear
spline啊，CE里面的音量包络线是spline啊。待会有空我再试试。多谢！ |
|
