confint() for difference of effect levels [🇷 for BE/BA]
Hi d_labes,
thanks for your input.
Except when we fit the model with intercept, treatment first. Now the SE for whichever treatment is not the intercept, is the SE of the difference (and the other treatment effect is unestimable because of the intercept). The funny logic of linear models. I need to get my head around that aspect.
When we look at treatment effect differences it does not matter how we define the effects by contrast coding. I don't absolutely need an LSMean; I can do with two model effects calcuated in any valid way as long as I can extract their difference.
Yes, I have heard of that big regulatory body. But there is a problem. All-at-once is the preferred method of another big regulatory body.



, and we are not throwing away treatments when fitting models, how would you approach the CI derivation in R?
Many thanks.
thanks for your input.
❝ No. confint()
gives you confidence intervals for the model parameters only.
Except when we fit the model with intercept, treatment first. Now the SE for whichever treatment is not the intercept, is the SE of the difference (and the other treatment effect is unestimable because of the intercept). The funny logic of linear models. I need to get my head around that aspect.
❝ Or if you prefer using such beasts like "least square means" have a look at package lsmeans.
When we look at treatment effect differences it does not matter how we define the effects by contrast coding. I don't absolutely need an LSMean; I can do with two model effects calcuated in any valid way as long as I can extract their difference.
❝ BTW: All-at-once is not the preferred method of a big regulatory body.
Yes, I have heard of that big regulatory body. But there is a problem. All-at-once is the preferred method of another big regulatory body.




❝ Pairwise lm objects (throwing away the data under those treatments not under consideration) will give you the opportunity to use confint()
to derive CI's.

Many thanks.
—
Pass or fail!
ElMaestro
Pass or fail!
ElMaestro
Complete thread:
- confint() for difference of effect levels ElMaestro 2017-08-29 10:20 [🇷 for BE/BA]
- confint() for difference of effect levels d_labes 2017-08-29 12:23
- confint() for difference of effect levelsElMaestro 2017-08-29 14:45
- Homecooked solution? ElMaestro 2017-08-29 23:18
- Restaurant solution mittyri 2017-08-30 08:11
- Great ElMaestro 2017-08-30 09:27
- Restaurant solution mittyri 2017-08-30 08:11
- All-at-once or Two-at-a-time d_labes 2017-08-30 19:55
- Homecooked solution? ElMaestro 2017-08-29 23:18
- confint() for difference of effect levelsElMaestro 2017-08-29 14:45
- confint() for difference of effect levels d_labes 2017-08-29 12:23