欠損データを含む時系列を操作する 

Mathematica は記号的であるため，時系列の値はどのような式でもよい．特に式Missing[…]は欠損観測値に注釈を付けるために使われ，Mathematica は特別な方法でそれを扱う．

イリノイ州シャンぺーンの2014年5月14日の気温の時系列．

In[1]:=

X wd = WeatherData[{"Champaign", "IL"}, "Temperature", {2014, 5, 14}]

Out[1]=

In[2]:=

X Cases[wd["Path"], {t_, mv_Missing} :> {DateObject[t], mv}]

Out[2]=

Out[940]=

欠損値は記号式として扱われる．

In[3]:=

X Mean[wd]

Out[3]=

MissingDataMethodオプションを使って，欠損値を充填する方法を指定する．

In[4]:=

X ts1 = TimeSeries[wd, MissingDataMethod -> "Interpolation"];

In[5]:=

X ts2 = TimeSeries[wd, MissingDataMethod -> {"Constant", Quantity[10.7, "DegreeCelsius"]}];

平均を比較する．

In[6]:=

X {Mean[ts1], Mean[ts2]}

Out[6]=