[R] error with effects package.

[Frédéric Hérault]

Dear John,

Thanks for the reply. I will look for another way to compute adjusted 
effect.

Frederic.

John Fox a écrit :
> Dear Frederic,
>
> The effect function won't handle a model with a nested effect.
>
> Sorry,
>  John
>
>   
>> -----Original Message-----
>> From: r-help-bounces at r-project.org [mailto:r-help-bounces at r-project.org]
>>     
> On
>   
>> Behalf Of Frédéric Hérault
>> Sent: March-18-09 9:57 AM
>> To: r-help at r-project.org
>> Subject: [R] error with effects package.
>>
>> Dear R helpers,
>>
>> I have the following model
>>
>> model_1<-
>>
>>     
> glm(y~A+B+C+E+A:D,contrasts=list(D=contrasts_D),data=mydata,na.action=na.omi
> t
>   
>> )
>>  with: options(contrasts=c("contr.sum", "contr.poly"))
>>
>> A,B and E are 2-levels factor,
>> C is covariate,
>> D is 20 levels factor  with 10  in relation with the first levels of
>> factor A, the other in relation with the second levels of factor A
>>
>>
>> Here is my data, and the contrast I use for factor D.
>> https://depot.jouy.inra.fr/get?k=xzsA9nLhQSNbzq55MXo
>>
>> I'm trying to use the effects package of Prof. J. Fox ,
>>
>> (eff_model_1<-effect("A",model1))
>>
>> but I get this error message:
>>  " Error in apply(mod.matrix[, components], 1, prod) :   subscript out
>> of bounds"
>>
>> I suppose that the problem comes from the interaction A:D, but I can't
>> find a solution to solve it.
>>
>> Does anyone know how to solve this problem??
>>
>> thanks.
>>
>> --------------------------------
>>
>> Frederic Herault
>> INRA
>> France
>>
>> ______________________________________________
>> R-help at r-project.org mailing list
>> https://stat.ethz.ch/mailman/listinfo/r-help
>> PLEASE do read the posting guide
>>     
> http://www.R-project.org/posting-guide.html
>   
>> and provide commented, minimal, self-contained, reproducible code.
>>     
>
> ______________________________________________
> R-help at r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
> and provide commented, minimal, self-contained, reproducible code.
>