I found a solution myself, but thanks for the answers. I solved it like this: D <- matrix(1:225*100,nrow=100,ncol=225) for(i in 1:225) D[,i] <- rt(100,df=225) end -- View this message in context: http://r.789695.n4.nabble.com/Storing-datasets-tp4632874p4633071.html Sent from the R help mailing list archive at Nabble.com.