# NAG CL Interfaceg10cac (data_​runningmedian)

## 1Purpose

g10cac computes a smoothed data sequence using running median smoothers.

## 2Specification

 #include
 void g10cac (Nag_Smooth_Type smoother, Integer n, const double y[], double smooth[], double rough[], NagError *fail)
The function may be called by the names: g10cac, nag_smooth_data_runningmedian or nag_running_median_smoother.

## 3Description

Given a sequence of $n$ observations recorded at equally spaced intervals, g10cac fits a smooth curve through the data using one of two smoothers. They are based on the use of running medians and averages to summarise the overlapping segments. The fit is called the smooth, the residuals the rough and they obey the following:
Data = Smooth + Rough
The two smoothers are :
1. 1.4253H, twice consisting of a running median of $4$, then $2$, then $5$, then 3 followed by Hanning. Hanning is a running weighted average, the weights being 1/4, 1/2 and 1/4. The result of this smoothing is then reroughed by computing residuals, applying the same smoother to them and adding the result to the smooth of the first pass.
2. 2.3RSSH, twice consisting of a running median of $3$, two splitting operations named S to improve the smooth sequence, each of which is followed by a running median of $3$, and finally Hanning. The end points are dealt with using the method described by Velleman and Hoaglin (1981). The full smoother 3RSSH, twice is produced by reroughing as described above.
The compound smoother 4253H, twice is recommended. The smoother 3RSSH, twice is popular when calculating by hand as it requires simpler computations and is included for comparison purposes.

## 4References

Tukey J W (1977) Exploratory Data Analysis Addison–Wesley
Velleman P F and Hoaglin D C (1981) Applications, Basics, and Computing of Exploratory Data Analysis Duxbury Press, Boston, MA

## 5Arguments

1: $\mathbf{smoother}$Nag_Smooth_Type Input
On entry: smoother must specify the method to be used.
${\mathbf{smoother}}=\mathrm{Nag_4253H}$
4253H, twice is used.
${\mathbf{smoother}}=\mathrm{Nag_3RSSH}$
Constraint: ${\mathbf{smoother}}=\mathrm{Nag_4253H}$ or $\mathrm{Nag_3RSSH}$.
2: $\mathbf{n}$Integer Input
On entry: the number, $n$, of the observations.
Constraint: ${\mathbf{n}}>6$.
If ${\mathbf{n}}\le 6$ then the sequence is not long enough to carry out smoothing.
3: $\mathbf{y}\left[{\mathbf{n}}\right]$const double Input
On entry: the sample observations.
4: $\mathbf{smooth}\left[{\mathbf{n}}\right]$double Output
On exit: contains the smooth.
5: $\mathbf{rough}\left[{\mathbf{n}}\right]$double Output
On exit: contains the rough.
6: $\mathbf{fail}$NagError * Input/Output
The NAG error argument (see Section 7 in the Introduction to the NAG Library CL Interface).

## 6Error Indicators and Warnings

On entry, argument smoother had an illegal value.
NE_INT_ARG_LE
On entry, ${\mathbf{n}}=〈\mathit{\text{value}}〉$.
Constraint: ${\mathbf{n}}>6$.

Not applicable.

## 8Parallelism and Performance

g10cac is not threaded in any implementation.

### 9.1Internal Changes

Internal changes have been made to this function as follows:
• At Mark 25: g10cac is a smoothing function with two possible smoothing methods. The function was previously using the incorrect method (i.e., if you asked for method A you would get method B, and vice versa).
For details of all known issues which have been reported for the NAG Library please refer to the Known Issues.

## 10Example

The example program reads in a sequence of 49 data taken from Tukey (1977), above. Results are obtained using the two smoothing methods described.

### 10.1Program Text

Program Text (g10cace.c)

### 10.2Program Data

Program Data (g10cace.d)

### 10.3Program Results

Program Results (g10cace.r)