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
1 2  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)

Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.