Fixed effects model: changing the F, p values [Two-Stage / GS Designs]
Dear mittyri,
for the results consistency (R versus Phoenix as you reported), there should be one more term of marvelous effect x) ... you posted it as number 7: (I resorted it for comparison of F and p values, note: int and Trt are missing in the results from R code)
If you change the muddle:
and run the mentioned code... Voilà!
The comparison is telling us that with groups it will be little bit more complicated.
» Now Group factor seems to be more suitable.
» What bothers me is that Group effect is not equal in R and PHX, cannot figure out the reason
Good luck!
zizou
for the results consistency (R versus Phoenix as you reported), there should be one more term of marvelous effect x) ... you posted it as number 7: (I resorted it for comparison of F and p values, note: int and Trt are missing in the results from R code)
» Dependent Hypothesis Numer_DF Denom_DF F_stat P_value
» 5 Ln(Var) Group*Per 2 14 3.2110481e+00 0.071153125
» 7 Ln(Var) Trt*Group 1 14 2.3942317e-01 0.632200219
» 2 Ln(Var) Group 1 14 9.9615453e-01 0.335181899
» 3 Ln(Var) Seq 1 14 9.7094813e-01 0.341167964
» 4 Ln(Var) Group*Seq 1 14 5.7364354e-02 0.814182229
» 1 Ln(Var) int 1 14 3.1047428e+03 0.000000000
» 6 Ln(Var) Trt 1 14 3.4404496e+00 0.084798088
If you change the muddle:
muddle <- lm(log(Var)~Group+Seq+Seq*Group+Subj%in%(Seq)+Per%in%Group+Trt+Group*Trt, data=data)
and run the mentioned code... Voilà!
Type III sum of squares:
(mimicing SAS/SPSS but with Seq tested against Subj)
log(Var) ~ Group + Seq + Seq * Group + Subj %in% (Seq * Group) + ...
Df Sum Sq Mean Sq F value Pr(>F)
Group:Per 2 0.04514 0.0225681 3.21105 0.071153 .
Group:Trt 1 0.00168 0.0016827 0.23942 0.632200
Group 1 0.22546 0.2254640 0.79946 0.386368
Seq 1 0.27383 0.2738263 0.97095 0.341168
Group:Seq 1 0.01618 0.0161779 0.05736 0.814182
Group:Seq:Subj 14 3.94827 0.2820195 40.12641 7.5802e-09 ***
Residuals 14 0.09840 0.0070283
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
The comparison is telling us that with groups it will be little bit more complicated.
» Now Group factor seems to be more suitable.
» What bothers me is that Group effect is not equal in R and PHX, cannot figure out the reason
Good luck!
zizou
Complete thread:
- Potvin C in the EU Helmut 2013-04-16 17:40
- Potvin C in the EU ElMaestro 2013-04-17 02:31
- Potvin C in the EU Helmut 2013-04-17 12:23
- 2 Groups model FDA d_labes 2013-04-17 13:10
- 2 Groups model FDA Helmut 2013-04-17 16:37
- Group effects obsolete? d_labes 2013-04-18 09:55
- Group effects FDA/EMA Helmut 2013-04-19 14:34
- Group effects FDA/EMA mittyri 2014-08-21 15:34
- Group effects FDA/EMA ElMaestro 2014-08-21 17:31
- Group effects EMA mittyri 2014-10-01 11:25
- Group effects EMA ElMaestro 2014-10-01 11:31
- Group effects EMA Helmut 2014-10-01 13:24
- Group effects EMA VStus 2016-10-06 22:06
- Group effects EMA ElMaestro 2016-10-07 12:45
- Group effects EMA mittyri 2016-10-07 15:06
- Fixed effects model with Group term mittyri 2016-10-09 13:27
- Fixed effects model: changing the F, p values mittyri 2016-10-10 18:33
- Fixed effects model: changing the F, p values ElMaestro 2016-10-10 20:10
- Fixed effects model: changing the F, p values Helmut 2016-10-11 00:17
- Holy War of type III d_labes 2016-10-11 13:17
- Significance of Group effect in Russia: why the type III is so 'important' mittyri 2016-10-11 15:07
- Fixed effects model: changing the F, p valueszizou 2016-10-11 01:06
- FDA group model in R mittyri 2016-10-11 13:43
- FDA group model in R ElMaestro 2016-10-12 01:07
- FDA group model in R VStus 2016-10-17 15:32
- FDA group model in R ElMaestro 2016-10-17 18:58
- FDA group model in R VStus 2016-10-20 13:57
- FDA group model in R ElMaestro 2016-10-17 18:58
- FDA group model in R mittyri 2016-10-11 13:43
- Fixed effects model: changing the F, p values ElMaestro 2016-10-10 20:10
- Fixed effects model with Group term VStus 2016-10-20 13:54
- Fixed effects model: changing the F, p values mittyri 2016-10-10 18:33
- Fixed effects model with Group term mittyri 2016-10-09 13:27
- Group effects EMA mittyri 2016-10-07 15:06
- Group effects EMA ElMaestro 2016-10-07 12:45
- Group effects EMA VStus 2016-10-06 22:06
- Group effects FDA/EMA zizou 2016-12-31 02:54
- Group effects FDA/EMA Helmut 2017-06-03 14:57
- Group effects EMA mittyri 2014-10-01 11:25
- Group effects FDA/EMA ElMaestro 2014-08-21 17:31
- Group effects FDA/EMA mittyri 2014-08-21 15:34
- Group effects FDA/EMA Helmut 2013-04-19 14:34
- 2 Groups model FDA hiren379 2013-08-22 14:32
- Homework Helmut 2013-08-23 11:53
- Group effects obsolete? d_labes 2013-04-18 09:55
- 2 Groups model FDA Helmut 2013-04-17 16:37
- 2 Groups model FDA d_labes 2013-04-17 13:10
- Potvin C in the EU Helmut 2013-04-17 12:23
- Potvin C in the EU ElMaestro 2013-04-17 02:31