MARIMA (многомерная модель ARIMA) Ошибка в optim(init[mask]
Всем привет,
В настоящее время я использую многомерную модель ARIMA с использованием многопараметрических данных. Я следовал некоторым учебникам на сайте, но до сих пор у меня все еще появляется ошибка в этом коде. Я пытаюсь оценить данные, используя несколько переменных в моих данных. Исходные данные, которые я хочу загрузить, выглядят следующим образом:
A B C
34.77452838 191.9425766 2
51.34502757 167.3464789 2
52.22462976 163.9371881 2
44.05790246 153.8406923 2
49.01569605 173.6158416 2
43.29661071 168.8358248 2
52.82309532 162.5718306 2
47.71516049 173.088253 2
40.46453333 173.1037056 2
43.96892905 156.5669003 2
50.53534484 168.1469926 2
54.29771912 191.7929009 2
49.84792352 153.2957975 2
47.29999697 157.3547558 2
50.74218512 167.5583718 2
51.04921973 169.2961421 2
42.94721591 161.0413471 2
46.857445 161.5720713 2
59.63383675 168.1864209 2
41.55891144 164.8500058 2
43.80523336 170.6241115 2
40.54447174 163.5741909 2
54.73254502 145.0296938 2
48.37875497 165.6939519 2
54.28330553 185.9884903 2
54.56595576 187.6435893 2
46.3341881 166.420677 2
48.64491069 197.8226879 2
52.16211736 156.9031211 2
49.69608331 174.5718085 2
44.3995806 165.870505 2
54.30303288 172.1445435 2
59.42004347 181.4705789 2
55.66045916 171.5254118 2
50.9424361 167.3649679 2
52.74883235 173.8478608 2
51.88352442 180.565498 2
48.82856619 181.4147966 2
52.91082788 166.5114416 2
61.81566179 186.9227225 2
59.18481278 171.9375519 2
53.73142791 165.4547203 2
45.90694845 167.5522616 2
53.20516884 177.0629353 2
58.42090189 179.9778995 2
48.20763016 181.7318594 2
64.8245225 182.4593924 2
53.4869312 178.0462498 2
50.98494649 198.5952871 2
57.18331504 180.2024605 2
56.07515264 184.3858003 2
53.02835584 119.6381436 2
52.83894384 190.4696428 2
64.60654461 190.7581833 2
57.41005075 187.114192 2
52.50923717 188.1313879 2
51.53594851 151.6875626 2
52.83034897 178.0217738 2
51.62433934 173.5677154 2
57.34844828 185.234996 2
56.90921938 197.0879734 2
51.47780192 201.7652231 2
53.35943985 184.3207373 2
55.12425864 202.156465 2
58.62819278 176.0629181 2
62.51672983 187.3645886 2
61.35114133 219.0208926 2
69.00147164 235.7991778 2
48.19376576 225.5629393 2
59.97428918 221.5355618 2
61.42004728 221.3157335 2
56.09382308 256.9632484 2
58.12355614 228.7762196 2
61.71775579 243.6576741 2
58.30987263 226.1145379 2
66.17636418 139.7724771 2
55.87241828 205.2602213 2
55.59852946 223.8182453 2
61.9962275 199.569491 2
61.64428604 209.228291 2
60.89585352 212.1291189 2
35.60671771 147.2426824 2
38.79108548 186.9673003 2
47.54260194 112.340896 2
55.62540507 178.5126261 2
53.74473393 211.8612091 2
46.5023706 196.6053041 2
56.1878252 216.0413246 2
60.25617886 213.8808684 2
55.06903768 191.5933198 2
46.60528803 193.228591 2
61.16524374 196.2754949 2
66.04131413 215.8749835 2
55.65532315 195.3569658 2
64.23627043 188.9093598 2
61.10553396 206.492755 2
47.89633095 124.2445289 2
58.98819304 124.9045309 2
57.35811675 189.5788539 2
55.1452173 196.5306579 2
55.75572991 221.4020661 2
63.98512149 232.496751 2
55.9264102 194.8117255 2
53.8304683 194.4430075 2
64.15636933 233.0761708 2
66.40452218 225.4115469 2
64.83224857 237.8602368 2
61.32524633 273.593856 2
59.92230034 255.1825145 2
55.16024733 233.8121403 2
48.41788888 213.9495593 2
45.87945926 182.9052714 2
37.12352931 214.9122403 2
36.23515606 184.767225 2
38.80360508 203.3612199 2
33.82899797 189.0318566 2
29.51428461 184.9648764 2
28.57971239 186.4950075 2
24.77523589 188.7776908 2
31.33283567 197.7679355 2
40.3275789 186.8421363 2
41.6556977 193.8692738 2
41.73315692 197.7107161 2
49.60428262 188.7511304 2
41.20251334 219.1196343 2
43.95248687 210.6029669 2
44.04459047 203.7891106 2
41.86283159 222.413766 2
44.25611329 219.7560903 2
53.94720805 244.3697403 2
47.52787864 233.2155174 2
59.11544025 83.04917738 2
51.94031441 221.0877964 2
54.44512963 180.7172578 2
74.13162887 141.8472405 2
58.41810071 52.32844725 2
34.09419501 137.1900364 2
44.34173286 204.0982826 2
52.73150897 204.2542315 2
52.38756931 215.1004646 2
54.43326163 198.5332844 2
52.59076023 213.0977851 2
56.87320662 200.364619 2
60.35769665 226.456495 2
55.56278813 183.1343606 2
51.77700472 204.6629449 2
56.84940159 220.0856916 2
50.40147328 234.4078768 2
57.58252537 227.4383631 2
54.13363087 200.6441259 2
64.28838575 187.2168003 2
50.38259399 229.655229 2
61.97002161 169.5985726 2
53.4693228 228.3208836 2
43.46435559 227.805616 2
50.76040995 221.5824154 2
37.74385309 243.5761809 2
44.50296307 260.2599051 2
52.93274188 208.5485301 2
48.94673824 229.0260541 2
44.69142139 235.2965628 2
45.87330532 230.5345631 2
50.15886772 231.5080565 2
50.46518803 233.3240541 2
39.64570737 252.3943345 2
50.31436396 257.1978301 2
45.72599959 255.5499076 2
35.90772772 235.547594 2
50.64034235 256.1243588 2
42.75332725 250.8706385 2
43.32600331 218.2298083 2
51.61082256 256.5452666 2
44.12634993 239.613958 2
46.19641519 238.9780285 2
51.83620155 215.1587094 2
43.97586417 219.6388278 2
48.92009437 251.0237063 2
49.06462419 71.92791063 2
55.18403387 262.2636031 2
32.57642794 236.0216013 2
43.32588315 217.0611338 2
47.7579248 244.944423 2
47.83448696 218.8873504 2
59.93806362 224.4503085 2
47.25084496 191.4388145 2
53.37664378 215.1277665 2
48.80090845 203.8133075 2
43.18675303 202.2225439 2
50.6951052 191.3762616 2
51.30918467 219.4169811 2
47.40389884 239.3385311 2
41.98253107 224.8199498 2
48.2856251 213.795973 2
44.51645684 218.5239546 2
47.280002 209.8027239 2
40.54046452 220.6853516 2
45.78141475 229.3941645 2
51.80096102 240.4784644 2
61.72670174 237.1822198 2
49.78334701 257.9867854 2
50.54248619 214.5982823 2
48.62856185 223.9721336 2
48.48320973 229.0910186 2
50.2515676 26.9172095 2
50.27965224 256.2915689 2
56.12183511 245.8979531 2
55.52810276 271.3756176 2
44.64924526 222.7377583 2
50.11476493 259.5392736 2
46.87115741 228.0808933 2
56.00424314 199.4617359 2
48.82094455 245.3930835 2
53.08217335 229.5334666 2
51.90354311 244.6746334 2
43.06275117 246.7410669 2
50.51734376 232.7116094 2
48.31047595 232.380875 2
49.67464817 217.3068271 2
57.93947721 231.5371941 2
56.18981779 248.290565 2
50.98111296 275.0819118 2
42.7040962 228.5930888 2
54.32273841 237.5194329 2
40.49863088 246.5566334 2
44.50957453 271.3497973 2
48.52619112 282.228831 2
55.82672036 271.9136146 2
49.16413927 279.3997284 2
57.59784043 281.7114483 2
52.18583846 272.4423745 2
49.37065029 267.0389461 2
47.12135267 236.7228498 2
52.00103676 248.8965285 2
51.79680061 279.691718 2
55.56773472 282.9113431 2
50.85593855 227.6402154 2
Вот код для вызова этих файлов и использования MARIMA:
# Remove All Variables rm(list=ls())
library(marima) # Calling the Library MARIMA library(tseries) # Calling the Library tseries library(ggplot2)
options(max.print = 1000000)
setwd("E:/LoRA Technology/myPaper")
mydata <- read.csv ("Копия combining_2.csv", заголовок = ИСТИНА)
Model5 <- define.model(kvar=3, ar=1, ma=1)
Marima5 <- marima(ts(mydata), Model5$ar.pattern, Model5$ma.pattern,
penalty=1) nstart <-418 nstep <-10
Тем не менее, я получил ошибку следующим образом, когда я пытался запустить Marima5
Error in optim(init[mask], armaCSS, method = optim.method, hessian = FALSE, :
non-finite value supplied by optim
Я просто не знаю, как решить это, у вас есть какое-то решение для этого?
С уважением,Идхам Хабиби