Lavoro per casa Espa 2018/2019 PDF

Preview

Summary

SCRIPT ls() TRUE) TRUE) TRUE) TRUE) mean, for(i in set(NULL) cumsum(table(x)) distribuzione di frequenza mean(scartz) numero casuale mean(scartv) OUTPUT ls() x x n f 1 86 8 2 73 7 3 100 10 4 93 9 5 179 18 6 115 11 7 109 11 8 113 11 9 50 5 10 58 5 976 Osservazioni mancanti: 36 TRUE) 5 TRUE) 5 TRUE) 6...

Description

SCRIPT load('~/Desktop/UNI/STATISTICA/SLIDES/homework E&M 2018-2019/WorldValueSurvey.rdata') ls() load('~/Desktop/UNI/STATISTICA/rmf_3.0-0.rdata') #1 frequenze(wvs$v120) mean(wvs$v120,na.rm = TRUE) median(wvs$v120,na.rm = TRUE) var(wvs$v120,na.rm = TRUE) sd(wvs$v120,na.rm = TRUE) #2 tapply(wvs$v120, wvs$etacat, mean, na.rm=TRUE) #3 boxplot(wvs$v120~wvs$etacat, horizzontal=TRUE) #4 cov(wvs$v208,wvs$v202,use="complete.obs") cor(wvs$v208,wvs$v202,use="complete.obs") plot(wvs$v208,wvs$v202) #5 dado sd(wvs$v120,na.rm = TRUE) [1] 2.531456 > #2 > tapply(wvs$v120, wvs$etacat, mean, na.rm=TRUE) [18,19] (19,29] (29,39] (39,49] (49,59] (59,69] (69,79] 5.173913 5.490798 5.260664 5.680723 4.946237 5.337500 5.149254

8

> #3

6 2

4

wvs$v120

> boxplot(wvs$v120~wvs$etacat, horizzontal=TRUE)

[18,19]

(19,29]

(29,39]

(39,49]

(49,59]

(59,69]

(69,79]

10

wvs$etacat

> #4

8

> cov(wvs$v208,wvs$v202,use="complete.obs") [1] 0.09038356 > cor(wvs$v208,wvs$v202,use="complete.obs") [1] 0.03862024

6 2

4

wvs$v202

> plot(wvs$v208,wvs$v202)

2

4

6 wvs$v208

> #5 > dado lanci for(i in 1:100000) {lanci[i] log=8 > sum(log==TRUE) [1] 88432 > sum(log==TRUE)/100000 [1] 0.88432 > #6 > ev Binomiale(n=262,p=0.07) ----------------------------------------------------------------------Distribuzione Binomiale ----------------------------------------------------------------------Numero delle prove

: 262

Probabilita' di successo: 0.07 Valore atteso (media) : 18.34 Varianza

: 17.0562

Somma delle probabilita': 1 x

f(x)

F(x)

[1,] 0 5.527548e-09 5.527548e-09 [2,] 1 1.090056e-07 1.145332e-07 [3,] 2 1.070717e-06 1.185250e-06 [4,] 3 6.984602e-06 8.169852e-06 [5,] 4 3.404055e-05 4.221040e-05 [6,] 5 1.322091e-04 1.744195e-04 [7,] 6 4.262440e-04 6.006635e-04 [8,] 7 1.173317e-03 1.773980e-03 [9,] 8 2.815014e-03 4.588995e-03 [10,] 9 5.979804e-03 1.056880e-02 [11,] 10 1.138735e-02 2.195614e-02 [12,] 11 1.963566e-02 4.159180e-02 [13,] 12 3.091385e-02 7.250565e-02 [14,] 13 4.474709e-02 1.172527e-01 [15,] 14 5.990336e-02 1.771561e-01 [16,] 15 7.454640e-02 2.517025e-01 [17,] 16 8.662012e-02 3.383226e-01 [18,] 17 9.434525e-02 4.326679e-01 [19,] 18 9.665598e-02 5.293238e-01 [20,] 19 9.342864e-02 6.227525e-01 [21,] 20 8.544200e-02 7.081945e-01 [22,] 21 7.411098e-02 7.823055e-01

[23,] 22 6.110715e-02 8.434126e-01 [24,] 23 4.799439e-02 8.914070e-01 [25,] 24 3.597429e-02 9.273813e-01 [26,] 25 2.577771e-02 9.531590e-01 [27,] 26 1.768619e-02 9.708452e-01 [28,] 27 1.163584e-02 9.824810e-01 [29,] 28 7.350598e-03 9.898316e-01 [30,] 29 4.464323e-03 9.942960e-01 [31,] 30 2.609789e-03 9.969058e-01 [32,] 31 1.470100e-03 9.983759e-01 [33,] 32 7.987738e-04 9.991746e-01 [34,] 33 4.190374e-04 9.995937e-01 [35,] 34 2.124342e-04 9.998061e-01 [36,] 35 1.041613e-04 9.999103e-01 [37,] 36 4.943617e-05 9.999597e-01 [38,] 37 2.272828e-05 9.999824e-01 [39,] 38 1.012933e-05 9.999926e-01 [40,] 39 4.379042e-06 9.999969e-01 [41,] 40 1.837550e-06 9.999988e-01 [42,] 41 7.488991e-07 9.999995e-01 [43,] 42 2.966070e-07 9.999998e-01 [44,] 43 1.142222e-07 9.999999e-01 [45,] 44 4.279147e-08 1.000000e+00 [46,] 45 1.560329e-08 1.000000e+00 [47,] 46 5.540300e-09 1.000000e+00 [48,] 47 1.916480e-09 1.000000e+00 [49,] 48 6.461251e-10 1.000000e+00 [50,] 49 2.123975e-10 1.000000e+00 [51,] 50 6.810423e-11 1.000000e+00 [52,] 51 2.130860e-11 1.000000e+00 [53,] 52 6.508023e-12 1.000000e+00 [54,] 53 1.940920e-12 1.000000e+00 [55,] 54 5.654252e-13 1.000000e+00 [56,] 55 1.609500e-13 1.000000e+00 [57,] 56 4.478044e-14 1.000000e+00 [58,] 57 1.218136e-14 1.000000e+00 [59,] 58 3.240684e-15 1.000000e+00 [60,] 59 8.433930e-16 1.000000e+00 [61,] 60 2.147780e-16 1.000000e+00 [62,] 61 5.353361e-17 1.000000e+00

[63,] 62 1.306309e-17 1.000000e+00 [64,] 63 3.121408e-18 1.000000e+00 [65,] 64 7.305311e-19 1.000000e+00 [66,] 65 1.674965e-19 1.000000e+00 [67,] 66 3.763076e-20 1.000000e+00 [68,] 67 8.285894e-21 1.000000e+00 [69,] 68 1.788464e-21 1.000000e+00 [70,] 69 3.784843e-22 1.000000e+00 [71,] 70 7.854567e-23 1.000000e+00 [72,] 71 1.598749e-23 1.000000e+00 [73,] 72 3.192245e-24 1.000000e+00 [74,] 73 6.253771e-25 1.000000e+00 [75,] 74 1.202229e-25 1.000000e+00 [76,] 75 2.268291e-26 1.000000e+00 [77,] 76 4.200896e-27 1.000000e+00 [78,] 77 7.637993e-28 1.000000e+00 [79,] 78 1.363551e-28 1.000000e+00 [80,] 79 2.390437e-29 1.000000e+00 [81,] 80 4.115793e-30 1.000000e+00 [82,] 81 6.960733e-31 1.000000e+00 [83,] 82 1.156471e-31 1.000000e+00 [84,] 83 1.887749e-32 1.000000e+00 [85,] 84 3.027841e-33 1.000000e+00 [86,] 85 4.772537e-34 1.000000e+00 [87,] 86 7.393315e-35 1.000000e+00 [88,] 87 1.125765e-35 1.000000e+00 [89,] 88 1.685071e-36 1.000000e+00 [90,] 89 2.479662e-37 1.000000e+00 [91,] 90 3.587659e-38 1.000000e+00 [92,] 91 5.104031e-39 1.000000e+00 [93,] 92 7.140633e-40 1.000000e+00 [94,] 93 9.824666e-41 1.000000e+00 [95,] 94 1.329510e-41 1.000000e+00 [96,] 95 1.769671e-42 1.000000e+00 [97,] 96 2.317143e-43 1.000000e+00 [98,] 97 2.984725e-44 1.000000e+00 [99,] 98 3.782485e-45 1.000000e+00 [100,] 99 4.716295e-46 1.000000e+00 [101,] 100 5.786336e-47 1.000000e+00 [102,] 101 6.985740e-48 1.000000e+00

[103,] 102 8.299524e-49 1.000000e+00 [104,] 103 9.704006e-50 1.000000e+00 [105,] 104 1.116683e-50 1.000000e+00 [106,] 105 1.264774e-51 1.000000e+00 [107,] 106 1.410008e-52 1.000000e+00 [108,] 107 1.547311e-53 1.000000e+00 [109,] 108 1.671478e-54 1.000000e+00 [110,] 109 1.777501e-55 1.000000e+00 [111,] 110 1.860903e-56 1.000000e+00 [112,] 111 1.918048e-57 1.000000e+00 [113,] 112 1.946406e-58 1.000000e+00 [114,] 113 1.944739e-59 1.000000e+00 [115,] 114 1.913189e-60 1.000000e+00 [116,] 115 1.853262e-61 1.000000e+00 [117,] 116 1.767711e-62 1.000000e+00 [118,] 117 1.660326e-63 1.000000e+00 [119,] 118 1.535658e-64 1.000000e+00 [120,] 119 1.398702e-65 1.000000e+00 [121,] 120 1.254570e-66 1.000000e+00 [122,] 121 1.108187e-67 1.000000e+00 [123,] 122 9.640231e-69 1.000000e+00 [124,] 123 8.258961e-70 1.000000e+00 [125,] 124 6.968409e-71 1.000000e+00 [126,] 125 5.790523e-72 1.000000e+00 [127,] 126 4.738959e-73 1.000000e+00 [128,] 127 3.819735e-74 1.000000e+00 [129,] 128 3.032299e-75 1.000000e+00 [130,] 129 2.370840e-76 1.000000e+00 [131,] 130 1.825684e-77 1.000000e+00 [132,] 131 1.384661e-78 1.000000e+00 [133,] 132 1.034322e-79 1.000000e+00 [134,] 133 7.609614e-81 1.000000e+00 [135,] 134 5.513948e-82 1.000000e+00 [136,] 135 3.935083e-83 1.000000e+00 [137,] 136 2.765883e-84 1.000000e+00 [138,] 137 1.914692e-85 1.000000e+00 [139,] 138 1.305404e-86 1.000000e+00 [140,] 139 8.765302e-88 1.000000e+00 [141,] 140 5.796410e-89 1.000000e+00 [142,] 141 3.774982e-90 1.000000e+00

[143,] 142 2.421180e-91 1.000000e+00 [144,] 143 1.529281e-92 1.000000e+00 [145,] 144 9.512329e-94 1.000000e+00 [146,] 145 5.826610e-95 1.000000e+00 [147,] 146 3.514504e-96 1.000000e+00 [148,] 147 2.087468e-97 1.000000e+00 [149,] 148 1.220874e-98 1.000000e+00 [150,] 149 7.030799e-100 1.000000e+00 [151,] 150 3.986639e-101 1.000000e+00 [152,] 151 2.225682e-102 1.000000e+00 [153,] 152 1.223369e-103 1.000000e+00 [154,] 153 6.620243e-105 1.000000e+00 [155,] 154 3.526914e-106 1.000000e+00 [156,] 155 1.849703e-107 1.000000e+00 [157,] 156 9.549404e-109 1.000000e+00 [158,] 157 4.852858e-110 1.000000e+00 [159,] 158 2.427420e-111 1.000000e+00 [160,] 159 1.195078e-112 1.000000e+00 [161,] 160 5.790666e-114 1.000000e+00 [162,] 161 2.761327e-115 1.000000e+00 [163,] 162 1.295804e-116 1.000000e+00 [164,] 163 5.983659e-118 1.000000e+00 [165,] 164 2.718775e-119 1.000000e+00 [166,] 165 1.215431e-120 1.000000e+00 [167,] 166 5.345757e-122 1.000000e+00 [168,] 167 2.313018e-123 1.000000e+00 [169,] 168 9.844836e-125 1.000000e+00 [170,] 169 4.121589e-126 1.000000e+00 [171,] 170 1.697125e-127 1.000000e+00 [172,] 171 6.872593e-129 1.000000e+00 [173,] 172 2.736835e-130 1.000000e+00 [174,] 173 1.071668e-131 1.000000e+00 [175,] 174 4.125875e-133 1.000000e+00 [176,] 175 1.561621e-134 1.000000e+00 [177,] 176 5.810285e-136 1.000000e+00 [178,] 177 2.124896e-137 1.000000e+00 [179,] 178 7.637509e-139 1.000000e+00 [180,] 179 2.697696e-140 1.000000e+00 [181,] 180 9.362972e-142 1.000000e+00 [182,] 181 3.192744e-143 1.000000e+00

[183,] 182 1.069530e-144 1.000000e+00 [184,] 183 3.519223e-146 1.000000e+00 [185,] 184 1.137290e-147 1.000000e+00 [186,] 185 3.609185e-149 1.000000e+00 [187,] 186 1.124610e-150 1.000000e+00 [188,] 187 3.440243e-152 1.000000e+00 [189,] 188 1.033017e-153 1.000000e+00 [190,] 189 3.044336e-155 1.000000e+00 [191,] 190 8.803938e-157 1.000000e+00 [192,] 191 2.497993e-158 1.000000e+00 [193,] 192 6.952859e-160 1.000000e+00 [194,] 193 1.898101e-161 1.000000e+00 [195,] 194 5.081380e-163 1.000000e+00 [196,] 195 1.333740e-164 1.000000e+00 [197,] 196 3.431665e-166 1.000000e+00 [198,] 197 8.653617e-168 1.000000e+00 [199,] 198 2.138262e-169 1.000000e+00 [200,] 199 5.176103e-171 1.000000e+00 [201,] 200 1.227237e-172 1.000000e+00 [202,] 201 2.849307e-174 1.000000e+00 [203,] 202 6.476388e-176 1.000000e+00 [204,] 203 1.440798e-177 1.000000e+00 [205,] 204 3.136463e-179 1.000000e+00 [206,] 205 6.679275e-181 1.000000e+00 [207,] 206 1.391080e-182 1.000000e+00 [208,] 207 2.832597e-184 1.000000e+00 [209,] 208 5.637665e-186 1.000000e+00 [210,] 209 1.096382e-187 1.000000e+00 [211,] 210 2.082733e-189 1.000000e+00 [212,] 211 3.863398e-191 1.000000e+00 [213,] 212 6.995502e-193 1.000000e+00 [214,] 213 1.236017e-194 1.000000e+00 [215,] 214 2.130207e-196 1.000000e+00 [216,] 215 3.579642e-198 1.000000e+00 [217,] 216 5.862716e-200 1.000000e+00 [218,] 217 9.354316e-202 1.000000e+00 [219,] 218 1.453393e-203 1.000000e+00 [220,] 219 2.197894e-205 1.000000e+00 [221,] 220 3.233461e-207 1.000000e+00 [222,] 221 4.625299e-209 1.000000e+00

[223,] 222 6.429627e-211 1.000000e+00 [224,] 223 8.680724e-213 1.000000e+00 [225,] 224 1.137595e-214 1.000000e+00 [226,] 225 1.446118e-216 1.000000e+00 [227,] 226 1.782018e-218 1.000000e+00 [228,] 227 2.127178e-220 1.000000e+00 [229,] 228 2.457832e-222 1.000000e+00 [230,] 229 2.746697e-224 1.000000e+00 [231,] 230 2.966279e-226 1.000000e+00 [232,] 231 3.092894e-228 1.000000e+00 [233,] 232 3.110669e-230 1.000000e+00 [234,] 233 3.014632e-232 1.000000e+00 [235,] 234 2.812105e-234 1.000000e+00 [236,] 235 2.521952e-236 1.000000e+00 [237,] 236 2.171719e-238 1.000000e+00 [238,] 237 1.793262e-240 1.000000e+00 [239,] 238 1.417822e-242 1.000000e+00 [240,] 239 1.071643e-244 1.000000e+00 [241,] 240 7.730044e-247 1.000000e+00 [242,] 241 5.311323e-249 1.000000e+00 [243,] 242 3.469139e-251 1.000000e+00 [244,] 243 2.149119e-253 1.000000e+00 [245,] 244 1.259620e-255 1.000000e+00 [246,] 245 6.965638e-258 1.000000e+00 [247,] 246 3.623179e-260 1.000000e+00 [248,] 247 1.766558e-262 1.000000e+00 [249,] 248 8.042343e-265 1.000000e+00 [250,] 249 3.403505e-267 1.000000e+00 [251,] 250 1.332124e-269 1.000000e+00 [252,] 251 4.793662e-272 1.000000e+00 [253,] 252 1.574979e-274 1.000000e+00 [254,] 253 4.685643e-277 1.000000e+00 [255,] 254 1.249664e-279 1.000000e+00 [256,] 255 2.950924e-282 1.000000e+00 [257,] 256 6.073390e-285 1.000000e+00 [258,] 257 1.067246e-287 1.000000e+00 [259,] 258 1.556789e-290 1.000000e+00 [260,] 259 1.809694e-293 1.000000e+00 [261,] 260 1.571694e-296 1.000000e+00 [262,] 261 9.065102e-300 1.000000e+00

0.08

[263,] 262 2.604273e-303 1.000000e+00 > pbinom(ev,size=262,prob=0.07,lower.tail=FALSE) [1] 0.4706762

0.06

> pbinom(19,size=262,prob=0.07,lower.tail=FALSE)

0.00

0.02

0.04

f(x)

[1] 0.3772475

0

50

100

150 x

> #7 > ProbNorm((11.1-23.5),30.2,11.1,23,5) La probabilita' di osservare un valore compreso fra -12.4 e 30.2 e' uguale a 0.643401303884524

200

250

0.010 0.005 0.000

dnorm(x, mean = mu, sd =

-50

0

50 x

> #8 > set.seed(200473); x cumsum(table(x)) -3 -2 -1 0 1 2 3 4 5 32 73 102 138 192 222 247 275 300 > z mediana scartz mean(scartz) [1] 2.03 > scartv mean(scartv) [1] 3.436667

COMMENTI: #1 Ho utilizzato le funzioni di R per calcolare la varianza e la deviazione standard #2 Si può notare la bassa variazione della media rispetto alle fasce di età #3

Si nota che non c'è un considerevole cambiamento di opinione al variare dell'età, infatti la mediana si trova sempre tra 5 e 6 #4 Si nota una correlazione quasi nulla tra le due variabili prese in considerazione. Per verificare questa affermazione graficamente, ho plottato lo scatter #5 Avendo ripetuto l’esperimento un numero considerevole di volte, il risultato della simulazione è un’ottima approssimazione della probabilità dell’evento #6 Si può notare la concentrazione della probabilità attorno al valore atteso, infatti nell'intervallo [18.34,19] è concentrato circa il 10% della probabilità totale...