On Jan 11, 2005, at 7:55 AM, Dimitris Rizopoulos wrote: > mat <- rnorm(9); dim(mat) <- c(3,3) > mat > ########### > cbind(i=rep(1:nrow(mat), each=ncol(mat)), > j=rep(1:ncol(mat), nrow(mat)), value=c(t(mat))) > That will do it. Thanks! Sean