Displaying 2 results from an estimated 2 matches for "rand_mat".
Did you mean:
rand_max
2007 May 25
0
Help with complex lme model fit
...GenStat...
The raw data are available on the web at
http://cmbeale.freehostia.com/OutData.txt and can be read directly
into R using:
gpdat <- read.table("http://cmbeale.freehostia.com/OutData.txt",
header = T)
gpdat$X7 <- as.factor(gpdat$X7)
gpdat$X4 <- as.factor(gpdat$X4)
rand_mat <- as.matrix(gpdat[,11:26])
gpdat <- groupedData(Y1 ~X1 + X2 + X3 + X4 + X5 + m_sum|.g, data =
gpdat)
the model fitted using:
library(Matrix)
library(nlme)
m_sum <- rowSums(gpdat[,11:27])
mod1 <- lme(fixed = Y1 ~ X1 + X2 + X3 + X4 + X5 + m_sum,...
2007 Jan 26
0
R crash with modified lmer code
...X1 = rnorm(400), X2 = rnorm(400),
F1 = as.factor(sample(1:4, 400, replace = T)))
forSm <- Matrix(c(runif(1200)), ncol = 3, sparse = T)
matDim <- dim(forSm)
smoother <- as.factor(sample(1:4, matDim[1], replace = T))
test <- lmer1 (Y ~ X1 + X2 + (1|F1) + (1|smoother), data = dat,
rand_mat = forSm)
the sessionInfo (called after sourcing the code but before the crash)
is:
R version 2.4.1 (2006-12-18)
i386-pc-mingw32
locale:
LC_COLLATE=English_United Kingdom.1252;LC_CTYPE=English_United
Kingdom.1252;LC_MONETARY=English_United
Kingdom.1252;LC_NUMERIC=C;LC_TIME=English_United Kingd...