library(dslabs)
library(dplyr)
## 
## Adjuntando el paquete: 'dplyr'
## The following objects are masked from 'package:stats':
## 
##     filter, lag
## The following objects are masked from 'package:base':
## 
##     intersect, setdiff, setequal, union
data("gapminder") 
library(RColorBrewer)
tabla_continente <- table(gapminder$continent)


texto <- "La base de datos de Gapminder incluye información sobre una variedad de indicadores sociales, economicos y de salud, como: Esperanza de vida ,Ingreso per cápita, Tasa de fertilidad, Mortalidad infantil, taza de natalidad 
conprendida en los años 1960 a 2016"

print(texto)
## [1] "La base de datos de Gapminder incluye información sobre una variedad de indicadores sociales, economicos y de salud, como: Esperanza de vida ,Ingreso per cápita, Tasa de fertilidad, Mortalidad infantil, taza de natalidad \nconprendida en los años 1960 a 2016"
head(gapminder, 7)
##               country year infant_mortality life_expectancy fertility
## 1             Albania 1960           115.40           62.87      6.19
## 2             Algeria 1960           148.20           47.50      7.65
## 3              Angola 1960           208.00           35.98      7.32
## 4 Antigua and Barbuda 1960               NA           62.97      4.43
## 5           Argentina 1960            59.87           65.39      3.11
## 6             Armenia 1960               NA           66.86      4.55
## 7               Aruba 1960               NA           65.66      4.82
##   population          gdp continent          region
## 1    1636054           NA    Europe Southern Europe
## 2   11124892  13828152297    Africa Northern Africa
## 3    5270844           NA    Africa   Middle Africa
## 4      54681           NA  Americas       Caribbean
## 5   20619075 108322326649  Americas   South America
## 6    1867396           NA      Asia    Western Asia
## 7      54208           NA  Americas       Caribbean
texto <- "country: El nombre del país.
continent: El continente al que pertenece el país.
year: El año del registro de datos.
lifeExp: Esperanza de vida  (en años).
gdp: Valor estimado del Producto Interno Bruto "
print(texto)
## [1] "country: El nombre del país.\ncontinent: El continente al que pertenece el país.\nyear: El año del registro de datos.\nlifeExp: Esperanza de vida  (en años).\ngdp: Valor estimado del Producto Interno Bruto "
barplot(tabla_continente,
        main = "Gráfico de barras",
        ylab = "Frecuencia",
        xlab = "continentes",
        ylim = c(0, 3000),
        col = c("blue", "green", "yellow"))


grid(nx = NA, ny = NULL, col = "lightgray", lty = "dotted", lwd = par("lwd"))

texto <- ("Este gráfico de barras muestra el número de países presentes en cada continente según el conjunto de dato")
print(texto)
## [1] "Este gráfico de barras muestra el número de países presentes en cada continente según el conjunto de dato"
barplot(table(gapminder$life_expectancy),
        main = "Gráfico de barras",
        ylab = "Frecuencia",
        xlab = "Esperanza de vida",
        ylim = c(0, 50),
       )


grid(nx = NA, ny = NULL, col = "lightgray", lty = "dotted", lwd = par("lwd"))

texto <- ("la esperanza de vida que mas se ha obtenido a travez der los años es alrededor de  74 años")
print(texto)
## [1] "la esperanza de vida que mas se ha obtenido a travez der los años es alrededor de  74 años"
hist(gapminder$fertility,
     main = "Histograma de Fertilidad",
      col = c("blue", "green", "yellow","green","brown","orange"),
     breaks = "sturges",
     probability = TRUE,
     xlab = "Tasa de fertilidad (hijos por mujer)"
)

texto <- "El histograma indica que la mayor parte de los países tienen tasas de fertilidad entre 2 y 3 hijos por mujer."
print(texto)
## [1] "El histograma indica que la mayor parte de los países tienen tasas de fertilidad entre 2 y 3 hijos por mujer."
gapminder_2000 <- gapminder %>%
  filter(year == 2000)

continente_2000 <-factor(gapminder_2000$continent,levels=c("Africa","Asia","Europe","Americas","Oceania"))
tabla1<-table(continente_2000)
barplot(tabla1,
        ylim = c(0, 60)
        )

texto <- ("orden desente de los paise con mas datos obtenidos en el año 2000")
print(texto)
## [1] "orden desente de los paise con mas datos obtenidos en el año 2000"
hist(table(gapminder_2000$infant_mortality),
     main = "Histograma de la tasa de mortalidad infantil en 2000",
     xlim = c(0, 5),
     ylim = c(0, 150),
     col = brewer.pal(3, "Paired"),
     xlab = "Tasa de mortalidad infantil ",
     labels = TRUE,
     border = brewer.pal(3, "Paired"))

hist(table(gapminder_2000$infant_mortality), 
     breaks = "Sturges",
     main = "Histograma de Mortalidad Infantil en el 2000",
     xlab = "Mortalidad infantil en el 2000",
     col = brewer.pal(9, "Set1"),
     ylim = c(0, 200),
     xlim = c(0, 5),  
     border = brewer.pal(9, "Set1"),
     labels = TRUE)

boxplot(gapminder_2000$life_expectancy,
        col="orange",ylab="eperanza de vida",
        outline=FALSE,main="Boxplot",ylim=c(0,100))
points(mean(gapminder_2000$life_expectancy),col="black",pch=20)
text(paste(" ", round(mean(gapminder_2000$life_expectancy), 2)),x=0.9,y=74) 

summary(gapminder_2000$life_expectancy)
##    Min. 1st Qu.  Median    Mean 3rd Qu.    Max. 
##   37.60   61.70   70.80   67.79   74.70   81.00
# Diagrama de dispersión PIB vs Esperanza de Vida

y <- gapminder_2000$life_expectancy
x <- gapminder_2000$gdp


plot(x, y, 
     main = "Esperanza de Vida vs PIB  ",
      xlim = c(0, 5000000000),#cincomilmillones
     col = rainbow(length(y)), 
     pch = 20, 
     xlab = "PIB (Producto Interno Bruto)", 
     ylab = "Esperanza de Vida")



legend("bottomright", 
       legend = paste("pib"),
       col = rainbow(length(y)), 
      
     )

texto <- ("muestra la esperanza de vida de acurdo al pib deacuerdo con la grafica los paise con pib de cinco mil miilones pueden tener la misma esperanza con paise de cien mil millones de us ")
print(texto)
## [1] "muestra la esperanza de vida de acurdo al pib deacuerdo con la grafica los paise con pib de cinco mil miilones pueden tener la misma esperanza con paise de cien mil millones de us "
boxplot(life_expectancy ~ continent, data = gapminder_2000, 
        col = brewer.pal(9, "Set1"), 
        outline = FALSE, 
        ylim = c(0, 100), 
        main = "Boxplot: Esperanza de vida vs Región")


abline(h = mean(gapminder_2000$life_expectancy), col = "blue", lwd = 2)


text(x = 1, y = mean(gapminder_2000$gdp), 
     labels = paste("Media PIB: ", round(mean(gapminder_2000$life_expectancy), 2)), 
     col = "blue", pos = 3)

par(mfrow=c(1,2))


# Boxplot PIB vs Región
boxplot(life_expectancy ~ continent, data = gapminder_2000, 
        col = brewer.pal(9, "Set1"), 
        outline = FALSE, 
        ylim = c(0, 100), 
        main = "Boxplot: Esperanza de vida vs Región")


abline(h = mean(gapminder_2000$life_expectancy), col = "blue", lwd = 2)


text(x = 1, y = mean(gapminder_2000$gdp), 
     labels = paste("Media PIB: ", round(mean(gapminder_2000$life_expectancy), 2)), 
     col = "blue", pos = 3)


y <- gapminder_2000$life_expectancy
x <- gapminder_2000$gdp


plot(x, y, 
     main = "PIB vs Esperanza de Vida",
      xlim = c(0, 5000000000),#cincomilmillones
     col = rainbow(length(y)), 
     pch = 20, 
     xlab = "PIB (Producto Interno Bruto)", 
     ylab = "Esperanza de Vida")



legend("bottomright", 
       legend = paste("pib"),
       col = rainbow(length(y)), 
      
     )