朱莉娅|如何对TimeArray数据集

时间:2016-08-15 00:40:24

标签: statistics time-series regression julia linear-regression

我有TimeArray类型的数据集,我想执行线性回归。但是,似乎julia目前不支持TimeArray类型的回归方法。

我可以将数据下载为DataFrame而不是TimeArray并使用GLM包,但TimeArray时间戳对以后的其他分析非常有用。我想直接在TimeArray数据集上执行线性回归。

编辑1:下面给出一个简单示例:

julia> using TimeSeries

       dates = collect(Date(1999,1,1):Date(1999,1,31))

       # Dependent variable
       y  = TimeArray(dates, rand(length(dates)))

       # Explanatory variables 
       x1 = TimeArray(dates, rand(length(dates))) # Explanatory variable 1
       x2 = TimeArray(dates, rand(length(dates))) # Explanatory variable 2
       x  = rename(merge(x1,x2), ["x1", "x2"])    # Merge x1 and x2 into a single TimeArray

       # Linear regression
       coefs = linreg(x, y) # Yields a method error since linreg does not support the TimeArray type. 

是否有人为此问题找到了解决方案或解决方法?

1 个答案:

答案 0 :(得分:4)

TimeArray类型似乎有一个.values字段,您可以使用该字段以正确的顺序获取与数组关联的值。因此,您可以使用以下方式执行线性回归:

coefs = linreg(x.values,y.values)