Hello Everyone,
I am using "tukeyhsd" a user written command (Philip Ender) to run pairwise comparisons. When I run my code accompanied with "if" statements the output doesn't include the values for HSD-test. Below I present the output with and without restrictions on sample. Any suggestions on how to deal with that? I am trying to identify statistically significant differences.
I am using "tukeyhsd" a user written command (Philip Ender) to run pairwise comparisons. When I run my code accompanied with "if" statements the output doesn't include the values for HSD-test. Below I present the output with and without restrictions on sample. Any suggestions on how to deal with that? I am trying to identify statistically significant differences.
Code:
. tukeyhsd grade_pa7
(2251 missing values generated)
Tukey HSD pairwise comparisons for variable grade_pa7
studentized range critical value(.05, 7, 1241) = 4.1762749
uses harmonic mean sample size = 74.930
mean
grp vs grp group means dif HSD-test
-------------------------------------------------------
1 vs 2 4446.2820 5149.5753 703.2933 4.4950*
1 vs 3 4446.2820 4993.6247 547.3426 3.4983
1 vs 4 4446.2820 4942.7287 496.4467 3.1730
1 vs 5 4446.2820 4342.6855 103.5965 0.6621
1 vs 6 4446.2820 4663.2158 216.9337 1.3865
1 vs 7 4446.2820 4083.5000 362.7820 2.3187
2 vs 3 5149.5753 4993.6247 155.9507 0.9967
2 vs 4 5149.5753 4942.7287 206.8467 1.3220
2 vs 5 5149.5753 4342.6855 806.8898 5.1571*
2 vs 6 5149.5753 4663.2158 486.3596 3.1085
2 vs 7 5149.5753 4083.5000 1066.0753 6.8137*
3 vs 4 4993.6247 4942.7287 50.8960 0.3253
3 vs 5 4993.6247 4342.6855 650.9391 4.1604
3 vs 6 4993.6247 4663.2158 330.4089 2.1118
3 vs 7 4993.6247 4083.5000 910.1247 5.8170*
4 vs 5 4942.7287 4342.6855 600.0431 3.8351
4 vs 6 4942.7287 4663.2158 279.5129 1.7865
4 vs 7 4942.7287 4083.5000 859.2287 5.4917*
5 vs 6 4342.6855 4663.2158 320.5302 2.0486
5 vs 7 4342.6855 4083.5000 259.1855 1.6566
6 vs 7 4663.2158 4083.5000 579.7158 3.7052
.
end of do-file
. do "/var/folders/jm/09q0bdrn4rx0_z64txpj22gm0000gn/T//SD24066.000000"
. tukeyhsd grade_pa7 if pa7==1 & weartimeacc7==1 ///
& !missing(chd_sex) & both_data7==3 & grade_pa7!=7 & !missing(chd_age_pa7) & ///
!missing(value7) & !missing(support_academic7)
(2251 missing values generated)
Tukey HSD pairwise comparisons for variable grade_pa7
studentized range critical value(.05, 7, 1241) = 4.1762749
uses harmonic mean sample size = .
mean
grp vs grp group means dif HSD-test
-------------------------------------------------------
1 vs 2 4413.2058 5241.2971 828.0913 .*
1 vs 3 4413.2058 5006.1609 592.9551 .*
1 vs 4 4413.2058 4976.0498 562.8440 .*
1 vs 5 4413.2058 4238.8231 174.3828 .*
1 vs 6 4413.2058 4670.8303 257.6245 .*
1 vs 7 4413.2058 . . .*
2 vs 3 5241.2971 5006.1609 235.1362 .*
2 vs 4 5241.2971 4976.0498 265.2473 .*
2 vs 5 5241.2971 4238.8231 1002.4741 .*
2 vs 6 5241.2971 4670.8303 570.4668 .*
2 vs 7 5241.2971 . . .*
3 vs 4 5006.1609 4976.0498 30.1111 .*
3 vs 5 5006.1609 4238.8231 767.3378 .*
3 vs 6 5006.1609 4670.8303 335.3306 .*
3 vs 7 5006.1609 . . .*
4 vs 5 4976.0498 4238.8231 737.2267 .*
4 vs 6 4976.0498 4670.8303 305.2195 .*
4 vs 7 4976.0498 . . .*
5 vs 6 4238.8231 4670.8303 432.0072 .*
5 vs 7 4238.8231 . . .*
6 vs 7 4670.8303 . . .*
.
end of do-file
