################################################################################### # # Empirical distribution funcion of a discrete variable # # Last review: 17/06/2018 # email: ja@fct.unl.pt # ################################################################################### ################################################################################### varlistindex <- 4 sortbyfreq <- FALSE colourbar <- rgb(1,0,0,0.5) colourline <- rgb(1,0,0,0.5) ################################################################################### par(mfrow=c(1,2),mar=c(3.5, 3.5, 2, 1), mgp=c(2.25, 0.8, 0)) # # Duplicate gtd -> gtdnew # gtdnew <- data.frame(gtd[,varlistindex]) names(gtdnew)[1]=names(gtd)[varlistindex] # # Calculations and graphic # freqabs <- table(gtdnew[,1]) if (sortbyfreq == TRUE) { freqabs <- sort(freqabs, decreasing = TRUE) } barplot(freqabs, col=colourbar, xlab=names(gtdnew)[1], ylab="f(x)", border="black", main = paste("Função de probabilidade empírica"), cex.names=0.8, cex.lab=1.0, cex.axis=1.0, cex.main=0.90, cex.sub=0.5) plot (cumsum(freqabs), type="s", col=colourline, lwd=2, xlab=names(gtdnew)[1], ylab="F(x)", main = "Função pseudo cumulativa empírica", cex.lab=1.0, cex.axis=1.0, cex.main=0.90, cex.sub=0.5)