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!

The code can easily be adjusted manually to plot density estimates of other variables (choose i and j).

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

; ----- Calculating the 3-Dimensional Grid------------------------------------------------------

i = 6                                                      ; choose variables 6 and 11 for the
j = 11                                                     ; 2-dimensional density estimation

o = min(data[,i|j])   ; origin of the grid

l = max(data[,i|j]) - o                                    ; length of the gridsides

n = 30*#(1,1)   ; number of gridpoints

w = l'./n   ; stepwidth

gridx = grid (o', w, n)

; ----- Compute Density Estimator --------------------------------------------------------------

h = 5*#(1,1)   ; bandwith of kernel

d2 = sker (sort(data[,i|j]), h', "qua", matrix(rows(data)), sort(gridx))

; ----- Choose Settings for Density Estimation -------------------------------------------------

cl=grcolorscheme("color",15)
createcolor(cl)

gr = grsurfacecol(sort(gridx)~d2, cl, 15)

d = createdisplay(1,1)

show(d,1,1,gr)
setmask(gr, "line", "blue")
setgopt(d, 1, 1, "title", "Surfaceplot of 'Tuberculosis' and 'Sanitation'")