A moving average is used to smooth out a
time series. Computing moving average is a typical case of ordered data computing.
Its basic computing method is to create a subset composed of N consecutive members
of a time series, compute the average of the set and shift the subset forward
one by one. The following example teaches you how to compute moving average in
R language.
Case description:
Data frame sales has two
fields: salesDate and Amount of this date. Requirement:
compute the moving average in three days. Computing steps include seeking sales
amount average of the previous day, the current day and the next day, and shift
forward along the dates. A part of the source data is as follows:
Code:
Computed result:
Code interpretation:
filter function can be used in R language to compute moving average, which
produces concise code. This method is quite convenient.
Despite the convenience of the filter function, it is difficult to
understand for beginners. For example, sales$Amount/3 means dividing the current value of field Amount by
three,but when it is used in filter
function, it may mean adding the three consecutive values together, then divide
the sum by three. [1,1,1] is the value of expression rep(1,3), which is used
here to specify the range of data fetching. In addition, because neither the
name nor the parameters of filter
function contain the words "average" and "moving", even many developers of R
language don't know its use for computing moving average.
In fact, filter function is a universal linear filter. Its use is more than
computing moving average. Its complete function reference is filter(x, filter, method =
c("convolution", "recursive"),sides = 2, circular = FALSE,
init).
Any modification of the requirement will
make the code more difficult to understand. For example, the code for computing
moving average of the current day and the previous two days cannot be written
as filter(sales$Amount/3, rep(0,2)),
it has to befilter(sales$Amount/3,
rep(1,3), sides = 1).
Summary:
R language can compute moving average, but
its code is rather elusive.
Third-party solutions
We
can also use Python, esProc and Perl to handle this case. As R language, all of
these languages can perform data statistics and analysis and compute moving
average. The following introduces solutions of Python and esProc briefly.
Python(pandas)
Pandas is Python's third-party library
function. It is powerful in processing structured data with basic data type
imitating R's dataframe. At present the latest version is 0.14. Its code for
handling this case is as follows:
pandas.stats.moments.rolling_mean(sales["Amount"],
3)
The name of rolling_mean function is clear, even a developer without experience
with pandas can understand it easily. The function’s usage is simple too. Its
first parameter is the sequence being computed and the second parameter is N, which
is the number of days in seeking moving average.
esProc
esProc is good at expressing business logic
freely with agile syntax. Its expressions for relative position can solve
computational problems of ordering data easily. The code is as follows:
sales.(Amount{-1,1}.avg())
{-1,1}
in the code represents a relative interval, that is, the three days of the
previous day, the current day and the next day. It can be seen that moving
average can be worked out clearly and flexibly by using a relative interval. If
it is required, for example, to compute the moving average of the current day
and the previous two days, we just need to change the interval to {-2,0}in esProc.
没有评论:
发表评论