----------------------------------------------------------------------------------------------------------------- name: log: C:\Users\IRILab001\Desktop\aula_lab1.log log type: text opened on: 29 Aug 2017, 20:13:02 . summarize Variable | Obs Mean Std. Dev. Min Max -------------+--------------------------------------------------------- patient | 120 60.5 34.78505 1 120 sex | 0 agegrp | 0 bp_before | 120 156.45 11.38985 138 185 bp_after | 120 151.3583 14.17762 125 185 . sum Variable | Obs Mean Std. Dev. Min Max -------------+--------------------------------------------------------- patient | 120 60.5 34.78505 1 120 sex | 0 agegrp | 0 bp_before | 120 156.45 11.38985 138 185 bp_after | 120 151.3583 14.17762 125 185 . summarize patient bp_before bp_after bp_after Variable | Obs Mean Std. Dev. Min Max -------------+--------------------------------------------------------- patient | 120 60.5 34.78505 1 120 bp_before | 120 156.45 11.38985 138 185 bp_after | 120 151.3583 14.17762 125 185 bp_after | 120 151.3583 14.17762 125 185 . summarize patient bp_before bp_after bp_after, detail patient ------------------------------------------------------------- Percentiles Smallest 1% 2 1 5% 6.5 2 10% 12.5 3 Obs 120 25% 30.5 4 Sum of Wgt. 120 50% 60.5 Mean 60.5 Largest Std. Dev. 34.78505 75% 90.5 117 90% 108.5 118 Variance 1210 95% 114.5 119 Skewness 0 99% 119 120 Kurtosis 1.799833 bp_before ------------------------------------------------------------- Percentiles Smallest 1% 140 138 5% 142 140 10% 142.5 140 Obs 120 25% 147 141 Sum of Wgt. 120 50% 154.5 Mean 156.45 Largest Std. Dev. 11.38985 75% 164 183 90% 172.5 184 Variance 129.7286 95% 175.5 185 Skewness .5472919 99% 185 185 Kurtosis 2.529885 bp_after ------------------------------------------------------------- Percentiles Smallest 1% 126 125 5% 131.5 126 10% 134 127 Obs 120 25% 140.5 131 Sum of Wgt. 120 50% 149.5 Mean 151.3583 Largest Std. Dev. 14.17762 75% 161 183 90% 170.5 184 Variance 201.005 95% 176.5 185 Skewness .3884312 99% 185 185 Kurtosis 2.466068 bp_after ------------------------------------------------------------- Percentiles Smallest 1% 126 125 5% 131.5 126 10% 134 127 Obs 120 25% 140.5 131 Sum of Wgt. 120 50% 149.5 Mean 151.3583 Largest Std. Dev. 14.17762 75% 161 183 90% 170.5 184 Variance 201.005 95% 176.5 185 Skewness .3884312 99% 185 185 Kurtosis 2.466068 . tabulate sex, sum(bp_before) | Summary of bp_before sex | Mean Std. Dev. Freq. ------------+------------------------------------ Female | 153.63333 10.7356 60 Male | 159.26667 11.413442 60 ------------+------------------------------------ Total | 156.45 11.389845 120 . tabulate sex, summarize(bp_before) | Summary of bp_before sex | Mean Std. Dev. Freq. ------------+------------------------------------ Female | 153.63333 10.7356 60 Male | 159.26667 11.413442 60 ------------+------------------------------------ Total | 156.45 11.389845 120 . histogram bp_before (bin=10, start=138, width=4.7) . histogram bp_before, discrete (start=138, width=1) . tabulate sex agegrp, summarize(bp_before) Means, Standard Deviations and Frequencies of bp_before | agegrp sex | 30-45 46-59 60+ | Total -----------+---------------------------------+---------- Female | 149.9 151.15 159.85 | 153.63333 | 8.378544 9.1379774 11.930831 | 10.7356 | 20 20 20 | 60 -----------+---------------------------------+---------- Male | 153.45 159.05 165.3 | 159.26667 | 9.9550305 12.377718 8.8442667 | 11.413442 | 20 20 20 | 60 -----------+---------------------------------+---------- Total | 151.675 155.1 162.575 | 156.45 | 9.2580872 11.459628 10.727122 | 11.389845 | 40 40 40 | 120 . histogram bp_before, discrete by(sex) . sum bp_before bp_after Variable | Obs Mean Std. Dev. Min Max -------------+--------------------------------------------------------- bp_before | 120 156.45 11.38985 138 185 bp_after | 120 151.3583 14.17762 125 185 . ztest2i 120 156.45 11.38985 120 151.3583 14.17762 x: Number of obs = 120 y: Number of obs = 120 ------------------------------------------------------------------------------ Variable | Mean Std. Err. z P>|z| [95% Conf. Interval] ---------+-------------------------------------------------------------------- x | 156.45 1.039746 150.469 0.0000 154.4121 158.4879 y | 151.3583 1.294234 116.948 0.0000 148.8216 153.895 ---------+-------------------------------------------------------------------- diff | 5.0917 1.660155 3.067 0.0022 1.837857 8.345543 ------------------------------------------------------------------------------ Ho: mean(x) - mean(y) = diff = 0 Ha: diff < 0 Ha: diff ~= 0 Ha: diff > 0 z = 3.0670 z = 3.0670 z = 3.0670 P < z = 0.9989 P > |z| = 0.0022 P > z = 0.0011 . display (156.45-151.3583)/(11.38985)^2/120+14.17762^2/120)^0.5 too many ')' or ']' r(132); . display (156.45-151.3583)/((11.38985)^2/120+14.17762^2/120)^0.5 3.0670035 . generate bpdiff= bp_before- bp_after . sum bpdiff Variable | Obs Mean Std. Dev. Min Max -------------+--------------------------------------------------------- bpdiff | 120 5.091667 16.7136 -41 45 . ztesti 120 5.091667 16.7136 0 One-sample z test ------------------------------------------------------------------------------ | Obs Mean Std. Err. Std. Dev. [95% Conf. Interval] ---------+-------------------------------------------------------------------- x | 120 5.091667 1.525736 16.7136 2.101279 8.082055 ------------------------------------------------------------------------------ mean = mean(x) z = 3.3372 Ho: mean = 0 Ha: mean < 0 Ha: mean != 0 Ha: mean > 0 Pr(Z < z) = 0.9996 Pr(|Z| > |z|) = 0.0008 Pr(Z > z) = 0.0004 . display 5.091667/(16.7136^2/120)^0.5 3.3371875 . ztesti 120 5.091667 16.7136 0, level(99) One-sample z test ------------------------------------------------------------------------------ | Obs Mean Std. Err. Std. Dev. [99% Conf. Interval] ---------+-------------------------------------------------------------------- x | 120 5.091667 1.525736 16.7136 1.161632 9.021702 ------------------------------------------------------------------------------ mean = mean(x) z = 3.3372 Ho: mean = 0 Ha: mean < 0 Ha: mean != 0 Ha: mean > 0 Pr(Z < z) = 0.9996 Pr(|Z| > |z|) = 0.0008 Pr(Z > z) = 0.0004 . . tab sex, sum( bp_before) | Summary of bp_before sex | Mean Std. Dev. Freq. ------------+------------------------------------ Female | 153.63333 10.7356 60 Male | 159.26667 11.413442 60 ------------+------------------------------------ Total | 156.45 11.389845 120 . tab sex agegrp, sum( bp_before) Means, Standard Deviations and Frequencies of bp_before | agegrp sex | 30-45 46-59 60+ | Total -----------+---------------------------------+---------- Female | 149.9 151.15 159.85 | 153.63333 | 8.378544 9.1379774 11.930831 | 10.7356 | 20 20 20 | 60 -----------+---------------------------------+---------- Male | 153.45 159.05 165.3 | 159.26667 | 9.9550305 12.377718 8.8442667 | 11.413442 | 20 20 20 | 60 -----------+---------------------------------+---------- Total | 151.675 155.1 162.575 | 156.45 | 9.2580872 11.459628 10.727122 | 11.389845 | 40 40 40 | 120 . ttest bp_before if sex=30-45 type mismatch r(109); . egen idage=group( agegrp) . ttest bp_before if agegrp=30-45 type mismatch r(109); . ttest bp_before if idage==1, by(sex) Two-sample t test with equal variances ------------------------------------------------------------------------------ Group | Obs Mean Std. Err. Std. Dev. [95% Conf. Interval] ---------+-------------------------------------------------------------------- Female | 20 149.9 1.873499 8.378544 145.9787 153.8213 Male | 20 153.45 2.226012 9.95503 148.7909 158.1091 ---------+-------------------------------------------------------------------- combined | 40 151.675 1.463832 9.258087 148.7141 154.6359 ---------+-------------------------------------------------------------------- diff | -3.55 2.90949 -9.439954 2.339954 ------------------------------------------------------------------------------ diff = mean(Female) - mean(Male) t = -1.2201 Ho: diff = 0 degrees of freedom = 38 Ha: diff < 0 Ha: diff != 0 Ha: diff > 0 Pr(T < t) = 0.1150 Pr(|T| > |t|) = 0.2299 Pr(T > t) = 0.8850 . egen idsex=group(sex) . ttest bp_before if idage==1, by(idsex) Two-sample t test with equal variances ------------------------------------------------------------------------------ Group | Obs Mean Std. Err. Std. Dev. [95% Conf. Interval] ---------+-------------------------------------------------------------------- 1 | 20 149.9 1.873499 8.378544 145.9787 153.8213 2 | 20 153.45 2.226012 9.95503 148.7909 158.1091 ---------+-------------------------------------------------------------------- combined | 40 151.675 1.463832 9.258087 148.7141 154.6359 ---------+-------------------------------------------------------------------- diff | -3.55 2.90949 -9.439954 2.339954 ------------------------------------------------------------------------------ diff = mean(1) - mean(2) t = -1.2201 Ho: diff = 0 degrees of freedom = 38 Ha: diff < 0 Ha: diff != 0 Ha: diff > 0 Pr(T < t) = 0.1150 Pr(|T| > |t|) = 0.2299 Pr(T > t) = 0.8850 . replace idsex=3 if idsex==1 (60 real changes made) . ttest bp_before if idage==1, by(idsex) Two-sample t test with equal variances ------------------------------------------------------------------------------ Group | Obs Mean Std. Err. Std. Dev. [95% Conf. Interval] ---------+-------------------------------------------------------------------- 2 | 20 153.45 2.226012 9.95503 148.7909 158.1091 3 | 20 149.9 1.873499 8.378544 145.9787 153.8213 ---------+-------------------------------------------------------------------- combined | 40 151.675 1.463832 9.258087 148.7141 154.6359 ---------+-------------------------------------------------------------------- diff | 3.55 2.90949 -2.339954 9.439954 ------------------------------------------------------------------------------ diff = mean(2) - mean(3) t = 1.2201 Ho: diff = 0 degrees of freedom = 38 Ha: diff < 0 Ha: diff != 0 Ha: diff > 0 Pr(T < t) = 0.8850 Pr(|T| > |t|) = 0.2299 Pr(T > t) = 0.1150 . replace idsex=1 if idsex==3 (60 real changes made) . ttest bpdiff=0 if idage==1 & idsex==1 One-sample t test ------------------------------------------------------------------------------ Variable | Obs Mean Std. Err. Std. Dev. [95% Conf. Interval] ---------+-------------------------------------------------------------------- bpdiff | 20 7.7 3.128898 13.99286 1.151142 14.24886 ------------------------------------------------------------------------------ mean = mean(bpdiff) t = 2.4609 Ho: mean = 0 degrees of freedom = 19 Ha: mean < 0 Ha: mean != 0 Ha: mean > 0 Pr(T < t) = 0.9882 Pr(|T| > |t|) = 0.0236 Pr(T > t) = 0.0118 . ttest bp_before= bp_after if idage==1 & idsex==1, paired option paired not allowed r(198); . ttest bp_before= bp_after if idage==1 & idsex==1, unpaired Two-sample t test with equal variances ------------------------------------------------------------------------------ Variable | Obs Mean Std. Err. Std. Dev. [95% Conf. Interval] ---------+-------------------------------------------------------------------- bp_bef~e | 20 149.9 1.873499 8.378544 145.9787 153.8213 bp_after | 20 142.2 2.039866 9.122557 137.9305 146.4695 ---------+-------------------------------------------------------------------- combined | 40 146.05 1.499551 9.483995 143.0169 149.0831 ---------+-------------------------------------------------------------------- diff | 7.7 2.769667 2.093103 13.3069 ------------------------------------------------------------------------------ diff = mean(bp_before) - mean(bp_after) t = 2.7801 Ho: diff = 0 degrees of freedom = 38 Ha: diff < 0 Ha: diff != 0 Ha: diff > 0 Pr(T < t) = 0.9958 Pr(|T| > |t|) = 0.0084 Pr(T > t) = 0.0042 . ttest bp_before= bp_after if idage==1 & idsex==1 Paired t test ------------------------------------------------------------------------------ Variable | Obs Mean Std. Err. Std. Dev. [95% Conf. Interval] ---------+-------------------------------------------------------------------- bp_bef~e | 20 149.9 1.873499 8.378544 145.9787 153.8213 bp_after | 20 142.2 2.039866 9.122557 137.9305 146.4695 ---------+-------------------------------------------------------------------- diff | 20 7.7 3.128898 13.99286 1.151142 14.24886 ------------------------------------------------------------------------------ mean(diff) = mean(bp_before - bp_after) t = 2.4609 Ho: mean(diff) = 0 degrees of freedom = 19 Ha: mean(diff) < 0 Ha: mean(diff) != 0 Ha: mean(diff) > 0 Pr(T < t) = 0.9882 Pr(|T| > |t|) = 0.0236 Pr(T > t) = 0.0118 . ttest bp_before= bp_after if idage==1 & idsex==1, unpaired Two-sample t test with equal variances ------------------------------------------------------------------------------ Variable | Obs Mean Std. Err. Std. Dev. [95% Conf. Interval] ---------+-------------------------------------------------------------------- bp_bef~e | 20 149.9 1.873499 8.378544 145.9787 153.8213 bp_after | 20 142.2 2.039866 9.122557 137.9305 146.4695 ---------+-------------------------------------------------------------------- combined | 40 146.05 1.499551 9.483995 143.0169 149.0831 ---------+-------------------------------------------------------------------- diff | 7.7 2.769667 2.093103 13.3069 ------------------------------------------------------------------------------ diff = mean(bp_before) - mean(bp_after) t = 2.7801 Ho: diff = 0 degrees of freedom = 38 Ha: diff < 0 Ha: diff != 0 Ha: diff > 0 Pr(T < t) = 0.9958 Pr(|T| > |t|) = 0.0084 Pr(T > t) = 0.0042 . ttesti 20 7.7 13.99286 0 One-sample t test ------------------------------------------------------------------------------ | Obs Mean Std. Err. Std. Dev. [95% Conf. Interval] ---------+-------------------------------------------------------------------- x | 20 7.7 3.128899 13.99286 1.15114 14.24886 ------------------------------------------------------------------------------ mean = mean(x) t = 2.4609 Ho: mean = 0 degrees of freedom = 19 Ha: mean < 0 Ha: mean != 0 Ha: mean > 0 Pr(T < t) = 0.9882 Pr(|T| > |t|) = 0.0236 Pr(T > t) = 0.0118 . tab idsex idage, sum( bp_before) Means, Standard Deviations and Frequencies of bp_before | group(agegrp) group(sex) | 1 2 3 | Total -----------+---------------------------------+---------- 1 | 149.9 151.15 159.85 | 153.63333 | 8.378544 9.1379774 11.930831 | 10.7356 | 20 20 20 | 60 -----------+---------------------------------+---------- 2 | 153.45 159.05 165.3 | 159.26667 | 9.9550305 12.377718 8.8442667 | 11.413442 | 20 20 20 | 60 -----------+---------------------------------+---------- Total | 151.675 155.1 162.575 | 156.45 | 9.2580872 11.459628 10.727122 | 11.389845 | 40 40 40 | 120 . tab idsex idage, sum( bp_after) Means, Standard Deviations and Frequencies of bp_after | group(agegrp) group(sex) | 1 2 3 | Total -----------+---------------------------------+---------- 1 | 142.2 144.3 155.1 | 147.2 | 9.1225574 10.063483 11.951217 | 11.742722 | 20 20 20 | 60 -----------+---------------------------------+---------- 2 | 146.45 157.25 162.85 | 155.51667 | 14.058056 15.562692 11.636037 | 15.243217 | 20 20 20 | 60 -----------+---------------------------------+---------- Total | 144.325 150.775 158.975 | 151.35833 | 11.893518 14.502851 12.286093 | 14.177622 | 40 40 40 | 120 . log close name: log: C:\Users\IRILab001\Desktop\aula_lab1.log log type: text closed on: 29 Aug 2017, 21:52:34 -----------------------------------------------------------------------------------------------------------------