Utilisation function sobol

Postez ici vos questions, réponses, commentaires ou suggestions - Les sujets seront ultérieurement répartis dans les archives par les modérateurs

Modérateur : Groupe des modérateurs

Charlotte Neau
Messages : 4
Enregistré le : 08 Juin 2018, 06:58

Utilisation function sobol

Messagepar Charlotte Neau » 05 Juil 2018, 08:48

Bonjour,
Je suis débutante sur R, et j'aimerais savoir comment utiliser la fonction sobol se trouvant dans le package 'Sensitivity'.

Voici mon jeu de données :

Code : Tout sélectionner

> Data_test
           V1        V2       V3
1   -3.00e+01 -3.00e+01 -0.82510
2   -2.70e+01 -3.00e+01 -0.44470
3   -2.40e+01 -3.00e+01  0.00969
4   -2.10e+01 -3.00e+01  0.62030
5   -1.80e+01 -3.00e+01  0.63180
6   -1.50e+01 -3.00e+01  0.04380
7   -1.20e+01 -3.00e+01 -0.43220
8   -9.00e+00 -3.00e+01 -0.76970
9   -6.00e+00 -3.00e+01 -1.08200
10  -3.00e+00 -3.00e+01 -1.37700
11   4.47e-08 -3.00e+01 -1.66700
12  -3.00e+01 -2.70e+01 -0.95950
13  -2.70e+01 -2.70e+01 -0.58140
14  -2.40e+01 -2.70e+01 -0.16720
15  -2.10e+01 -2.70e+01  0.56540
16  -1.80e+01 -2.70e+01  0.66190
17  -1.50e+01 -2.70e+01  0.50180
18  -1.20e+01 -2.70e+01 -0.26960
19  -9.00e+00 -2.70e+01 -0.63860
20  -6.00e+00 -2.70e+01 -0.95400
21  -3.00e+00 -2.70e+01 -1.25300
22   4.47e-08 -2.70e+01 -1.53900
23  -3.00e+01 -2.40e+01 -1.09300
24  -2.70e+01 -2.40e+01 -0.71960
25  -2.40e+01 -2.40e+01 -0.31740
26  -2.10e+01 -2.40e+01  0.45650
27  -1.80e+01 -2.40e+01  0.66300
28  -1.50e+01 -2.40e+01  0.61010
29  -1.20e+01 -2.40e+01 -0.06570
30  -9.00e+00 -2.40e+01 -0.49600
31  -6.00e+00 -2.40e+01 -0.82940
32  -3.00e+00 -2.40e+01 -1.13200
33   4.47e-08 -2.40e+01 -1.41800
34  -3.00e+01 -2.10e+01 -1.21600
35  -2.70e+01 -2.10e+01 -0.85690
36  -2.40e+01 -2.10e+01 -0.46820
37  -2.10e+01 -2.10e+01  0.05110
38  -1.80e+01 -2.10e+01  0.63990
39  -1.50e+01 -2.10e+01  0.65960
40  -1.20e+01 -2.10e+01  0.40340
41  -9.00e+00 -2.10e+01 -0.33630
42  -6.00e+00 -2.10e+01 -0.69930
43  -3.00e+00 -2.10e+01 -1.01200
44   4.47e-08 -2.10e+01 -1.30300
45  -3.00e+01 -1.80e+01 -1.33100
46  -2.70e+01 -1.80e+01 -0.98290
47  -2.40e+01 -1.80e+01 -0.61170
48  -2.10e+01 -1.80e+01 -0.18100
49  -1.80e+01 -1.80e+01  0.59240
50  -1.50e+01 -1.80e+01  0.68230
51  -1.20e+01 -1.80e+01  0.55680
52  -9.00e+00 -1.80e+01 -0.17560
53  -6.00e+00 -1.80e+01 -0.55760
54  -3.00e+00 -1.80e+01 -0.88740
55   4.47e-08 -1.80e+01 -1.18800
56  -3.00e+01 -1.50e+01 -1.44700
57  -2.70e+01 -1.50e+01 -1.09900
58  -2.40e+01 -1.50e+01 -0.74260
59  -2.10e+01 -1.50e+01 -0.34840
60  -1.80e+01 -1.50e+01  0.46130
61  -1.50e+01 -1.50e+01  0.68640
62  -1.20e+01 -1.50e+01  0.63460
63  -9.00e+00 -1.50e+01  0.02080
64  -6.00e+00 -1.50e+01 -0.42350
65  -3.00e+00 -1.50e+01 -0.75850
66   4.47e-08 -1.50e+01 -1.06800
67  -3.00e+01 -1.20e+01 -1.55300
68  -2.70e+01 -1.20e+01 -1.22100
69  -2.40e+01 -1.20e+01 -0.86220
70  -2.10e+01 -1.20e+01 -0.49010
71  -1.80e+01 -1.20e+01 -0.00469
72  -1.50e+01 -1.20e+01  0.65050
73  -1.20e+01 -1.20e+01  0.67880
74  -9.00e+00 -1.20e+01  0.45220
75  -6.00e+00 -1.20e+01 -0.29170
76  -3.00e+00 -1.20e+01 -0.63730
77   4.47e-08 -1.20e+01 -0.94760
78  -3.00e+01 -9.00e+00 -1.66800
79  -2.70e+01 -9.00e+00 -1.33100
80  -2.40e+01 -9.00e+00 -0.98740
81  -2.10e+01 -9.00e+00 -0.61770
82  -1.80e+01 -9.00e+00 -0.20530
83  -1.50e+01 -9.00e+00  0.57500
84  -1.20e+01 -9.00e+00  0.69500
85  -9.00e+00 -9.00e+00  0.57710
86  -6.00e+00 -9.00e+00 -0.14670
87  -3.00e+00 -9.00e+00 -0.52130
88   4.47e-08 -9.00e+00 -0.83590
89  -3.00e+01 -6.00e+00 -1.77800
90  -2.70e+01 -6.00e+00 -1.44500
91  -2.40e+01 -6.00e+00 -1.10300
92  -2.10e+01 -6.00e+00 -0.74620
93  -1.80e+01 -6.00e+00 -0.35400
94  -1.50e+01 -6.00e+00  0.43500
95  -1.20e+01 -6.00e+00  0.68210
96  -9.00e+00 -6.00e+00  0.64730
97  -6.00e+00 -6.00e+00  0.05320
98  -3.00e+00 -6.00e+00 -0.40390
99   4.47e-08 -6.00e+00 -0.72800
100 -3.00e+01 -3.00e+00 -1.87700
101 -2.70e+01 -3.00e+00 -1.56000
102 -2.40e+01 -3.00e+00 -1.21800
103 -2.10e+01 -3.00e+00 -0.87010
104 -1.80e+01 -3.00e+00 -0.49280
105 -1.50e+01 -3.00e+00 -0.01310
106 -1.20e+01 -3.00e+00  0.64180
107 -9.00e+00 -3.00e+00  0.68390
108 -6.00e+00 -3.00e+00  0.46570
109 -3.00e+00 -3.00e+00 -0.28160
110  4.47e-08 -3.00e+00 -0.62090
111 -3.00e+01  4.47e-08 -1.97200
112 -2.70e+01  4.47e-08 -1.66800
113 -2.40e+01  4.47e-08 -1.33500
114 -2.10e+01  4.47e-08 -0.98480
115 -1.80e+01  4.47e-08 -0.62550
116 -1.50e+01  4.47e-08 -0.20570
117 -1.20e+01  4.47e-08  0.56930
118 -9.00e+00  4.47e-08  0.69260
119 -6.00e+00  4.47e-08  0.58810
120 -3.00e+00  4.47e-08 -0.14830
121  4.47e-08  4.47e-08 -0.51540


Et lorsque j'utilise la fonction 'sobol', le message d'erreur suivant d'affiche :

Code : Tout sélectionner

> sobol(model = NULL, Data_test$V1, Data_test$V2, order = 1)
Error in if ((ncol(X1) != ncol(X2)) | (nrow(X1) != nrow(X2))) { :
  argument is of length zero


Pourriez-vous m'aider svp ??

Merci

Eric Casellas
Messages : 767
Enregistré le : 06 Jan 2009, 14:59

Re: Utilisation function sobol

Messagepar Eric Casellas » 05 Juil 2018, 09:43

Bonjour,

D'après l'aide la la fonction

Code : Tout sélectionner

?sobol
, X1 et X2 correspondent à tes Data_test$V1 et Data_test$V2 respectivement.
Et donc si tu regarde ce que donne

Code : Tout sélectionner

((ncol(X1) != ncol(X2)) | (nrow(X1) != nrow(X2)))
, ça donne un

Code : Tout sélectionner

logical(0)
car tes X1 et X2 sont des vecteurs et non pas des data.frame/matrix comme attendus et c'est ce que te dis le message d'erreur...

Eric
Eric


Retourner vers « Questions en cours »

Qui est en ligne

Utilisateurs parcourant ce forum : Google [Bot] et 1 invité