---
title: "Analysis of Statistics Major students"
author: "Shawn Chen"
output: 
  flexdashboard::flex_dashboard:
    vertical_layout: fill
    source_code: embed
---

```{r setup, include=FALSE}
#load packages
library(flexdashboard)
library(readr)
library(tidyverse)
library(ggthemes)
library(plotly)
library(knitr)
library(kableExtra)
library(leaflet)
set.seed(1992-01-16)
```

```{r, message=FALSE, warning=FALSE, include=FALSE}
#load data and data clean
gre <- read_csv("~/Downloads/Analysis of Statistics Major/gre.csv")
count = read_csv("~/Downloads/Analysis of Statistics Major/count.csv")
glimpse(gre)
names(gre)[7] = paste("Enroll_year")
names(gre)[9] = paste("Rank")


#create ID
gre$ID = seq(1, 544, 1)

# fill information
gre$College_ranking[gre$College_ranking == "A"] = "USA"


#handle mising values of Toefl variable
sum(gre$TOEFL == "WAIVED")
gre$TOEFL[which(is.na(gre$TOEFL))] = median(gre$TOEFL,na.rm = T)

# missing ratio
129 / 544
gre$TOEFL[gre$TOEFL == "WAIVED"] = NA
class(gre$TOEFL)
median(gre$TOEFL,na.rm = T) #mdian is 104
gre$TOEFL[which(is.na(gre$TOEFL))] = 104
gre$TOEFL = as.numeric(as.character(gre$TOEFL))


#missing value of GRE writng score
gre %>%
  select(College_ranking,GRE_wirting) %>%
  group_by(College_ranking) %>%
  summarize(median(GRE_wirting,na.rm = T))

class(gre$GRE_wirting)

#211
college_211 = gre %>%
  filter(College_ranking == "211") %>%
  select(GRE_wirting, ID)

college_211$GRE_wirting[which(is.na(college_211$GRE_wirting))] = 3.50

#985
college_985 = gre %>%
  filter(College_ranking == "985") %>%
  select(GRE_wirting, ID)

college_985$GRE_wirting[which(is.na(college_985$GRE_wirting))] = 3.00

#USA
college_USA = gre %>%
  filter(College_ranking == "USA") %>%
  select(GRE_wirting, ID)

college_USA$GRE_wirting[which(is.na(college_USA$GRE_wirting))] = 3.50

#C9
college_C9 = gre %>%
  filter(College_ranking == "C9") %>%
  select(GRE_wirting, ID)

college_C9$GRE_wirting[which(is.na(college_C9$GRE_wirting))] = 3.50

#OTHER
college_OTHER = gre %>%
  filter(College_ranking == "OTHER") %>%
  select(GRE_wirting, ID)

college_OTHER$GRE_wirting[which(is.na(college_OTHER$GRE_wirting))] = 3.00

gre_wirthing_nomissing = college_211 %>%
  rbind(college_985) %>%
  rbind(college_USA) %>%
  rbind(college_C9) %>%
  rbind(college_OTHER)

gre_wirthing_nomissing = gre_wirthing_nomissing %>% arrange(ID)

#replacing missing value
gre$GRE_wirting = gre_wirthing_nomissing$GRE_wirting

# missing value on Statistics_or_Not
gre$Statistics_or_Not[which(is.na(gre$Statistics_or_Not))] = 0
sum(is.na(gre$Statistics_or_Not))

gre$Statistics_or_Not[gre$Statistics_or_Not == "1"] = "STAT"
gre$Statistics_or_Not[gre$Statistics_or_Not == "0"] = "NOT_STAT"

# Enroll year
table(gre$Enroll_year)
gre$Enroll_year[gre$Enroll_year == "1"] = "2015"
gre$Enroll_year[gre$Enroll_year == "2"] = "2016"
gre$Enroll_year[gre$Enroll_year == "3"] = "2017"
```


```{r, message=FALSE, warning=FALSE, include=FALSE}
#select subset to show
show = gre
colnames(show)
show2 = show[c("College_ranking" ,"GPA", "TOEFL", "GRE_total", "Statistics_or_Not", "Enroll_year", "ID")]
colnames(show2) = c("UNIVERSITY", "GPA", "TOEFL", "GRE", "MAJOR", "ACADEMIC YEAR", "ID")
```

Data Summary
=======================================================================
Column 
-----------------------------------------------------------------------

### Data information

```{r}
kable(show2, "html") %>%
  kable_styling() %>%
  scroll_box()
```

Column 
-----------------------------------------------------------------------

### The count information for each College Type

```{r}
p1 = ggplot(data = gre)  + 
  geom_bar(aes(x = College_ranking, fill = College_ranking,y = ..count.., alpha =0.6)) +
  coord_flip() +
  theme_hc() + 
  theme(plot.title = element_text(hjust = 0.5)) +
  labs(
    title = paste("The number of students in each university"),
    x = paste(" "),
    y = paste("Count")) + 
  scale_fill_discrete(name = " ")

ggplotly(p1)
```

### Chart The count information for student's major

```{r}
p10 = ggplot(data = gre) + 
  geom_bar(aes(x = factor(Statistics_or_Not), fill = Statistics_or_Not, y = ..count.., alpha =0.6)) +
  coord_flip() + 
  theme_hc() + 
  theme(plot.title = element_text(hjust = 0.5)) +
  labs(
    title = paste("The number of students STAT major and non-STAT major"),
    x = paste(" "),
    y = paste(" ")) + 
  scale_fill_discrete(name = paste(" "))

ggplotly(p10) 
```

### The count information for each academic year
```{r}
p15 = ggplot(data = gre) + 
  geom_bar(aes(x = factor(Enroll_year), fill = Enroll_year, y = ..count.., alpha =0.6)) +
  coord_flip() + 
  theme_hc() + 
  theme(plot.title = element_text(hjust = 0.5)) +
  labs(
    title = paste("The number of students in each academic year"),
    x = paste(" "),
    y = paste(" ")) + 
  scale_fill_discrete(name = paste(" "))

ggplotly(p15) 

```

GPA Information
===========================================================================

Column
---------------------------------------------------------------------------
###Distribution plot for GPA
```{r}
p2 = ggplot(data = gre) + 
  geom_density(aes(x = GPA),fill = "#b4dfdc", alpha = 0.6) +
  theme_hc() + 
  theme(plot.title = element_text(hjust = 0.5)) +
  labs(
    title = paste("GPA score's distribution"),
    x = paste("GPA"),
    y = paste(" "))

ggplotly(p2)
```

###The Boxplot for GPA score with different University Type
```{r}
plot_ly(data = gre, x = ~ GPA, y = ~ factor(College_ranking), color = ~College_ranking, type = "box") %>%
  layout(
    title = "GPA score in each university",
    xaxis = list(title = "GPA"),
    yaxis = list(title = "")
  )
```

Column
--------------------------------------------------------------------------
###The Boxplot for GPA score with different Major Type
```{r}
plot_ly(data = gre, x = ~ GPA, color = ~factor(Statistics_or_Not), type = "box") %>%
  layout(
    title = "GPA score in STAT and non-STAT major",
    xaxis = list(title = "GPA"),
    yaxis = list(title = "")
  )
```

###The Boxplot for GPA score in different academic year
```{r}
plot_ly(data = gre, x = ~GPA, color = ~Enroll_year, type = "box") %>%
  layout(
    title = "GPA Total score in each academic year",
    xaxis = list(title = "GPA"),
    yaxis = list(title = "")
  )
```



TOEFL Information
===========================================================================

Row
---------------------------------------------------------------------------
###The distribution for TOEFL score
```{r}
pp = ggplot(data = gre) + 
  geom_density(aes(x = TOEFL),fill = "#b4dfdc", alpha =0.6) +
  theme_hc() + 
  theme(plot.title = element_text(hjust = 0.5)) +
  labs(
    title = paste("TOEFL score's distribution"),
    x = paste("TOEFL SCORE"),
    y = paste(" "))
ggplotly(pp)
```

###The Boxplot for TOEFL score with different University Type
```{r}
plot_ly(data = gre, x = ~ TOEFL, y = ~ factor(College_ranking), color = ~College_ranking, type = "box") %>%
  layout(
    title = "TOEFL score in each university",
    xaxis = list(title = "TOEFL"),
    yaxis = list(title = "")
  )
```

Row
---
###The Boxplot for TOEFL score with different Major Type
```{r}
plot_ly(data = gre, x = ~ TOEFL, color = ~factor(Statistics_or_Not), type = "box") %>%
  layout(
    title = "TOEFL score in STAT and non-STAT major",
    xaxis = list(title = "TOEFL"),
    yaxis = list(title = "")
  )
```

###The Boxplot for TOEFL score in different Enroll year
```{r, echo=FALSE, message=FALSE, warning=FALSE}
plot_ly(data = gre, x = ~TOEFL, color = ~Enroll_year, type = "box") %>%
  layout(
    title = "TOEFL Total score in each academic year",
    xaxis = list(title = "TOEFL"),
    yaxis = list(title = "")
  )
```


GRE
==================================================================================
Row
----------------------------------------------------------------------------------
###The distribution for GRE score
```{r}
p6 = ggplot(data = gre) + 
  geom_density(aes(x = GRE_total),fill = "#b4dfdc", alpha =0.6) +
  theme_hc() + 
  theme(plot.title = element_text(hjust = 0.5)) +
  labs(
    title = paste("GRE score's distribution"),
    x = paste("GRE SCORE"),
    y = paste(" "))
ggplotly(p6)
```

###The Boxplot for GRE score with different University Type
```{r}
plot_ly(data = gre, x = ~GRE_total, y = ~ factor(College_ranking), color = ~College_ranking, type = "box") %>%
  layout(
    title = "GRE score in each university",
    xaxis = list(title = "GRE"),
    yaxis = list(title = "")
  )
```

Row
--------------------------------------------------------------------
###The Boxplot for GRE score with different Major Type
```{r}
###
plot_ly(data = gre, x = ~GRE_total, color = ~Statistics_or_Not, type = "box") %>%
  layout(
    title = "GRE score in STAT and non-STAT major",
    xaxis = list(title = "GRE"),
    yaxis = list(title = "")
  )
```

###The Boxplot for GRE Total score in different academic year
```{r}
plot_ly(data = gre, x = ~GRE_total, color = ~Enroll_year, type = "box") %>%
  layout(
    title = "GRE Total score in each academic year",
    xaxis = list(title = "GRE"),
    yaxis = list(title = "")
  )
```

GIS
=============================================================
Row{.tabset}
-------------------------------------------------------------

### Map
```{r}
latitude = c( 40.8075, 40.116421, 42.2780,38.8997, 41.3163, 47.6553, 42.4534, 40.4982, 37.8719, 41.8077,
              42.3505, 33.7925, 33.6405, 36.0014, 38.0336, 38.6488, 39.3299,
              41.7886, 29.7174, 42.3770, 41.8268, 44.9740, 30.4419, 34.0689,
              35.9049, 32.8801, 38.9076, 40.6069, 40.0142, 39.9522, 35.7847,
              43.0392, 40.4444, 32.9858, 38.5382, 42.055984, 40.7295, 29.9403,
              42.7018, 33.7756, 40.4428, 34.4140, 32.8412, 43.0783, 40.4237,
              40.4444, 41.5043, 44.5638, 39.6780, 42.2746, 29.6436, 37.4275,
              29.7034, 36.1352, 38.9404, 40.9124, 33.424564, 43.1306, 36.9914,
              30.6185, 41.6627, 39.5105, 40.7448, 34.0224, 49.2606, 39.1761,
              40.5734)

longitude =c(-73.9626, -88.243385, -83.7382, -77.0486, -72.9223, -122.3035, -76.4735, -74.4468, -122.2585, - 72.2540,
             -71.1054, -84.3240, -117.8443, -78.9382, -78.5080, -90.3108, -76.6205,
             -87.5987, -95.4018, -71.1167, -71.4025, -93.2277, -84.2985, -118.4452,
             -79.0469, -117.2340, -77.0723, -75.3783, -83.0309, -75.1932, -78.6821,
             -76.1351, -79.9608, -96.7501, -121.7617, -87.675171, -73.9965, -90.1207,
             -84.4822, -84.3963, -79.9430, -119.8489, -96.7845, -87.8820, -86.9212,
             -79.9608, -81.6084, -123.2794, -75.7506, -71.8063, -82.3549, -122.1697,
             -95.4030, -80.2763, -92.3277, -73.1234, -111.928001, -77.6260, -122.0609,
             -96.3365, -91.5549, -84.7309, -74.0257, -118.2851, -123.2460, -86.5131,
             -105.0865)

COLUMBIA = makeIcon(
  iconUrl = "/Users/xiaochen/Downloads/Analysis of Statistics Major/logo/Columbia University.png", 
  iconWidth = 60, iconHeight = 60,
  iconAnchorX = 31*215/230/2, iconAnchorY = 16
)

UIUC <- makeIcon(
  iconUrl = "/Users/xiaochen/Downloads/Analysis of Statistics Major/logo/UIUC.png", 
  iconWidth = 60, iconHeight = 60,
  iconAnchorX = 31*215/230/2, iconAnchorY = 16
)

UMICH <- makeIcon(
  iconUrl = "/Users/xiaochen/Downloads/Analysis of Statistics Major/logo/UMICH.png", 
  iconWidth = 60, iconHeight = 60,
  iconAnchorX = 31*215/230/2, iconAnchorY = 16
)

GWU <- makeIcon(
  iconUrl = "/Users/xiaochen/Downloads/Analysis of Statistics Major/logo/GWU.png", 
  iconWidth = 60, iconHeight = 60,
  iconAnchorX = 31*215/230/2, iconAnchorY = 16
)

YALE <- makeIcon(
  iconUrl = "/Users/xiaochen/Downloads/Analysis of Statistics Major/logo/Yale.png", 
  iconWidth = 60, iconHeight = 60,
  iconAnchorX = 31*215/230/2, iconAnchorY = 16
)

UW <- makeIcon(
  iconUrl = "/Users/xiaochen/Downloads/Analysis of Statistics Major/logo/UW.png", 
  iconWidth = 60, iconHeight = 60,
  iconAnchorX = 31*215/230/2, iconAnchorY = 16
)

CORNELL <- makeIcon(
  iconUrl = "/Users/xiaochen/Downloads/Analysis of Statistics Major/logo/CORNELL.png", 
  iconWidth = 60, iconHeight = 60,
  iconAnchorX = 31*215/230/2, iconAnchorY = 16
)

RUTGERS <- makeIcon(
  iconUrl = "/Users/xiaochen/Downloads/Analysis of Statistics Major/logo/RUTGERS.png", 
  iconWidth = 60, iconHeight = 60,
  iconAnchorX = 31*215/230/2, iconAnchorY = 16
)

UCB <- makeIcon(
  iconUrl = "/Users/xiaochen/Downloads/Analysis of Statistics Major/logo/UCB.png", 
  iconWidth = 100, iconHeight = 60,
  iconAnchorX = 31*215/230/2, iconAnchorY = 16
)

UCONN <- makeIcon(
  iconUrl = "/Users/xiaochen/Downloads/Analysis of Statistics Major/logo/UCONN.png", 
  iconWidth = 60, iconHeight = 60,
  iconAnchorX = 31*215/230/2, iconAnchorY = 16
)

BU <- makeIcon(
  iconUrl = "/Users/xiaochen/Downloads/Analysis of Statistics Major/logo/BU.png", 
  iconWidth = 60, iconHeight = 60,
  iconAnchorX = 31*215/230/2, iconAnchorY = 16
)

EMORY <- makeIcon(
  iconUrl = "/Users/xiaochen/Downloads/Analysis of Statistics Major/logo/EMORY.png", 
  iconWidth = 60, iconHeight = 80,
  iconAnchorX = 31*215/230/2, iconAnchorY = 16
)

UCI <- makeIcon(
  iconUrl = "/Users/xiaochen/Downloads/Analysis of Statistics Major/logo/UCI.png", 
  iconWidth = 60, iconHeight = 80,
  iconAnchorX = 31*215/230/2, iconAnchorY = 16
)

DUKE <- makeIcon(
  iconUrl = "/Users/xiaochen/Downloads/Analysis of Statistics Major/logo/DUKE.png", 
  iconWidth = 60, iconHeight = 80,
  iconAnchorX = 31*215/230/2, iconAnchorY = 16
)

UVA <- makeIcon(
  iconUrl = "/Users/xiaochen/Downloads/Analysis of Statistics Major/logo/UVA.png", 
  iconWidth = 60, iconHeight = 60,
  iconAnchorX = 31*215/230/2, iconAnchorY = 16
)

WUSTL <- makeIcon(
  iconUrl = "/Users/xiaochen/Downloads/Analysis of Statistics Major/logo/WUSTL.png", 
  iconWidth = 60, iconHeight = 60,
  iconAnchorX = 31*215/230/2, iconAnchorY = 16
)

JHU <- makeIcon(
  iconUrl = "/Users/xiaochen/Downloads/Analysis of Statistics Major/logo/JHU.png", 
  iconWidth = 60, iconHeight = 60,
  iconAnchorX = 31*215/230/2, iconAnchorY = 16
)

UCHICAGO <- makeIcon(
  iconUrl = "/Users/xiaochen/Downloads/Analysis of Statistics Major/logo/UCHICAGO.png", 
  iconWidth = 60, iconHeight = 60,
  iconAnchorX = 31*215/230/2, iconAnchorY = 16
)

RICE <- makeIcon(
  iconUrl = "/Users/xiaochen/Downloads/Analysis of Statistics Major/logo/RICE.png", 
  iconWidth = 60, iconHeight = 60,
  iconAnchorX = 31*215/230/2, iconAnchorY = 16
)

HARVARD <- makeIcon(
  iconUrl = "/Users/xiaochen/Downloads/Analysis of Statistics Major/logo/HARVARD.png", 
  iconWidth = 60, iconHeight = 60,
  iconAnchorX = 31*215/230/2, iconAnchorY = 16
)

BROWN <- makeIcon(
  iconUrl = "/Users/xiaochen/Downloads/Analysis of Statistics Major/logo/Brown.png", 
  iconWidth = 150, iconHeight = 60,
  iconAnchorX = 20, iconAnchorY = 16
)

UMN <- makeIcon(
  iconUrl = "/Users/xiaochen/Downloads/Analysis of Statistics Major/logo/UMN.png", 
  iconWidth = 60, iconHeight = 60,
  iconAnchorX = 31*215/230/2, iconAnchorY = 16
)


FSU <- makeIcon(
  iconUrl = "/Users/xiaochen/Downloads/Analysis of Statistics Major/logo/FSU.png", 
  iconWidth = 60, iconHeight = 60,
  iconAnchorX = 31*215/230/2, iconAnchorY = 16
)


UCLA <- makeIcon(
  iconUrl = "/Users/xiaochen/Downloads/Analysis of Statistics Major/logo/UCLA.png", 
  iconWidth = 60, iconHeight = 60,
  iconAnchorX = 31*215/230/2, iconAnchorY = 16
)


UNC <- makeIcon(
  iconUrl = "/Users/xiaochen/Downloads/Analysis of Statistics Major/logo/UNC.png", 
  iconWidth = 60, iconHeight = 60,
  iconAnchorX = 31*215/230/2, iconAnchorY = 16
)


UCSD <- makeIcon(
  iconUrl = "/Users/xiaochen/Downloads/Analysis of Statistics Major/logo/UCSD.png", 
  iconWidth = 60, iconHeight = 60,
  iconAnchorX = 31*215/230/2, iconAnchorY = 16
)


GERGETOWN <- makeIcon(
  iconUrl = "/Users/xiaochen/Downloads/Analysis of Statistics Major/logo/GERGETOWN.png", 
  iconWidth = 60, iconHeight = 60,
  iconAnchorX = 31*215/230/2, iconAnchorY = 16
)


LEHIGH <- makeIcon(
  iconUrl = "/Users/xiaochen/Downloads/Analysis of Statistics Major/logo/LEHIGH.png", 
  iconWidth = 60, iconHeight = 80,
  iconAnchorX = 31*215/230/2, iconAnchorY = 16
)


OSU <- makeIcon(
  iconUrl = "/Users/xiaochen/Downloads/Analysis of Statistics Major/logo/OSU.png", 
  iconWidth = 60, iconHeight = 60,
  iconAnchorX = 31*215/230/2, iconAnchorY = 16
)


UPEN <- makeIcon(
  iconUrl = "/Users/xiaochen/Downloads/Analysis of Statistics Major/logo/upenn.png", 
  iconWidth = 100, iconHeight = 60,
  iconAnchorX = 31*215/230/2, iconAnchorY = 16
)


NCSU <- makeIcon(
  iconUrl = "/Users/xiaochen/Downloads/Analysis of Statistics Major/logo/NCSU.png", 
  iconWidth = 60, iconHeight = 60,
  iconAnchorX = 31*215/230/2, iconAnchorY = 16
)


SYRACUSE <- makeIcon(
  iconUrl = "/Users/xiaochen/Downloads/Analysis of Statistics Major/logo/SYRACUSE.png", 
  iconWidth = 60, iconHeight = 60,
  iconAnchorX = 31*215/230/2, iconAnchorY = 16
)


UPIT <- makeIcon(
  iconUrl = "/Users/xiaochen/Downloads/Analysis of Statistics Major/logo/UPIT.png", 
  iconWidth = 60, iconHeight = 60,
  iconAnchorX = 31*215/230/2, iconAnchorY = 16
)


UTD<- makeIcon(
  iconUrl = "/Users/xiaochen/Downloads/Analysis of Statistics Major/logo/UTD.png", 
  iconWidth = 100, iconHeight = 60,
  iconAnchorX = 31*215/230/2, iconAnchorY = 16
)


UCD <- makeIcon(
  iconUrl = "/Users/xiaochen/Downloads/Analysis of Statistics Major/logo/UCD.png", 
  iconWidth = 60, iconHeight = 60,
  iconAnchorX = 31*215/230/2, iconAnchorY = 16
)


NORTHWESTERN <- makeIcon(
  iconUrl = "/Users/xiaochen/Downloads/Analysis of Statistics Major/logo/NORTHWESTERN.png", 
  iconWidth = 60, iconHeight = 60,
  iconAnchorX = 31*215/230/2, iconAnchorY = 16
)


NYU <- makeIcon(
  iconUrl = "/Users/xiaochen/Downloads/Analysis of Statistics Major/logo/NYU.png", 
  iconWidth = 60, iconHeight = 60,
  iconAnchorX = 31*215/230/2, iconAnchorY = 16
)


TULANE <- makeIcon(
  iconUrl = "/Users/xiaochen/Downloads/Analysis of Statistics Major/logo/TULANE.png", 
  iconWidth = 100, iconHeight = 60,
  iconAnchorX = 31*215/230/2, iconAnchorY = 16
)


MSU <- makeIcon(
  iconUrl = "/Users/xiaochen/Downloads/Analysis of Statistics Major/logo/MSU.png", 
  iconWidth = 60, iconHeight = 60,
  iconAnchorX = 31*215/230/2, iconAnchorY = 16
)


GATECH <- makeIcon(
  iconUrl = "/Users/xiaochen/Downloads/Analysis of Statistics Major/logo/GATECH.png", 
  iconWidth = 60, iconHeight = 60,
  iconAnchorX = 31*215/230/2, iconAnchorY = 16
)


CMU <- makeIcon(
  iconUrl = "/Users/xiaochen/Downloads/Analysis of Statistics Major/logo/CMU.png", 
  iconWidth = 60, iconHeight = 60,
  iconAnchorX = 31*215/230/2, iconAnchorY = 16
)


UCSB <- makeIcon(
  iconUrl = "/Users/xiaochen/Downloads/Analysis of Statistics Major/logo/UCSB.png", 
  iconWidth = 100, iconHeight = 60,
  iconAnchorX = 31*215/230/2, iconAnchorY = 16
)


SMU <- makeIcon(
  iconUrl = "/Users/xiaochen/Downloads/Analysis of Statistics Major/logo/SMU.png", 
  iconWidth = 60, iconHeight = 60,
  iconAnchorX = 31*215/230/2, iconAnchorY = 16
)


UWM <- makeIcon(
  iconUrl = "/Users/xiaochen/Downloads/Analysis of Statistics Major/logo/UWM.png", 
  iconWidth = 60, iconHeight = 80,
  iconAnchorX = 31*215/230/2, iconAnchorY = 16
)


PURDUE <- makeIcon(
  iconUrl = "/Users/xiaochen/Downloads/Analysis of Statistics Major/logo/PURDUE.png", 
  iconWidth = 60, iconHeight = 60,
  iconAnchorX = 31*215/230/2, iconAnchorY = 16
)


PITT <- makeIcon(
  iconUrl = "/Users/xiaochen/Downloads/Analysis of Statistics Major/logo/PITT.png", 
  iconWidth = 60, iconHeight = 60,
  iconAnchorX = 31*215/230/2, iconAnchorY = 16
)


CWRU <- makeIcon(
  iconUrl = "/Users/xiaochen/Downloads/Analysis of Statistics Major/logo/CWRU.png", 
  iconWidth = 60, iconHeight = 80,
  iconAnchorX = 31*215/230/2, iconAnchorY = 16
)


OREGONSTATE <- makeIcon(
  iconUrl = "/Users/xiaochen/Downloads/Analysis of Statistics Major/logo/OREGONSTATE.png", 
  iconWidth = 80, iconHeight = 60,
  iconAnchorX = 31*215/230/2, iconAnchorY = 16
)


UDELAWARE <- makeIcon(
  iconUrl = "/Users/xiaochen/Downloads/Analysis of Statistics Major/logo/UDELAWARE.png", 
  iconWidth = 60, iconHeight = 60,
  iconAnchorX = 31*215/230/2, iconAnchorY = 16
)


WPI <- makeIcon(
  iconUrl = "/Users/xiaochen/Downloads/Analysis of Statistics Major/logo/WPI.png", 
  iconWidth = 60, iconHeight = 60,
  iconAnchorX = 31*215/230/2, iconAnchorY = 16
)


UFL <- makeIcon(
  iconUrl = "/Users/xiaochen/Downloads/Analysis of Statistics Major/logo/UFL.png", 
  iconWidth = 100, iconHeight = 60,
  iconAnchorX = 31*215/230/2, iconAnchorY = 16
)


STANFORD <- makeIcon(
  iconUrl = "/Users/xiaochen/Downloads/Analysis of Statistics Major/logo/STANFORD.png", 
  iconWidth = 60, iconHeight = 80,
  iconAnchorX = 31*215/230/2, iconAnchorY = 16
)


UTH <- makeIcon(
  iconUrl = "/Users/xiaochen/Downloads/Analysis of Statistics Major/logo/UTH.png", 
  iconWidth = 60, iconHeight = 60,
  iconAnchorX = 31*215/230/2, iconAnchorY = 16
)


WFU <- makeIcon(
  iconUrl = "/Users/xiaochen/Downloads/Analysis of Statistics Major/logo/WFU.png", 
  iconWidth = 60, iconHeight = 60,
  iconAnchorX = 31*215/230/2, iconAnchorY = 16
)



MIZZOU <- makeIcon(
  iconUrl = "/Users/xiaochen/Downloads/Analysis of Statistics Major/logo/Mizzou.png", 
  iconWidth = 60, iconHeight = 80,
  iconAnchorX = 31*215/230/2, iconAnchorY = 16
)

SBU <- makeIcon(
  iconUrl = "/Users/xiaochen/Downloads/Analysis of Statistics Major/logo/SBU.png", 
  iconWidth = 60, iconHeight = 60,
  iconAnchorX = 31*215/230/2, iconAnchorY = 16
)


ASU <- makeIcon(
  iconUrl = "/Users/xiaochen/Downloads/Analysis of Statistics Major/logo/ASU.png", 
  iconWidth = 60, iconHeight = 60,
  iconAnchorX = 31*215/230/2, iconAnchorY = 16
)


ROCHESTER <- makeIcon(
  iconUrl = "/Users/xiaochen/Downloads/Analysis of Statistics Major/logo/ROCHESTER.png", 
  iconWidth = 100, iconHeight = 60,
  iconAnchorX = 31*215/230/2, iconAnchorY = 16
)


UCSC <- makeIcon(
  iconUrl = "/Users/xiaochen/Downloads/Analysis of Statistics Major/logo/UCSC.png", 
  iconWidth = 60, iconHeight = 60,
  iconAnchorX = 31*215/230/2, iconAnchorY = 16
)


TAMU <- makeIcon(
  iconUrl = "/Users/xiaochen/Downloads/Analysis of Statistics Major/logo/TAMU.png", 
  iconWidth = 60, iconHeight = 60,
  iconAnchorX = 31*215/230/2, iconAnchorY = 16
)


UIOWA <- makeIcon(
  iconUrl = "/Users/xiaochen/Downloads/Analysis of Statistics Major/logo/UIOWA.png", 
  iconWidth = 60, iconHeight = 60,
  iconAnchorX = 31*215/230/2, iconAnchorY = 16
)


MIAMI<- makeIcon(
  iconUrl = "/Users/xiaochen/Downloads/Analysis of Statistics Major/logo/MIAMI.png", 
  iconWidth = 60, iconHeight = 60,
  iconAnchorX = 31*215/230/2, iconAnchorY = 16
)


SIT <- makeIcon(
  iconUrl = "/Users/xiaochen/Downloads/Analysis of Statistics Major/logo/SIT.png", 
  iconWidth = 60, iconHeight = 60,
  iconAnchorX = 31*215/230/2, iconAnchorY = 16
)


USC <- makeIcon(
  iconUrl = "/Users/xiaochen/Downloads/Analysis of Statistics Major/logo/USC.png", 
  iconWidth = 60, iconHeight = 60,
  iconAnchorX = 31*215/230/2, iconAnchorY = 16
)

UBC <- makeIcon(
  iconUrl = "/Users/xiaochen/Downloads/Analysis of Statistics Major/logo/UBC.png", 
  iconWidth = 60, iconHeight = 80,
  iconAnchorX = 31*215/230/2, iconAnchorY = 16
)

IUB <- makeIcon(
  iconUrl = "/Users/xiaochen/Downloads/Analysis of Statistics Major/logo/IUB.png", 
  iconWidth = 60, iconHeight = 60,
  iconAnchorX = 31*215/230/2, iconAnchorY = 16
)

CSU <- makeIcon(
  iconUrl = "/Users/xiaochen/Downloads/Analysis of Statistics Major/logo/CSU.png", 
  iconWidth = 60, iconHeight = 80,
  iconAnchorX = 31*215/230/2, iconAnchorY = 16
)





iconlist = c(
             COLUMBIA,
             UIUC,
             UMICH,
             GWU,
             YALE,
             UW,
             CORNELL,
             RUTGERS,
             UCB,
             UCONN,
             BU,
             EMORY,
             UCI,
             DUKE,
             UVA,
             WUSTL,
             JHU,
             UCHICAGO,
             RICE,
             HARVARD,
             BROWN,
             UMN,
             FSU,
             UCLA,
             UNC,
             UCSD,
             GERGETOWN,
             LEHIGH,
             OSU,
             UPEN,
             NCSU,
             SYRACUSE,
             UPIT,
             UTD,
             UCD,
             NORTHWESTERN,
             NYU,
             TULANE,
             MSU,
             GATECH,
             CMU,
             UCSB,
             SMU,
             UWM,
             PURDUE,
             PITT,
             CWRU,
             OREGONSTATE,
             UDELAWARE,
             WPI,
             UFL,
             STANFORD,
             UTH,
             WFU,
             MIZZOU,
             SBU,
             ASU,
             ROCHESTER,
             UCSC,
             TAMU,
             UIOWA,
             MIAMI,
             SIT,
             USC,
             UBC,
             IUB,
             CSU)



leaflet() %>%
  addTiles() %>%
  addMarkers(lat = latitude[1] , lng = longitude[1],icon = iconlist[1:5],
             popup = paste("Columbia University")) %>%
  addMarkers(lat = latitude[2] , lng = longitude[2],icon = iconlist[6:10],
             popup = paste("University of Illinois at Urbana–Champaign")) %>%
  addMarkers(lat = latitude[3] , lng = longitude[3],icon = iconlist[11:15],
             popup = paste("University of Michigan")) %>%
  addMarkers(lat = latitude[4] , lng = longitude[4],icon = iconlist[16:20],
             popup = paste("George Washington University")) %>%
  addMarkers(lat = latitude[5] , lng = longitude[5],icon = iconlist[21:25],
             popup = paste("Yale University")) %>%
  addMarkers(lat = latitude[6] , lng = longitude[6],icon = iconlist[26:30],
             popup = paste("University of Washington")) %>%
  addMarkers(lat = latitude[7] , lng = longitude[7],icon = iconlist[31:35],
             popup = paste("Cornell University")) %>%
  addMarkers(lat = latitude[8] , lng = longitude[8],icon = iconlist[36:40],
             popup = paste("Rutgers University")) %>%
  addMarkers(lat = latitude[9] , lng = longitude[9],icon = iconlist[41:45],
             popup = paste("University of California, Berkeley")) %>%
  addMarkers(lat = latitude[10] , lng = longitude[10],icon = iconlist[46:50],
             popup = paste("University of Connecticut")) %>%
  addMarkers(lat = latitude[11] , lng = longitude[11],icon = iconlist[51:55],
             popup = paste("Boston University")) %>%
  addMarkers(lat = latitude[12] , lng = longitude[12],icon = iconlist[56:60],
             popup = paste("Emory University")) %>%
  addMarkers(lat = latitude[13] , lng = longitude[13],icon = iconlist[61:65],
             popup = paste("University of California, Irvine")) %>%
  addMarkers(lat = latitude[14] , lng = longitude[14],icon = iconlist[66:70],
             popup = paste("Duke University")) %>%
  addMarkers(lat = latitude[15] , lng = longitude[15],icon = iconlist[71:75],
             popup = paste("University of Virginia")) %>%
  addMarkers(lat = latitude[16] , lng = longitude[16],icon = iconlist[76:80],
             popup = paste("Washington University in St. Louis")) %>%
  addMarkers(lat = latitude[17] , lng = longitude[17],icon = iconlist[81:85],
             popup = paste("Johns Hopkins University")) %>%
  addMarkers(lat = latitude[18] , lng = longitude[18],icon = iconlist[86:90],
             popup = paste("University of Chicago")) %>%
  addMarkers(lat = latitude[19] , lng = longitude[19],icon = iconlist[91:95],
             popup = paste("Rice University")) %>%
  addMarkers(lat = latitude[20] , lng = longitude[20],icon = iconlist[96:100],
             popup = paste("Harvard University")) %>%
  addMarkers(lat = latitude[21] , lng = longitude[21],icon = iconlist[101:105],
             popup = paste("Brown University")) %>%
  addMarkers(lat = latitude[22] , lng = longitude[22],icon = iconlist[106:110],
             popup = paste("University of Wisconsin–Milwaukee")) %>%
  addMarkers(lat = latitude[23] , lng = longitude[23],icon = iconlist[111:115],
             popup = paste("Florida State University")) %>%
  addMarkers(lat = latitude[24] , lng = longitude[24],icon = iconlist[116:120],
             popup = paste("University of California, Los Angeles")) %>%
  addMarkers(lat = latitude[25] , lng = longitude[25],icon = iconlist[121:125],
             popup = paste("University of North Carolina at Chapel Hill")) %>%
  addMarkers(lat = latitude[26] , lng = longitude[26],icon = iconlist[126:130],
             popup = paste("University of California, San Diego")) %>%
  addMarkers(lat = latitude[27] , lng = longitude[27],icon = iconlist[131:135],
             popup = paste("Georgetown University")) %>%
  addMarkers(lat = latitude[28] , lng = longitude[28],icon = iconlist[136:140],
             popup = paste("Lehigh University University")) %>%
  addMarkers(lat = latitude[29] , lng = longitude[29],icon = iconlist[141:145],
             popup = paste("Ohio State University")) %>%
  addMarkers(lat = latitude[30] , lng = longitude[30],icon = iconlist[146:150],
             popup = paste("University of Pennsylvania")) %>%
  addMarkers(lat = latitude[31] , lng = longitude[31],icon = iconlist[151:155],
             popup = paste("North Carolina State University")) %>%
  addMarkers(lat = latitude[32] , lng = longitude[32],icon = iconlist[156:160],
             popup = paste("Syracuse University")) %>%
  addMarkers(lat = latitude[33] , lng = longitude[33],icon = iconlist[161:165],
             popup = paste("University of Pitești")) %>%
  addMarkers(lat = latitude[34] , lng = longitude[34],icon = iconlist[166:170],
             popup = paste("University of Texas at Dallas")) %>%
  addMarkers(lat = latitude[35] , lng = longitude[35],icon = iconlist[171:175],
             popup = paste("University College Dublin")) %>%
  addMarkers(lat = latitude[36] , lng = longitude[36],icon = iconlist[176:180],
             popup = paste("Northwestern University")) %>%
  addMarkers(lat = latitude[37] , lng = longitude[37],icon = iconlist[181:185],
             popup = paste("New York University")) %>%
  addMarkers(lat = latitude[38] , lng = longitude[38],icon = iconlist[186:190],
             popup = paste("Tulane University")) %>%
  addMarkers(lat = latitude[39] , lng = longitude[39],icon = iconlist[191:195],
             popup = paste("Michigan State University")) %>%
  addMarkers(lat = latitude[40] , lng = longitude[40],icon = iconlist[196:200],
             popup = paste("Georgia Institute of Technology")) %>%
  addMarkers(lat = latitude[41] , lng = longitude[41],icon = iconlist[201:205],
             popup = paste("Carnegie Mellon University")) %>%
  addMarkers(lat = latitude[42] , lng = longitude[42],icon = iconlist[206:210],
             popup = paste("University of California, Santa Barbara")) %>%
  addMarkers(lat = latitude[43] , lng = longitude[43],icon = iconlist[211:215],
             popup = paste("Southern Methodist University")) %>%
  addMarkers(lat = latitude[44] , lng = longitude[44],icon = iconlist[216:220],
             popup = paste("University of Wisconsin–Milwaukee")) %>%
  addMarkers(lat = latitude[45] , lng = longitude[45],icon = iconlist[221:225],
             popup = paste("Purdue University")) %>%
  addMarkers(lat = latitude[46] , lng = longitude[46],icon = iconlist[226:230],
             popup = paste("University of Pittsburgh")) %>%
  addMarkers(lat = latitude[47] , lng = longitude[47],icon = iconlist[231:235],
             popup = paste("Case Western Reserve University")) %>%
  addMarkers(lat = latitude[48] , lng = longitude[48],icon = iconlist[236:240],
             popup = paste("Oregon State University")) %>%
  addMarkers(lat = latitude[49] , lng = longitude[49],icon = iconlist[241:245],
             popup = paste("University of Delaware")) %>%
  addMarkers(lat = latitude[50] , lng = longitude[50],icon = iconlist[246:250],
             popup = paste("Worcester Polytechnic Institute")) %>%
  addMarkers(lat = latitude[51] , lng = longitude[51],icon = iconlist[251:255],
             popup = paste("University of Florida")) %>%
  addMarkers(lat = latitude[52] , lng = longitude[52],icon = iconlist[256:260],
             popup = paste("Stanford University")) %>%
  addMarkers(lat = latitude[53] , lng = longitude[53],icon = iconlist[261:265],
             popup = paste("University of Texas Health Science Center at Houston")) %>%
  addMarkers(lat = latitude[54] , lng = longitude[54],icon = iconlist[266:270],
             popup = paste("Wake Forest University")) %>%
  addMarkers(lat = latitude[55] , lng = longitude[55],icon = iconlist[271:275],
             popup = paste("University of MIssouri")) %>%
  addMarkers(lat = latitude[56] , lng = longitude[56],icon = iconlist[276:280],
             popup = paste("Stony Brook University")) %>%
  addMarkers(lat = latitude[57] , lng = longitude[57],icon = iconlist[281:285],
             popup = paste("Arizona State University")) %>%
  addMarkers(lat = latitude[58] , lng = longitude[58],icon = iconlist[286:290],
             popup = paste("University of Rochester")) %>%
  addMarkers(lat = latitude[59] , lng = longitude[59],icon = iconlist[291:300],
             popup = paste("University of California, Santa Cruz")) %>%
  addMarkers(lat = latitude[60] , lng = longitude[60],icon = iconlist[296:305],
             popup = paste("Texas A&M University")) %>%
  addMarkers(lat = latitude[61] , lng = longitude[61],icon = iconlist[301:310],
             popup = paste("University of Iowa")) %>%
  addMarkers(lat = latitude[62] , lng = longitude[62],icon = iconlist[306:315],
             popup = paste("Miami University")) %>%
  addMarkers(lat = latitude[63] , lng = longitude[63],icon = iconlist[311:320],
             popup = paste("Stevens Institute of Technology")) %>%
  addMarkers(lat = latitude[64] , lng = longitude[64],icon = iconlist[316:325],
             popup = paste("University of Southern California")) %>%
  addMarkers(lat = latitude[65] , lng = longitude[65],icon = iconlist[321:330],
             popup = paste("University of British Columbia")) %>%
  addMarkers(lat = latitude[66] , lng = longitude[66],icon = iconlist[326:335],
             popup = paste("Indiana University Bloomington")) %>%
  addMarkers(lat = latitude[67] , lng = longitude[67],icon = iconlist[331:340],
             popup = paste("Colorado State University")) 

```

### Bar plot

```{r echo=FALSE, message=FALSE, warning=FALSE}
count$School = factor(count$School, levels = unique(count$School[order(count$Count,decreasing = F)]))
plot_ly(data = count, x = ~Count, y = ~School, type = 'bar')
```