@@ -38,26 +38,32 @@ for(i in 1:length(mods))
38
38
library(jointest )
39
39
res = join_flipscores(mods ,n_flips = 5000 , seed = 1 , tested_coeffs = " X" )
40
40
summary(res )
41
- # > Model Coeff Estimate Score Std. Error z value Part. Cor Pr(>|t|)
42
- # > 1 mod1 X 1.1254046 19.78375 6.309305 3.135646 0.7605058 0.01739652
43
- # > 2 mod2 X 0.9552644 15.38036 5.053879 3.043278 0.7608195 0.01519696
44
- # > 3 mod3 X 1.0121921 14.31866 4.836669 2.960437 0.7643816 0.01859628
45
- # > 4 mod4 X 1.1696906 17.62979 5.875380 3.000621 0.7501552 0.00799840
41
+ # > Model Coeff Estimate Score Std. Error z value Part. Cor p
42
+ # > 1 mod1 X 1.1254046 19.78375 6.309305 3.135646 0.7605058 0.01499700
43
+ # > 2 mod2 X 0.9552644 15.38036 5.053879 3.043278 0.7608195 0.01339732
44
+ # > 3 mod3 X 1.0121921 14.31866 4.836669 2.960437 0.7643816 0.01279744
45
+ # > 4 mod4 X 1.1696906 17.62979 5.875380 3.000621 0.7501552 0.00979804
46
+ ```
47
+
48
+ ``` r
46
49
summary(combine(res ))
47
- # > Coeff Stat nMods S p
48
- # > 1 X maxT 4 19.78375 0.01879624
49
- res = p.adjust.fwer(res )
50
+ # > Model Coeff Stat nMods S p
51
+ # > 1 Overall X maxT 4 19.78375 0.00739852
52
+ ```
53
+
54
+ ``` r
55
+ res = jointest ::: p.adjust.jointest(res )
50
56
summary(res )
51
- # > Model Coeff Estimate Score Std. Error z value Part. Cor Pr(>|t|)
52
- # > 1 mod1 X 1.1254046 19.78375 6.309305 3.135646 0.7605058 0.01739652
53
- # > 2 mod2 X 0.9552644 15.38036 5.053879 3.043278 0.7608195 0.01519696
54
- # > 3 mod3 X 1.0121921 14.31866 4.836669 2.960437 0.7643816 0.01859628
55
- # > 4 mod4 X 1.1696906 17.62979 5.875380 3.000621 0.7501552 0.00799840
56
- # > p.adj
57
- # > 1 0.01879624
58
- # > 2 0.02019596
59
- # > 3 0.02019596
60
- # > 4 0.01879624
57
+ # > Model Coeff Estimate Score Std. Error z value Part. Cor p
58
+ # > 1 mod1 X 1.1254046 19.78375 6.309305 3.135646 0.7605058 0.01499700
59
+ # > 2 mod2 X 0.9552644 15.38036 5.053879 3.043278 0.7608195 0.01339732
60
+ # > 3 mod3 X 1.0121921 14.31866 4.836669 2.960437 0.7643816 0.01279744
61
+ # > 4 mod4 X 1.1696906 17.62979 5.875380 3.000621 0.7501552 0.00979804
62
+ # > p.adj
63
+ # > 1 0.019996
64
+ # > 2 0.019996
65
+ # > 3 0.019996
66
+ # > 4 0.019996
61
67
```
62
68
63
69
## References
0 commit comments