Après vous avoir contacté par mail, je fais part de mon problème d'ajustement par mail afin que cette question soit accessible à toute la communauté RGeoS.
Je résume le problème : je lance de nombreux ajustement de variogrammes afin de kriger des données de pollution atmosphérique. Une carte est ainsi produite chaque jour. Les structures variographiques utilisées sont les suivantes (j'ai corrigé l'ordre sur vos conseils) :
1) Nugget Effect
2) Exponential
3) Spherical
4) Linear
Il arrive que la procédure d'ajustement plante et me renvoie le type de message suivant :
- Code: Select all
Error in parameter 2: lower bound (0.000009) > default value (0.000009)
Error in model.auto(data.vario, struct = struct_) :
Error in R_model_auto_f
Le message d'erreur fait apparaître un pb avec le paramètre 2 mais le même message apparaît parfois avec le paramètre 4.
Suite à vos conseils pour obtenir plus d'informations permettant de démarrer un debug a distance, j'ai placé l'ordre suivant avant le lancement:
debug.define("converge",1).
Je fournis la sortie obtenue avec cette commande sur l'un des variogrammes qui pose problème :
Je précise que je travaille sous une architecture SUSE Linux.
- Code: Select all
Optimization Conditions
=======================
- Number of variables = 1
- Space dimension = 2
- Number of experiments = 20
- Number of parameters = 6
- Constrained Minimization (OFF)
- Anisotropy (OFF)
- Anisotropy Rotation (OFF)
- Global Rotation (OFF)
- Rotation around Z only (OFF)
- Lock third dimension (OFF)
- Consider Isotropic 2-D (OFF)
- Use the Goulard option (OFF)
Model
-----
Structure : Nugget Effect
1 - Sill : 17.431838 [0.000000,NA[
Structure : Exponential
2 - Sill : 17.431838 [0.000000,NA[
3 - Range U : 2.371696 [0.000009,NA[
Structure : Spherical
4 - Sill : 17.431838 [0.000000,NA[
5 - Range U : 3.557544 [0.000009,NA[
Structure : Linear
6 - Sill : 17.431838 [0.000000,NA[
Trajectory of parameters in Foxleg Algorithm
--------------------------------------------
Iteration Score Delta Par-1 Par-2 Par-3 Par-4 Par-5 Par-6
1 0.00897 1 17.432 17.432 2.372 17.432 3.558 17.432
2 0.00236 3 34.864 26.534 3.558 0.000 2.372 34.864
3 0.00116 3 57.408 0.000 7.115 7.427 0.000 21.163
4 0.000236 3 54.407 0.000 7.240 0.000 0.000 34.943
5 0.000236 1.5 2.111 0.000 10.798 52.296 0.000 34.943
Convergence has not been reached
Statistics for the Minimization procedure
-----------------------------------------
Number of experimental values = 20
Number of parameters = 6
Number of iterations = 5/1000
Value of minimized function = 0.000236251
Initial increment value = 1
Current increment value = 0.75
Increment Stopping Criterion = 1e-05
Stopping Value = 1.72881e-08
Stopping Criterion = 1e-05
Model
-----
Structure : Nugget Effect
1 - Sill : 2.111069 [0.000000,NA[
Structure : Exponential
2 - Sill : 0.000000 [0.000000,NA[
3 - Range U : 10.797805 [0.000009,NA[
Structure : Spherical
4 - Sill : 52.295515 [0.000000,NA[
5 - Range U : 0.000009 [0.000009,NA[
Structure : Linear
6 - Sill : 34.943023 [0.000000,NA[
Suppressing the unnecessary basic structures
============================================
Structure 'Nugget Effect' in model #1 is suppressed
Structure 'Exponential' in model #1 is suppressed
Optimization Conditions
=======================
- Number of variables = 1
- Space dimension = 2
- Number of experiments = 20
- Number of parameters = 3
- Constrained Minimization (OFF)
- Anisotropy (OFF)
- Anisotropy Rotation (OFF)
- Global Rotation (OFF)
- Rotation around Z only (OFF)
- Lock third dimension (OFF)
- Consider Isotropic 2-D (OFF)
- Use the Goulard option (OFF)
Model
-----
Structure : Spherical
1 - Sill : 52.295515 [0.000000,NA[
2 - Range U : 0.000009 [0.000009,NA[
Structure : Linear
3 - Sill : 34.943023 [0.000000,NA[
Error in parameter 2: lower bound (0.000009) > default value (0.000009)
Error in model.auto(data.vario, struct = struct_) :
Error in R_model_auto_fit
En espérant faire avance le pb,
Maxime