Re: [R] lme nesting/interaction advice

2008-05-12 Thread Federico Calboli
On 11 May 2008, at 22:45, Andrew Robinson wrote: lme(y ~ selection * males, random = ~1|replica/selection/males, mydata) forgive me, but I seem to see nesting in the random statement. That is what happens when we separate factors with a '/'; they are nested. We would expect that

Re: [R] lme nesting/interaction advice

2008-05-12 Thread Federico Calboli
On 12 May 2008, at 01:05, Andrew Robinson wrote: On Mon, May 12, 2008 at 10:34:40AM +1200, Rolf Turner wrote: On 12/05/2008, at 9:45 AM, Andrew Robinson wrote: On Sun, May 11, 2008 at 07:52:50PM +0100, Federico Calboli wrote: The main point of my question is, having a 3 way anova (or

Re: [R] lme nesting/interaction advice

2008-05-12 Thread Federico Calboli
On 12 May 2008, at 09:29, Dieter Menne wrote: Federico: First, mixed models are different from standard 101 Anova, and quite a lot of the nesting stuff I used to ponder about 30 year ago when I started teaching this is no longer relevant and works implicitely when you code the parameters

Re: [R] lme nesting/interaction advice

2008-05-12 Thread Federico Calboli
On 11 May 2008, at 23:34, Rolf Turner wrote: It doesn't seem to me to be a complaint as such. It is a request for insight. I too would like some insight as to what on earth is going on. And why do you say Federico shows no evidence of having searched the

Re: [R] lme nesting/interaction advice

2008-05-12 Thread Andrew Robinson
On Mon, May 12, 2008 at 10:50:03AM +0100, Federico Calboli wrote: On 12 May 2008, at 01:05, Andrew Robinson wrote: On Mon, May 12, 2008 at 10:34:40AM +1200, Rolf Turner wrote: On 12/05/2008, at 9:45 AM, Andrew Robinson wrote: On Sun, May 11, 2008 at 07:52:50PM +0100, Federico Calboli

Re: [R] lme nesting/interaction advice

2008-05-12 Thread Federico Calboli
On 12 May 2008, at 11:16, Andrew Robinson wrote: Well. I have documentation relevant to nlme that goes back about 10 years. I don't know when it was first added to S-plus, but I assume that it was about then. Now, do you think that if the thing that you want to do was really bog standard,

Re: [R] lme nesting/interaction advice

2008-05-11 Thread Federico Calboli
On 10 May 2008, at 07:36, Kingsford Jones wrote: Federico, I think you'll be more likely to receive the type of response you're looking for if you formulate your question more clearly. The inclusion of commented, minimal, self-contained, reproducible code (as is requested at the bottom of

Re: [R] lme nesting/interaction advice

2008-05-11 Thread Andrew Robinson
On Sun, May 11, 2008 at 07:52:50PM +0100, Federico Calboli wrote: The main point of my question is, having a 3 way anova (or ancova, if you prefer), with *no* nesting, 2 fixed effects and 1 random effect, why is it so boneheaded difficult to specify a bog standard fully crossed model?

Re: [R] lme nesting/interaction advice

2008-05-11 Thread Rolf Turner
On 12/05/2008, at 9:45 AM, Andrew Robinson wrote: On Sun, May 11, 2008 at 07:52:50PM +0100, Federico Calboli wrote: The main point of my question is, having a 3 way anova (or ancova, if you prefer), with *no* nesting, 2 fixed effects and 1 random effect, why is it so boneheaded difficult to

Re: [R] lme nesting/interaction advice

2008-05-11 Thread Andrew Robinson
On Mon, May 12, 2008 at 10:34:40AM +1200, Rolf Turner wrote: On 12/05/2008, at 9:45 AM, Andrew Robinson wrote: On Sun, May 11, 2008 at 07:52:50PM +0100, Federico Calboli wrote: The main point of my question is, having a 3 way anova (or ancova, if you prefer), with *no* nesting, 2 fixed

Re: [R] lme nesting/interaction advice

2008-05-10 Thread Kingsford Jones
On Fri, May 9, 2008 at 4:04 AM, Federico Calboli [EMAIL PROTECTED] wrote: Note that random can be a list: a one-sided formula of the form ~x1+...+xn, or a pdMat object with a formula (i.e. a non-NULL value for formula(object)), or a list of such formulas or pdMat objects. If you can

Re: [R] lme nesting/interaction advice

2008-05-09 Thread Dieter Menne
Federico Calboli f.calboli at imperial.ac.uk writes: Hi everyone, I am confused on how to specify some nesting and interaction terma with lme(). lme(y ~ selection * males, random = ~1|replica/selection/males, mydata) Note that random can be a list: a one-sided formula of the form

Re: [R] lme nesting/interaction advice

2008-05-09 Thread Federico Calboli
Note that random can be a list: a one-sided formula of the form ~x1+...+xn, or a pdMat object with a formula (i.e. a non-NULL value for formula(object)), or a list of such formulas or pdMat objects. If you can translate that into *informative* English I'd be grateful. I have the Pinheiro

[R] lme nesting/interaction advice

2008-05-08 Thread Federico Calboli
Hi everyone, I am confused on how to specify some nesting and interaction terma with lme(). I have a dataset where some flies where selected for accessory gland size, made to mate in presence/absence of another male and the level of some protein measured. Now the complex stuff. The