NAG Library Routine Document

G10ZAF

1  Purpose

2  Specification

3  Description

4  References

5  Parameters

1:     $\mathrm{WEIGHT}$ – CHARACTER(1)Input

On entry: indicates whether user-defined weights are to be used.

• If $\mathbf{WEIGHT}=\text{'W'}$, user-defined weights are to be used and must be supplied in WT.
• If $\mathbf{WEIGHT}=\text{'U'}$, the data is treated as unweighted.

Constraint: $\mathbf{WEIGHT}=\text{'W'}$ or $\text{'U'}$.

2:     $\mathrm{N}$ – INTEGERInput

On entry: $n$, the number of observations.

Constraint: $\mathbf{N}\ge 1$.

3:     $\mathrm{X}\left(\mathbf{N}\right)$ – REAL (KIND=nag_wp) arrayInput

On entry: the values, $x_{\mathit{i}}$, for $\mathit{i}=1,2,\dots ,n$.

4:     $\mathrm{Y}\left(\mathbf{N}\right)$ – REAL (KIND=nag_wp) arrayInput

On entry: the values $y_{\mathit{i}}$, for $\mathit{i}=1,2,\dots ,n$.

5:     $\mathrm{WT}\left(*\right)$ – REAL (KIND=nag_wp) arrayInput

Note: the dimension of the array WT must be at least $\mathbf{N}$ if $\mathbf{WEIGHT}=\text{'W'}$.

On entry: if $\mathbf{WEIGHT}=\text{'W'}$, WT must contain the $n$ weights. Otherwise WT is not referenced and unit weights are assumed.

Constraints:

• if $\mathbf{WEIGHT}=\text{'W'}$, $\mathbf{WT}\left(\mathit{i}\right)>0.0$, for $\mathit{i}=1,2,\dots ,n$;
• if $\mathbf{WEIGHT}=\text{'W'}$, ${\sum }_{i=1}^n\mathbf{WT}\left(i\right)>0$.

6:     $\mathrm{NORD}$ – INTEGEROutput

On exit: the number of distinct observations.

7:     $\mathrm{XORD}\left(\mathbf{N}\right)$ – REAL (KIND=nag_wp) arrayOutput

On exit: the first NORD elements contain the ordered and distinct $x_i$.

8:     $\mathrm{YORD}\left(\mathbf{N}\right)$ – REAL (KIND=nag_wp) arrayOutput

On exit: the first NORD elements contain the values $y^{\prime }$ corresponding to the values in XORD.

9:     $\mathrm{WTORD}\left(\mathbf{N}\right)$ – REAL (KIND=nag_wp) arrayOutput

On exit: the first NORD elements contain the values $w^{\prime }$ corresponding to the values of XORD and YORD.

10:   $\mathrm{RSS}$ – REAL (KIND=nag_wp)Output

On exit: the within group sum of squares for tied observations.

11:   $\mathrm{IWRK}\left(\mathbf{N}\right)$ – INTEGER arrayWorkspace

12:   $\mathrm{IFAIL}$ – INTEGERInput/Output

On entry: IFAIL must be set to $0$, $-1\text{ or }1$. If you are unfamiliar with this parameter you should refer to Section 3.3 in the Essential Introduction for details.

For environments where it might be inappropriate to halt program execution when an error is detected, the value $-1\text{ or }1$ is recommended. If the output of error messages is undesirable, then the value $1$ is recommended. Otherwise, if you are not familiar with this parameter, the recommended value is $0$. When the value $-\mathbf{1}\text{ or }\mathbf{1}$ is used it is essential to test the value of IFAIL on exit.

On exit: $\mathbf{IFAIL}=\mathbf{0}$ unless the routine detects an error or a warning has been flagged (see Section 6).

6  Error Indicators and Warnings

$\mathbf{IFAIL}=1$

On entry,	$\mathbf{WEIGHT}\ne \text{'W'}$ or $\text{'U'}$,
or	$\mathbf{N}<1$.

$\mathbf{IFAIL}=2$

On entry,

$\mathbf{WEIGHT}=\text{'W'}$ and at least one element of WT is $<0.0$, or all elements of WT are $0.0$.

$\mathbf{IFAIL}=-99$

An unexpected error has been triggered by this routine. Please contact NAG.

See Section 3.8 in the Essential Introduction for further information.

$\mathbf{IFAIL}=-399$

Your licence key may have expired or may not have been installed correctly.

See Section 3.7 in the Essential Introduction for further information.

$\mathbf{IFAIL}=-999$

Dynamic memory allocation failed.

See Section 3.6 in the Essential Introduction for further information.

7  Accuracy

8  Parallelism and Performance

9  Further Comments

10  Example

10.1  Program Text

Program Text (g10zafe.f90)

10.2  Program Data

Program Data (g10zafe.d)

10.3  Program Results

Program Results (g10zafe.r)