全部版块 我的主页
论坛 数据科学与人工智能 数据分析与数据科学 R语言论坛
6070 12
2007-11-19

不同于平时的指数平滑,系数为一个阿尔法,我的是平滑系数为logistic funtion,这种情况一般用Newton-Raphson算法,但是我在实际操作中发现设定初值后,结构非常复杂,没法进行下去,不知道说清楚没有

哪位大虾知道有什么解决的方法 用R语言怎么实现呀

[此贴子已经被作者于2007-11-19 20:24:08编辑过]

二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

全部回复
2007-11-19 21:13:00

你还真说对了,我还不太清楚题意.

不过不管求算几个参数,假设五个,

只要你能明确定出functionsderivative functions

#function values

fun = function (x)

{

f = array(0,c(…))

f[1] =…

f[2] =…

f[3] =…

f[4] =…

f[5] =…

f

}

#derivative values,Jacobian Matrix 5x5

Jac = function (x)

{

j = Matrix 5x5

j

}

给个初始值,就能用Newton.Raphson求解

result=Newton.Raphson (fun, Jac, x0=c(…), tol=1.0e-06, max.it=50)

二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

2007-11-19 22:07:00

[灌水]顶 顶

我说的结构复杂就是指导数函数及其复杂。平滑系数是logistic function 

也就是说参数在分母上,求导后变得极其复杂。不知道大虾的邮箱我可以把我看的论文发给您

二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

2007-11-20 08:15:00

似乎有拉近距离,较了解题意了.

在作STAR model,

有用到smooth transition function,

常用的有logistic,exponential,cumulative dist.function

你的是不是类似这个?

二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

2007-11-20 09:30:00

[灌水]hehe

hehe,差不多吧

我在做稍微复杂一点的波动率模型,因此比较复杂

怎么办呀,大虾 呵呵 超郁闷

二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

2007-11-20 10:39:00

建议你先安装 Package "tsDyn",

 Time series analysis based on dynamical systems theory

先行参考.

而且根据你的描述似乎应该用optim()

optim(par, fn, gr = NULL, ...,

      method = c("Nelder-Mead", "BFGS", "CG", "L-BFGS-B", "SANN"),

      lower = -Inf, upper = Inf,

      control = list(), hessian = FALSE)

二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

点击查看更多内容…
相关推荐
栏目导航
热门文章
推荐文章

说点什么

分享

扫码加好友,拉您进群
各岗位、行业、专业交流群