XMOVVAR

Purpose:

Calculates the centered moving variance of a series given a duration.

Syntax:

XMOVVAR(series, xdur, dof, "naflag", "edgeflag")

series

-

A series or array, the input series

xdur

-

A real, the X duration used to compute the integer number of samples for the moving interval.

dof

-

Optional. An integer, the degrees of freedom normalization mode.

0

:

sample variance, normalize by 1/(N-1) (default)

1

:

population variance, normalize by 1/N

"naflag"

-

Optional. A string, the NA handling method.

"omitnan"

:

ignore NA values (default)

"includenan"

:

include NA values

"edgeflag"

-

Optional. A string, the edge processing flag.

"full"

:

process all segments (default)

"center"

:

center the result within the overlapping N point segments

"discard"

:

ignore partial segments

"shrink"

:

shrink the segment size to actual number of samples at the edge segments (same as rampflag == 1)  

Returns:

A series or table.

Example:

W1: {1, 2, 4, 7};setdeltax(1/100)

W2: movvar(w1, 3)

W3: xmovvar(w1, 0.02)

 

W1 contains the input data with a sample rate of 100 Hz.

W2 computes a 3 point centered moving variance.

W3 computes the same moving variance except the block size is specified in terms of seconds. Since a 0.02 interval equates to 3 samples, both functions return the series:

{0.5, 2.333, 6.333, 4.5}

Example:

W1: integ(gnorm(10000, 1/10000)) * 1000

W2: movvar(w1, 51)

W3: xmovvar(w1, 0.005)

W4: W2 - W3

 

W1 synthesizes 10000 samples of data with a sample rate of 10kHz.

W2 computes a 51 point centered moving variance.

W3 computes a centered moving variance with a moving interval of 0.005 seconds.

Since a 0.005 interval equates to 51 samples, the two moving variance results are identical.

Remarks:

The interval duration xdur is converted to the integer sample interval N with:

 

N = xdur / deltax(s) + 1

 

where xdur >= deltax(s)

 

NaN values are included by default. Set naflag to "omitnan" to ignore NaN values.

 

See MOVVAR for more details on naflag and edgeflag.

 

See XMOVVARIANCE to compute the non-centered moving variance given a duration.

See Also:

CUMVAR

MOVAVG

MOVRMS

MOVSTD

MOVVAR

STDEV

VAR

XMOVVARIANCE