Hello-

I'm working on an R script. I have this so far:

> curGLM <- glm(QIDSSR16 ~ TreatmentGroupDescription + baseline, family = gaussian, data=AllData);
> summary(curGLM);

Which gives a nice output, and does the calculation I want; but I can't figure out how to get all the pairwise p-Values out.

The "Summary(curGLM);" statement gives:
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) -2.6379 2.7638 -0.954 0.341983
TreatmentGroupDescription_TheFirstGroup 2.0502 1.1303 1.814 0.072473 .
TreatmentGroupDescriptionPlacebo 1.0421 1.2753 0.817 0.415656
baseline 0.5369 0.1508 3.560 0.000554 ***

However, there are 3 treatment groups, and I am only getting the p-Value for the diff of means between 1&2 and 1&3. I need 2&3 also. Any ideas?

Furthermore, I would like to combine two of my groups and then compare against the third (i.e. 1+2 vs 3), using the same model (curGLM) that I already have. Help on this also appreciated.

Thanks,
Leo