支持向量机 - 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的标题来判断可能是不适合的时间序列是一个好主意。