Re: [R] Post-hoc tests in MASS using glm.nb

2011-05-17 Thread Bill.Venables
you can mostly use data= arguments on the individual functions. Bill Venables. -Original Message- From: Venables, Bill (CMIS, Dutton Park) Sent: Wednesday, 18 May 2011 9:07 AM To: 'Bert Gunter'; 'Peter Ehlers' Cc: 'R list' Subject: RE: [R] Post-hoc tests in MASS

Re: [R] Post-hoc tests in MASS using glm.nb

2011-05-17 Thread Bill.Venables
03 1.72e-01 1.50e-02 5.60e-06 ... $ x2 : num 0.268811 0.034224 0.309802 0.060844 0.000315 ... > Bill Venables. -Original Message- From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On Behalf Of Bert Gunter Sent: Wednesday, 18 May 2011 12:08 AM To: Peter Ehle

Re: [R] Post-hoc tests in MASS using glm.nb

2011-05-17 Thread Bert Gunter
Folks: > Only if the user hasn't yet been introduced to the with() function, > which is linked to on the ?attach page. > > Note also this sentence from the ?attach page: >  " attach can lead to confusion." > > I can't remember the last time I needed attach(). > > Peter Ehlers Yes. But perhaps

Re: [R] Post-hoc tests in MASS using glm.nb

2011-05-17 Thread Peter Ehlers
On 2011-05-17 02:22, Timothy Bates wrote: Dear Bryony: the suggestion was not to change the name of the data object, but to explicitly tell glm.nb what dataset it should look in to find the variables you mention in the formula. so the salient difference is: m1<- glm.nb(Cells ~ Cryogel*Day, da

Re: [R] Post-hoc tests in MASS using glm.nb

2011-05-17 Thread David Winsemius
011 03:21 To: Bryony Tolhurst; r-help@r-project.org Subject: RE: [R] Post-hoc tests in MASS using glm.nb ?relevel Also, you might want to fit the models as follows Model1 <- glm.nb(Cells ~ Cryogel*Day, data = myData) myData2 <- within(myData, Cryogel <- relevel(Cryogel, ref = "2"

Re: [R] Post-hoc tests in MASS using glm.nb

2011-05-17 Thread Timothy Bates
h(data) > > But obviously I am still using the attach function, if not with 'data'!! > > Thanks again > > Bryony Tolhurst > > -Original Message- > From: bill.venab...@csiro.au [mailto:bill.venab...@csiro.au] > Sent: 17 May 2011 03:21 > To: Bryony Tolhurst

Re: [R] Post-hoc tests in MASS using glm.nb

2011-05-17 Thread Bryony Tolhurst
!! Thanks again Bryony Tolhurst -Original Message- From: bill.venab...@csiro.au [mailto:bill.venab...@csiro.au] Sent: 17 May 2011 03:21 To: Bryony Tolhurst; r-help@r-project.org Subject: RE: [R] Post-hoc tests in MASS using glm.nb ?relevel Also, you might want to fit the models as foll

Re: [R] Post-hoc tests in MASS using glm.nb

2011-05-16 Thread Bill.Venables
elp-boun...@r-project.org] On Behalf Of bryony Sent: Tuesday, 17 May 2011 3:46 AM To: r-help@r-project.org Subject: [R] Post-hoc tests in MASS using glm.nb I am struggling to generate p values for comparisons of levels (post-hoc tests) in a glm with a negative binomial distribution I am try

[R] Post-hoc tests in MASS using glm.nb

2011-05-16 Thread bryony
I am struggling to generate p values for comparisons of levels (post-hoc tests) in a glm with a negative binomial distribution I am trying to compare cell counts on different days as grown on different media (e.g. types of cryogel) so I have 2 explanatory variables (Day and Cryogel), which are bot