search for: fjohannes

Dear all, I have three concerns: 1) I am running models with the lme4 package. I cannot find a way to pull out a vector of the fitted values and the residuals. Does anybody know how to do it? 2) How can I nest a random effect variable into a "two-level" fixed effect variable? 3) Suppose I have the following model: y = a + b|c + d + error, where 'a' is a fixed effect, 'c'

Dear all, I am analyzing some data that requires a mixed model. I have been reading Pinheiro and Bates' book, but cannot find the notation to fit the following model: Suppose I have the dataset below. Here I am fitting variable p as a fixed effect, variable h as a random effect and variable t as nested within h. I would like to include the variable j as well as an independent (non-nested)

I have a large matrix and want to create a 3d surface of it. Suppose the matrix looks something like Matrix K: a<-c(1:1200) b<-c(rep(1:30,40)) c<-c(a+b^2) K<-data.frame(a,b,c) The vector values are not ordered (and repeat themselves as in b). Whenever I try commands like "persp" it wants ordered values. The only command that approximates what I

HI, Suppose I have the following data structure. LRT tp 1 1.50654010 522 2 0.51793929 522 3 0.90340299 522 4 1.20293325 522 5 1.05578774 523 6 0.01617942 523 7 0.68183543 523 8 0.43820244 523 9 1.14123995 524 10 0.05809550 524 11 0.93061597 524 12 1.39739700 524 13 1.05220953 525 14 0.03471461 525 15 0.63168798 525 16 1.40592603 525 17 1.41884492 526 18