MARS：如何使用地球包装拟合时间序列并在R中找到结

Question

我们正在开发一个单变量的时间序列数据。

我们怎样才能将火星模型与地球包裹相匹配？我想我也能做到这一点。 在拟合如何在拟合模型中找到结之后。

请找到我们正在使用的代码

library(tsoutliers)
library(expsmooth)
library(fma)
library(tseries)
library(forecast)
library(earth)

setwd("./files/")
examples <- list.files("./")
print(examples)
length(examples)

labels <- c('Raw', 'Smoothed', 'second time smoothed')
cols <- c('red', 'yellow', 'blue')
par(new=TRUE)
for (i in c(2:10)){ # length(examples)
  # par(mfrow=2:1)
  dataseries <- read.csv(examples[i])
  series <- ts(dataseries[2], start = 1, end = length(dataseries[,2]),na.omit(daily_data$cnt_ma), frequency = 1)
  plot(series, col = cols[1], main=examples[i])
  smoothedseries <- filter(series, filter=rep(1/140, 140),
                           method= 'convolution', sides=2)
  lines(smoothedseries, col=cols[2], lwd=2)

  second <- filter(smoothedseries, filter=rep(1/140, 140),
                   method= 'convolution', sides=2)
  lines(second, col=cols[3], lwd=2)
  legend("bottomleft", labels, col=cols, lwd=2)
  df <- data.frame(smoothedseries)
  df[is.na(df)] <- 0
  df1 <- cbind(dataseries[1], df)
  mars <-  earth(x = df1[2], y = df1[1])
  plot(mars)
  #plot.earth.models(mars)
  #earth_plotmodsel(mars)
  #print(mars$grsq)
  # cp <- ts(dataseries[4], start = 1, end = length(dataseries[,4]),na.omit(daily_data$cnt_ma), frequency = 1)
  # plot(cp,  main=examples[i])
}

Answer 1

我能够找到MARS适合的结。通过使用模型$ cuts [model $ selected_cuts]参数