
1 - Problemes ressolts d'estadística descriptiva
Enunciat 1.1
Es disposa de les dades de l’alçada en centímetres de 50 persones, agafats a l’atzar:

- Agrupis les dades en intervals d’igual longitud —tancats per l’esquerra i oberts per la dreta—. Com afecta la distribució dels intervals a la precisió de les dades obtingudes? Calculis els valors d’interès descriptiu.
- Representis gràficament.
Solució
a) Agrupis les dades en intervals d’igual longitud —tancats per l’esquerra i oberts per la dreta—. Com afecta la distribució dels intervals a la precisió de les dades obtingudes? Calculis els valors d’interès descriptiu.
Amb R:

Al agrupar-se les dades en intervals, hi ha una pèrdua d’informació. Pèrdua que augmenta amb l’amplitud dels intervals.
Calculem ara els valors descriptius d'interès.
Mitjana aritmètica \(\bar x\) : mean(alçades) \(\rightarrow\) 174.1
Mitjana geomètrica : geometric(alçades) \(\rightarrow\) 173.9395
Mediana (\(M_e\)) : median(alçades) \(\rightarrow\) 173.5
Moda (\(M_o\)) : mfv(alçades) \(\rightarrow\) 176
Variància mostral (\(S_N^2\)) : variancia(alçades) \(\rightarrow\) 56.29
Variància mostral corregida (\(S^2_{(N-1)}\))[La quasivariància. Per defecte amb R.] : var(alçades) \(\rightarrow\) 47.43878
Desviació estàndard (\(S_N\)) : sd(alçades) \(\rightarrow\) 7.578837
Rang(\(R\)) : range(alçades) \(\rightarrow\) 160 \(\quad\) 192
Quantils (a la imatge anterior).
Si es vol calcular —per exemple— únicament el quantil del 25 %: quantile(alçades, c(.25)) \(\rightarrow\) 25% \(\quad\) 168.25
Rang interquartíl·lic (\(IQR\)) : IQR(alçades) \(\rightarrow\) 10.75
Skewness : skewness(alçades) \(\rightarrow\) 0.3157094
Kurtosis : kurtosis(alçades) \(\rightarrow\) 2.459508
El coeficient és positiu. La distribució s’anomena leptocúrtica. Hi ha una major concentració de dades en torn a la mitjana.
a) Agrupis les dades en intervals d’igual longitud —tancats per l’esquerra i oberts per la dreta—. Com afecta la distribució dels intervals a la precisió de les dades obtingudes? Calculis els valors d’interès descriptiu.
Es representen les dades amb un histograma. S’ajusten les dades a una funció de densitat normal —de mitjana i desviació estàndard calculades anteriorment—.
Amb R:

Enunciat 1.2
Una bàscula registra les següents mesures (en Kg):

Calculis la mitjana aritmètica \(\bar x\), la desviació estàndard \(S_N\), la mediana \(M_e\) i els quartils.
Solució
Amb R:

Enunciat 1.3
Es pretén estudiar la diversitat comparativa entre una zona boscosa en recuperació després d’un incendi i altra que no ha estat cremada.
L’estudi recull les dades de l’Índex de Comparació Seqüencial (ICS).Un bon indicador de la diversitat específica. A més diversitat major és l’ICS.
S’obtenen les següents dades estadístiques:

Extreguis la informació pertinent de les dades estadístiques recollides.
Solució
S’observa que a la part cremada (A) l’ICS és major. Hi ha més diversitat específica —per terme mitjà—.
Per tal de comparar la variabilitat entre A i B, es calcula el coeficient de variació: \[\begin{align*} &V_A = \frac{\sigma_A}{ICS_A} = \frac{0.74749}{1.3403330} \cdot 100 \approx 55.8\% \\ &V_B = \frac{\sigma_B}{ICS_B} = \frac{0.62651}{1.00684147} \cdot 100 \approx 62.2\% \end{align*}\]
A la part cremada les dades són més homogènies (\(V_A < V_B\)). En tot dos casos el coeficient de variació és inferior al 80 %, per tant la mitjana és representativa de les dades —més en la part cremada—.
Respecte la forma, l’índex de la zona cremada és quasi leptocúrtic. La curtosis és - 0.88523.
L’índex a la zona no cremada (B) mostra asimetria positiva i curtosis negativa —menor que a la zona cremada—.
Enunciat 1.4
Es llança una moneda 50 vegades i se n’obtenen 22 cares i 28 creus. Quina és la distribució de freqüències i quina la distribució de probabilitat?
Solució
Distribució de freqüències o empírica: 22 cares i 28 creus.
Distribució de probabilitat o teòrica: 25 cares i 25 creus.
S’observa que la moneda no és ideal —encara que són poques dades…— la distribució empírica és diferent de la teòrica.
Amb R:

Enunciat 1.5
Es disposa d’una mostra \(X = \{4, 7, 6, 7, 5, 5, 4\}\). Calculis la seva mitjana aritmètica, \(\overline X\).
Solució
dades <- c(4,7,6,7,5,5,4); mean(dades) \(\rightarrow\) 5.428571
Calculat com: \[\frac{4 + 7 + 6 + 7 + 5 + 5 + 4}{7} \approx 5.43\] o bé com: \[ \frac{(4 \cdot 2) + (7 \cdot 2) + (5 \cdot 2) + 6}{7} \approx 5.43\]
S’observa a més que \[\displaystyle\sum_{i=1}^{n} \left( x_i - \overline{X}_n \right) = 0\]
Enunciat 1.6
Es disposa de la mostra \(X = \{4, 7, 6, 7, 5, 8, 9\}\). Calculis la seva mediana.
Solució
dades <- c(4,7,6,7,5,8,9); median(dades) \(\rightarrow\) 7
Enunciat 1.7
Calculis la mediana dels elements registrats a la taula següent:

Solució
Els elements són agrupats en classes. Com \(n\) és senar, la mitjana serà el valor que ocupa la posició 31 —dins la classe \([11,13)\) — \[\tilde{x} = 11 + \left( \frac{\frac{61}{2} - (5 + 11 + 10)}{7} \right) \cdot 2 = \frac{86}{7} \approx 12.29\]
Enunciat 1.8
Amb les dades de l’exercici anterior, calculis la seva moda.
Solució
Es té que \[D = 7 + \left( \frac{11 - 5}{2 \cdot 11 - 5 - 10} \right) \cdot 2 = 12\]
S’observa que \(D = 12, \tilde x = 12.29, \bar x = 12.2951\). S’acompleix així la inequació: \(D < \tilde x < \bar x\)
La moda és sempre menor que la mediana i la mediana és sempre menor que la mitjana.
Enunciat 1.9
Es disposa de la mostra ordenada \(X = \{4, 5, 6, 7, 7, 8, 9\}\). Calculis la seva moda.
Solució
Amb R:

La moda de la mostra és 7, el valor central.
Enunciat 1.10
Un empleat obté uns augments de sou en tres anys consecutius, del 6, 10 i 12 % respectivament, sobre el sou de l’any anterior. Calculis l’augment mitjà del sou.
Solució
Amb R:

L’augment mitjà aproximat ha estat del 9.31 % —donat que \(1 + 0.0931 = 1.0931\)—.
Enunciat 1.11
En un cultiu bacterià, el nombre de bactèries per unitat de superfície s’incrementa en tres dies de 200 a 5000. Calculis l’increment mitjà diari en %.
Solució
Sigui \(x\) l’increment buscat, l’increment serà igual a \(200 \cdot (1 + x)^i\), on \(i\) és el número del dia.
Aleshores, segons l’enunciat, l’increment és \(200 \cdot (1 + x)^3 = 5000 \implies x \approx 1.924 = 192.4 \text{ %}\).
Enunciat 1.12
Es disposa de la mostra \(X = \{4, 7, 6, 7, 5, 8, 9\}\). Calculis la seva mitjana geomètrica.
Solució
Amb R:

La mitjana geomètrica de la mostra és 6.37.
O també, amb igual resultat:
Amb R:

Enunciat 1.13
En un determinat comerç es ven un article a diferents preus segons la quantitat: 10 articles a 1 €; 5 articles a 0.75 €; 1 article 0.20 €. Quin és el nombre mitjà d’articles per cada cèntim d’euro?
Solució
\[\overline{x}_H = \frac{3}{\frac{100}{10} + \frac{75}{5} + \frac{20}{1}} = \frac{1}{15} \approx 0.07\]
El nombre mitjà d’articles per cada cèntim d’euro és 0.07.
També es pot fer el següent raonament: Respectivament, els articles surten a 10, 15 i 20 cèntims d’euro. Per tant 3 articles a un preu mitjà de \[\frac{3}{10 + 15 + 20} = \frac{1}{15} \approx 0.07\]
Enunciat 1.14
Un mòbil es trasllada de \(A\) fins a \(B\) amb una velocitat mitjana de 20 km/h; per tornar de \(B\) fins \(A\) ho fa en canvi a una velocitat mitjana de 50 km/h. Calculis la velocitat mitjana global.
Solució
La mitjana global és \[\overline{x}_G = \frac{1 + 1}{\frac{1}{20} + \frac{1}{50}} = \frac{200}{7} \approx 28.57 \text{ km/h}\]
Enunciat 1.15
Es disposa de la mostra \(X = \{4, 7, 6, 7, 5, 8, 9\}\). Calculis la seva mitjana harmònica.
Solució
Amb R:

La mitjana harmònica de la mostra és 6.15.
Enunciat 1.16
Es disposa de la mostra \(X = \{4, 7, 6, 7, 5, 8, 9\}\). Calculis la seva mitjana quadràtica.
Solució
Amb R:

La mitjana quadràtica de la mostra és 6.76.
Enunciat 1.17
Es disposa de la mostra \(X = \{4, 7, 6, 7, 5, 8, 9\}\). Calculis la seva mitjana ponderada.
Solució
Amb R:

La mitjana ponderada de la mostra és 6.44.
Com més gran és el pes d’un element, més importància adquireix en el càlcul. Si tots els elements tenen el mateix pes la mitjana ponderada coincideix amb la mitjana —en aquest cas igual a 6.571429—.
Amb R:

Enunciat 1.18
Generis un vector de dades amb els nombres naturals del 0 al 100 i calculis els seus quantils.
Solució
Amb R:

Enunciat 1.19
Es disposa d’una mostra \(M\) = “alçades de la població \(X\)” amb els valors següents:

Calculis:
- Les freqüències i la distribució per quartils.
- El percentil al 43 %.
- El rang de distribució.
- El rang interquartíl·lic.
- La variància i la quasivariància i la desviació típica.
Solució
a) Les freqüències i la distribució per quartils.
Amb R:

b) El percentil al 43 %.
quantile(alçades, c(.43)) \(\rightarrow\) 172.07
El 43 % de les dades són inferiors a 172.07.
c) El rang de distribució.
range(alçades) \(\rightarrow\) 160 \(\quad\) 188
En aquest exemple —aproximat a una distribució normal— es podria estimar la desviació tipus de la població a partir de \[\hat{s} \approx \frac{R}{6} = \frac{188 - 160}{6} = \frac{28}{6} = 4.\overline{6}\]
Això és així donat que en una distribució normal el 99.7 % dels seus valors són dins l’interval \(6 \sigma\).
d) El rang interquartíl·lic.
Q1 <- quantile(alçades, c(.25)) \(\rightarrow\) 25% \(\quad\) 168
Q3 <- quantile(alçades, c(.75)) \(\rightarrow\) 75% \(\quad\) 177.75
IQR(alçades) \(\rightarrow\) 9.75
Q3 - Q1 \(\rightarrow\) 75%
e) La variància i la quasivariància i la desviació típica.
Amb R:

La variància de la mostra és 46.04 i la quasivariància 46.98.
La desviació típica de la mostra, calculadad de dues maneres, és 6.85.
Enunciat 1.20
Es vol concedir una beca d’estudis i n’hi ha dos candidats possibles. Un alumne amb 9 punts, que prové d’un centre de mitjana 6 i desviació típica 1. Un altre amb 10 punts, que prové d’un centre de mitjana 7 i desviació típica 2. Trobis quin alumne la mereix més.
Solució
Al normalitzar les dades es poden comparar millor els dos alumnes: \[\begin{align*} &z_1 = \frac{9 - 6}{1} = 3 \\ &z_2 = \frac{10 - 7}{2} = 1.5 \end{align*}\]
Ara queda evident quin és el millor candidat. El primer alumne té millors notes —malgrat tenir menys punts—.
Enunciat 1.21
Calculis el coeficient de variació relativa d’una mostra amb: \(n = 80; \quad S = 2; \quad \bar x = 10.\)
Solució
Es té \[\begin{align*} &V(x) = r = \frac{S}{|\overline{x}|} = \frac{2}{10} = 0.2 \\ &V_r = \frac{r}{\sqrt{n-1}} \cdot 100 = \frac{0.2}{\sqrt{80 - 1}} \cdot 100 \approx 2.25\% \end{align*}\]
El coeficient de variació relativa de la mostra és 2.25 %.
Enunciat 1.22
Es disposa de cinc mostres, \(M_x\) = “alçades de la població \(X\)”:

Trobis en quina de les mostres hi ha més variabilitat.
Solució
Amb R:

S’observa més variabilitat a les mostres 1 i 3.
Enunciat 1.23
Es disposa de la mostra:
\(X = \{ 3, 4, 4, 5, 5, 5, 5, 5, 5, 6, 7, 7, 7, 7, 7, 8, 8, 9, 9, 9, 10, 11, 12, 14 \}\)
Trobis entre quin interval de valors es troba la massa central de les dades de la mostra.
Solució
Amb R:

S’observa que la massa central de les dades es troba entre \(6.693/1.463 \approx 4.58\) i \(6.693 \cdot 1.463 \approx 9.79\).
Enunciat 1.24
Es disposa d’una mostra \(M\) = “alçades de la població \(X\)” amb:

Calculis el coeficient d’asimetria de la mostra \(M\).
Solució
Amb R:

El coeficient d’asimetria de Fisher de la mostra és 0.045. El que ens indica una distribució quasi normal.
Enunciat 1.25
R disposa de la base de dades de 272 observacions d’erupcions volcàniques del Guèiser “Old Faithful” al parc nacional de Yellowstone. Amb dues variables: duració de l’erupció i temps d’espera fins la propera erupció —ambdues en minuts—. Analitzis la distribució de les dades.
Solució
Amb R:

La curtosis \(-1.51\) ens indica una distribució platicúrtica. Les erupcions s’agrupen en una distribució bimodal —respecte a durada i temps d’espera—. Hi ha una relació directament proporcional entre la duració de l’erupció i el temps d’espera fins la propera erupció. Les erupcions més freqüents són les de més durada.
Enunciat 1.26
Es disposa de les notes d’un examen de filosofia de 40 estudiants d’un Institut de Secundària:
- Trobis la seva distribució per freqüències i calculis la variància i desviació típiques.
- Representis les dades amb un diagrama de barres i en un diagrama de sectors.
Solució
Amb R:

La desviació típica de les notes és 2.22.
S’observen, al diagrama de caixa, els bigotis i els intervals de confiança per a la mediana.
Enunciat 1.27
L’estació meteorològica de l’Institut Europa ha recollit les dades de la temperatura màxima corresponents al temporal meteorològic Filomena durant la primera setmana de gener de 2021. —en °C—:

Representis les dades normalitzades en un histograma.
Solució
Amb R:

S’agrupa la variable \(T_max\) en intervals de 2 °C i es compta el nombre de dies dins cada interval —les freqüències absolutes— i es dibuixa l’histograma i la normalització de les dades.
Es poden calcular també les freqüències relatives —el quocient entre la freqüència absoluta i el nombre total de dades—.
Enunciat 1.28
En un estudi sobre paràsits en 44 tortugues marines es recullen dades sobre la variable aleatòria \(X = \{ \text{nombre de paràsits/individu} \}\)
\(X = \left\{ \begin{aligned} &0, 2, 0, 0, 2, 2, 0, 0, 1, 1, 3, 0, 0, 1, 0, 0, 1, 0, 1, 4, 0, 0, 0, 1, 4, 2, 0, 0, 1, 0, \\ &0, 2, 2, 1, 1, 0, 0, 0, 5, 1, 3, 0, 1, 0, 1, 0, 2, 2, 0, 1, 1, 3, 0, 0, 1, 0, 0, 0, \\ &1, 4, 0, 1, 4, 2, 0, 2, 0, 1, 0, 0, 1, 1, 3, 0, 0, 5 \end{aligned} \right\}\)
Representis les dades, en freqüències i gràficament.
Solució
La mida de la mostra és \(n = 44\). El vector d’observacions és \(X = (x_1, x_2, \cdots ,x_44)\). La variable \(k = 6\) pren els diferents valors \(x = (0, 1, 2, 3, 4, 5)\).
Amb R:

Enunciat 1.29
Es disposa de les edats de 27 mestres d’una escola de primària.
E = {23, 24, 24, 25, 25, 32, 37, 38, 38, 39, 40, 41, 42, 43, 44, 45, 47, 48, 48, 48, 49, 49, 51, 55, 57, 59, 63}
Agrupis les dades per intervals.
Solució
Amb R:

Si es disposa de les dades agrupades es pot calcular la mitjana com: \[\overline{x}_{27} = \frac{(25 \cdot 5) + (35 \cdot 5) + (45 \cdot 12) + (55 \cdot 4) + (65 \cdot 1)}{27} \approx 41.667\]
La desviació típica és: \(s_{27} \approx 10.541\).
La classe de la mediana és \((40,50)\) ja que \(F_3 = 0.815\ (18.52 \text{ %})\ \text{ i } F_2 = 0.370\ (37.04 \text{ %})\), per tant \[M = 40 + \frac{50 - 40}{12} \left( \frac{27}{2} - 10 \right) \approx 42.91\]
i el primer quartil resulta \[Q_1 = 30 + \frac{30 - 20}{5} \left( \frac{27}{24} - 5 \right) \approx 33.5\]
Enunciat 1.30
Es volen generar 75 dades normalitzades, situant la seva distribució sobre un diagrama de caixa i assenyalant la mitjana de la distribució amb una línia discontínua.
Solució
Amb R:

Enunciat 1.31
Es disposa dels resultats del % de ciment en unes perforacions geològiques.
C = {10, 20, 31, 37, 14, 13, 21, 13, 18, 16, 49, 32, 12, 24, 17, 32, 25, 27, 12, 36, 16, 13, 19, 25 }
Representis \(C\) en un diagrama de tija i fulles.
Solució
Amb R:

Enunciat 1.32
Es disposa de les següents magnituds d’un terratrèmol en l’escala de Richter.
\(R = \left\{ \begin{aligned} &1.0, 8.3, 3.1, 1.1, 5.1, 1.2, 1.0, 4.1, 1.1, 4.0, \\ &2.0, 1.9, 6.3, 1.4, 1.3, 3.3, 2.2, 2.3, 2.1, 2.1, \\ &1.4, 2.7, 2.4, 3.0, 4.1, 5.0, 2.2, 1.2, 7.7, 1.5 \end{aligned} \right\}\)
Representis \(R\) en un diagrama de tija i fulles.
Solució
Amb R:

El primer valor \((1.0)\) ens indica la magnitud mínima i l’últim \((8.3)\) la magnitud màxima. A l’esquerra de la barra hi tenim les parts enteres de la variable i a la dreta els valors del primer decimal.
Per exemple, hi ha tres terratrèmols de magnituds 4 (4.0, 4.1, 4.1).
S’observa asimetria positiva —la freqüència disminueix a mesura que augmenta la magnitud—.
Enunciat 1.33
Es disposa dels resultats acadèmics recollits en la taula següent.

Mostris els resultats en un diagrama de coordenades polars.
Solució
Amb R:

Observis que el gràfic de barres horitzontals és més senzill i més clar.
Enunciat 1.34
Generis una simulació de dades normalitzades i representades en un gràfic de coordenades polars
[Adaptació de l’exemple, Ross Ihaka, Department of Statistics, University of Auckland. Private Bag 92019, Auckland New Zealand.]
Solució
Amb R:

Observis que el gràfic ens permet analitzar la dispersió i identificar-hi els diferents clusters.
Enunciat 1.35
Es disposa de les dades de la temperatura mitjana màxima i mínima (en °C), de la precipitació total acumulada (en mm), de la velocitat mitjana del vent (en km/h), dels índex HDD i CDD [34] dels anys 2011 al 2020, per tal de comparar-les entre sí.
Les dades han estat recollides per l’estació meteorològica de l’Institut Europa de L’Hospitalet de Llobregat.
[Observació: Heating Degree Day (HDD) i Cooling Degree Day (CDD) són índexs quantitatius que indiquen la demanda d’energia per escalfar o refredar edificis. La base per determinar si un dia específic és un dia de grau de calefacció (HDD) o un dia de grau de refrigeració (CDD) és la temperatura mitjana diària: T ̅=[(T_max+T_min)/2]. Cada dia amb una temperatura mitjana inferior a 18 °C es compta com un dia de calefacció. El dia del grau de refrigeració és la contrapart estival del dia del grau de calefacció.

Solució
Amb R:

S’observa que les diferències es donen principalment en la precipitació total acumulada i en la velocitat del vent. L’any més regular ha estat l’any 2020.
Enunciat 1.36
Es disposa dels valors de la pluviometria (en mm) registrada a l’estació meteorològica de l’Institut Europa de l’Hospitalet de Llobregat, entre els anys 2008 i 2020.
any = 2008, 2009, 2010, 2011, 2012, 2013, 2014, 2015, 2016, 2017, 2018, 2019, 2020
mm = 707.8, 614.2, 747.5, 875.4, 415.7, 432.0, 519.1, 329.4, 450.5, 435.4, 86.1, 503.2, 899.1
Es desitja saber quina és la tendència de la pluviometria local.
Solució
Amb R:

Enunciat 1.37
Es volen comparar les mesures del pètal, comparant les espècies Iris setosa, Iris versicolor i Iris virginica del dataset “iris” de R, amb un gràfic de dispersió.
Solució
Amb R:

Enunciat 1.38
Es disposa de les dades de la direcció predominant del vent, durant 4555 dies, registrada a l’estació meteorològica de l’Institut Europa de l’Hospitalet de Llobregat, entre els anys 2008 i 2020.
Es desitja saber, utilitzant un diagrama de Pareto, quines són les direccions locals predominants del vent.
Solució
Amb R:

S’observa que els vents predominants són del sud-oest (entre WSW i SSW).
Jordi Coll Vera - 2022