******************************************************************************************** TABLE 2 summary of summaries. see run _13 for other table 1's. ******************************************************************************************** TABLE 2. ******************************************************************************************** TABLE 2. ******************************************************************************************** TABLE 2. ******************************. ******************************.******************************. get FILE ='C:\weight\OVERWT_yhl\data_OVERWT_yhl\OVERWT_yhl_BOTH_05.sav'. var labels lostweight 'Lost 10+ lbs'. frequencies smoke. compute cursmoke = smoke. recode cursmoke (3=1)(1 2 = 0). var labels cursmoke 'Current Smoker'. frequencies diabada. compute diabetes = diabada. recode diabetes (3=1)(1 2 4 = 0). var labels diabetes 'Diabetes'. frequencies diabetes. *** I took these vars out of the tables because they made the tables too long. but they can be put back in. frequencies bmi_1990. value labels bmi_1990 2 'Low' 3 'Medium' 4 'High'. var labels bmi_1990 'BMI Category *'. VAR LABELS yhl_ALIVE 'Alive' yhl_BD 'BED' yhl_BK 'BLOCK' yhl_CP 'SOC' yhl_CG 'COG' yhl_DP 'DEP' yhl_FH 'FLUSH' yhl_FLW 'FLW' yhl_HP 'HOSP' yhl_LV 'LEV' yhl_SPL 'SPL' yhl_TW 'TWLK' yhl_VG 'EVGG' yhl_XD 'ADL' yhl_XS 'EXSTR' yhl_YD 'IADL'. compute cursmoke = cursmoke*100. compute lostweight = 100*lostweight. * General Tables. TABLES /FORMAT BLANK MISSING('.') /OBSERVATION= N yhl_ALIVE yhl_BD yhl_BK yhl_CP yhl_CG yhl_DP yhl_FH yhl_FLW yhl_HP yhl_LV yhl_SPL yhl_TW yhl_VG yhl_XD yhl_XS yhl_YD age1990 cursmoke lostweight diabetes /GBASE=CASES /TABLE=BLACK > ( N +age1990 +cursmoke+lostweight+ yhl_ALIVE+ yhl_SPL+ yhl_FLW +yhl_CP + yhl_DP + yhl_BD +yhl_FH + yhl_CG +yhl_HP +yhl_VG +yhl_LV + yhl_XS + yhl_YD + yhl_BK + yhl_XD + yhl_TW ) BY gend01 > BMI_1990 /STATISTICS sum( N) mean( yhl_SPL) mean( yhl_FLW) mean( yhl_CP) mean( yhl_ALIVE) mean( yhl_BD) mean( yhl_BK) mean( yhl_CG) mean( yhl_DP) mean( yhl_FH) mean( yhl_HP) mean( yhl_LV) mean( yhl_TW) mean( yhl_VG) mean( yhl_XD) mean( yhl_XS) mean( yhl_YD) mean(age1990) mean(cursmoke) mean(lostweight) / title 'Table 2, Mean YOL and YHL by BMI *, Gender and Race' /caption '* Low =18.5-24.9 (normal), Medium=25-29.9 (overweight), High = 30+ (obese).' 'tables_04.sps 9-24-2011' . var labels waistcat 'Waist Circumference *'. * General Tables. TABLES /FORMAT BLANK MISSING('.') /OBSERVATION= N yhl_ALIVE yhl_BD yhl_BK yhl_CP yhl_CG yhl_DP yhl_FH yhl_FLW yhl_HP yhl_LV yhl_SPL yhl_TW yhl_VG yhl_XD yhl_XS yhl_YD age1990 cursmoke lostweight diabetes /GBASE=CASES /TABLE=BLACK > ( N +age1990 +cursmoke+lostweight + yhl_ALIVE + yhl_SPL+ yhl_FLW +yhl_CP + yhl_DP + yhl_BD +yhl_FH + yhl_CG +yhl_HP +yhl_VG +yhl_LV + yhl_XS + yhl_YD + yhl_BK + yhl_XD + yhl_TW ) BY gend01 > waistcat /STATISTICS sum( N) mean( yhl_SPL) mean( yhl_FLW) mean( yhl_CP) mean( yhl_ALIVE) mean( yhl_BD) mean( yhl_BK) mean( yhl_CG) mean( yhl_DP) mean( yhl_FH) mean( yhl_HP) mean( yhl_LV) mean( yhl_TW) mean( yhl_VG) mean( yhl_XD) mean( yhl_XS) mean( yhl_YD) mean(age1990) mean(cursmoke) mean(lostweight) / title 'Table 3, Mean YOL and YHL by Waist Circumference, Gender and Race' /caption '* Medium = 84.5-96.5 for White Female, 94.0-107.5 for Black Female, 93.0-101.5 for Male' 'tables_04.sps 9-24-2011' . *** appendix table 11. crosstabs bmi_1990 by waistcat. ***appendix table 12. means age1990 cursmoke lostweight diabetes by bmi_1990 . means age1990 cursmoke lostweight diabetes by waistcat. ***make a table of yhl/yol . temporary. compute yhl_spl = yhl_spl /yhl_alive. * General Tables. TABLES /FORMAT BLANK MISSING('.') /OBSERVATION= N yhl_ALIVE yhl_BD yhl_BK yhl_CP yhl_CG yhl_DP yhl_FH yhl_FLW yhl_HP yhl_LV yhl_SPL yhl_TW yhl_VG yhl_XD yhl_XS yhl_YD age1990 cursmoke lostweight diabetes /GBASE=CASES /TABLE=BLACK > ( N +age1990 +cursmoke+lostweight+ yhl_ALIVE+ yhl_SPL+ yhl_FLW +yhl_CP + yhl_DP + yhl_BD +yhl_FH + yhl_CG +yhl_HP +yhl_VG +yhl_LV + yhl_XS + yhl_YD + yhl_BK + yhl_XD + yhl_TW ) BY gend01 > waistcat /STATISTICS sum( N) mean( yhl_SPL) mean( yhl_FLW) mean( yhl_CP) mean( yhl_ALIVE) mean( yhl_BD) mean( yhl_BK) mean( yhl_CG) mean( yhl_DP) mean( yhl_FH) mean( yhl_HP) mean( yhl_LV) mean( yhl_TW) mean( yhl_VG) mean( yhl_XD) mean( yhl_XS) mean( yhl_YD) mean(age1990) mean(cursmoke) mean(lostweight) / title 'Table 2, Mean YOL and YHL by BMI *, Gender and Race' /caption '* Low =18.5-24.9 (normal), Medium=25-29.9 (overweight), High = 30+ (obese).' 'tables_04.sps 9-24-2011' . ***70.11 69.14 67.96. ******************************************************************************************** TABLE 2_A table 3_a table 4_a. ******************************************************************************************** TABLE 2_A. ******************************************************************************************** TABLE 2_A. ******************************************************************************************** TABLE 2_A. *get file = 'C:\weight\overwt_YHL\data_overwt_yhl\runs_long_01.sav' . GET file = 'C:\weight\overwt_YHL\data_overwt_yhl\RESULTS_COMBINED_01.SAV'. FREQUENCIES Y. ***THIS IS FOR THE YHL OUTPUT. SELECT IF (Y = 1). ***CHANGE, PUT ALIVE AT THE TOP OF THE LIST. (NOT DONE YET). compute vnum3 = vnum. recode vnum3 (1=0)(2=6)(3=14)(4=8)(5=5)(6=7) (7=2)(8=9) (9=1) (10=16)(11=10) (12=15)(13=12)(14=13)(15=3)(16=11). frequencies vnum3. means b_a by vnum3. var labels vnum2 'vnum3 sorted by b_a'. value labels vnum3 0 'ALIVE' 6 'BED' 14 'BLOCK' 8 'COG' 5 'DEP' 7 'FLUSH' 2 'FLW' 9 'HOSP' 1 'SPL' 16 'TWLK' 10 'EVGG' 15 'ADL' 12 'EXSTR' 13 'IADL' 3 'SOC' 11 'LEV' . means b_a by vnum3. crosstabs vnum2 by vnum3. frequencies x2. value labels x2 1 'Med BMI' 2 'Med WC' 3 'Hi BMI' 4 'Hi WC'. ************************ MODEL A age only. Table 4, confounders only. TABLES /FORMAT BLANK MISSING('.') /OBSERVATION= b_a_sig /GBASE=CASES /TABLE=black > vnum3 > b_a_sig BY male > (STATISTICS) > X2 > (STATISTICS) /title 'Table 4 Significant (p<.10) effects of fatness on YHL, in years ' /caption 'Adjusted for age, smoking, and weight loss tables_19.sps 9-24-2011' . **************************************************************************new version, plump only, p-values. frequencies x x2. var labels b_a 'Difference' . var labels p_a 'Significance'. descriptives b_a p_a. compute x3 = x2. recode x3 (1=1)(2=2)(3=sysmis) (4=sysmis). value labels x3 1 'Med BMI' 2 'Med WC'. crosstabs x2 by x3. frequencies x y x3. * General Tables. TABLES /FORMAT BLANK MISSING('.') /OBSERVATION= b_a /GBASE=CASES /TABLE=black > vnum3 BY male > (STATISTICS) > x3 > (STATISTICS) > b_a > (STATISTICS)/ title 'coefficients only' . * General Tables. TABLES /FORMAT BLANK MISSING('.') /OBSERVATION= p_a /GBASE=CASES /TABLE=black > vnum3 BY male > (STATISTICS) > x3 > (STATISTICS) > p_a > (STATISTICS) / title 'p-values only' . **************************************************************************new version, HIGH only, p-values. frequencies x x2. var labels b_a 'Difference' . var labels p_a 'Significance'. descriptives b_a p_a. compute x4 = x2. recode x4 (1=sysmis)(2=sysmis)(3=3) (4=4). value labels x4 3 'High BMI' 4 'High WC'. crosstabs x4 by x2. frequencies x4. * General Tables. TABLES /FORMAT BLANK MISSING('.') /OBSERVATION= b_a /GBASE=CASES /TABLE=black > vnum3 BY male > (STATISTICS) > x4 > (STATISTICS) > b_a > (STATISTICS)/ title 'Appendix Table 4 HIGH coefficients only' . * General Tables. TABLES /FORMAT BLANK MISSING('.') /OBSERVATION= p_a /GBASE=CASES /TABLE=black > vnum3 BY male > (STATISTICS) > x4 > (STATISTICS) > p_a > (STATISTICS) / title 'Appendix Table 4 HIGH p-values only' . ************************ MODEL A BONFERRONI age only. not fixed yet. TABLES /FORMAT BLANK MISSING('.') /OBSERVATION= b_a_BON /GBASE=CASES /TABLE=black > vnum3 > b_a_BON BY male > (STATISTICS) > X2 > (STATISTICS) /title 'Significant (p<.10) effects of fatness on YHL, in years (Model A)' /caption 'MODEL A regression coefficient controlling for age only tables_19.sps 9-24-2011' . *************************************************************************************************** TABLE 2_B table 3_b table 4_b . *************************************************************************************************** TABLE 2_B. *************************************************************************************************** TABLE 2. *************************************************************************************************** TABLE 2. *get file = 'C:\weight\overwt_YHL\data_overwt_yhl\runs_long_01_B.sav' . GET file = 'C:\weight\overwt_YHL\data_overwt_yhl\RESULTS_COMBINED_01.SAV'. DESCRIPTIVES ALL. FREQUENCIES Y. compute vnum3 = vnum. recode vnum3 (1=0)(2=6)(3=14)(4=8)(5=5)(6=7) (7=2)(8=9) (9=1) (10=16)(11=10) (12=15)(13=12)(14=13)(15=3)(16=11). frequencies vnum3. means b_a by vnum3. var labels vnum2 'vnum3 sorted by b_a'. value labels vnum3 0 'ALIVE' 6 'BED' 14 'BLOCK' 8 'COG' 5 'DEP' 7 'FLUSH' 2 'FLW' 9 'HOSP' 1 'SPL' 16 'TWLK' 10 'EVGG' 15 'ADL' 12 'EXSTR' 13 'IADL' 3 'SOC' 11 'LEV' . means b_a by vnum3. crosstabs vnum2 by vnum3. frequencies x2. value labels x2 1 'Med BMI' 2 'Med WC' 3 'Hi BMI' 4 'Hi WC'. ***THIS IS FOR THE YHL/YHL OUTPUT. SELECT IF (Y = 2). ***these tables do not have "alive" in them, so i don't need to change to vnum3. **********************************. *** Model A. for yhl/yol. temporary. select if (vnum ne 1). TABLES /FORMAT BLANK MISSING('.') /OBSERVATION= b_a_sig /GBASE=CASES /TABLE=black > vnum3 > b_a_sig BY male > (STATISTICS) > X2 > (STATISTICS) /title 'Table 5 Significant (p<.10) effects of fatness on YHL/YOL, in % points' /caption 'Adjusted for age, smoking, and recent weight loss. tables_19.sps 9-24-2011' . ****************************************************** WHAT IS SIGNIFICANT?. COUNT THE NIUMBER SIGNIFICANT IN DIFFERENT CATEGORIES. ***see run _13 for more info in this area. ************************************************************************************** THIS IS THE ACTUAL RUN. appendix table 10. GET file = 'C:\weight\overwt_YHL\data_overwt_yhl\RESULTS_COMBINED_01.SAV'. descriptives all. ******* get rid of yol, then do a count. temporary. select if (vnum > 1). frequencies black male x y sig05a sig05b sig05c. ************************** from version 26 on, I keep death in. ***SELECT IF (VNUM > 1). frequencies black male x y sig05a sig05b sig05c. *** I would life signifi pos and signif neg. do repeat sigall = sig05a sig05b sig05c/ beta = b_a b_b b_c/ sigpos = sigposa sigposb sigposc/ signeg = signega signegb signegc. compute sigpos = sigall. if (beta le 0) sigpos = 0. compute signeg = sigall. if (beta ge 0) signeg = 0. end repeat. crosstabs sigposa by x by y by sexrace. ***but I no longer let Y = 2. frequencies y. temporary. select if (y = 1). crosstabs sigposa by x by y by sexrace. temporary. select if (y = 1). crosstabs signega by x by y by sexrace. descriptives sig05a sig05b sig05c sigposa to signegc. crosstabs sigposa signega by sig05a. crosstabs sigposb signegb by sig05b. crosstabs sigposc signegc by sig05c. *how many are significant for each variable?.... frequencies sigposa signega. means sigposa signega by vnum2. ***goes in table A10, number sig for each variable. temporary. select if (y = 1). MEANS TABLES=signega sigposa BY vnum2 /CELLS COUNT SUM . temporary. select if (y = 1). crosstabs sigposa by x by y by sexrace. temporary. select if (y = 1). crosstabs signega by x by y by sexrace. ***HOW MANY ARE SIGNIFICANT IN EACH QUADRANT OF TABLE 2A?.. CROSSTABS SIG05A BY MALE BY X BY BLACK BY Y . CROSSTABS SIG05A BY MALE BY BLACK BY Y . CROSSTABS SIG05B BY MALE BY BLACK BY Y . CROSSTABS SIG05C BY MALE BY BLACK BY Y . ***POSITIVE ONLY. CROSSTABS SIGposA BY MALE BY BLACK BY Y . CROSSTABS SIGposB BY MALE BY BLACK BY Y . CROSSTABS SIGposC BY MALE BY BLACK BY Y . ***negative ONLY. CROSSTABS SIGnegA BY MALE BY BLACK BY Y . CROSSTABS SIGnegB BY MALE BY BLACK BY Y . CROSSTABS SIGnegC BY MALE BY BLACK BY Y . frequencies y. *2a. temporary. select if (sig05a = 1 and y = 1 ). list sigposa vnum2 b_a p_a sig05a sigposa signega x y. *2b. temporary. select if (sig05a = 1 and y = 2 ). list sigposa vnum2 b_a p_a sig05a sigposa signega x y. *3a. temporary. select if (sig05b = 1 and y = 1 ). list sigposa vnum2 b_a p_a sig05b sigposb signegb x y. *3b. temporary. select if (sig05b = 1 and y = 2 ). list sigposa vnum2 b_a p_a sig05b sigposb signegb x y. *4a. temporary. select if (sig05c = 1 and y = 1 ). list sigposa vnum2 b_a p_a sig05c sigposc signegc x y. *4b. temporary. select if (sig05c = 1 and y = 2 ). list sigposa vnum2 b_a p_a sig05c sigposc signegc x y. ************************ MODEL A age only. do significant positive only. temporary. select if (sigposa = 1 and y = 1). TABLES /FORMAT BLANK MISSING('.') /OBSERVATION= b_a_sig /GBASE=CASES /TABLE=black > vnum2 > b_a_sig BY male > (STATISTICS) > X2 > (STATISTICS) /title 'Table 2A Significant (p<.10) effects of fatness on YHL, in years (Model A)' /caption ' MODEL A regression coefficient controlling for age only tables_19.sps 9-24-2011' . ***4 entries from the real table 2a that have a positive coeff are not in this table. temporary. select if (sigposb = 1 and y = 1). TABLES /FORMAT BLANK MISSING('.') /OBSERVATION= b_b_sig /GBASE=CASES /TABLE=black > vnum2 > b_b_sig BY male > (STATISTICS) > X2 > (STATISTICS) /title 'Table 3a Significant (p<.10) effects of fatness on YHL, in years (Model A)' /caption ' MODEL A regression coefficient controlling for age only tables_19.sps 9-24-2011' . temporary. select if (sigposc = 1 and y = 1). TABLES /FORMAT BLANK MISSING('.') /OBSERVATION= b_c_sig /GBASE=CASES /TABLE=black > vnum2 > b_c_sig BY male > (STATISTICS) > X2 > (STATISTICS) /title 'Table 4a Significant (p<.10) effects of fatness on YHL, in years (Model A)' /caption ' MODEL A regression coefficient controlling for age only tables_19.sps 9-24-2011' . temporary. select if (vnum2 = 5 and male = 1 and x=1 and black = 0 and y = 1). descriptives all. frequencies x y. **x is 4 categories med bmi, med wc, hi bmi, hi wc.. **y = is yhl or yhl/yol. temporary. select if (sig05a = 1 and x = 2). list sigposa vnum2 b_a p_a sig05a sigposa signega x y. temporary. select if (sigposa = 1 and x = 1). list sigposa vnum b_a p_a sig05a sigposa signega x y. temporary. select if (sigposa = 1 and x = 2). list sigposa vnum2 b_a p_a sig05a sigposa signega x y. means sigposa sigposb sigposc by sexrace. means sigposa sigposb sigposc by y. means signega signegb signegc by sexrace. means signega signegb signegc by y. ***************************************************** but i now want only model a, which is the confounders model. crosstabs sigposa by x by y by sexrace . crosstabs signega by x by y by sexrace . ******************************************************************************************************************* list the p-values and the coeffs?. ******************************************************************************************************************* list the p-values and the coeffs?. ******************************************************************************************************************* list the p-values and the coeffs?. GET file = 'C:\weight\overwt_YHL\data_overwt_yhl\RESULTS_COMBINED_01.SAV'. descriptives all. compute vnum3 = vnum. recode vnum3 (1=0)(2=6)(3=14)(4=8)(5=5)(6=7) (7=2)(8=9) (9=1) (10=16)(11=10) (12=15)(13=12)(14=13)(15=3)(16=11). frequencies vnum3. means b_a by vnum3. var labels vnum2 'vnum3 sorted by b_a'. value labels vnum3 0 'ALIVE' 6 'BED' 14 'BLOCK' 8 'COG' 5 'DEP' 7 'FLUSH' 2 'FLW' 9 'HOSP' 1 'SPL' 16 'TWLK' 10 'EVGG' 15 'ADL' 12 'EXSTR' 13 'IADL' 3 'SOC' 11 'LEV' . means b_a by vnum3. crosstabs vnum2 by vnum3. frequencies x2. value labels x2 1 'Med BMI' 2 'Med WC' 3 'Hi BMI' 4 'Hi WC'. ************************ MODEL A age only. temporary. select if ( y = 1). TABLES /FORMAT BLANK MISSING('.') /OBSERVATION= b_a /GBASE=CASES /TABLE=black > vnum3 > b_a BY male > (STATISTICS) > X2 > (STATISTICS) /title 'Appendix Table A4 effects of fatness on YHL, in years ' /caption ' tables_19.sps 9-24-2011' . temporary. select if ( y = 2). TABLES /FORMAT BLANK MISSING('.') /OBSERVATION= b_a /GBASE=CASES /TABLE=black > vnum3 > b_a BY male > (STATISTICS) > X2 > (STATISTICS) /title 'Appendix Table A5 effects of fatness on YHL/YOL, in % points ' /caption ' tables_19.sps 9-24-2011' .