#####################################################################################
##
##  Intro STATS : Summary Statistics
##  JDRS
##
#####################################################################################


#Get WB Data

WB <- read.csv("http://reuningscherer.net/stat10x/data/WorldBank2013.csv")

#see variables available
names(WB)

#summary stats for fertility rates - includes mean, median, min, max, quartiles
summary(WB$Fertility)

#variance and standard deviation
var(WB$Fertility)  #notice this doesn't work due to missing values

#fix missing value problem
var(WB$Fertility, na.rm=T)

#standard deviation
sd(WB$Fertility, na.rm=T)

#percentiles
quantile(WB$Fertility, c(.25, .50,  .75, .90, .99), na.rm=T)