library(bigmemory)
nrows <- 50000
ncols <- 15
bkFile <- "bigmat4.bk"
descFile <- "bigmatk4.desc"
bigmat <- filebacked.big.matrix(nrow=nrows, ncol=ncols, type="double",
backingfile=bkFile, backingpath=".",
descriptorfile=descFile,
dimnames=c(NULL,NULL))
# Each column value with be the column number multiplied by
# samples from a standard normal distribution.
set.seed(123)
for (i in 1:ncols) bigmat[,i] = rnorm(nrows)*i
y <- rnorm(nrows) + bigmat[,1] - bigmat[,2]
system.time(lmr1 <- bigLm(bigmat, y))
summary(lmr1)
Run the code above in your browser using DataLab