https://github.com/cran/robCompositions
Tip revision: e168b07923b144993e99ea7b65180aed9fe3efa8 authored by Matthias Templ on 11 December 2009, 00:00:00 UTC
version 1.3.3
version 1.3.3
Tip revision: e168b07
compLW.R
x <- read.table("lw09.csv", sep=";", dec=".", header=TRUE)
# nur Gemeinden:
gem <- x[substr(x[,"Nr."],5,5) != 0,]
head(gem)
# nur Anteile:
nchar(colnames(x))
w <- which(substr(colnames(gem),nchar(colnames(gem)), nchar(colnames(gem))) == "p")
anteile <- gem[, c(2,w)]
head(anteile)
# biplots:
library(compositions)
library(robCompositions)
library(robustbase)
# standard biplot for raw data
x <- anteile[,5:8]
rownames(x) <- anteile[,1]
## Problem mit rownames noch nicht geloest...:
rn <- ifelse(duplicated(substr(rownames(x),1,3)), paste(substr(rownames(x),1,3),1, sep=""), substr(rownames(x),1,3))
j=0
for(i in 1:nrow(x)){
TR <- ifelse(substr(rn[i], 4,4) == "1", TRUE, FALSE)
if(TR) j <- j + 1
if(j == 9) j <- 3
rn[i] <- ifelse(substr(rn[i], 4,4) == "1", paste(substr(rn[i],1,3), j, sep=""), rn[i])
}
rownames(x) <- rn
##
pca <- princomp(x)
summary(pca)
biplot(pca,main="Election - standard biplot",cex=1)
# robust biplot for raw data
x11()
pcarob <- princomp(x,covmat=covMcd(x))
summary(pcarob)
biplot(pcarob,main="Election - robust biplot",cex=1)
#Standard compositional biplot
x11()
x <- acomp(x)
y <- clr(x)
cpca <- princomp(x)
summary(cpca)
biplot(cpca,main="Election - standard CoDa biplot",cex=1)
# robust compositional biplot
x11()
pdf("ooeWahl.pdf")
pca <- pcaCoDa(x)
summary(princomp(ilr(x),covmat <- covMcd(ilr(x))))
rownames(pca$scores) <- rownames(x)
biplot(-pca$scores[,1:2],-pca$loadings[,1:2],main="Election - robust CoDa biplot",cex=1)
dev.off()