ETS と季節ナイーブの比較

この演習では、cement の20年分（四半期ごとのセメント生産量）のデータに対して、4期先の時系列クロスバリデーションを用い、ETS 予測と季節ナイーブ予測を比較します。実行に時間がかかるため、作業スペースには短縮版の cement 系列が用意されています。

tsCV() の2番目の引数は forecast オブジェクトを返す必要があるため、モデルを当てはめて予測を返す関数が必要です。思い出してください:

> args(tsCV)
function (y, forecastfunction, h = 1, ...)

この演習では、既存の予測関数と、あらかじめ用意された関数の両方を使います。ときには、シンプルな手法のほうが高度な手法より良い結果を出すこともあります。

ETS の予測を返す関数 fets() は用意されています。
ETS と季節ナイーブの両手法に対して、cement データに tsCV() を適用し、予測ホライズンを 4 にします。予測関数の引数には、新しく作成された fets と既存の snaive をそれぞれ e1 と e2 に使ってください。
得られた 4 期先誤差の MSE を計算し、欠損値を除外します。MSE を計算する式は用意されていますが、2つ目のオプション引数は未指定です（以前に使いましたね）。
最良の MSE を bestmse として保存します。前の指示で最良の MSE を生成するコード行をそのままコピーして構いません。