- html - 出于某种原因,IE8 对我的 Sass 文件中继承的 html5 CSS 不友好?
- JMeter 在响应断言中使用 span 标签的问题
- html - 在 :hover and :active? 上具有不同效果的 CSS 动画
- html - 相对于居中的 html 内容固定的 CSS 重复背景?
我曾尝试在 python
中同时在 R
上使用 auto arima
获取相同的数据,但得到了不同的 ARIMA
模型选择是具有不同 AIC
的最佳模型。你能告诉我为什么我从两种语言中得到具有不同 AIC
的不同最佳模型吗?
R 的数据和代码
wineind <- c(15136., 16733., 20016., 17708., 18019., 19227., 22893., 23739.,
21133., 22591., 26786., 29740., 15028., 17977., 20008., 21354.,
19498., 22125., 25817., 28779., 20960., 22254., 27392., 29945.,
16933., 17892., 20533., 23569., 22417., 22084., 26580., 27454.,
24081., 23451., 28991., 31386., 16896., 20045., 23471., 21747.,
25621., 23859., 25500., 30998., 24475., 23145., 29701., 34365.,
17556., 22077., 25702., 22214., 26886., 23191., 27831., 35406.,
23195., 25110., 30009., 36242., 18450., 21845., 26488., 22394.,
28057., 25451., 24872., 33424., 24052., 28449., 33533., 37351.,
19969., 21701., 26249., 24493., 24603., 26485., 30723., 34569.,
26689., 26157., 32064., 38870., 21337., 19419., 23166., 28286.,
24570., 24001., 33151., 24878., 26804., 28967., 33311., 40226.,
20504., 23060., 23562., 27562., 23940., 24584., 34303., 25517.,
23494., 29095., 32903., 34379., 16991., 21109., 23740., 25552.,
21752., 20294., 29009., 25500., 24166., 26960., 31222., 38641.,
14672., 17543., 25453., 32683., 22449., 22316., 27595., 25451.,
25421., 25288., 32568., 35110., 16052., 22146., 21198., 19543.,
22084., 23816., 29961., 26773., 26635., 26972., 30207., 38687.,
16974., 21697., 24179., 23757., 25013., 24019., 30345., 24488.,
25156., 25650., 30923., 37240., 17466., 19463., 24352., 26805.,
25236., 24735., 29356., 31234., 22724., 28496., 32857., 37198.,
13652., 22784., 23565., 26323., 23779., 27549., 29660., 23356.)
tswineind<-ts(wineind, start=c(1985,1), frequency=12)
library(forecast)
tswineindbest<-auto.arima(tswineind,approximation = FALSE)
tswineindbest
R 的结果
ARIMA(0,1,3)(0,1,1)[12]
Python 的数据和代码
import numpy as np
import pmdarima as pm
from pmdarima.datasets import load_wineind
# this is a dataset from R
wineind = load_wineind().astype(np.float64)
# fit stepwise auto-ARIMA
stepwise_fit = pm.auto_arima(wineind, start_p=1, start_q=1,
max_p=3, max_q=3, m=12,
start_P=0, seasonal=True,
d=1, D=1, trace=True,
error_action='ignore', # don't want to know if an order does not work
suppress_warnings=True, # don't want convergence warnings
stepwise=True) # set to stepwise
stepwise_fit.summary()
Python 结果
SARIMAX(1, 1, 2)x(0, 1, 1, 12) AIC 3066.742
我希望 R
和 Python
具有相同的最佳模型和相同的 AIC
。
最佳答案
我在网上转了一圈,发现这个 python 代码非常有用
# import package
import itertools
# Define the p, d and q parameters to take any value between 0 and 2
p = d = q = range(0, 3)
# Generate all different combinations of p, q and q triplets
pdq = list(itertools.product(p, d, q))
# Generate all different combinations of seasonal p, q and q
triplets
seasonal_pdq = [(x[0], x[1], x[2], 12) for x in
list(itertools.product(p, d, q))]
print('Examples of parameter combinations for Seasonal ARIMA...')
print('SARIMAX: {} x {}'.format(pdq[1], seasonal_pdq[1]))
print('SARIMAX: {} x {}'.format(pdq[1], seasonal_pdq[2]))
print('SARIMAX: {} x {}'.format(pdq[2], seasonal_pdq[3]))
print('SARIMAX: {} x {}'.format(pdq[2], seasonal_pdq[4]))
和
warnings.filterwarnings("ignore") # specify to ignore warning messages
for param in pdq:
for param_seasonal in seasonal_pdq:
try:
mod = sm.tsa.statespace.SARIMAX(ts,
order=param,
seasonal_order=param_seasonal,
enforce_stationarity=False,
enforce_invertibility=False)
results = mod.fit()
print('ARIMA{}x{}12 - AIC:{}'.format(param, param_seasonal, results.aic))
except:
continue
这里输入的变量tp是我在python代码中用ts
表示的单维时间序列数据。结果与 R
中的 auto.arima
相同。
关于r - auto arima : r and python suggest different arima models for same data, 为什么?,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/56802974/
已结束。此问题正在寻求书籍、工具、软件库等的推荐。它不满足Stack Overflow guidelines 。目前不接受答案。 我们不允许提出寻求书籍、工具、软件库等推荐的问题。您可以编辑问题,以便
我的计算机上安装了 SE JDK 1.8.0_45。我希望编译器的源版本和目标版本都是 1.7。 我在 gradle 中做到了这一点,但我不知道如何在 IntelliJ IDEA 14.1.3(社区版
我正致力于根据以前成功搜索中使用的术语在我们的搜索应用程序(使用 Solr)中提供自动建议功能。在 Solr 建议文档 ( http://wiki.apache.org/solr/Suggester
我想过滤来自我的建议者的自动完成结果 假设我有一张书 table Table (Id Guid, BookName String, BookOwner id) 我希望每个用户都能从自己的书中获得一个自
你好 Stack Overflow 的人。我想就以下问题提出一些建议。我正在使用 Java。 我有一个包含多个字符串的数组#1。例如,其中两个字符串可能是:“一个苹果落在牛顿的头上”和“苹果长在树上”
当有人使用我的应用程序时,我想关闭在软/虚拟键盘上显示“建议的单词”(仅在某些 Activity 上)。对于默认的 Android 键盘,可以在“设置”下(在 Word Suggestion Sett
我正在尝试在 Solr 上使用自动完成功能,我找到的方法是使用 Solr Suggester .但是,它没有按预期工作。如果有人调用 http://localhost:8983/solr/techpr
Google建议如何工作?根据遥远的Google数据库中的信息,它如何快速地更新客户端上的网页?如果网页经常更新,为什么页面看起来不“跳动”? 最佳答案 它使用AJAX。 当您编写查询时,它会搜索与您
每当我运行 Build & Reload来自 RStudio 的 Build Pane ,我收到消息 ==> Rcmd.exe INSTALL --no-multiarch rwiots Err
几天来,我正在使用 NetBeans 6.8 进行 PHP 工作。 但是即使包含一个类文件并且方法是公共(public)的并且使用了 phpDoc,NetBeans 每次都会在窗口中显示“No Sug
我不知道我在这里使用的术语是否正确。然而,这就是我想要实现的目标,我想就如何实现这一目标提出一些建议。我想要一个可见边框的圆圈。现在这是困难的部分,我什至不知道如何开始。我想以这样一种方式来操纵圆圈,
我即将开发一个处理极其有值(value)的数据的应用程序。如果用户丢失这些数据,代价将非常高昂,因此我有兴趣了解更多有关满足我们需求的最佳架构设计的信息。 用户每天都会在 iPhone 中输入这些数据
wiki page of the Solr Suggester component没有提到如何搜索提供的字段?仅是前缀,还是也可以进行中缀搜索? 最佳答案 是的,支持。编辑你的 solrconfig.
我创建了一个新的 Azure 搜索建议器,但使用以下代码实现了模糊搜索: ISearchIndexClient indexClient = CreateSearchIndexClient(); var
是否可以让 Elasticsearch 完成建议程序按排序顺序返回结果?我正在遵循 this Elasticsearch blog post 中的示例并增加了额外的酒店: {"name": "Merc
我无法借助 search:suggest 函数提供这个简单的自动完成功能。 基于 Marklogic 演示数据中的奥斯卡语料库,我尝试提供一个建议查询,例如,即使用户当前正在写入“Robert Lo”
这是用 C 语言编写的 CGI 程序的一部分。当客户端单击链接时,我希望开始下载文件,并使用建议的默认文件名。 我知道规范明确指出 Content-disposition header 中指定的文件名
我正在尝试将谷歌建议链接到我的网站。我认为我正在使用的链接或我尝试读取文件的方式可能有问题。无论如何,这是我在我的 php 文件中使用的链接和代码 $filehandle=fopen("http://
我正在寻找一个不错的模板引擎或一小段代码来扩展 Java 字符串中类似 Ant 的变量。示例: String result = expand ("${firstName} ${familyName}"
关闭。这个问题是off-topic .它目前不接受答案。 想改进这个问题? Update the question所以它是on-topic对于堆栈溢出。 8年前关闭。 Improve this que
我是一名优秀的程序员,十分优秀!