来自网络资源的数据
下载过去一年中各大湖的每日平均气温.
In[1]:=

rawdata =
Import["http://coastwatch.glerl.noaa.gov/ftp/glsea/avgtemps/glsea-\
temps_1024.dat", "Data"];
观测是从第 11 行开始的.
In[2]:=

rawdata[[;; 11]]
Out[2]=

In[3]:=

obs = rawdata[[11 ;;]];
观测次数.
In[4]:=

Length[obs]
Out[4]=

定义湖泊列表.
In[5]:=

lakes = {"Superior", "Michigan", "Huron", "Erie", "Ontario",
"St. Clair"};
显示完整的 Wolfram 语言输入
Out[6]=

创建日期用作 TimeSeries 的时间戳.
In[7]:=

dates = Table[
DatePlus[{obs[[k, 1]], 1, 1}, Quantity[obs[[k, 2]], "Days"]], {k,
1, Length[obs]}];
用 QuantityArray 来存储温度和单位.
In[8]:=

values = QuantityArray[obs[[All, 3 ;;]], "DegreesCelsius"]
Out[8]=

以 MetaInformation 存储湖泊列表以便以后提取使用,如用于产生图例.
In[9]:=

ts = TimeSeries[values, {dates},
MetaInformation -> {"Lakes" -> lakes}]
Out[9]=

In[10]:=

DateListPlot[ts, PlotLegends -> ts["Lakes"], PlotTheme -> "Detailed",
FrameLabel -> Automatic]
Out[10]=

研究每个湖泊的描述性统计量.
In[11]:=

stats = {Min, Max, Mean, Median, StandardDeviation};
summary[lake_] := Map[#[lake] &, stats]
In[12]:=

TableForm[
Table[summary[ts["PathComponent", k]], {k, 1, Length[lakes]}],
TableHeadings -> {lakes, stats}]
Out[12]//TableForm=
