search for: xvar2

Hi all, I try to make a split-plot with poisson errors using glmmPQL, but I have some doubts about the model simplification. Look my system: Block = 3 blocks Xvar1 = 2 levels Xvar2 = 13 levels Yvar = Count data Response I need know about the behaviour of Var1, Var2 and interaction Var1:Var2. Look the levels: > levels(Xvar1) [1] "A" "B" > levels(Xvar2) [1] "a" "b" "c" "d" "e" "f" "g&...

...(0,length(data.mult))), intercept3=c(rep(0,length(data.mult)),rep(0,length(data.mult)),rep(0,length(data.mult)),rep(1,length(data.mult))), xvar1=c(rep(0,length(data.mult)),rep(1,length(data.mult)),rep(0,length(data.mult)),rep(0,length(data.mult)))*rep(xvar,4), xvar2=c(rep(0,length(data.mult)),rep(0,length(data.mult)),rep(1,length(data.mult)),rep(0,length(data.mult)))*rep(xvar,4), xvar3=c(rep(0,length(data.mult)),rep(0,length(data.mult)),rep(0,length(data.mult)),rep(1,length(data.mult)))*rep(xvar,4), stratum=rep(1:length(data.mult),4))...

...efficients, digits=digits), ", ", sep="", collapse="")) cat(", \\\\ \n \\text{where }", substr(est, 1, (nchar(est)-2)), file=file) } # Write the closing math mode tag cat(math.env[2], "\n", file=file) } # END latex.lm Xvar1 <- rnorm(20) Xvar2 <- rnorm(20) Xvar3 <- factor(rep(c("A","B"),10)) Y.var <- rnorm(20) D <- data.frame(Xvar1, Xvar2, Xvar3, Y.var) x1 <- lm(Y.var ~ pol(Xvar1, 3) + Xvar2*Xvar3, data=D) latex.lm(x1) -- Ista Zahn Graduate student University of Rochester Department of Clinical and S...