qat_analyse_histogram_test_jsd_2d: Perform a histogram test with the metric JSD

Description Usage Arguments Details Value Author(s) References See Also Examples

View source: R/qat_analyse_histogram_test_jsd_2d.R

Description

This check divides the data into blocks, estimates their probability density functions by histograms and compares them by using the Jenson-Shannon Divergence.

Usage

1
2
qat_analyse_histogram_test_jsd_2d(measurement_vector, blocksize, numofbars, 
factorofbar)

Arguments

measurement_vector

The measurement vector (2d array), which should be tested

blocksize

Number of elements in the first dimension, which should be used for each block

numofbars

Number of bins of the histogram

factorofbar

Correction factor for non-value bins

Details

The field will be divided into blocks in the first dimension, with a length given by the parameter blocksize. From these blocks histograms are computed and afterwards compared. As a metric for the comparison the Jenson-Shannon Divergence is used. As a result a field is generated, which includes the results of the comparison between every combination of blocks.

Value

It returns a list with the following entries:

field

Result matrix of the comparison.

blocksize

Size of blocks in the first dimension.

numofbars

Number of bins of the used histograms.

factorofbar

Correction factor used for the used histograms.

metric

Used metric in the comparisons.

runs

Number of blocks, which are compared.

Author(s)

Andre Duesterhus

References

Duesterhus, A., Hense, A. (2012) Advanced Information Criterion for Environmental Data Quality Assurance, \_Advances in Science and Research\_, *8*, 99-104.

See Also

qat_analyse_histogram_test_jsd_1d, qat_analyse_histogram_test_kld_2d, qat_analyse_histogram_test_rms_2d, qat_analyse_histogram_test_ms_2d, qat_analyse_histogram_test_emd_2d

Examples

1
2
3
4
5
vec <- array(rnorm(1000), c(100, 20))
vec[51:100, ] <- round(vec[51:100, ])
result <- qat_analyse_histogram_test_jsd_2d(vec, 4, 65, 100)
qat_plot_histogram_test(result$field, "test_jsd_2d", result$blocksize, 
result$numofbars, result$factorofbar, "jsd", result$runs)

qat documentation built on May 30, 2017, 8:27 a.m.