Re: [R] why NA coefficients

2011-11-08 Thread William Dunlap
From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On > Behalf Of array chip > Sent: Tuesday, November 08, 2011 10:57 AM > To: Dennis Murphy > Cc: r-help@r-project.org > Subject: Re: [R] why NA coefficients > > Hi Dennis, Thanks very much for the details. Al

Re: [R] why NA coefficients

2011-11-08 Thread array chip
and should not be used? Thanks John From: Dennis Murphy Cc: "r-help@r-project.org" Sent: Tuesday, November 8, 2011 10:22 AM Subject: Re: [R] why NA coefficients The cell mean mu_{12} is non-estimable because it has no data in the cell. How can yo

Re: [R] why NA coefficients

2011-11-08 Thread Dennis Murphy
Because "treat 1_group 2" has no cases, so why it is still estimable while > on the contrary, "treat 7_group 2" which has 2 cases is not? > Thanks > John > > > > > > From: Dennis Murphy > To: array chip > Sent: Monday, November

Re: [R] why NA coefficients

2011-11-08 Thread David Winsemius
vid. Thanks John From: David Winsemius To: David Winsemius Cc: array chip ; "r-help@r-project.org" > Sent: Monday, November 7, 2011 10:27 PM Subject: Re: [R] why NA coefficients But this output suggests there may be alligators in the swamp: > predict(lmod, newdata=data.

Re: [R] why NA coefficients

2011-11-08 Thread array chip
true, why it has to omit "treat 7-group 2" Thanks again From: David Winsemius Cc: "r-help@r-project.org" Sent: Monday, November 7, 2011 10:19 PM Subject: Re: [R] why NA coefficients On Nov 7, 2011, at 10:07 PM, array chip wrote:

Re: [R] why NA coefficients

2011-11-08 Thread array chip
Sure it does, but still struggling with what is going on... Thanks John From:David Winsemius To:David Winsemius oject.org> Sent:Monday, November 7, 2011 10:27 PM Subject:Re: [R] why NA coefficients But this output suggests there may be alligators in

Re: [R] why NA coefficients

2011-11-07 Thread David Winsemius
nday, November 7, 2011 5:13 PM Subject: Re: [R] why NA coefficients On Nov 7, 2011, at 7:33 PM, array chip wrote: > Hi, I am trying to run ANOVA with an interaction term on 2 factors (treat has 7 levels, group has 2 levels). I found the coefficient for the last interaction term is al

Re: [R] why NA coefficients

2011-11-07 Thread David Winsemius
reat 7-group 2" is not estimable? Well, it had to omit one of them didn't it? (But I don't know why that level was chosen.) -- David. Thanks John From: David Winsemius To: array chip Cc: "r-help@r-project.org" Sent: Monday, November 7, 2011 5:13 PM Subject: Re: [

Re: [R] why NA coefficients

2011-11-07 Thread array chip
hich has 2 cases is not? Thanks John From: Dennis Murphy Sent: Monday, November 7, 2011 9:29 PM Subject: Re: [R] why NA coefficients Hi John: What is the estimate of the cell mean \mu_{12}? Which model effects involve that cell mean? With this data ar

Re: [R] why NA coefficients

2011-11-07 Thread array chip
imable?   Thanks John From: David Winsemius Cc: "r-help@r-project.org" Sent: Monday, November 7, 2011 5:13 PM Subject: Re: [R] why NA coefficients On Nov 7, 2011, at 7:33 PM, array chip wrote: > Hi, I am trying to run ANOVA with an interaction term on 2 factors

Re: [R] why NA coefficients

2011-11-07 Thread David Winsemius
On Nov 7, 2011, at 7:33 PM, array chip wrote: Hi, I am trying to run ANOVA with an interaction term on 2 factors (treat has 7 levels, group has 2 levels). I found the coefficient for the last interaction term is always 0, see attached dataset and the code below: test<-read.table("test.t

[R] why NA coefficients

2011-11-07 Thread array chip
Hi, I am trying to run ANOVA with an interaction term on 2 factors (treat has 7 levels, group has 2 levels). I found the coefficient for the last interaction term is always 0, see attached dataset and the code below: > test<-read.table("test.txt",sep='\t',header=T,row.names=NULL) > lm(y~factor(t