library(MASS)
load("mult.RData")
## or load from the web by
load("http://www.stats.ox.ac.uk/~ripley/MultAnal_HT2008/mult.RData")
## or browse for it by
load(choose.files())

ft2 <- ft1$FTSc
names(ft2) <- row.names(ft1)
FTScore <- sort(ft2)
nm <- names(FTScore); nm <- factor(nm, levels=nm)

windows(pointsize=7)
library(lattice)
dotplot(nm ~ FTScore, xlab="Financial Times overall Score", xlim=c(0, 100), cex=0.7)

Times$type <- c(rep(1, 5), rep(2, 51), rep(3, 40)); Times$type[52] <- 3

library(rggobi)
gg <- ggobi(Times[1:8])
d <- displays(gg)[[1]]
pmode(d) <- "2D Tour"
variables(d) <- list(X=0:7) # 0-based
glyph_colour(gg[[1]]) <- Times$type
colorscheme(gg) <- "Paired 4"


FTtype <- c(rep(1, 5), rep(2, 51), rep(3, 41))
gg$FT <- ft
glyph_colour(gg$FT) <- FTtype
d2 <- display(gg$FT, "2D Tour")
variables(d2) <- list(X=0:13) # 0-based

close(gg)

T2 <- Times[1:8]
T2[] <- lapply(T2, function(x) (x-median(x))/IQR(x))
biplot(princomp(Times[,1:8]),cex=0.7); title("Times: Raw")
biplot(princomp(T2), expand=0.8, cex=0.7); title("Times: IQR scaled")

Times.sam <- sammon(dist(T2))
eqscplot(Times.sam$points, type="n", xlab="", main="Sammon: Times, IQR")
text(Times.sam$points, row.names(Times), cex=0.7)
Times.iso <- isoMDS(dist(T2))
eqscplot(Times.iso$points, type="n", xlab="", main="isoMDS: Times, IQR")
text(Times.iso$points, row.names(Times), cex=0.7)