Displaying 3 results from an estimated 3 matches for "factora2".
Did you mean:
factor2
2011 Jan 21
1
TRADUCING lmer() syntax into lme()
...r Rsociety,
I'd like to kingly ask to anyone is willing to answer me how to implement a
NON NESTED random effects structure in lme()
In particular I've tried the following translation from lmer to lme, as
suggested from some web example
mod1<-lmer(y~x*z+(x*z|factorA1/factorB)+(x*z|factorA2/factorB)) # y,x,z
continuous
mod2<-lme(y~x*z, random= pdBlocked(list(pdIdent(~1|factorA1/factorB
),pdIdent(~1|factorA2/factorB))))
In detail check how I've tried to state in mod1 that Iwant to evaluate
randomness in the interaction x*z (i.e intercept, slope, interaction)
grouped by by a...
2008 Mar 25
0
Behaviour of interactions in glm
...ctorA+FactorB+FactorA:FactorB, FactorA:FactorB -
but whatever I do I always get 6 coefficients.
Really what I would like to do is to ask for FactorA*FactorB less the
entries in the design matrix that I get from FactorA and FactorB. This
would leave me with the design matrix being:
Record Mean FactorA2:FactorB1 FactorA2:FactorB3
1 1 0 0
2 1 0 0
3 1 0 0
4 1 1 0
5 1 0 0
6 1 0 1
If anyone has any advice on how I could make this happen, I'd...
2011 Jan 22
0
how to call BayesX in R to see the graph
...r Rsociety,
I'd like to kingly ask to anyone is willing to answer me how to implement a
NON NESTED random effects structure in lme()
In particular I've tried the following translation from lmer to lme, as
suggested from some web example
mod1<-lmer(y~x*z+(x*z|factorA1/factorB)+(x*z|factorA2/factorB)) # y,x,z
continuous
mod2<-lme(y~x*z, random= pdBlocked(list(pdIdent(~1|factorA1/factorB
),pdIdent(~1|factorA2/factorB))))
In detail check how I've tried to state in mod1 that Iwant to evaluate
randomness in the interaction x*z (i.e intercept, slope, interaction)
grouped by by a...