Program Code

Attention! For repeating the computation a transformed dataset is needed! If you have not yet computated and saved the transformation, run the program for transformation on the wikipage en:Analysis of Tuberculosis first! Additionally, after copying this program code to your XploRe input window, the arguments of "graxes(argument)" have to be rewritten into one line. Otherwise, one gets an error message like "One parenthesis not closed".

library("xplore")
library("stats")

; ----- Reading Data --------------------------------------------------------------------------

choose = "Read data from:"

defaults = "C:\Dokumente und Einstellungen\All Users\Desktop\UN_data_ordered.csv"

v1 = readvalue(choose, defaults)

x = readcsvm(v1)

data = x.double
country = x.text

; ----- Standardization to 0-1-Scale ----------------------------------------------------------

dats = (data - min(data))/(max(data)-min(data))

; ----- Create Colorvector for Variables -------------------------------------------------------

randomize(25)
bunt = uniform(10, 3)*255
createcolor(bunt)

; ----- Generate, Rotate and Move Boxplots -----------------------------------------------------

gr8 = grbox(data[,8], bunt[2])
gr9 = grbox(data[,9], bunt[3])
gr10 = grbox(data[,10], bunt[4])
gr11 = grbox(data[,11], bunt[5])
gr12 = grbox(data[,12], bunt[6])

gr8 = grrot(gr8, 1)
gr9 = grrot(gr9, 1)
gr9 = grmove(gr9, #(1.5, 0))
gr10 = grrot(gr10, 1)
gr10 = grmove(gr10, #(3, 0))
gr11 = grrot(gr11, 1)
gr11 = grmove(gr11, #(4.5, 0))
gr12 = grrot(gr12, 1)
gr12 = grmove(gr12, #(6, 0))

gr7 = grbox(dats[,7], bunt[1])
gr13 = grbox(dats[,13], bunt[7])
gr14 = grbox(dats[,14], bunt[8])
gr15 = grbox(dats[,15], bunt[9])
gr16 = grbox(dats[,16], bunt[10])

gr7 = grrot(gr7, 1)
gr13 = grrot(gr13, 1)
gr13 = grmove(gr13, #(1.5, 0))
gr14 = grrot(gr14, 1)
gr14 = grmove(gr14, #(3, 0))
gr15 = grrot(gr15, 1)
gr15 = grmove(gr15, #(4.5, 0))
gr16 = grrot(gr16, 1)
gr16 = grmove(gr16, #(6, 0))

; ----- Graphical Settings --------------------------------------------------------------------

; ---------- Generate New Axes and Descriptions -----------

axesoff()

xtext1 = "Drugs"|"Education"|"Literacy"|"Sanitation"|"Water"
xtext2 = "Condom"|"CO2"|"Internet"|"PC"|"Telephone"

axes1 = graxes((-0.5|7.5)~(-5|115), "origin", 7.5, "ytextpos", 9, "xtextpos", 6,
               "xticks", 0.5|2|3.5|5|6.5, "xtext", xtext1,
               "yticks", 0|25|50|75|100, "ytext", "0%"|"25%"|"50%"|"75%"|"100%",
               "xtextsize", 16, "ytextsize", 16)

axes2 = graxes((-0.5|7.5)~(-0.1|1.1), "origin", 7.5, "ytextpos", 9, "xtextpos", 6,
               "xticks", 0.5|2|3.5|5|6.5, "xtext", xtext2,
               "yticks", 0|0.5|1, "ytext", "0"|"standardized"|"1",
               "xtextsize", 16, "ytextsize", 16)

; ---------- Label Certain Data Points -----------

ptt1 = "37 countries"
pt1 = #(0.5,50)'
setmaskt(pt1, ptt1, bunt[2], 6, 16)

ptt2 = max(data[,7|13|14|15|16])'
pt2 = (0.5|2|3.5|5|6.5)~matrix(5,1)
fmt2 = "%8.5g"
ptt2 = string(fmt2, ptt2)
setmaskt(pt2, ptt2, bunt[1|7:10], 3, 16)
setmaskp(pt2, 0, 0, 0)

; ------------------------------------------------

setsize(900, 600)
d2 = createdisplay(2, 1)

show(d2, 1, 1, pt1, gr8, gr9, gr10, gr11, gr12, axes1)
show(d2, 2, 1, pt2, gr7, gr13, gr14, gr15, gr16, axes2)

setgopt(d2, 1, 1, "title", "Boxplots", "yoffset", 20|5)
setgopt(d2, 2, 1, "ylabel", "standardized")