[R] pairwise combinatons of variables
Christos Hatzis
christos at silicoinsights.com
Sat Mar 25 20:26:25 CET 2006
You can try 'outer'.
Here is a way to get the index combinations that you want in a character
vector:
as.vector(outer(1:4, 1:4, paste, sep=","))
-Christos
-----Original Message-----
From: r-help-bounces at stat.math.ethz.ch
[mailto:r-help-bounces at stat.math.ethz.ch] On Behalf Of Mark Hempelmann
Sent: Saturday, March 25, 2006 12:50 PM
To: r-help at stat.math.ethz.ch
Subject: [R] pairwise combinatons of variables
Dear WizaRds,
although this might be a trivial question to the community, I was
unable to find anything solving my problem in the help files on CRAN. Please
help.
Suppose I have 4 variables and want to use all possible combinations:
1,2
1,3
1,4
2,3
2,4
3,4
for a further kmeans partitioning.
I tried permutations() of package e1071, but this is not what I need. Thank
you for your help and support.
mark.
---------------------------------------------------------------------
Additionally: For anybody who is willing to offer some advise, here is my
complete approach:
library(e1071)
mat <- matrix( c(6,7,8,2,3,4,12,14,14, 14,15,13,3,1,2,3,4,2,
15,3,10,5,11,7,13,6,1, 15,4,10,6,12,8,12,7,1), ncol=9, byrow=T )
rownames(mat) <- paste("v", 1:4, sep="" )
tmat <- t(mat)
cluster <- c(1, 2, 1, 3, 3, 3, 1, 2, 2)
centroids <- matrix( 0, ncol=3, nrow=4 )
obj <- vector(mode="list", length=3)
for (j in 1:4){
for (i in 1:3){
where <- which(cluster==i) # which obj is in which class?
centroids[j,i] <- mean( tmat[ where,j ] )
obj[[i]] <- tmat[where,]
}
}
colnames(centroids) <- paste( c("Cluster"), 1:3)
rownames(centroids) <- rownames(mat)
centroids
obj
## now I want to do kmeans of all possible variable pairs, e.g. v1 and
v3
## automization in a second step later
wjk <- kmeans(tmat[,c(1,3)], centers=centroids[c(1,3),], iter.max=10,
algorithm="MacQueen") ## obviously wrong
______________________________________________
R-help at stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide!
http://www.R-project.org/posting-guide.html
More information about the R-help
mailing list