已经求出了时间序列yt的fit后的残差的平方,residuals(fit).并且对其平方,然后求得了acf,可是
a=acf((resduals(fit))^2,lag.max=100)完了出来的数据没法直接求mse=sum((a-mean(a))^2)/length(a).
a
0 1 2 3 4 5 6 7 8 9 10
1.000 0.297 0.210 0.221 0.061 0.074 -0.001 0.053 0.017 0.001 0.018
11 12 13 14 15 16 17 18 19 20 21
0.047 0.067 -0.009 0.012 -0.013 -0.032 -0.044 -0.017 -0.035 0.004 0.006
22 23 24 25 26 27 28 29 30 31 32
-0.002 0.046 -0.017 0.094 0.137 -0.019 0.024 0.005 0.013 0.015 -0.007
33 34 35 36 37 38 39 40 41 42 43
0.096 0.041 0.095 0.160 0.158 0.171 0.030 0.060 0.037 0.028 -0.016
44 45 46 47 48 49 50 51 52 53 54
0.010 0.025 0.017 0.030 -0.008 -0.004 -0.021 -0.015 -0.005 -0.013 -0.005
55 56 57 58 59 60 61 62 63 64 65
-0.021 -0.002 -0.024 -0.019 -0.015 -0.005 0.012 -0.011 0.003 0.006 0.003
66 67 68 69 70 71 72 73 74 75 76
0.020 0.088 0.015 -0.001 0.003 -0.007 -0.019 -0.018 0.002 -0.018 -0.011
77 78 79 80 81 82 83 84 85 86 87
0.029 -0.001 -0.010 -0.026 -0.013 -0.007 -0.021 0.005 0.068 0.009 0.015
88 89 90 91 92 93 94 95 96 97 98
0.025 0.017 0.001 -0.005 0.016 0.008 0.011 -0.015 -0.010 0.003 0.017
99 100
0.020 0.011
length(a)得到的是6,而且mse公式也没法算,是a的数据类型的原因好像,我的想法是转换成可以算的类型,可是试了下没成功,求朋友们帮帮我,非常感谢