pomdev {Pomic}

POMDEV calculation for numerical patterns

Description

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

Usage

pomdev(object1, object2, eps=10^-30, nrange=1000)

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

Details

This function returns the POMDEV value as computed in Piou et al. It corresponds to twice the sum of the log of an approximated likelihood given by the approximating function of density from the simulation results applied on the field data.

Value

Return a value corresponding to the pomdev score of model results to reproduce the pattern.

Author(s)

Cyril Piou
[e-mail]

See Also

pomdev.extra pomdev.corrected

Examples

pomdev(rnorm(100,10,5),rnorm(100*5,10.5,6),eps=10^-20)