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