pomdev.extra: POMDEV calculation for numerical patterns as well as other goodness of fit indicators

Description

This function calculate the POMDEV value of simulations results compared to a vector taken as a field pattern and compute also other goodness of fit indicators.

Usage

pomdev.extra(object1, object2, eps=10^-30, nrange=1000, 
              fullmsd=FALSE,  plotting=FALSE,...)

Arguments

object1

Numerical, vector of the field pattern

object2

Numerical, vector of simulations results

eps

Numerical, value to use when the likelihood = 0 for mathematical stability

nrange

Numerical, number of interval to consider for the probability density function creation and integration

fullmsd

Boolean, should the full cross MSD be calculated as well

plotting

Boolean, should a plot illustrating the calculation be created

…

eventual information to give to the last plot (if plotting=T)...

Value

Return a list containing:

field_data

the field data

sim_data

the simulation data

kernel_estimator

information about the kernel estimator used

result

a table of 4 or 5 rows (depending if fullmsd=TRUE or not) with:

- POMDEV

pomdev score of model results to reproduce the pattern

- overlap

indicates if the model results overlap the range of the field data

- KLdiv

KL divergence of the model results from the field pattern

- MSD

\((\sum(field)/n_{field}-\sum(simulation)/n_{simulation})^2\)

- CrossMSD

if fullmsd is true, calculate the mean square deviation among each pair of simulation and field data possible

Details

The KL divergence is only calculated comparing the simulations to the field vector, as if it is the truth (different from the flexmix KLdiv calculations)

Examples

Run this code

# NOT RUN {
pomdev.extra(rnorm(100,10,5),rnorm(100*5,10.5,6),eps=10^-20)
pomdev.extra(rnorm(100,10,5),rnorm(100*5,10.5,6),eps=10^-20,fullmsd=TRUE,plotting=TRUE)
# }