pomdev.extra {Pomic}R Documentation

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)...

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)

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

Author(s)

Cyril Piou
cyril.piou@yahoo.fr

See Also

pomdev pomdev.corrected

Examples

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)

[Package Pomic version 1.0.1 Index]