setwd("~/UCE/III/estadistica")
#Carga de Datos
datos <- read.csv("Conjunto_Datos_Minerales.csv", header = TRUE, sep = ";",
dec = ".")
PrimerAnio <- na.omit(datos$yr_fst_prd)
# Extraer variable
PrimerAnio <- datos$yr_fst_prd
# Tabla de frecuencia
TDFprimerAnio <- table(PrimerAnio)
TablaPrimerAnio <- as.data.frame(TDFprimerAnio)
# Porcentaje (hi)
hi <- (TablaPrimerAnio$Freq / sum(TablaPrimerAnio$Freq)) * 100
TablaPrimerAnio <- data.frame(TablaPrimerAnio, round(hi, 2))
# Frecuencias acumuladas
Niasc <- cumsum(TablaPrimerAnio$Freq)
Hiasc <- cumsum(TablaPrimerAnio$round.hi)
Nidsc <- rev(cumsum(rev(TablaPrimerAnio$Freq)))
Hidsc <- rev(cumsum(rev(TablaPrimerAnio$round.hi)))
# Tabla final con nombres ni y hi
TablaPrimerAnioFinal <- data.frame(
PrimerAnio = TablaPrimerAnio$PrimerAnio,
ni = TablaPrimerAnio$Freq,
hi = round(TablaPrimerAnio$round.hi, 2),
Niasc, Hiasc,
Nidsc, Hidsc
)
#Visualizar la tabla
print(TablaPrimerAnioFinal)
## PrimerAnio ni hi Niasc Hiasc Nidsc Hidsc
## 1 1 2 0.01 2 0.01 16024 100.11
## 2 225 1 0.01 3 0.02 16022 100.10
## 3 770 1 0.01 4 0.03 16021 100.09
## 4 950 1 0.01 5 0.04 16020 100.08
## 5 968 1 0.01 6 0.05 16019 100.07
## 6 1000 2 0.01 8 0.06 16018 100.06
## 7 1036 1 0.01 9 0.07 16016 100.05
## 8 1050 1 0.01 10 0.08 16015 100.04
## 9 1100 1 0.01 11 0.09 16014 100.03
## 10 1250 1 0.01 12 0.10 16013 100.02
## 11 1311 1 0.01 13 0.11 16012 100.01
## 12 1318 1 0.01 14 0.12 16011 100.00
## 13 1320 1 0.01 15 0.13 16010 99.99
## 14 1400 1 0.01 16 0.14 16009 99.98
## 15 1446 1 0.01 17 0.15 16008 99.97
## 16 1492 1 0.01 18 0.16 16007 99.96
## 17 1500 4 0.02 22 0.18 16006 99.95
## 18 1505 1 0.01 23 0.19 16002 99.93
## 19 1542 1 0.01 24 0.20 16001 99.92
## 20 1545 1 0.01 25 0.21 16000 99.91
## 21 1546 1 0.01 26 0.22 15999 99.90
## 22 1549 1 0.01 27 0.23 15998 99.89
## 23 1550 4 0.02 31 0.25 15997 99.88
## 24 1554 1 0.01 32 0.26 15993 99.86
## 25 1555 1 0.01 33 0.27 15992 99.85
## 26 1560 1 0.01 34 0.28 15991 99.84
## 27 1568 1 0.01 35 0.29 15990 99.83
## 28 1572 1 0.01 36 0.30 15989 99.82
## 29 1575 2 0.01 38 0.31 15988 99.81
## 30 1591 1 0.01 39 0.32 15986 99.80
## 31 1595 2 0.01 41 0.33 15985 99.79
## 32 1600 4 0.02 45 0.35 15983 99.78
## 33 1607 1 0.01 46 0.36 15979 99.76
## 34 1623 2 0.01 48 0.37 15978 99.75
## 35 1625 1 0.01 49 0.38 15976 99.74
## 36 1630 1 0.01 50 0.39 15975 99.73
## 37 1634 1 0.01 51 0.40 15974 99.72
## 38 1640 3 0.02 54 0.42 15973 99.71
## 39 1650 5 0.03 59 0.45 15970 99.69
## 40 1665 1 0.01 60 0.46 15965 99.66
## 41 1672 1 0.01 61 0.47 15964 99.65
## 42 1675 2 0.01 63 0.48 15963 99.64
## 43 1677 4 0.02 67 0.50 15961 99.63
## 44 1680 2 0.01 69 0.51 15957 99.61
## 45 1687 1 0.01 70 0.52 15955 99.60
## 46 1692 1 0.01 71 0.53 15954 99.59
## 47 1700 3 0.02 74 0.55 15953 99.58
## 48 1701 1 0.01 75 0.56 15950 99.56
## 49 1707 3 0.02 78 0.58 15949 99.55
## 50 1710 8 0.05 86 0.63 15946 99.53
## 51 1713 2 0.01 88 0.64 15938 99.48
## 52 1715 1 0.01 89 0.65 15936 99.47
## 53 1720 3 0.02 92 0.67 15935 99.46
## 54 1721 1 0.01 93 0.68 15932 99.44
## 55 1722 2 0.01 95 0.69 15931 99.43
## 56 1727 3 0.02 98 0.71 15929 99.42
## 57 1733 2 0.01 100 0.72 15926 99.40
## 58 1734 2 0.01 102 0.73 15924 99.39
## 59 1737 1 0.01 103 0.74 15922 99.38
## 60 1742 1 0.01 104 0.75 15921 99.37
## 61 1746 1 0.01 105 0.76 15920 99.36
## 62 1749 1 0.01 106 0.77 15919 99.35
## 63 1750 5 0.03 111 0.80 15918 99.34
## 64 1751 1 0.01 112 0.81 15913 99.31
## 65 1753 2 0.01 114 0.82 15912 99.30
## 66 1755 3 0.02 117 0.84 15910 99.29
## 67 1756 2 0.01 119 0.85 15907 99.27
## 68 1757 1 0.01 120 0.86 15905 99.26
## 69 1758 1 0.01 121 0.87 15904 99.25
## 70 1760 4 0.02 125 0.89 15903 99.24
## 71 1764 1 0.01 126 0.90 15899 99.22
## 72 1765 3 0.02 129 0.92 15898 99.21
## 73 1767 1 0.01 130 0.93 15895 99.19
## 74 1770 2 0.01 132 0.94 15894 99.18
## 75 1771 2 0.01 134 0.95 15892 99.17
## 76 1772 2 0.01 136 0.96 15890 99.16
## 77 1773 1 0.01 137 0.97 15888 99.15
## 78 1774 2 0.01 139 0.98 15887 99.14
## 79 1775 11 0.07 150 1.05 15885 99.13
## 80 1776 8 0.05 158 1.10 15874 99.06
## 81 1777 3 0.02 161 1.12 15866 99.01
## 82 1778 3 0.02 164 1.14 15863 98.99
## 83 1779 1 0.01 165 1.15 15860 98.97
## 84 1780 2 0.01 167 1.16 15859 98.96
## 85 1784 1 0.01 168 1.17 15857 98.95
## 86 1785 4 0.02 172 1.19 15856 98.94
## 87 1786 2 0.01 174 1.20 15852 98.92
## 88 1787 1 0.01 175 1.21 15850 98.91
## 89 1789 4 0.02 179 1.23 15849 98.90
## 90 1790 3 0.02 182 1.25 15845 98.88
## 91 1791 1 0.01 183 1.26 15842 98.86
## 92 1792 2 0.01 185 1.27 15841 98.85
## 93 1795 3 0.02 188 1.29 15839 98.84
## 94 1797 1 0.01 189 1.30 15836 98.82
## 95 1800 37 0.23 226 1.53 15835 98.81
## 96 1802 1 0.01 227 1.54 15798 98.58
## 97 1804 4 0.02 231 1.56 15797 98.57
## 98 1805 2 0.01 233 1.57 15793 98.55
## 99 1807 1 0.01 234 1.58 15791 98.54
## 100 1808 1 0.01 235 1.59 15790 98.53
## 101 1809 2 0.01 237 1.60 15789 98.52
## 102 1810 6 0.04 243 1.64 15787 98.51
## 103 1812 3 0.02 246 1.66 15781 98.47
## 104 1813 2 0.01 248 1.67 15778 98.45
## 105 1814 1 0.01 249 1.68 15776 98.44
## 106 1815 5 0.03 254 1.71 15775 98.43
## 107 1816 1 0.01 255 1.72 15770 98.40
## 108 1817 2 0.01 257 1.73 15769 98.39
## 109 1818 4 0.02 261 1.75 15767 98.38
## 110 1820 3 0.02 264 1.77 15763 98.36
## 111 1821 4 0.02 268 1.79 15760 98.34
## 112 1823 2 0.01 270 1.80 15756 98.32
## 113 1824 3 0.02 273 1.82 15754 98.31
## 114 1825 53 0.33 326 2.15 15751 98.29
## 115 1826 2 0.01 328 2.16 15698 97.96
## 116 1827 4 0.02 332 2.18 15696 97.95
## 117 1828 27 0.17 359 2.35 15692 97.93
## 118 1829 5 0.03 364 2.38 15665 97.76
## 119 1830 19 0.12 383 2.50 15660 97.73
## 120 1831 15 0.09 398 2.59 15641 97.61
## 121 1832 6 0.04 404 2.63 15626 97.52
## 122 1833 5 0.03 409 2.66 15620 97.48
## 123 1834 9 0.06 418 2.72 15615 97.45
## 124 1835 32 0.20 450 2.92 15606 97.39
## 125 1836 14 0.09 464 3.01 15574 97.19
## 126 1837 7 0.04 471 3.05 15560 97.10
## 127 1838 8 0.05 479 3.10 15553 97.06
## 128 1839 6 0.04 485 3.14 15545 97.01
## 129 1840 19 0.12 504 3.26 15539 96.97
## 130 1841 6 0.04 510 3.30 15520 96.85
## 131 1842 10 0.06 520 3.36 15514 96.81
## 132 1843 4 0.02 524 3.38 15504 96.75
## 133 1844 11 0.07 535 3.45 15500 96.73
## 134 1845 29 0.18 564 3.63 15489 96.66
## 135 1846 13 0.08 577 3.71 15460 96.48
## 136 1847 13 0.08 590 3.79 15447 96.40
## 137 1848 22 0.14 612 3.93 15434 96.32
## 138 1849 18 0.11 630 4.04 15412 96.18
## 139 1850 143 0.89 773 4.93 15394 96.07
## 140 1851 15 0.09 788 5.02 15251 95.18
## 141 1852 25 0.16 813 5.18 15236 95.09
## 142 1853 150 0.94 963 6.12 15211 94.93
## 143 1854 38 0.24 1001 6.36 15061 93.99
## 144 1855 60 0.37 1061 6.73 15023 93.75
## 145 1856 26 0.16 1087 6.89 14963 93.38
## 146 1857 19 0.12 1106 7.01 14937 93.22
## 147 1858 17 0.11 1123 7.12 14918 93.10
## 148 1859 16 0.10 1139 7.22 14901 92.99
## 149 1860 62 0.39 1201 7.61 14885 92.89
## 150 1861 51 0.32 1252 7.93 14823 92.50
## 151 1862 27 0.17 1279 8.10 14772 92.18
## 152 1863 28 0.17 1307 8.27 14745 92.01
## 153 1864 65 0.41 1372 8.68 14717 91.84
## 154 1865 92 0.57 1464 9.25 14652 91.43
## 155 1866 45 0.28 1509 9.53 14560 90.86
## 156 1867 39 0.24 1548 9.77 14515 90.58
## 157 1868 91 0.57 1639 10.34 14476 90.34
## 158 1869 44 0.27 1683 10.61 14385 89.77
## 159 1870 135 0.84 1818 11.45 14341 89.50
## 160 1871 63 0.39 1881 11.84 14206 88.66
## 161 1872 101 0.63 1982 12.47 14143 88.27
## 162 1873 109 0.68 2091 13.15 14042 87.64
## 163 1874 42 0.26 2133 13.41 13933 86.96
## 164 1875 224 1.40 2357 14.81 13891 86.70
## 165 1876 34 0.21 2391 15.02 13667 85.30
## 166 1877 52 0.32 2443 15.34 13633 85.09
## 167 1878 44 0.27 2487 15.61 13581 84.77
## 168 1879 84 0.52 2571 16.13 13537 84.50
## 169 1880 234 1.46 2805 17.59 13453 83.98
## 170 1881 65 0.41 2870 18.00 13219 82.52
## 171 1882 151 0.94 3021 18.94 13154 82.11
## 172 1883 63 0.39 3084 19.33 13003 81.17
## 173 1884 69 0.43 3153 19.76 12940 80.78
## 174 1885 272 1.70 3425 21.46 12871 80.35
## 175 1886 58 0.36 3483 21.82 12599 78.65
## 176 1887 87 0.54 3570 22.36 12541 78.29
## 177 1888 70 0.44 3640 22.80 12454 77.75
## 178 1889 48 0.30 3688 23.10 12384 77.31
## 179 1890 176 1.10 3864 24.20 12336 77.01
## 180 1891 56 0.35 3920 24.55 12160 75.91
## 181 1892 87 0.54 4007 25.09 12104 75.56
## 182 1893 70 0.44 4077 25.53 12017 75.02
## 183 1894 52 0.32 4129 25.85 11947 74.58
## 184 1895 147 0.92 4276 26.77 11895 74.26
## 185 1896 82 0.51 4358 27.28 11748 73.34
## 186 1897 82 0.51 4440 27.79 11666 72.83
## 187 1898 79 0.49 4519 28.28 11584 72.32
## 188 1899 82 0.51 4601 28.79 11505 71.83
## 189 1900 293 1.83 4894 30.62 11423 71.32
## 190 1901 110 0.69 5004 31.31 11130 69.49
## 191 1902 157 0.98 5161 32.29 11020 68.80
## 192 1903 156 0.97 5317 33.26 10863 67.82
## 193 1904 142 0.89 5459 34.15 10707 66.85
## 194 1905 162 1.01 5621 35.16 10565 65.96
## 195 1906 191 1.19 5812 36.35 10403 64.95
## 196 1907 197 1.23 6009 37.58 10212 63.76
## 197 1908 244 1.52 6253 39.10 10015 62.53
## 198 1909 118 0.74 6371 39.84 9771 61.01
## 199 1910 235 1.47 6606 41.31 9653 60.27
## 200 1911 131 0.82 6737 42.13 9418 58.80
## 201 1912 116 0.72 6853 42.85 9287 57.98
## 202 1913 127 0.79 6980 43.64 9171 57.26
## 203 1914 154 0.96 7134 44.60 9044 56.47
## 204 1915 265 1.65 7399 46.25 8890 55.51
## 205 1916 277 1.73 7676 47.98 8625 53.86
## 206 1917 284 1.77 7960 49.75 8348 52.13
## 207 1918 227 1.42 8187 51.17 8064 50.36
## 208 1919 123 0.77 8310 51.94 7837 48.94
## 209 1920 114 0.71 8424 52.65 7714 48.17
## 210 1921 106 0.66 8530 53.31 7600 47.46
## 211 1922 94 0.59 8624 53.90 7494 46.80
## 212 1923 76 0.47 8700 54.37 7400 46.21
## 213 1924 90 0.56 8790 54.93 7324 45.74
## 214 1925 179 1.12 8969 56.05 7234 45.18
## 215 1926 92 0.57 9061 56.62 7055 44.06
## 216 1927 89 0.56 9150 57.18 6963 43.49
## 217 1928 81 0.51 9231 57.69 6874 42.93
## 218 1929 90 0.56 9321 58.25 6793 42.42
## 219 1930 121 0.76 9442 59.01 6703 41.86
## 220 1931 70 0.44 9512 59.45 6582 41.10
## 221 1932 113 0.71 9625 60.16 6512 40.66
## 222 1933 75 0.47 9700 60.63 6399 39.95
## 223 1934 102 0.64 9802 61.27 6324 39.48
## 224 1935 134 0.84 9936 62.11 6222 38.84
## 225 1936 96 0.60 10032 62.71 6088 38.00
## 226 1937 141 0.88 10173 63.59 5992 37.40
## 227 1938 109 0.68 10282 64.27 5851 36.52
## 228 1939 119 0.74 10401 65.01 5742 35.84
## 229 1940 157 0.98 10558 65.99 5623 35.10
## 230 1941 186 1.16 10744 67.15 5466 34.12
## 231 1942 259 1.62 11003 68.77 5280 32.96
## 232 1943 154 0.96 11157 69.73 5021 31.34
## 233 1944 80 0.50 11237 70.23 4867 30.38
## 234 1945 97 0.61 11334 70.84 4787 29.88
## 235 1946 58 0.36 11392 71.20 4690 29.27
## 236 1947 61 0.38 11453 71.58 4632 28.91
## 237 1948 90 0.56 11543 72.14 4571 28.53
## 238 1949 81 0.51 11624 72.65 4481 27.97
## 239 1950 174 1.09 11798 73.74 4400 27.46
## 240 1951 119 0.74 11917 74.48 4226 26.37
## 241 1952 127 0.79 12044 75.27 4107 25.63
## 242 1953 142 0.89 12186 76.16 3980 24.84
## 243 1954 143 0.89 12329 77.05 3838 23.95
## 244 1955 201 1.25 12530 78.30 3695 23.06
## 245 1956 127 0.79 12657 79.09 3494 21.81
## 246 1957 126 0.79 12783 79.88 3367 21.02
## 247 1958 97 0.61 12880 80.49 3241 20.23
## 248 1959 99 0.62 12979 81.11 3144 19.62
## 249 1960 107 0.67 13086 81.78 3045 19.00
## 250 1961 90 0.56 13176 82.34 2938 18.33
## 251 1962 88 0.55 13264 82.89 2848 17.77
## 252 1963 83 0.52 13347 83.41 2760 17.22
## 253 1964 98 0.61 13445 84.02 2677 16.70
## 254 1965 140 0.87 13585 84.89 2579 16.09
## 255 1966 95 0.59 13680 85.48 2439 15.22
## 256 1967 123 0.77 13803 86.25 2344 14.63
## 257 1968 131 0.82 13934 87.07 2221 13.86
## 258 1969 119 0.74 14053 87.81 2090 13.04
## 259 1970 107 0.67 14160 88.48 1971 12.30
## 260 1971 112 0.70 14272 89.18 1864 11.63
## 261 1972 99 0.62 14371 89.80 1752 10.93
## 262 1973 84 0.52 14455 90.32 1653 10.31
## 263 1974 101 0.63 14556 90.95 1569 9.79
## 264 1975 128 0.80 14684 91.75 1468 9.16
## 265 1976 89 0.56 14773 92.31 1340 8.36
## 266 1977 99 0.62 14872 92.93 1251 7.80
## 267 1978 68 0.42 14940 93.35 1152 7.18
## 268 1979 89 0.56 15029 93.91 1084 6.76
## 269 1980 114 0.71 15143 94.62 995 6.20
## 270 1981 77 0.48 15220 95.10 881 5.49
## 271 1982 84 0.52 15304 95.62 804 5.01
## 272 1983 66 0.41 15370 96.03 720 4.49
## 273 1984 75 0.47 15445 96.50 654 4.08
## 274 1985 63 0.39 15508 96.89 579 3.61
## 275 1986 79 0.49 15587 97.38 516 3.22
## 276 1987 66 0.41 15653 97.79 437 2.73
## 277 1988 68 0.42 15721 98.21 371 2.32
## 278 1989 81 0.51 15802 98.72 303 1.90
## 279 1990 52 0.32 15854 99.04 222 1.39
## 280 1991 42 0.26 15896 99.30 170 1.07
## 281 1992 32 0.20 15928 99.50 128 0.81
## 282 1993 19 0.12 15947 99.62 96 0.61
## 283 1994 21 0.13 15968 99.75 77 0.49
## 284 1995 15 0.09 15983 99.84 56 0.36
## 285 1996 29 0.18 16012 100.02 41 0.27
## 286 1997 8 0.05 16020 100.07 12 0.09
## 287 2003 1 0.01 16021 100.08 4 0.04
## 288 2005 1 0.01 16022 100.09 3 0.03
## 289 2006 1 0.01 16023 100.10 2 0.02
## 290 2007 1 0.01 16024 100.11 1 0.01
# Gráfico de barras: FRECUENCIA LOCAL
barplot(TablaPrimerAnioFinal$ni,
main = "Gráfico No.1: Distribución de cantidad del primer año de producción",
xlab = "Año", ylab = "cantidad",
col = "red",
names.arg = TablaPrimerAnioFinal$PrimerAnio,
las=1,,
cex.names = 0.5
)
# Gráfico de barras: FRECUENCIA GLOBAL
barplot(TablaPrimerAnioFinal$ni,
main = "Gráfico No.2: Distribución de cantidad del primer año de producción",
xlab = "Año", ylab = "cantidad",
col = "red",
names.arg = TablaPrimerAnioFinal$PrimerAnio,
ylim = c(0, length(PrimerAnio)),
las=1,
cex.names = 0.5,
cex.axis = 0.5)
# Gráfico de barras: PORCENTAJE LOCAL
barplot(TablaPrimerAnioFinal$hi,
main = "Gráfico No.3: Distribución de porcentaje del ultimo año de producción",
xlab = "Año", ylab = "Porcentaje",
col = "blue",
names.arg = TablaPrimerAnioFinal$PrimerAnio,
las=1,
cex.names = 0.5)
# Gráfico de barras: PORCENTAJE GLOBAL
barplot(TablaPrimerAnioFinal$hi,
main = "Gráfico No.4: Distribución de porcentaje del ultimo año de producción",
xlab = "Año", ylab = "Porcentaje",
col = "blue",
names.arg = TablaPrimerAnioFinal$PrimerAnio,
ylim = c(0, 100),
las=1,
cex.names = 0.5,
cex.axis = 0.5
)
# Boxplot
boxplot(PrimerAnio, horizontal = TRUE, col = "green",
main = "Gráfico No.5: Distribución del Primer Año de Producción")
# Ojiva Ni ascendente
x <- as.numeric(as.character(TablaPrimerAnioFinal$PrimerAnio))
y <- TablaPrimerAnioFinal$Niasc
plot(x, y, type = "p", main = "Gráfico No.6: Ojiva Ascendente - Primer Año de Producción",
xlab = "Año", ylab = "Frecuencia acumulada", col = "blue")
# Ojiva Ni descendente
y <- TablaPrimerAnioFinal$Nidsc
plot(x, y, type = "p", main = "Gráfico No.7: Ojiva Descendente - Primer Año de Producción",
xlab = "Año", ylab = "Frecuencia acumulada", col = "blue")
# Ojiva Hi ascendente
y <- TablaPrimerAnioFinal$Hiasc
plot(x, y, type = "p", main = "Gráfico No.8: Ojiva Ascendente - Primer Año de Producción",
xlab = "Año", ylab = "Porcentaje acumulado (%)", col = "blue")
# Ojiva Hi descendente
y <- TablaPrimerAnioFinal$Hidsc
plot(x, y, type = "p", main = "Gráfico No.9: Ojiva Descendente - Primer Año de Producción",
xlab = "Año", ylab = "Porcentaje acumulado (%)", col = "blue")
# Intersección Ni (Ojivas)
plot(x, TablaPrimerAnioFinal$Niasc, type = "p",
main = "Gráfico No 10 : Ojiva Asc. y Desc. (Ni) - Primer Año de Producción",
xlab = "Año", ylab = "Frecuencia acumulada", col = "green")
lines(x, TablaPrimerAnioFinal$Nidsc, type = "p", col = "blue")
# Intersección Hi (Ojivas)
plot(x, TablaPrimerAnioFinal$Hiasc, type = "p",
main = "Gráfico No 11: Ojiva Asc. y Desc. (Hi%) - Primer Año de Producción",
xlab = "Año", ylab = "Porcentaje acumulado", col = "green")
lines(x, TablaPrimerAnioFinal$Hidsc, type = "p", col = "blue")
