pomic {Pomic}R Documentation

Pomic calculation for numerical patterns

Description

This function calculate the POMIC value of simulations results compared to a vector taken as a field pattern.

Usage

pomic(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:
- POMIC pomic 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

pomic.simple pomic.corrected

Examples

pomic(rnorm(100,10,5),rnorm(100*5,10.5,6),eps=10^-20)
pomic(rnorm(100,10,5),rnorm(100*5,10.5,6),eps=10^-20,fullmsd=TRUE,plotting=TRUE)

[Package Pomic version 0.9.1 Index]