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) > > > > > Nina Hobbhahn, PhD > Claude Leon post-doctoral fellow > Lab of Prof. S. D. Johnson > School of Life Sciences > University of KwaZulu-Natal > Private Bag X01 > Scottsville, Pietermaritzburg, 3201 > South Africa > > _______________________________________________ > R-sig-phylo mailing list - R-sig-phylo@r-project.org > https://stat.ethz.ch/mailman/listinfo/r-sig-phylo > Searchable archive at http://www.mail-archive.com/r-sig-phylo@r-project.org/ > _______________________________________________ R-sig-phylo mailing list - R-sig-phylo@r-project.org https://stat.ethz.ch/mailman/listinfo/r-sig-phylo Searchable archive at http://www.mail-archive.com/r-sig-phylo@r-project.org/