Assingment 1-4 Post Module

First we run de package readr and run our database.

library(readr)
internet_readr <- read_csv("API_IT.NET.USER.ZS_DS2_en_csv_v2_213153.csv", skip=4)

## New names:
## Rows: 266 Columns: 69
## ── Column specification
## ──────────────────────────────────────────────────────── Delimiter: "," chr
## (4): Country Name, Country Code, Indicator Name, Indicator Code dbl (51): 1960,
## 1965, 1970, 1975, 1976, 1977, 1978, 1979, 1980, 1981, 1982, ... lgl (14): 1961,
## 1962, 1963, 1964, 1966, 1967, 1968, 1969, 1971, 1972, 1973, ...
## ℹ Use `spec()` to retrieve the full column specification for this data. ℹ
## Specify the column types or set `show_col_types = FALSE` to quiet this message.
## • `` -> `...69`

#Then we clean some columns

library(dplyr)

## 
## Attaching package: 'dplyr'

## The following objects are masked from 'package:stats':
## 
##     filter, lag

## The following objects are masked from 'package:base':
## 
##     intersect, setdiff, setequal, union

internet_readr <- internet_readr %>% select(`Country Code`, `Country Name`, `1992`:`2022`) 

#Then we group the data that is not a country.

non_countries <- read_csv("non_countries.csv")

## Rows: 79 Columns: 2
## ── Column specification ────────────────────────────────────────────────────────
## Delimiter: ","
## chr (2): code, region
## 
## ℹ Use `spec()` to retrieve the full column specification for this data.
## ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.

#And filter out the non countries from the data base

internet_readr <-internet_readr[!internet_readr$`Country Code` %in% non_countries$code, ]

#Q1:let’s graph the histograms from 1992 to 2022.

par(mfrow=c(1,3))
hist(internet_readr$`1992`, breaks=8, 
     main= "Internet usage for 1992", 
     col="blue", xlab=" ")
hist(internet_readr$`1993`, breaks=8, 
     main= "Internet usage for 1993", 
     col="blue", xlab=" ")
hist(internet_readr$`1994`, breaks=8, 
     main= "Internet usage for 1994", 
     col="blue", xlab=" ")

hist(internet_readr$`1995`, breaks=8, 
     main= "Internet usage for 1995", 
     col="blue", xlab=" ")
hist(internet_readr$`1996`, breaks=8, 
     main= "Internet usage for 1996", 
     col="blue", xlab=" ")
hist(internet_readr$`1997`, breaks=8, 
     main= "Internet usage for 1997", 
     col="blue", xlab=" ")

hist(internet_readr$`1998`, breaks=8, 
     main= "Internet usage for 1998", 
     col="blue", xlab=" ")
hist(internet_readr$`1999`, breaks=8, 
     main= "Internet usage for 1999", 
     col="blue", xlab=" ")
hist(internet_readr$`2000`, breaks=8, 
     main= "Internet usage for 2000", 
     col="blue", xlab=" ")

hist(internet_readr$`2001`, breaks=8, 
     main= "Internet usage for 2001", 
     col="blue", xlab=" ")
hist(internet_readr$`2002`, breaks=8, 
     main= "Internet usage for 2002", 
     col="blue", xlab=" ")
hist(internet_readr$`2003`, breaks=8, 
     main= "Internet usage for 2003", 
     col="blue", xlab=" ")

hist(internet_readr$`2004`, breaks=8, 
     main= "Internet usage for 2004", 
     col="blue", xlab=" ")
hist(internet_readr$`2005`, breaks=8, 
     main= "Internet usage for 2005", 
     col="blue", xlab=" ")
hist(internet_readr$`2006`, breaks=8, 
     main= "Internet usage for 2006", 
     col="blue", xlab=" ")

hist(internet_readr$`2007`, breaks=8, 
     main= "Internet usage for 2007", 
     col="blue", xlab=" ")
hist(internet_readr$`2008`, breaks=8, 
     main= "Internet usage for 2008", 
     col="blue", xlab=" ")
hist(internet_readr$`2009`, breaks=8, 
     main= "Internet usage for 2009", 
     col="blue", xlab=" ")

hist(internet_readr$`2010`, breaks=8, 
     main= "Internet usage for 2010", 
     col="blue", xlab=" ")
hist(internet_readr$`2011`, breaks=8, 
     main= "Internet usage for 2011", 
     col="blue", xlab=" ")
hist(internet_readr$`2012`, breaks=8, 
     main= "Internet usage for 2012", 
     col="blue", xlab=" ")

hist(internet_readr$`2013`, breaks=8, 
     main= "Internet usage for 2013", 
     col="blue", xlab=" ")
hist(internet_readr$`2014`, breaks=8, 
     main= "Internet usage for 2014", 
     col="blue", xlab=" ")
hist(internet_readr$`2015`, breaks=8, 
     main= "Internet usage for 2015", 
     col="blue", xlab=" ")

hist(internet_readr$`2016`, breaks=8, 
     main= "Internet usage for 2016", 
     col="blue", xlab=" ")
hist(internet_readr$`2017`, breaks=8, 
     main= "Internet usage for 2017", 
     col="blue", xlab=" ")
hist(internet_readr$`2018`, breaks=8, 
     main= "Internet usage for 2018", 
     col="blue", xlab=" ")

hist(internet_readr$`2019`, breaks=8, 
     main= "Internet usage for 2019", 
     col="blue", xlab=" ")
hist(internet_readr$`2020`, breaks=8, 
     main= "Internet usage for 2020", 
     col="blue", xlab=" ")
hist(internet_readr$`2021`, breaks=8, 
     main= "Internet usage for 2021", 
     col="blue", xlab=" ")

hist(internet_readr$`2022`, breaks=8, 
     main= "Internet usage for 2022", 
     col="blue", xlab=" ")

#Before advancing into Q2, let’s see how many countries have data for each year

total_countries_1992<- length(internet_readr$`1992`) - sum(is.na(internet_readr$`1992`))
total_countries_1993<- length(internet_readr$`1993`) - sum(is.na(internet_readr$`1993`))
total_countries_1994<- length(internet_readr$`1994`) - sum(is.na(internet_readr$`1994`))
total_countries_1995<- length(internet_readr$`1995`) - sum(is.na(internet_readr$`1995`))
total_countries_1996<- length(internet_readr$`1996`) - sum(is.na(internet_readr$`1996`))
total_countries_1997<- length(internet_readr$`1997`) - sum(is.na(internet_readr$`1997`))
total_countries_1998<- length(internet_readr$`1998`) - sum(is.na(internet_readr$`1998`))
total_countries_1999<- length(internet_readr$`1999`) - sum(is.na(internet_readr$`1999`))
total_countries_2000<- length(internet_readr$`2000`) - sum(is.na(internet_readr$`2000`))
total_countries_2001<- length(internet_readr$`2001`) - sum(is.na(internet_readr$`2001`))
total_countries_2002<- length(internet_readr$`2002`) - sum(is.na(internet_readr$`2002`))
total_countries_2003<- length(internet_readr$`2003`) - sum(is.na(internet_readr$`2003`))
total_countries_2004<- length(internet_readr$`2004`) - sum(is.na(internet_readr$`2004`))
total_countries_2005<- length(internet_readr$`2005`) - sum(is.na(internet_readr$`2005`))
total_countries_2006<- length(internet_readr$`2006`) - sum(is.na(internet_readr$`2006`))
total_countries_2007<- length(internet_readr$`2007`) - sum(is.na(internet_readr$`2007`))
total_countries_2008<- length(internet_readr$`2008`) - sum(is.na(internet_readr$`2008`))
total_countries_2009<- length(internet_readr$`2009`) - sum(is.na(internet_readr$`2009`))
total_countries_2010<- length(internet_readr$`2010`) - sum(is.na(internet_readr$`2010`))
total_countries_2011<- length(internet_readr$`2011`) - sum(is.na(internet_readr$`2011`))
total_countries_2012<- length(internet_readr$`2012`) - sum(is.na(internet_readr$`2012`))
total_countries_2012<- length(internet_readr$`2012`) - sum(is.na(internet_readr$`2012`))
total_countries_2013<- length(internet_readr$`2013`) - sum(is.na(internet_readr$`2013`))
total_countries_2014<- length(internet_readr$`2014`) - sum(is.na(internet_readr$`2014`))
total_countries_2015<- length(internet_readr$`2015`) - sum(is.na(internet_readr$`2015`))
total_countries_2016<- length(internet_readr$`2016`) - sum(is.na(internet_readr$`2016`))
total_countries_2017<- length(internet_readr$`2017`) - sum(is.na(internet_readr$`2017`))
total_countries_2018<- length(internet_readr$`2018`) - sum(is.na(internet_readr$`2018`))
total_countries_2019<- length(internet_readr$`2019`) - sum(is.na(internet_readr$`2019`))
total_countries_2020<- length(internet_readr$`2020`) - sum(is.na(internet_readr$`2020`))
total_countries_2021<- length(internet_readr$`2021`) - sum(is.na(internet_readr$`2021`))
total_countries_2022<- length(internet_readr$`2022`) - sum(is.na(internet_readr$`2022`))

#Q2: Now let’s graph the boxplots.

par(mfrow=c(1,2))
boxplot(internet_readr$`1992`, 
     main= "Int. usage for 1992 208 Countries", 
     col="blue", xlab=paste("The median is:",median(internet_readr$`1992`)),horizontal = TRUE,na.rm=TRUE)

boxplot(internet_readr$`1993`, 
     main= "Int. usage for 1993 208 Countries", 
     col="blue", xlab=paste("The median is:",median(internet_readr$`1993`)),horizontal = TRUE,na.rm=TRUE)

boxplot(internet_readr$`1994`, 
     main= "Int. usage for 1994 208 Countries", 
     col="blue", xlab=paste("The median is:",median(internet_readr$`1994`)),horizontal = TRUE,na.rm=TRUE)

boxplot(internet_readr$`1995`, 
     main= "Int. usage for 1995 208 Countries", 
     col="blue", xlab=paste("The median is:",median(internet_readr$`1995`)),horizontal = TRUE,na.rm=TRUE)

boxplot(internet_readr$`1996`, 
     main= "Int. usage for 1996 171 Countries", 
     col="blue", xlab=paste("The median is:",median(internet_readr$`1996`)),horizontal = TRUE,na.rm=TRUE)

boxplot(internet_readr$`1997`,
     main= "Int. usage for 1997 181  Countries", 
     col="blue", xlab=paste("The median is:",median(internet_readr$`1997`)),horizontal = TRUE,na.rm=TRUE)

boxplot(internet_readr$`1998`, 
     main= "Int. usage for 1998 187 Countries", 
     col="blue", xlab=paste("The median is:",median(internet_readr$`1998`)),horizontal = TRUE,na.rm=TRUE)

boxplot(internet_readr$`1999`,
     main= "Int. usage for 1999 194 Countries", 
     col="blue", xlab=paste("The median is:",median(internet_readr$`1999`)),horizontal = TRUE,na.rm=TRUE)

boxplot(internet_readr$`2000`, 
     main= "Int. usage for 2000 196 Countries", 
     col="blue", xlab=paste("The median is:",median(internet_readr$`2000`)),horizontal = TRUE,na.rm=TRUE)

boxplot(internet_readr$`2001`, 
     main= "Int. usage for 2001 199 Countries", 
     col="blue", xlab=paste("The median is:",median(internet_readr$`2001`)),horizontal = TRUE,na.rm=TRUE)

boxplot(internet_readr$`2002`, 
     main= "Int. usage for 2002 202 Countries", 
     col="blue", xlab=paste("The median is:",median(internet_readr$`2002`)),horizontal = TRUE,na.rm=TRUE)

boxplot(internet_readr$`2003`, 
     main= "Int. usage for 2003 196 Countries", 
     col="blue", xlab=paste("The median is:",median(internet_readr$`2003`)),horizontal = TRUE,na.rm=TRUE)

boxplot(internet_readr$`2004`, 
     main= "Int. usage for 2004 199 Countries", 
     col="blue", xlab=paste("The median is:",median(internet_readr$`2004`)),horizontal = TRUE,na.rm=TRUE)

boxplot(internet_readr$`2005`, 
     main= "Int. usage for 2005 200 Countries", 
     col="blue", xlab=paste("The median is:",median(internet_readr$`2005`)),horizontal = TRUE,na.rm=TRUE)

boxplot(internet_readr$`2006`, 
     main= "Int. usage for 2006 199 Countries", 
     col="blue", xlab=paste("The median is:",median(internet_readr$`2006`)),horizontal = TRUE,na.rm=TRUE)

boxplot(internet_readr$`2007`, 
     main= "Int. usage for 2007 204 Countries", 
     col="blue", xlab=paste("The median is:",median(internet_readr$`2007`)),horizontal = TRUE,na.rm=TRUE)

boxplot(internet_readr$`2008`, 
     main= "Int. usage for 2008 202 Countries", 
     col="blue", xlab=paste("The median is:",median(internet_readr$`2008`)),horizontal = TRUE,na.rm=TRUE)

boxplot(internet_readr$`2009`, 
     main= "Int. usage for 2009 202 Countries", 
     col="blue", xlab=paste("The median is:",median(internet_readr$`2009`)),horizontal = TRUE,na.rm=TRUE)

boxplot(internet_readr$`2010`, 
     main= "Int. usage for 2010 201 Countries", 
     col="blue", xlab=paste("The median is:",median(internet_readr$`2010`)),horizontal = TRUE,na.rm=TRUE)

boxplot(internet_readr$`2011`, 
     main= "Int. usage for 2011 204 Countries", 
     col="blue", xlab=paste("The median is:",median(internet_readr$`2011`)),horizontal = TRUE,na.rm=TRUE)

boxplot(internet_readr$`2012`, 
     main= "Int. usage for 2012 202 Countries", 
     col="blue", xlab=paste("The median is:",median(internet_readr$`2012`)),horizontal = TRUE,na.rm=TRUE)

boxplot(internet_readr$`2013`, 
     main= "Int. usage for 2013 201 Countries", 
     col="blue", xlab=paste("The median is:",median(internet_readr$`2013`)),horizontal = TRUE,na.rm=TRUE)

boxplot(internet_readr$`2014`, 
     main= "Int. usage for 2014 201 Countries", 
     col="blue", xlab=paste("The median is:",median(internet_readr$`2014`)),horizontal = TRUE,na.rm=TRUE)

boxplot(internet_readr$`2015`, 
     main= "Int. usage for 2015 200 Countries", 
     col="blue", xlab=paste("The median is:",median(internet_readr$`2015`)),horizontal = TRUE,na.rm=TRUE)

boxplot(internet_readr$`2016`, 
     main= "Int. usage for 2016 203 Countries", 
     col="blue", xlab=paste("The median is:",median(internet_readr$`2016`)),horizontal = TRUE,na.rm=TRUE)

boxplot(internet_readr$`2017`, 
     main= "Int. usage for 2017 205 Countries", 
     col="blue", xlab=paste("The median is:",median(internet_readr$`2017`)),horizontal = TRUE,na.rm=TRUE)

boxplot(internet_readr$`2018`, 
     main= "Int. usage for 2018 169 Countries", 
     col="blue", xlab=paste("The median is:",median(internet_readr$`2018`)),horizontal = TRUE,na.rm=TRUE)

boxplot(internet_readr$`2019`, 
     main= "Int. usage for 2019 186 Countries", 
     col="blue", xlab=paste("The median is:",median(internet_readr$`2019`)),horizontal = TRUE,na.rm=TRUE)

boxplot(internet_readr$`2020`, 
     main= "Int. usage for 2020 185 Countries", 
     col="blue", xlab=paste("The median is:",median(internet_readr$`2020`)),horizontal = TRUE,na.rm=TRUE)

boxplot(internet_readr$`2021`, 
     main= "Int. usage for 2021 182 Countries", 
     col="blue", xlab=paste("The median is:",median(internet_readr$`2021`)),horizontal = TRUE,na.rm=TRUE)

boxplot(internet_readr$`2022`, 
     main= "Int usage for 2022 57 Countries", 
     col="blue", xlab=paste("The median is:",median(internet_readr$`2022`)),horizontal = TRUE,na.rm=TRUE)

#Now let’s move to the exercise 3.

library(readr)
library(dplyr)
bikeshare <- read_csv("bikesharedailydata.csv")

## Rows: 731 Columns: 16
## ── Column specification ────────────────────────────────────────────────────────
## Delimiter: ","
## chr  (1): dteday
## dbl (15): instant, season, yr, mnth, holiday, weekday, workingday, weathersi...
## 
## ℹ Use `spec()` to retrieve the full column specification for this data.
## ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.

Average_temp<-function(bikeshare){
  mean(bikeshare$temp,na.rm=TRUE)
}
Mean_temp<-Average_temp(bikeshare)
print(Mean_temp)

## [1] 0.4953848