根据R/S分析法算法编写的rsFit函数来求Hurst指数-经管之家官网！

坛友互助群

扫码加入各岗位、行业、专业交流群

rsFit

function (x, levels = 50, minnpts = 3, cut.off = 10^c(0.7, 2.5),

doplot = FALSE, trace = FALSE, title = NULL, description = NULL)

{

call = match.call()

data = list(x = x)

x = as.vector(x)

n = length(x)

increment = (log10(n/minnpts))/levels

M = floor(10^((1:levels) * increment))

M = M[M > 1]

Y = cumsum(x)

Y2 = cumsum(x * x)

RS = NULL

for (m in M) {

S = sqrt(Y2[m]/m - (Y[m]/m)^2)

Z = Y[1:m] - (1:m) * Y[m]/m

STATS = (max(Z) - min(Z))/S

RS = c(RS, STATS)

if (trace)

cat("\n\tm = \t", m, "\tR/S = \t", STATS)

}

if (trace)

cat("\n")

这是重标极差法（R/S分析法）的算法程序，下图为所依据的算法步骤的一部分：

http://ask.pinggu.org/data/attach/1512/kZy7rS9t.png首先，我对照算法来看程序，不知求所有子区间a的极差的均值这一步在程序哪里体现；

其次，RS = c(RS, STATS)和if (trace)以及cat函数这两步是什么意思；

还有，我用数据运行过，发现S,Z显示都是null，且报错如下，

Error in if (trace) cat("\n\tm = \t", m, "\tR/S = \t", STATS) :

argument is not interpretable as logical；Error in if (trace) cat("\n") : argument is not interpretable as logical

本人刚接触编程，还想根据经典R/S法来自己编写修正的R/S法和V/S法用于论文，但一上来就被卡住了，希望各位高手帮忙解答下，如能直接提供以上两个程序，必有重谢，非常感谢！

扫码或添加微信号：坛友素质互助

「经管之家」APP：经管人学习、答疑、交友，就上经管之家！
免流量费下载资料----在经管之家app可以下载论坛上的所有资源，并且不额外收取下载高峰期的论坛币。
涵盖所有经管领域的优秀内容----覆盖经济、管理、金融投资、计量统计、数据分析、国贸、财会等专业的学习宝库，各类资料应有尽有。
来自五湖四海的经管达人----已经有上千万的经管人来到这里，你可以找到任何学科方向、有共同话题的朋友。
经管之家（原人大经济论坛），跨越高校的围墙，带你走进经管知识的新世界。
扫描下方二维码下载并注册APP