支持向量機 - R的代碼 - 預測時間序列

Question

我試圖預測，使用R代碼的時間序列殘差的殘差。我的數據集有以下兩列（我會放一個樣品與第10行）：

Observation Residuals 
1 -0,087527458 
2 -0,06907199 
3 -0,066604145 
4 -0,07796713 
5 -0,081723932 
6 -0,094046868 
7 -0,101535816 
8 -0,101884203 
9 -0,11131246 
10 -0,092548176 

對於我建立一個支持向量機使用R中的預測：

# Load the data from the csv file 
dataDirectory <- "C://" 
data <- read.csv(paste(dataDirectory, "Data_SVM_Test.csv", sep=""),sep=";", header = TRUE) 
head(data) 
# Plot the data 
plot(data, pch=16) 

# Create a linear regression model 
model <- lm(Residuals ~ Observation, data) 

# Add the fitted line 
abline(model) 

predictedY <- predict(model, data) 

# display the predictions 
points(data$Observation, predictedY, col = "blue", pch=4) 

# This function will compute the RMSE 
rmse <- function(error) 
{ 
    sqrt(mean(error^2)) 
} 

error <- model$residuals # same as data$Y - predictedY 
predictionRMSE <- rmse(error) # 5.70377 

plot(data, pch=16) 

plot.new() 
# svr model ============================================== 
if(require(e1071)){ 
    model <- svm(Residuals ~ Observation , data) 

    predictedY <- predict(model, data) 

    points(data$Observation, predictedY, col = "red", pch=4) 

    # /!\ this time svrModel$residuals is not the same as data$Y - predictedY 
    # so we compute the error like this 
    error <- data$Residuals - predictedY 
    svrPredictionRMSE <- rmse(error) # 3.157061 
} 

當我執行上面的代碼我收到以下錯誤消息，並沒有任何輸出：

Warning message: 
In Ops.factor(data$Residuals, predictedY) : ‘-’ not meaningful for factors 

任何人都有一個想法如何解決這個錯誤？

非常感謝！

Answer 1

當使用svm進行分類，輸出類型因子。這是從文檔：

Output of svm: A vector of predicted values (for classification: a vector of labels, for density estimation: a logical vector).

這可以從下面的例子可以看出：

library(e1071) 
model <- svm(Species ~ ., data = iris) 
> str(predict(model, iris)) 
Factor w/ 3 levels "setosa","versicolor",..: 1 1 1 1 1 1 1 1 1 1 ... 
- attr(*, "names")= chr [1:150] "1" "2" "3" "4" ... 

這是你的數據相同。水平表明PredictedY是一個因素：

> predictedY <- predict(model, df) 
> predictedY 
      1   2   3   4   5   6   7   8   9   10 
-0,087527458 -0,06907199 -0,066604145 -0,07796713 -0,081723932 -0,094046868 -0,101535816 -0,101884203 -0,11131246 -0,092548176 
Levels: -0,066604145 -0,06907199 -0,07796713 -0,081723932 -0,087527458 -0,092548176 -0,094046868 -0,101535816 -0,101884203 -0,11131246 

在你的代碼predictedY <- predict(model, data)線，predictedY的類型的因素。如果你試圖從一個因子（反之亦然）扣除一個數字，你讓你的錯誤：

> 1:10 - as.factor(1:10) 
[1] NA NA NA NA NA NA NA NA NA NA 
Warning message: 
In Ops.factor(1:10, as.factor(1:10)) : ‘-’ not meaningful for factors 

如果要使其工作，你需要使用as.numeric因素轉換爲數字。 1:10 - as.numeric(as.factor(1:10))。

我不知道你的數據是什麼樣子，但是我從這個問題svm的標題來判斷可能是不適合的時間序列是一個好主意。