Aplique uma transformação Box–Cox a uma série temporal 

O TimeSeriesMap pode ser usado para transformar os valores das séries temporais.

Os dados contém a taxa de desemprego civil nos Estados Unidos de 1948 a 1993.

mostre o input completo de Wolfram Languageoculte o input

In[1]:=

X timeFrame = Interval[{DateObject[{1948, 1, 1}], DateObject[{1993, 12, 1}]}];

In[2]:=

X data = EntityValue[Entity["Country", "UnitedStates"], EntityProperty["Country", "CivilianUnemploymentRate", {"Date" -> timeFrame, "Frequency" -> "Monthly"}]]

Out[2]=

In[3]:=

X DateListPlot[data, ImageSize -> Medium]

Out[3]=

Definir uma transformação Box - Cox.

In[4]:=

X BoxCox[x_, \[Lambda]_, \[Alpha]_, \[Sigma]_] := With[{y = QuantityMagnitude[x]}, Piecewise[{{(((y - \[Alpha])/\[Sigma])^\[Lambda] - 1)/\[Lambda], \[Lambda] > 0 || \[Lambda] < 0}}, Log[(y - \[Alpha])/\[Sigma]]]]

Use o TimeSeriesMap para aplicar uma Box - Cox aos dados com poucas opções de λ.

In[5]:=

X lambdas = {1/2, 0, -1/2}; newdata = Table[TimeSeriesMap[BoxCox[#, \[Lambda], .02, .016] &, data], {\[Lambda], lambdas}];

mostre o input completo de Wolfram Languageoculte o input

In[6]:=

X Grid[{MapThread[ DateListPlot[#1, ImageSize -> 180, PlotRange -> All, FrameTicks -> {{Table[{x, NumberForm[N[x], {5, 2}]}, {x, Min[#1], Max[#1], (Max[#1] - Min[#1])/5}], Automatic}, {Automatic, Automatic}}, PlotLabel -> Row[{"\[Lambda] = ", #2}]] &, {newdata, lambdas}]}]

Out[6]=