ARIMA錯誤

Question

我正在嘗試在.csv文件中的時間數據集上運行ARIMA。 到目前為止，這是我的代碼：

Oil_all <- read.delim("/Users/Jkels/Documents/Introduction to Computational
Statistics/Oil production.csv",sep="\t",header=TRUE,stringsAsFactors=FALSE)
Oil_all

該文件如下所示：

year.mbbl
1     1880,30
2     1890,77
3    1900,149
4    1905,215
5    1910,328
6    1915,432
7    1920,689
8   1925,1069
9   1930,1412
10  1935,1655
11  1940,2150
12  1945,2595
13  1950,3803
14  1955,5626
15  1960,7674
16  1962,8882
17 1964,10310
18 1966,12016
19 1968,14104
20 1970,16690
21 1972,18584
22 1974,20389
23 1976,20188
24 1978,21922
25 1980,21732
26 1982,19403
27 1984,19608

碼：

apply(Oil_all,1,function(x) sum(is.na(x)))

結果：

[1] 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0

當我運行ARIMA時：

library(forecast)
auto.arima(Oil_all,xreg=year)

這是錯誤：

Error in lm.fit(x, y, offset = offset, singular.ok = singular.ok, ...) : 
  0 (non-NA) cases
In addition: Warning message:
In data.matrix(data) : NAs introduced by coercion

因此，我能夠調用數據集並進行打印。 但是，當我去檢查值是否與apply函數一起出現時，我看到的全為0，所以我知道出了點問題，這可能就是為什么我得到此錯誤。 我只是不確定錯誤的含義或如何在代碼中修復它。

有什么建議嗎？

Answer 1

如果我對您的問題正確，應該是這樣的：

Oil_all <- read.csv("myfolder/myfile.csv",header=TRUE) 
## I don't have your source data, so I tried to reproduce it with the data you printed
Oil_all
   year value
1  1880    30
2  1890    77
3  1900   149
4  1905   215
5  1910   328
6  1915   432
7  1920   689
8  1925  1069
9  1930  1412
10 1935  1655
11 1940  2150
12 1945  2595
13 1950  3803
14 1955  5626
15 1960  7674
16 1962  8882
17 1964 10310
18 1966 12016
19 1968 14104
20 1970 16690
21 1972 18584
22 1974 20389
23 1976 20188
24 1978 21922
25 1980 21732
26 1982 19403
27 1984 19608


library(forecast)
auto.arima(Oil_all$value,xreg=Oil_all$year)

Series: Oil_all$value 
ARIMA(3,0,0) with non-zero mean 

Coefficients:
         ar1     ar2      ar3  intercept  Oil_all$year
      1.2877  0.0902  -0.4619  -271708.4      144.2727
s.e.  0.1972  0.3897   0.2275   107344.4       55.2108

sigma^2 estimated as 642315:  log likelihood=-221.07
AIC=454.15   AICc=458.35   BIC=461.92

Answer 2

您的導入應為

Oil_all<-read.csv("/Users/Jkels/Documents/Introduction to Computational Statistics/Oil production.csv")

這就是為什么您的數據很奇怪的原因。 抱歉，我沒有發表評論的名聲。我和Nemesi一樣，后來奏效了。 我認為您正在嘗試將csv導入為制表符分隔的文件。