[R-sig-ME] lme interaction result strange

Ista Zahn istazahn at gmail.com
Tue May 1 20:09:23 CEST 2012


Hi Charles,

You can print and assign contrasts in R using the contrasts function.
See ?contrasts for details. You may also wish to look at ?contr.SAS

HTH,
Ista

On Tue, May 1, 2012 at 1:53 PM, Charles Determan Jr <deter088 at umn.edu> wrote:
> Thank you for your response Ista.  Do you have a recommendations how I might
> address variation in the contrasts?
>
> Regards,
> Charles
>
>
> On Tue, May 1, 2012 at 12:51 PM, Ista Zahn <istazahn at gmail.com> wrote:
>>
>> Hi Charles,
>>
>> My first guess is that you have (a) categorical variable(s) in your
>> predictors, and that the contrasts in SAS are different than those in
>> R.
>>
>> Best,
>> Ista
>>
>> On Tue, May 1, 2012 at 1:39 PM, Charles Determan Jr <deter088 at umn.edu>
>> wrote:
>> > Dear R users,
>> >
>> > I have been working on replicating some linear mixed models from SAS.
>> >  The
>> > first one matches perfectly when the SAS model is simple with the three
>> > separate factors.
>> >
>> > fit=lme(var~group+Event_name+Died,
>> >    data=liv34,
>> >    random=~1|ID)
>> > anova.lme(fit, type="marginal", adjustSigma=F)
>> >
>> > However, once I put an interaction into the formula the values don't
>> > match.
>> >
>> > fit=lme(var~group+Event_name+Died+Event_name*Died,
>> >    data=liv34,
>> >    random=~1|ID)
>> > anova.lme(fit, type="marginal", adjustSigma=F)
>> >
>> >                          numDF denDF     F-value      p-value
>> > (Intercept)               1        91       111.20483  <.0001
>> > group                      1        23        0.46632     0.5015
>> > Event_name            5        91        1.14042     *0.3449*
>> > Died                       1        23        0.50989    * 0.4824*
>> > Event_name:Died     5       91        1.10436     0.3637
>> > Done.
>> >
>> > The numbers *bold* don't match up.  They should be approximately .0290
>> > and
>> > .1318 respectively.  The other two are still exact matches.  I know
>> > looking
>> > for exact matches is ambitious but the numbers should be at least
>> > similar
>> > that the conclusions don't change so drastically.
>> >
>> > Any thoughts as to why this discrepancy is happening would be most
>> > appreciated.
>> >
>> > Regards,
>> >
>> > Charles
>> >
>> >        [[alternative HTML version deleted]]
>> >
>> > _______________________________________________
>> > R-sig-mixed-models at r-project.org mailing list
>> > https://stat.ethz.ch/mailman/listinfo/r-sig-mixed-models
>
>



More information about the R-sig-mixed-models mailing list