#Gr?ficos num?ricas
idh<-read.csv(“IDH.csv”) names(idh);attach(idh)
#El comando “hist” genera en el eje de las Y las frecuencias #y las cajas est?n en intervalos iguales. En lugar de frecuencias #alquien usa densidad, la suma de todas las cajas es uno y cada caja #es una probabilidad.
par(mfrow=c(1,1)) hist(IDH2015) hist(IDH2016)
#Vamos a generar el histograma con 50 quiebres y para que #se genere la densidad.
hist(IDH2016,breaks=40, freq=T)
p<-na.omit(IDH2016) hist(p,breaks=40, freq=T, main=“Histograma del IDH 2016”, xlim=c(0,2), ylim=c(0,6),col=“yellow”,border=F) curve(dnorm(x,mean=mean(p), sd=sd(p)), add = TRUE,col=“red”) rug(jitter(IDH2016)) #mejor densidad de los datos
#Curva de la densidad de Kernel
d<-density(IDH2016, na.rm=TRUE) plot(d, main=“Densidad de Kernel del IDH 2016”) polygon(d,col=“ivory2”,border=“blue”)
e<-density(IDH2015, na.rm=TRUE) plot(e, main=“Densidad de Kernel del IDH 2015”) polygon(e, col=“white”,border=“red”)
#BOXPLOTS (se usa el paquete base)
par(mfrow=c(1,1)) boxplot(IDH2015) boxplot(IDH2016) boxplot(IDH2016, plot=F)
#vector con cinco valores: el primer valor extremo, #el primer cuartil, el segundo y tercer cuartil y el ?ltimo #valor que cae en los bigotes. El n es igual a 169, #No hay outliers.
boxplot(IDH2016, notch=T)
par(mfrow=c(1,3)) boxplot(IDH2016,ylim=c(0,1),col=“green”,range=1.5) boxplot(IDH2016,ylim=c(0,1),col=“red”,range=1) boxplot(IDH2016,ylim=c(0,1),col=“yellow”,range=0.5)
#Los graficamos en la misma escala. De forma predeterminada #son 1.5 veces del rango intercuartilico.
par(mfrow=c(2,1),mar=c(3,5,3,5)) hist(IDH2016,breaks=40, freq=F, main=“Histograma del IDH 2016”, xlim=c(0,1), ylim=c(0,5),col=“blue”,border=F ) boxplot(IDH2016,ylim=c(0,1),col=“blue”,horizontal=T)
par(mfrow=c(1,1) boxplot(IDH2016,IDH2015,ylim=c(0,1),col=c(“lightblue”,“gray”),horizontal=F, main=“Boxplots del IDH”, names= c(“IDH 2015”,“IDH 2016”), notch=T)
#Boxplots Usando otros paquetes
install.packages(“gplots”) library(gplots) boxplot(IDH2016)
install.packages(“Rlab”) library(Rlab) bplot(IDH2015, ylim=c(0,1))
install.packages(“ggplot2”) library(ggplot2)
nume<-read.csv(“PAISES.csv”) names(nume) str(nume)
library(Hmisc) describe(nume$Elrec)
#Histrograma Usando comando qplot par(mfrow=c(1,2)) qplot(IPC, data=nume, geom=“histogram”)
qplot(IDH, data=nume, geom=“histogram”)
#Histograma Usando comando ggplot p <- ggplot(nume, aes(x=IDH)) p + geom_histogram() p + geom_histogram(bindwith=0.1) p + geom_histogram() + geom_density() p + stat_bin(geom=“area”) p + stat_bin(geom=“point”) p + stat_bin(geom=“line”) p + geom_histogram(aes(fill =EL)) p + geom_histogram(aes(fill = ..count.. )) + scale_fill_gradient(“Count”, low = “green”, high = “red”) p + geom_histogram(aes(y = ..density..)) p + geom_histogram(aes(y = ..density..)) + geom_density()
p + geom_histogram(colour=“darkgreen”, fill = “white”, bindwith = 0.5)
#Faceting (mostrando histogramas por categor?as) p <- p + geom_histogram(binwidth = 0.1) #Definiendo p + facet_grid(Elrec ~ IPCrec) #muestra el histograma por categor?as de 2 variables p + facet_grid(~ IPCrec) #muestra el histograma por cada categor?a de 1 variable
#Histograma de IDH rellenado por Estado de Libertad q <- ggplot(nume, aes(x=IDH, fill=EL)) q + geom_bar() # defaults to stacking es para variables categóricas q + geom_bar(position=“fill”) q + geom_bar(position=“dodge”)
#Gr?fico de densidad (sin definir) ggplot(nume, aes(IPC, ..density.., colour = EL)) + geom_freqpoly(binwidth = 0.5) ggplot(nume, aes(IPC, colour = EL)) + geom_density() ggplot(nume, aes(IPC, fill = EL)) + geom_density(alpha = 0.2)
#Gr?fico de barras a<- ggplot(nume, aes(factor(EL))) a + geom_bar(fill=“darkblue”, colour=“gray”) a + geom_bar(fill=“darkblue”, colour=“gray”) + coord_flip()
#Gráfico de barras rellenado por otras categóricas qplot(factor(EL), data=nume, geom=“bar”, fill=factor(IPC2)) qplot(factor(IPC2), data=nume, geom=“bar”, fill=factor(EL))
#Boxplots h <- ggplot(nume, aes(factor(EL), IDH)) h + geom_boxplot() h + geom_boxplot() + geom_jitter() h + geom_boxplot() + coord_flip() h + geom_boxplot(outlier.colour = “green”, outlier.size = 3) h + geom_boxplot(aes(fill = factor(EL))) h + geom_boxplot(aes(fill = factor(IPC2)))
########******************************************####################