Question

我完全没有运气通过RSRuby和R使用HoltWinters功能。你究竟如何1）通过RSRuby创建一个时间序列对象，2）成功地将它传递给HoltWinters以获得输出？

示例：

@r = RSRuby.instance
=> #<RSRuby:0x106bfe6c0 @proc_table={}, @class_table={}, @default_mode=-1, @cache={"get"=>#<RObj:0x106bfe580>, "helpfun"=>#<RObj:0x106bfd3d8>, "help"=>#<RObj:0x106bfd3d8>, "NaN"=>NaN, "FALSE"=>false, "TRUE"=>true, "F"=>false, "NA"=>-2147483648, "eval"=>#<RObj:0x106bfdf18>, "T"=>true, "parse"=>#<RObj:0x106bfe0d0>}, @caching=true>
@r.assign('mytime',@r.ts(:data => [1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29,30,31,32,33,34], :frequency => 12, :start => [1993,3], :end => [1995,3]))
=> [1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25]
@r.HoltWinters(@r.mytime)
    RException: Error in decompose(ts(x[1L:wind], start = start(x), frequency = f), seasonal) : 
    time series has no or less than 2 periods

rsruby（0.5.1.1）

R版本2.12.2（2011-02-25）

平台：x86_64-apple-darwin9.8.0 / x86_64（64位）

：编辑：只有R内部的类似示例...如果我可以通过RSRuby从HoltWinters获得任何输出（除了错误）我会非常高兴

> z <- ts(1:34, frequency = 12, start = c(1993,3), end = c(1995,3))
> z
     Jan Feb Mar Apr May Jun Jul Aug Sep Oct Nov Dec
1993           1   2   3   4   5   6   7   8   9  10
1994  11  12  13  14  15  16  17  18  19  20  21  22
1995  23  24  25                                    
> HoltWinters(z)
Holt-Winters exponential smoothing with trend and additive seasonal component.

Call:
 HoltWinters(x = z) 

Smoothing parameters:
 alpha:  1 
 beta :  0 
 gamma:  0 

Coefficients:
             [,1]
a    2.500000e+01
b    1.000000e+00
s1  -8.141636e-16
s2  -8.141636e-16
s3   9.621933e-16
s4   2.738550e-15
s5  -8.141636e-16
s6  -8.141636e-16
s7   7.401487e-17
s8  -8.141636e-16
s9   9.621933e-16
s10 -8.141636e-16
s11 -8.141636e-16
s12  9.621933e-16

Answer 1

如果我理解你的问题，你真的在寻找使用特定界面的指导。在前面，让我告诉你我不使用RSRuby，但我确实使用不同的实用程序进行R和Ruby集成。

经过进一步检查，我发现你真的想要一个时间序列选项...我相信......是的，可以使用as.ts(data_frame) fxn从数据框中轻松生成。 http://stat.ethz.ch/R-manual/R-patched/library/stats/html/ts.html

我用过的东西很简单。如果这是一个主要问题，也许这对您有用。

require 'rserve/simpler'
r_object = Rserve::Simpler.new

接下来我要做的是获取一个Hash，其中键对应于dataframe列，值是数组，我在它们上运行Rserve::Simpler fxn Hash.to_dataframe，以便它们可以准备好转换。

data = Hash.new()#insert data here
datafr = data.to_dataframe
r_object.converse(df: datafr) do 
  %Q{df$time <- strptime(as.character(df$time), "%Y-%m-%d %X")
     df$name <- factor(df$name)
  }
end

这在converse块中使用直接R代码并很好地处理所有内容。我没有尝试过你的特定问题，但我知道这对我来说是有用的，可以将数组中的时间数据（Ruby Datetime类）导入到R中，这样我就可以绘制它。祝你好运！

如何通过RSRuby将ts对象传递给R.

1 个答案: