Hi Ben, thanks - we tried that, but it still doesn't give us all combinations.
Nina On 2013-12-02, at 4:05 PM, Ben Bolker wrote: > On 13-12-02 07:00 AM, Nina Hobbhahn wrote: >> Dear fellow list users, >> >> how do I extract estimates for all possible combinations of effects in >> significant interactions from compar.gee? >> >> To illustrate: I'm interested in flower-head size in males and females of >> wind- or insect-pollinated species. My model looks like this: >> >> compar.gee(FHsize~poll+sex+poll*sex, phy=phy, data=data, family=gaussian) >> >> The output looks like this: >> >> Coefficients: >> Estimate S.E. t Pr(T > |t|) >> (Intercept) 244.51760 1.256536e+02 1.945966e+00 0.1285677 >> pollW 54.77936 5.384660e+01 1.017323e+00 0.3703143 >> sexM -29.03269 8.148837e-05 -3.562802e+05 0.0000000 >> pollW:sexM -134.18408 1.502553e-04 -8.930409e+05 0.0000000 >> >> >> However, to illustrate this significant interaction I need estimates for >> >> pollW:sexM >> pollW:sexF >> pollI:sexM >> pollI:sexF >> >> How do I make compar.gee return these estimates? >> >> Many thanks, >> >> Nina > > Try compar.gee(FHsize~poll:sex-1, phy=phy, data=data, family=gaussian) _______________________________________________ R-sig-phylo mailing list - [email protected] https://stat.ethz.ch/mailman/listinfo/r-sig-phylo Searchable archive at http://www.mail-archive.com/[email protected]/
