pomdev.extra: POMDEV calculation for numerical patterns as well as other...
In Pomic: Pattern Oriented Modelling Information Criterion

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

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.

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

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

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)

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`	(∑(field)/n_{field}-∑(simulation)/n_{simulation})^2
`- CrossMSD`	if fullmsd is true, calculate the mean square deviation among each pair of simulation and field data possible

Cyril Piou
cyril.piou@cirad.fr

pomdev pomdev.corrected

1 2	pomdev.extra(rnorm(100,10,5),rnorm(1005,10.5,6),eps=10^-20) pomdev.extra(rnorm(100,10,5),rnorm(1005,10.5,6),eps=10^-20,fullmsd=TRUE,plotting=TRUE)