******************************************************************* ******************************************************************* * This is the do file for * * Schipper, B.C. (2022/11). Sex hormones and choice under risk, * University of California, Davis * ******************************************************************* ******************************************************************* * Notice to users ******************************************************************* * Thanks for your interest. * By using the data the user agrees to cite * * Schipper (2022) "Sex hormones and choice of risk", University of * California, Davis * * Schipper (2015) "Sex hormones and competitive bidding", Management * Science 61, 249--266 * * Pearson and Schipper (2013) "Menstrual cycle and competitive bidding", * Games and Economic Behavior 78, 1-20 * * Pearson and Schipper (2012) "The Visible Hand: Finger ratio (2D:4D) * and competitive bidding", Experimental Economics 15, 510-529, * * in any work by using this data or this STATA do file or any derivative * thereof. The data has been used in above papers. * You also agree to email me any paper in which you use the data. ******************************************************************* * The data can be downloaded from * http://www.econ.ucdavis.edu/faculty/schipper/sexhormones_mens_dr_data.zip * Be careful, the file is huge (about 100 MB). * Email: bcschipper@ucdavis.edu * I am happy if the data are used again. Have fun. ******************************************************************* clear clear matrix set more off set mem 50m ******************************************* ******************************************* *COMBINE DATA ******************************************* ******************************************* * mensdr study data set (2007 data set) use "C:\bcschipper\burkhard\research\mensdr\stats\mensdr_quest_20110626.dta" sort session subject save "C:\bcschipper\burkhard\research\mensdr\stats\mensdr_quest_20110626.dta", replace clear use "C:\bcschipper\burkhard\research\mensdr\stats\mensdr_dr_20110626.dta" sort session subject period save "C:\bcschipper\burkhard\research\mensdr\stats\mensdr_dr_20110626.dta", replace clear use "C:\bcschipper\burkhard\research\mensdr\stats\mensdr_auction_20110626.dta" sort session subject period save "C:\bcschipper\burkhard\research\mensdr\stats\mensdr_auction_20110626.dta", replace merge session subject using "C:\bcschipper\burkhard\research\mensdr\stats\mensdr_quest_20110626.dta", uniqusing sort session subject period drop _merge merge session subject period using "C:\bcschipper\burkhard\research\mensdr\stats\mensdr_dr_20110626.dta", uniqusing sort session subject period drop _merge gen new = 0 * generates variable indicating 2007 data set save "C:\bcschipper\burkhard\research\mensdr\stats\mensdr_data_20110626.dta", replace clear use "C:\bcschipper\burkhard\research\mensdr\stats\mensdr_data_20110626.dta", sort session subject period save "C:\bcschipper\burkhard\research\mensdr\stats\mensdr_data_20110626.dta", replace clear * hormone study data set (2010 data set) use "C:\bcschipper\burkhard\research\hormone\stats\hormone_quest_20100824.dta" sort session subject save "C:\bcschipper\burkhard\research\hormone\stats\hormone_quest_20100824.dta", replace clear use "C:\bcschipper\burkhard\research\hormone\stats\hormone_saliva_20100824.dta" *egen sdtestf = sd(testf) *label variable sdtestf "Testosterone" *egen sdestf = sd(estf) *label variable sdestf "Estradiol" *egen sdprogf = sd(progf) *label variable sdprogf "Progesterone" *egen sdcortf = sd(cortf) *label variable sdcortf "Cortisol" *egen sdtests = sd(tests) *egen sdests = sd(ests) *egen sdprogs = sd(progs) *egen sdcorts = sd(progs) *egen meantestf = mean(testf) *egen meanestf = mean(estf) *egen meanprogf = mean(progf) *egen meancortf = mean(cortf) *egen meantests = mean(tests) *egen meanests = mean(ests) *egen meanprogs = mean(progs) *egen meancorts = mean(corts) *gen testavg = (testf + tests)/2 *gen estavg = (estf + ests)/2 *gen progavg = (progf + progs)/2 *gen cortavg = (cortf + corts)/2 *egen sdtestavg = sd(testavg) *egen sdestavg = sd(estavg) *egen sdprogavg = sd(progavg) *egen sdcortavg = sd(cortavg) sort session subject save "C:\bcschipper\burkhard\research\hormone\stats\hormone_saliva_20100824.dta", replace clear use "C:\bcschipper\burkhard\research\hormone\stats\hormone_risk_20100824.dta" sort session subject save "C:\bcschipper\burkhard\research\hormone\stats\hormone_risk_20100824.dta", replace clear use "C:\bcschipper\burkhard\research\hormone\stats\hormone_dr_20110728.dta" sort session subject save "C:\bcschipper\burkhard\research\hormone\stats\hormone_dr_20110728.dta", replace clear use "C:\bcschipper\burkhard\research\hormone\stats\hormone_auction_20100824.dta" sort session subject period save "C:\bcschipper\burkhard\research\hormone\stats\hormone_auction_20100824.dta", replace merge session subject using "C:\bcschipper\burkhard\research\hormone\stats\hormone_quest_20100824.dta", uniqusing sort session subject period drop _merge merge session subject using "C:\bcschipper\burkhard\research\hormone\stats\hormone_saliva_20100824.dta", uniqusing sort session subject period drop _merge merge session subject using "C:\bcschipper\burkhard\research\hormone\stats\hormone_risk_20100824.dta", uniqusing sort session subject period drop _merge merge session subject using "C:\bcschipper\burkhard\research\hormone\stats\hormone_dr_20110728.dta", uniqusing sort session subject period drop _merge gen new = 1 * indicates 2010 data set save "C:\bcschipper\burkhard\research\hormone\stats\hormone_data_20120906.dta", replace clear use "C:\bcschipper\burkhard\research\hormone\stats\hormone_data_20120906.dta" sort session subject period append using "C:\bcschipper\burkhard\research\mensdr\stats\mensdr_data_20110626.dta" *./. ****************************************** ****************************************** * Generate labels and covariates ********* ****************************************** ****************************************** ******************************** * GENERATE ARBITRARY ID VARIABLE ******************************** gen test = substr(session, 3, 7) gen test2 = regexr(test,"_16","1") gen test3 = regexr(test2,"_17","2") gen test4 = regexr(test3,"_15","3") destring test4, gen(test5) egen id=concat(test5 subject) destring id, replace drop test test2 test3 test4 test5 label variable id "Arbitrary ID number" ********************** * Demographics ********************** label variable siblings "Num. of Siblings" label variable age "Age" label variable height "Height" label variable weight "Weight" * One subject states a weight of 60 pounds. Probably mixed it up with 60kg. * We replace the weight with 132 pounds. replace weight = 132 if weight == 60 gen bmi = (weight * 703)/(height^2) if new == 1 label variable bmi "BMI" gen dbmi = bmi - 21.75 if new == 1 label variable dbmi "Difference Standard BMI" gen obese = 0 if new == 1 replace obese = 1 if bmi > 30 label variable obese "Obese" gen age2 = age^2 label variable age2 "Age Squared" gen ysiblings = youngersib label variable ysiblings "Num. of Younger Siblings" gen osiblings = oldersib label variable osiblings "Num. of Older Siblings" gen totalsiblings = ysiblings + osiblings gen consistentsib = totalsiblings - siblings egen meanage = mean(age) egen meansib = mean(siblings) gen cage = age - meanage gen csib = siblings - meansib gen female = regexm(gender,"Female") gen female1 = regexm(gender,"female") replace female = 1 if female == 1 | female1 == 1 drop female1 label variable female "Female" gen gender1 = "Female" if female == 1 replace gender1 = "Male" if female == 0 *(This is for consistent labels in graphs.) * Ethnicities ************* gen white = regexm(ethnic,"White") label variable white "White" gen black = regexm(ethnic,"African") label variable black "Black" gen asian = regexm(ethnic,"Asian") label variable asian "Asian" gen hispanic = regexm(ethnic,"Hispanic") label variable hispanic "Hispanic" gen other = regexm(ethnic,"Other") label variable other "Other" * Generate variable comparable to 'Other Ethnicities' in CKO 2005 gen otherethn = 0 replace otherethn = 1 if hispanic==1 | other==1 | black==1 label variable otherethn "Other" * Academic variables ******************** gen science = regexm(major,"Chemistry") label variable science "Chemistry" * category is "Chemistry/Physics/Biology" gen othersci = regexm(major,"Natural Science") label variable othersci "Natural Science" gen humanities = regexm(major,"Humanities") label variable humanities "Humanities" gen econ = regexm(major,"Economics") label variable econ "Economics" gen socsci = regexm(major,"Social Science") label variable socsci "Social Science" gen eng = regexm(major,"Engineering") label variable eng "Engineering" gen math = regexm(major,"Mathematics") label variable math "Mathematics" gen english = regexm(major,"English") label variable english "English" replace humanities = 1 if english == 1 * Generate variable comparable to majors in CKO 2005 gen sciothersci=0 replace sciothersci=1 if science==1 | othersci==1 | eng==1 label variable sciothersci "Science \& Engineering" * Generate more aggregate variables of majors gen msci = 0 replace msci = 1 if sciothersci == 1 | math == 1 label variable msci "Science, Eng., Math." gen sochum = 0 replace sochum = 1 if socsci == 1 | humanities == 1 | english == 1 label variable sochum "Social Sci., Humanities" * Further academic variables label variable gpa "GPA" label variable quarters "Quarters Studied" label variable sat "SAT" label variable satyear "SAT Year" gen sat2005 = 1 if satyear < 2005 & new == 1 label variable sat2005 "SAT before 2005" * Sexual orientation, dating and relationship status **************************************************** gen hetero = regexm(sex_orient, "Heterosexual") if new == 1 label variable hetero "heterosexual" gen homo = regexm(sex_orient, "Homosexual") if new == 1 label variable homo "homosexual" gen bi = regexm(sex_orient, "Bisexual") if new == 1 label variable bi "bisexual" gen homobi = 0 if new == 1 replace homobi = 1 if homo == 1 | bi == 1 label variable homobi "Homo/Bisexual" gen bgfriend = regexm(relationship, "Boyfriend/Girlfriend") if new == 1 label variable bgfriend "Boy/Girlfriend" gen marry = regexm(relationship, "Married") if new == 1 label variable marry "Married" gen relation = 0 if new == 1 replace relation = 1 if bgfriend == 1 | marry == 1 label variable relation "In relationship" gen date1 = regexm(date, "1 person") if new == 1 gen date2 = regexm(date, "2 persons") if new == 1 gen date3 = regexm(date, "3 persons") if new == 1 gen date4 = regexm(date, "4 persons") if new == 1 gen date5 = regexm(date, "5 persons") if new == 1 gen date6 = regexm(date, "6 persons") if new == 1 gen date7 = regexm(date, "7 persons") if new == 1 gen date8 = regexm(date, "8 persons") if new == 1 gen date9 = regexm(date, "9 persons") if new == 1 gen date10 = regexm(date, "10 persons") if new == 1 gen date11 = regexm(date, "11 persons") if new == 1 gen date12 = regexm(date, "12 persons") if new == 1 gen dating = 0 if new == 1 replace dating = date1 + 2*date2 + 3*date3 + 4*date4 + 5*date5 + 6*date6 + 7*date7 + 8*date8 + 9*date9 + 10*date10 + 11*date11 + 12*date12 drop date1 date2 date3 date4 date5 date6 date7 date8 date9 date10 date11 date12 label variable dating "No. of Person's Dating" gen kids = regexm(num_kids, "1 child") if new == 1 * Remark: No one in the sample has more than one child label variable kids "Child" gen preg = regexm(pregnant, "Yes") if new == 1 label variable preg "Pregnant" gen roommate = regexm(living, "roommate") if new == 1 label variable roommate "Living with roommate" gen partner = regexm(living, "partner") if new == 1 label variable partner "Living with boyfriend/girlfriend/partner/spouse" gen parents = regexm(living, "parents") if new == 1 label variable parents "Living with parents" gen livalone = regexm(living, "alone") if new == 1 label variable livalone "Living alone" * Dietary preferences, smoking, live style ****************************************** gen smokeO = regexm(smoke, "Occasionally") if new == 1 gen smokeD = regexm(smoke, "Daily") if new == 1 gen smoking = 0 if new == 1 replace smoking = 1 if smokeO == 1 replace smoking = 2 if smokeD == 1 label variable smoking "Smoking" gen supplement = regexm(supplements, "Yes") if new == 1 label variable supplement "Supplements" gen vegetarian = regexm(veg, "Yes") if new == 1 label variable vegetarian "Vegetarian" gen soyN = regexm(soy, "Not very often") if new == 1 label variable soyN "Occational soybean food" gen soyW = regexm(soy, "several times a week") if new == 1 label variable soyW "Weekly soybean food" gen soyD = regexm(soy, "daily") if new == 1 label variable soyD "Daily soybean food" gen soybean = 0 if new == 1 replace soybean = soyN + 2*soyW + 3*soyD label variable soybean "Soy consum." * Exercising behavior *********************** gen workout1 = regexm(exercise, "once") if new == 1 gen workout2 = regexm(exercise, "twice") if new == 1 gen workout3 = regexm(exercise, "three") if new == 1 gen workout4 = regexm(exercise, "Four") if new == 1 gen workout = 0 if new == 1 replace workout = workout1 + 2*workout2 + 3*workout3 + 4*workout4 drop workout1 workout2 workout3 workout4 label variable workout "Workout freq." * Daylight Savings Time Change March 14, 2010 *********************************************** gen dst2010 = 0 if new == 1 replace dst2010 = 1 if session == "100315_1550" | session == "100316_1538" label variable dst2010 "Daylight Sav. Time" * Psychological variables ************************* * Personality replace personality = "optimistic" if personality == "Optimistic" gen optimistic = regexm(personality, "optimistic") label variable optimistic "Optimistic Personality" replace personality = "pessimistic" if personality == "Pessimistic" gen pessimistic = regexm(personality, "pessimistic") label variable optimistic "Pessimistic Personality" replace personality = "neither" if personality == "Neither" gen neither = regexm(personality, "neither") label variable neither "Neither Optimistic/Pessimist" * Emotions gen anger = regexm(emotions, "anger") gen anger1 = regexm(emotions, "Anger") replace anger = 1 if anger1 == 1 drop anger1 label variable anger "Anger" gen anxiety = regexm(emotions, "anxiety") gen anxiety1 = regexm(emotions, "Anxiety") replace anxiety = 1 if anxiety1 == 1 drop anxiety1 label variable anxiety "Anxiety" gen confusion = regexm(emotions, "confusion") gen confusion1 = regexm(emotions, "Confusion") replace confusion = 1 if confusion1 == 1 drop confusion1 label variable confusion "Confusion" gen contentment = regexm(emotions, "contentment") gen contentment1 = regexm(emotions, "Contentment") replace contentment = 1 if contentment1 == 1 drop contentment1 label variable contentment "Contentment" gen fatigue = regexm(emotions, "fatigue") gen fatigue1 = regexm(emotions, "Fatigue") replace fatigue = 1 if fatigue1 == 1 drop fatigue1 label variable fatigue "Fatique" gen happiness = regexm(emotions, "happiness") gen happiness1 = regexm(emotions, "Happiness") replace happiness = 1 if happiness1 == 1 drop happiness1 label variable happiness "Happiness" gen irritation = regexm(emotions, "irritation") gen irritation1 = regexm(emotions, "Irritation") replace irritation = 1 if irritation1 == 1 drop irritation1 label variable irritation "Irritation" gen moodswings = regexm(emotions, "mood swings") gen moodswings1 = regexm(emotions, "Mood swings") replace moodswings = 1 if moodswings1 == 1 drop moodswings1 label variable moodswings "Mood Swings" gen withdrawal = regexm(emotions, "withdrawl") gen withdrawal1 = regexm(emotions, "Withdrawl") replace withdrawal = 1 if withdrawal1 == 1 drop withdrawal1 label variable withdrawal "Withdrawal" *************************** * Auction related variables *************************** label variable v "Value" label variable b "Bid" label variable win "Won round" label variable b_adj "Bid + U(0,1) to break ties" label variable bh "The Winning Bid" label variable disttype "1=high 2=low distribution" drop v21 v22 v23 v24 v25 v26 gen v2 = v^2 label variable v2 "Value Squared" gen v3 = v^3 label variable v3 "Value Cubed" gen period2 = period^2 label variable period2 "Period Squared" gen period3 = period^2 label variable period3 "Period Cubed" gen learning = 0 replace learning = 1 / ln(1 + period) if period > 0 gen dom = 0 replace dom = 1 if b > v * weakly dominated bid * PERIOD DUMMIES tab period, gen(pdum) * MOMENTS OF INDIVIDUAL EMPIRICAL VALUE DISTRIBUTIONS egen meanv = mean(v) if period > 0, by(id) egen sdv = sd(v) if period > 0, by(id) egen skewv = skew(v) if period > 0, by(id) * GENERATE DOLLAR PROFITS gen dollarprofit = totalprofit/13 label variable dollarprofit "Total Profit in US Dollars" * Markup (markdown) gen markup = b - v egen meanmarkup = mean(markup), by(id) * Overbidding Risk Neutral Nash Equilibrium gen ornne = b - (v/2) * High value gen hv = 0 replace hv = 1 if v > 50 * Low value gen lv = 1 - hv * Early bidding periods gen fperiod = 0 replace fperiod = 1 if period > 0 & period < 16 * Later bidding periods gen lperiod = 0 replace lperiod = 1 if fperiod == 0 & period > 0 gen datum = substr(session, 1, 6) destring(datum), replace label variable datum "Session Date" * generate session dummies tab datum, gen(sessiondum) ************************* * 2D4D related variables ************************* * mensdr sample (2007 data set) * GENERATE DIGIRATIO AVERAGE replace fourd_mp = fourd_mp/11.811 label variable fourd_mp "Fourth Digit Length in mm, Matthew" replace twod_mp = twod_mp/11.811 label variable twod_mp "Second Digit Length in mm, Matthew" replace fourd_bcs = fourd_bcs/11.811 label variable fourd_bcs "Fourth Digit Length in mm, Burkhard" replace twod_bcs = twod_bcs/11.811 label variable twod_bcs "Second Digit Length in mm, Burkhard" gen dr_mp = twod_mp/fourd_mp label variable dr_mp "Matthew's DR measurement" gen dr_bcs = twod_bcs/fourd_bcs label variable dr_bcs "Burkhard's DR Measurement" gen avdr = (dr_bcs + dr_mp)/2 label variable avdr "Average digit ratio" gen avdrfem = avdr*female label variable avdrfem "Av. Digitratio*Female" gen twod = (twod_mp + twod_bcs)/2 gen fourd = (fourd_mp + fourd_bcs)/2 label variable twod "Second Digit, mm" label variable fourd "Fourth Digit, mm" gen dr = twod/fourd label variable dr "Digit Ratio" gen drfem = dr*female la var drfem "Digit Ratio*Female" gen dr2 = dr^2 label variable dr2 "Digit Ratio Squared" * only 2010 data set gen broken = regexm(brokenfinger, "Yes") if new == 1 label variable brokenf "Broken Finger" gen broken2d4d = regexm(_4_digit, "Yes") if new == 1 label variable broken2d4d "Broken 2D:4D" gen unclear2 = 0 replace unclear2 = 1 if unclear2_bcs == 1 | unclear2_mp == 1 gen unclear4 = 0 replace unclear4 = 1 if unclear4_bcs == 1 | unclear4_mp == 1 gen uncleardr = 0 replace uncleardr = 1 if unclear2 == 1 | unclear4 == 1 label variable uncleardr "DR difficult to measure" gen drfperiod = dr*fperiod label variable drfperiod "Digit Ratio and Early Bids" gen drlperiod = dr*lperiod label variable drlperiod "Digit Ratio and Late Bids" gen cagedr = cage*dr gen csibdr = csib*dr gen drhv = dr*hv label var drhv "Digit Ratio and High Values" gen drlv = dr*lv label var drlv "Digit Ratio and Low Values" gen femaledr = dr*female label variable femaledr "Digit Ratio and Female" gen maledr = dr*(1 - female) label variable maledr "Digit Ratio and Male" gen mathdr = math*dr gen sciotherscidr = sciothersci*dr gen econdr = econ*dr gen socscidr = socsci*dr gen humanitiesdr = humanities*dr gen asiandr = asian*dr gen otherethndr = otherethn*dr ****************************** * Menstrual cycle & fertility ****************************** * Contraceptives gen pill = regexm(contraceptive,"Yes") *(2009 data) gen pill1 = regexm(on_pill,"Yes") *(2007 data) replace pill = 1 if pill == 1 | pill1 == 1 drop pill1 label variable pill "Contraceptives" * PMS replace symptoms_pms = "none" if symptoms_pms == "None" gen none = regexm(symptoms_pms,"none") label variable none "No PMS" replace symptoms_pms = "mild" if symptoms_pms == "Mild" gen mild = regexm(symptoms_pms,"mild") label variable mild "Mild PMS" replace symptoms_pms = "severe" if symptoms_pms == "Severe" gen severe = regexm(symptoms_pms,"severe") label variable severe "Severe PMS" replace usual_pms = "none" if usual_pms == "None" gen usual_none = regexm(usual_pms, "none") if new == 1 label variable none "Usually no PMS" replace usual_pms = "mild" if usual_pms == "Mild" gen usual_mild = regexm(usual_pms, "mild") if new == 1 label variable mild "Usually mild PMS" replace usual_pms = "severe" if usual_pms == "Severe" gen usual_severe = regexm(usual_pms, "none") if new == 1 label variable severe "Usually severe PMS" * there is only one participant with severe PMS: combine with mild gen pms = 0 replace pms = 1 if severe == 1 | mild == 1 label variable pms "Currently PMS" gen upms = 0 replace upms = 1 if usual_mild == 1 | usual_severe == 1 label variable upms "Usually PMS" *Generate standardized phase dummies of the menstrual cycle rename men_awaydays mens *label variable mens "Days from 1st day of last period" la var mens "Day of Menstrual Cycle" gen menst=0 gen folic=0 gen peri=0 gen lut=0 gen premens=0 replace menst=1 if mens<6 & female==1 replace folic=1 if 535" replace duration = "24" if duration=="less than 25" | duration=="< 25" destring duration, replace label variable duration "Length Menstrual Cycle" gen mens_adj = 0 label variable mens_adj "Mens/Duration" replace mens_adj = mens/duration if female==1 gen mensttest = 5.5/28 gen folictest = 13.5/28 gen peritest = 16.5/28 gen luttest = 23.5/28 replace menst_adj = 1 if mens_adj<=mensttest & female==1 replace folic_adj =1 if mens_adj>mensttest & mens_adj<=folictest replace peri_adj =1 if mens_adj>folictest & mens_adj<=peritest replace lut_adj =1 if mens_adj>peritest & mens_adj<=luttest replace premens_adj = 1 if mens_adj>luttest drop mensttest folictest peritest luttest *GENERATE folicular phase adjusted menstrual phases gen premensfoltest = 0 gen lutfoltest = 0 gen perifoltest = 0 replace premensfoltest = duration - 5 replace lutfoltest = duration - 13 replace perifoltest = duration - 16 gen menst_fol = 0 label variable menst_fol "Fol. Adj. Menstrual Phase" gen folic_fol = 0 label variable folic_fol "Fol. Adj. Follicular Phase" gen peri_fol = 0 label variable peri_fol "Fol. Adj. Peri-ovular Phase" gen lut_fol = 0 label variable lut_fol "Fol. Adj. Luteal Phase" gen premens_fol = 0 label variable premens_fol "Fol. Adj. Premenstrual Phase" replace premens_fol = 1 if mens > premensfoltest & female==1 replace lut_fol = 1 if mens > lutfoltest & premens_fol == 0 & female==1 replace peri_fol = 1 if mens > perifoltest & lut_fol == 0 & premens_fol == 0 & female==1 replace menst_fol = 1 if menst ==1 & female==1 replace folic_fol = 1 if menst_fol==0 & premens_fol==0 & lut_fol==0 & peri_fol==0 & female==1 drop premensfoltest lutfoltest perifoltest * Women on the pill don't have a natural menstrual cycle. Generate phase dummies for naturally * cycling women gen nmenst = 0 replace nmenst = menst if pill == 0 label variable nmenst "Menstrual Phase" gen nfolic = 0 replace nfolic = folic if pill == 0 label variable nfolic "Follicular Phase" gen nperi = 0 replace nperi = peri if pill == 0 label variable nperi "Peri-ovular Phase" gen nlut = 0 replace nlut = lut if pill == 0 label variable nlut "Luteal Phase" gen npremens = 0 replace npremens = premens if pill == 0 label variable npremens "Premenstrual Phase" gen nmenst_adj = 0 replace nmenst_adj = menst_adj if pill == 0 label variable nmenst_adj "Uni. Adj. Menstrual Phase" gen nfolic_adj = 0 replace nfolic_adj = folic_adj if pill == 0 label variable nfolic_adj "Uni. Adj. Follicular Phase" gen nperi_adj = 0 replace nperi_adj = peri_adj if pill == 0 label variable nperi_adj "Uni. Adj. Peri-ovular Phase" gen nlut_adj = 0 replace nlut_adj = lut_adj if pill == 0 label variable nlut_adj "Uni. Adj. Luteal Phase" gen npremens_adj = 0 replace npremens_adj = premens_adj if pill == 0 label variable npremens_adj "Uni. Adj. Premenstrual Phase" gen nmenst_fol = 0 replace nmenst_fol = menst_fol if pill == 0 label variable nmenst_fol "Fol. Adj. Menstrual Phase" gen nfolic_fol = 0 replace nfolic_fol = folic_fol if pill == 0 label variable nfolic_fol "Fol. Adj. Follicular Phase" gen nperi_fol = 0 replace nperi_fol = peri_fol if pill == 0 label variable nperi_fol "Fol. Adj. Peri-ovular Phase" gen nlut_fol = 0 replace nlut_fol = lut_fol if pill == 0 label variable nlut_fol "Fol. Adj. Luteal Phase" gen npremens_fol = 0 replace npremens_fol = premens_fol if pill == 0 label variable npremens_fol "Fol. Adj. Premenstrual Phase" gen nduration = 0 replace nduration = duration if pill == 0 & female == 1 label variable nduration "Length Menstrual Cycle" * High progesterone phases of the menstrual cycle (second half) gen shalf = 0 replace shalf = 1 if nperi == 1 | nlut == 1 | npremens == 1 gen shalf_adj = 0 replace shalf_adj = 1 if nperi_adj == 1 | nlut_adj == 1 | npremens_adj == 1 gen shalf_fol = 0 replace shalf_fol = 1 if nperi_fol == 1 | nlut_fol == 1 | npremens_fol == 1 * Generate another type high progesterone phase dummy of the menstrual cycle gen hprog = 0 replace hprog = 1 if lut==1 | premens==1 gen hprog_adj = 0 replace hprog_adj = 1 if lut_adj == 1 | premens_adj == 1 gen hprog_fol = 0 replace hprog_fol = 1 if lut == 1 | premens_fol == 1 * Generate "high" estradiol phase dummy of the menstrual cycle gen hest = 0 replace hest = 1 if folic == 1 | peri == 1 | lut == 1 * Alternative definition of follicular phase gen follic = 0 replace follic = 1 if lut == 0 & premens == 0 & female == 1 label variable follic "Follicular Phase" gen nfollic_fol = 0 replace nfollic_fol = 1 if nlut_fol == 0 & npremens_fol == 0 & female == 1 label variable nfollic_fol "Follicular Phase" gen pillorfol = 0 replace pillorfol = 1 if (pill == 1 | nfollic_fol == 1) & new == 1 * Regularity of menstrual cycle gen irregular = regexm(men_changes, "irregular") if new == 1 label variable irregular "Irregular menses" gen regular = regexm(men_changes, "same number") if new == 1 label variable regular "Regular menses" * Two females with natural cycles have duration missing (IDs 113013, 120516) gen durationmissing = 0 replace durationmissing = 1 if duration == . & female == 1 & pill == 0 * One female with natural cycle reports mens = 160 (ID 31138) gen longmens = 0 replace longmens = 1 if mens ~= . & mens > 100 & female == 1 * Calendar keeping gen calendar = regexm(men_calendar, "Yes") if new == 1 label variable calendar "Menstrual Calendar" gen ncalendar = regexm(men_calendar, "No") if new == 1 label variable ncalendar "No Menstrual Calendar" * normalized days in the cycle gen nmens_adj = mens_adj*28 label variable nmens_adj "28-Days Normalized Day of the Menstrual Cycle" * To do: Consider alternative asking for menstrual cycle info * generate specific contraceptives variables ********************************************** gen tmp1 = regexm(bcontrol_type, "DepoProvera") gen tmp2 = regexm(pill_name, "DepoProvera") gen tmp3 = regexm(bcontrol_type, "Seasonale") gen tmp4 = regexm(pill_name, "Seasonale") gen tmp5 = regexm(bcontrol_type, "Seasonellw") gen tmp6 = regexm(pill_name, "Seasonellw") gen threemonthpill = 0 replace threemonthpill = 1 if tmp1 == 1 | tmp2 == 1 | tmp3 == 1 | tmp4 == 1 | tmp5 == 1 | tmp6 == 1 drop tmp1 tmp2 tmp3 tmp4 tmp5 tmp6 * 24 days pills gen yaz1 = regexm(bcontrol_type, "Yaz") gen yaz2 = regexm(pill_name, "yaz") gen yaz = 0 replace yaz = 1 if yaz1 == 1 | yaz2 == 2 drop yaz1 yaz2 * reported name of pill gen rep01 = regexm(bcontrol_type, "DepoProvera") gen rep02 = regexm(pill_name, "DepoProvera") gen rep03 = regexm(bcontrol_type, "Seasonale") gen rep04 = regexm(pill_name, "Seasonale") gen rep05 = regexm(bcontrol_type, "Seasonellw") gen rep06 = regexm(pill_name, "Seasonellw") gen rep07 = regexm(bcontrol_type, "Yaz") gen rep08 = regexm(pill_name, "yaz") gen rep09 = regexm(bcontrol_type, "Microgestin") gen rep10 = regexm(bcontrol_type, "Nuvaring") gen rep11 = regexm(bcontrol_type, "Junel") gen rep12 = regexm(bcontrol_type, "sprintec") gen rep13 = regexm(bcontrol_type, "YAZ") gen rep14 = regexm(bcontrol_type, "Desogen") gen rep15 = regexm(pill_name, "Ortho Tri cyclen Lo") gen rep16 = regexm(pill_name, "Yasmin") gen rep17 = regexm(pill_name, "Avian?") gen rep18 = regexm(pill_name, "Ortho-Tricylcen Lo") gen rep19 = regexm(pill_name, "ortho lo") gen rep20 = regexm(pill_name, "Ortho-tri-cyclen") gen rep21 = regexm(pill_name, "yazmine") gen rep22 = regexm(pill_name, "Yasmin") gen rep23 = regexm(pill_name, "Orthotricyclen") gen reppill = 0 replace reppill = 1 if rep01 + rep02 + rep03 + rep04 + rep05 + rep06 + rep07 + rep08 + rep09 + rep10 + rep11 + rep12 + rep13 + rep14 + rep15 + rep16 + rep17 + rep18 + rep19 + rep20 + rep21 + rep22 + rep23 > 0 drop rep01 rep02 rep03 rep04 rep05 rep06 rep07 rep08 rep09 rep10 rep11 rep12 rep13 rep14 rep15 rep16 rep17 rep18 rep19 rep20 rep21 rep22 rep23 label variable reppill "Name of Pill reported" * we assume that women on 3-month contraceptives are not in the pill break since we have no way of measuring it * and because of the long intake/effect it is more likely that they are not in a break gen pillbreak1 = 1 if yaz == 1 & mens > 24 & mens < 29 & female == 1 gen pillbreak2 = 1 if pill == 1 & mens > 21 & mens < 29 & female == 1 gen pillbreak = 0 replace pillbreak = 1 if pillbreak1 == 1 | pillbreak2 == 1 drop pillbreak1 pillbreak2 label variable pillbreak "Placebo/Break" * doesn't make much sense in our samples given that only a small number of women are on pillbreak gen pilluse = 0 replace pilluse = 1 if pill == 1 & pillbreak == 0 label variable pilluse "Active Ingredient" *GENERATE Probability of Conception (use data from Wilcox et al., 2001) ***************************************************************************** * Average probabilities ("all women" in Wilcox et al., 2001, Table 1, p. 213) gen all_prob = 0.000 replace all_prob = 0.000 if mens==0 | mens==1 | mens==2 replace all_prob = 0.001 if mens==3 replace all_prob = 0.002 if mens==4 replace all_prob = 0.004 if mens==5 replace all_prob = 0.009 if mens==6 replace all_prob = 0.017 if mens==7 replace all_prob = 0.029 if mens==8 replace all_prob = 0.044 if mens==9 replace all_prob = 0.061 if mens==10 replace all_prob = 0.075 if mens==11 replace all_prob = 0.084 if mens==12 replace all_prob = 0.086 if mens==13 replace all_prob = 0.081 if mens==14 replace all_prob = 0.072 if mens==15 replace all_prob = 0.061 if mens==16 replace all_prob = 0.050 if mens==17 replace all_prob = 0.040 if mens==18 replace all_prob = 0.032 if mens==19 replace all_prob = 0.025 if mens==20 replace all_prob = 0.020 if mens==21 replace all_prob = 0.016 if mens==22 replace all_prob = 0.013 if mens==23 | mens==38 replace all_prob = 0.011 if mens==24 | mens==36 replace all_prob = 0.009 if mens==25 | mens==33 | mens==34 replace all_prob = 0.008 if mens==26 | mens==31 | mens==32 replace all_prob = 0.007 if mens==27 | mens==28 | mens==29 | mens==30 replace all_prob = 0.010 if mens==35 replace all_prob = 0.012 if mens==37 replace all_prob = 0.014 if mens==39 replace all_prob = 0.015 if mens==40 | mens==41 | mens==42 | mens==43 | mens==44 | mens==45 | mens==46 | mens==47 | mens==48 | mens==49 | mens==50 * women with regular cycles gen r_prob = 0.000 replace r_prob = 0.001 if mens==3 replace r_prob = 0.002 if mens==4 replace r_prob = 0.004 if mens==5 replace r_prob = 0.009 if mens==6 replace r_prob = 0.018 if mens==7 replace r_prob = 0.032 if mens==8 replace r_prob = 0.050 if mens==9 replace r_prob = 0.069 if mens==10 replace r_prob = 0.085 if mens==11 replace r_prob = 0.094 if mens==12 replace r_prob = 0.093 if mens==13 replace r_prob = 0.085 if mens==14 replace r_prob = 0.073 if mens==15 replace r_prob = 0.059 if mens==16 replace r_prob = 0.047 if mens==17 replace r_prob = 0.036 if mens==18 replace r_prob = 0.028 if mens==19 replace r_prob = 0.021 if mens==20 replace r_prob = 0.016 if mens==21 replace r_prob = 0.013 if mens==22 replace r_prob = 0.010 if mens==23 | mens==36 | mens==37 replace r_prob = 0.008 if mens==24 | mens==34 replace r_prob = 0.007 if mens==25 | mens==32 | mens==33 replace r_prob = 0.006 if mens==26 | mens==30 | mens==31 replace r_prob = 0.005 if mens==27 | mens==28 | mens==29 replace r_prob = 0.009 if mens==35 replace r_prob = 0.011 if mens==38 replace r_prob = 0.012 if mens==39 replace r_prob = 0.013 if mens==40 | mens==41 | mens==42 | mens==43 | mens==44 | mens==45 | mens==46 | mens==47 | mens==48 | mens==49 | mens==50 * women with irregular cycles gen ir_prob = 0.000 replace ir_prob = 0.001 if mens==4 | mens==5 replace ir_prob = 0.002 if mens==6 replace ir_prob = 0.004 if mens==7 replace ir_prob = 0.007 if mens==8 replace ir_prob = 0.011 if mens==9 replace ir_prob = 0.018 if mens==10 replace ir_prob = 0.026 if mens==11 replace ir_prob = 0.035 if mens==12 replace ir_prob = 0.045 if mens==13 replace ir_prob = 0.055 if mens==14 replace ir_prob = 0.061 if mens==15 replace ir_prob = 0.065 if mens==16 | mens==17 replace ir_prob = 0.062 if mens==18 replace ir_prob = 0.057 if mens==19 replace ir_prob = 0.051 if mens==20 replace ir_prob = 0.044 if mens==21 replace ir_prob = 0.038 if mens==22 replace ir_prob = 0.032 if mens==23 replace ir_prob = 0.027 if mens==24 replace ir_prob = 0.023 if mens==25 replace ir_prob = 0.020 if mens==26 replace ir_prob = 0.017 if mens==27 replace ir_prob = 0.015 if mens==28 replace ir_prob = 0.013 if mens==29 | mens==36 | mens==37 replace ir_prob = 0.012 if mens==30 | mens==35 replace ir_prob = 0.011 if mens==31 | mens==32 | mens==33 | mens==34 replace ir_prob = 0.014 if mens==38 | mens==39 | mens==40 | mens==41 | mens==42 | mens==43 | mens==44 | mens==45 | mens==46 | mens==47 | mens==48 | mens==49 | mens==50 * OVERALL PROBABILITY OF CONCEPTION DEPENDING ON THE REGULARITY OF THE CYCLES (as in mens paper) gen conception = 0 label variable conception "Prob. of Conception" replace conception = r_prob replace conception = ir_prob if mens > duration | mens > 40 | duration > 40 gen conceptionpill = conception*pill * conception prob for nat. cycling women only gen nconception = 0 replace nconception = conception if pill == 0 *********************** * Hormones Data *********************** * Label hormone variables label variable testf "Ex ante Testosterone" label variable estf "Ex ante Estradiol" label variable progf "Ex ante Progesterone" label variable cortf "Ex ante Cortisol" label variable tests "Ex post Testosterone" label variable ests "Ex post Estradiol" label variable progs "Ex post Progesterone" label variable corts "Ex post Cortisol" label variable testavg "Avg. Testosterone" label variable estavg "Avg. Estradiol" label variable progavg "Avg. Progesterone" label variable cortavg "Avg. Cortisol" label variable sdtestf "SD EA Testosterone" label variable sdestf "SD EA Estradiol" label variable sdprogf "SD EA Progesterone" label variable sdcortf "SD EA Cortisol" label variable sdtests "SD EP Testosterone" label variable sdests "SD EP Estradiol" label variable sdprogs "SD EP Progesterone" label variable sdcorts "SD EP Cortisol" label variable sdtestavg "SD Testosterone" label variable sdestavg "SD Estradiol" label variable sdprogavg "SD Progesterone" label variable sdcortavg "SD Cortisol" * Normalized hormones (by standard deviation) gen ntestf = testf/sdtestf label variable ntestf "Testosterone" gen nestf = estf/sdestf label variable nestf "Estradiol" gen nprogf = progf/sdprogf label variable nprogf "Progesterone" gen ncortf = cortf/sdcortf label variable ncortf "Cortisol" gen ntests = tests/sdtests label variable ntests "Testosterone" gen nests = ests/sdests label variable nests "Estradiol" gen nprogs = progs/sdprogs label variable nprogs "Progesterone" gen ncorts = corts/sdcorts label variable ncorts "Cortisol" gen ntestavg = testavg/sdtestavg label variable ntestavg "Testosterone" gen nestavg = estavg/sdestavg label variable nestavg "Estradiol" gen nprogavg = progavg/sdprogavg label variable nprogavg "Progesterone" gen ncortavg = cortavg/sdcortavg label variable ncortavg "Cortisol" * Interaction terms gen femtestf = female*ntestf gen femestf = female*nestf gen femprogf = female*nprogf gen femcortf = female*ncortf gen maletestf = (1 - female)*ntestf gen maleestf = (1 - female)*nestf gen maleprogf = (1 - female)*nprogf gen malecortf = (1 - female)*ncortf gen menstprogf = menst*progf gen folicprogf = folic*progf gen periprogf = peri*progf gen lutprogf = lut*progf gen premensprogf = premens*progf gen progmensprogf = hprog*progf gen progmensprogf_adj = hprog_adj*progf gen progmensprogf_fol = hprog_fol*progf * Assay quality variables * See original assay results from CNPRC for brief explanation of what is LDD. * Essentially LDD refers to score at the lower limit of the assay. label variable limittestf "LDD Ex Ante Testosterone" label variable limitestf "LDD Ex Ante Estradiol" label variable limitprogf "LDD Ex Ante Progesterone" label variable limitcortf "LDD Ex Ante Cortisol" label variable limittests "LDD Ex Post Testosterone" label variable limitests "LDD Ex Post Estradiol" label variable limitprogs "LDD Ex Post Progesterone" label variable limitcorts "LDD Ex Post Cortisol" * Interaction terms gen limtestf = limittestf*ntestf if new == 1 label variable limtestf "LDD x Testost." gen limestf = limitestf*nestf if new == 1 label variable limestf "LDD x Estrad." gen limprogf = limitprogf*nprogf if new == 1 label variable limprogf "LDD x Progest." gen limcortf = limitcortf*ncortf if new == 1 label variable limcortf "LDD x Cort." gen follicprogf = follic*nprogf if new == 1 label variable follicprogf "Foll. Ph. x Progest." gen dbmiprogf = dbmi*nprogf if new == 1 label variable dbmiprogf "Diff. Standard BMI x Progest." gen obeseprogf = obese*nprogf if new == 1 label variable obeseprogf "Obese x Progest." gen pilltestf = pill*ntestf if new == 1 label variable pilltestf "Contracept. x Testost." gen pillestf = pill*nestf if new == 1 label variable pillestf "Contracept. x Estradiol" gen pillprogf = pill*nprogf if new == 1 label variable pillprogf "Contracept. x Progest." gen smokeprogf = smoking*nprogf if new == 1 label variable smokeprogf "Smoking x Progest." gen soyprogf = soybean*nprogf if new == 1 label variable soyprogf "Soybean x Progest." gen workoutprogf = workout*nprogf if new == 1 label variable workoutprogf "Workout x Progest." gen homoprogf = homobi*nprogf if new == 1 label variable homoprogf "Homosex. x Progest." * Lunch before saliva collection gen skipped = regexm(lunch, "skipped") if new == 1 label variable skipped "Skipped lunch" gen lunch11 = regexm(lunch, "11.00 am") if new == 1 gen lunch12 = regexm(lunch, "12.00 pm") if new == 1 gen lunch13 = regexm(lunch, "1.00 pm") if new == 1 gen lunch14 = regexm(lunch, "2.00 pm") if new == 1 gen lunch15 = regexm(lunch, "3.00 pm") if new == 1 gen luncht = lunch11 + 2*lunch12 + 3*lunch13 + 4*lunch14 + 5*lunch15 if new == 1 * measures the contamination of saliva with food label variable luncht "When lunch today" gen lunchtcortf = luncht*ncortf label variable lunchtcortf "Lunch today x Cort." gen lunchcortf = (1 - skipped)*ncortf label variable lunchcortf "Lunch today x Cort." gen lunchtestf = (1 - skipped)*ntestf label variable lunchtestf "Lunch today x Testost." * Last food intake gen food30 = regexm(time_last_ate, "30 min") if new == 1 gen food60 = regexm(time_last_ate, "1") if new == 1 gen food120 = regexm(time_last_ate, "2") if new == 1 gen food180 = regexm(time_last_ate, "3 hours") if new == 1 gen food240b = regexm(time_last_ate, "4") if new == 1 gen food300 = regexm(time_last_ate, "ago") if new == 1 gen food240 = 0 if new == 1 replace food240 = food240b if food300 == 0 & new == 1 gen foodtime = 30*food30 + 60*food60 + 120*food120 + 180*food180 + 240*food240 + 300*food300 if new == 1 * measures the minutes ago when last eaten label variable foodtime "Time last eaten" gen drink30 = regexm(drinks_when, "30 min") if new == 1 gen drink60 = regexm(drinks_when, "1") if new == 1 gen drink90 = regexm(drinks_when, "1.5") if new == 1 gen drink120 = regexm(drinks_when, "2 hours") if new == 1 * 3 subjects stated drink60 AND drink90. In this case we set drink90 to zero since * the lowest matters for the quality of the assays. replace drink90 = 0 if drink60 == 1 gen drinktime = 30*drink30 + 60*drink60 + 90*drink90 + 120*drink120 if new == 1 * measures the minutes ago when last eaten label variable drinktime "Time last drank" * Mouth injuries gen mouth = regexm(mouth_injuries, "Yes") if new == 1 label variable mouth "Mouth injuries" * Teeth brushing gen brush1 = regexm(brushperday, "Once") if new == 1 gen brush2 = regexm(brushperday, "Twice") if new == 1 gen brush3 = regexm(brushperday, "Three") if new == 1 gen brush4 = regexm(brushperday, "More") if new == 1 gen brush = brush1 + 2*brush2 + 3*brush3 + 4*brush4 if new == 1 label variable brush "Freq. toothbrush." gen blast30 = regexm(brush_last, "30") if new == 1 gen blast60 = regexm(brush_last, "1 hour") if new == 1 gen blast120 = regexm(brush_last, "2 hours") if new == 1 gen blast180 = regexm(brush_last, "3 hours") if new == 1 gen blast240 = regexm(brush_last, "More") if new == 1 * Some subjects brushed several times. We look for the most recent time since * this is what matters for the assays. replace blast60 = 0 if blast30 ==1 replace blast120 = 0 if blast60 == 1 | blast30 == 1 replace blast180 = 0 if blast120 == 1 | blast60 == 1 | blast30 == 1 replace blast240 = 0 if blast180 == 1 | blast120 == 1 | blast60 == 1 | blast30 == 1 gen blast = 30*blast30 + 60*blast60 + 120*blast120 + 180*blast180 + 240*blast240 if new == 1 label variable blast "Time last toothbrush." gen pillorfolprogf = pillorfol*nprogf label variable pillorfolprogf "Contraception or Follicular x Prog." * Generate quasi z scores ****************************** * 1. Since standard deviation differs for males and females, we create gender specific quasi z scores. * 2. Because of well-known problems with collinearity in interaction terms (because there may be no observation with z score of null), * we use the observation (of the same gender) closest to the mean instead of the mean in defining the gender specific z score. * males sum testf if period == 30 & female == 0 gen meantestf_m = r(mean) gen sdtestf_m = r(sd) gen dtestf_m = meantestf_m - testf if female == 0 gen absdtestf_m = abs(dtestf_m) if female == 0 egen minabsdtestf_m = min(absdtestf_m) gen help1 = 0 if female == 0 & new == 1 replace help1 = 1 if testf + minabsdtestf_m == meantestf_m | testf - minabsdtestf_m == meantestf_m gen help2 = 0 if female == 0 & new == 1 replace help2 = testf if help1 == 1 sum help2 if period == 30 & female == 0 gen nbymeantestf_m = r(max) gen testfz_m = (testf - nbymeantestf_m)/sdtestf_m drop absdtestf_m minabsdtestf_m help1 help2 sum testfz_m if female == 0 & period == 30 label variable testfz_m "Testosterone" sum estf if period == 30 & female == 0 gen meanestf_m = r(mean) gen sdestf_m = r(sd) gen destf_m = meanestf_m - estf if female == 0 gen absdestf_m = abs(destf_m) if female == 0 egen minabsdestf_m = min(absdestf_m) gen help1 = 0 if female == 0 & new == 1 replace help1 = 1 if estf + minabsdestf_m == meanestf_m | estf - minabsdestf_m == meanestf_m gen help2 = 0 if female == 0 & new == 1 replace help2 = estf if help1 == 1 sum help2 if period == 30 & female == 0 gen nbymeanestf_m = r(max) gen estfz_m = (estf - nbymeanestf_m)/sdestf_m drop absdestf_m minabsdestf_m help1 help2 sum estfz_m if female == 0 & period == 30 label variable estfz_m "Estradiol" sum progf if period == 30 & female == 0 gen meanprogf_m = r(mean) gen sdprogf_m = r(sd) gen dprogf_m = meanprogf_m - progf if female == 0 gen absdprogf_m = abs(dprogf_m) if female == 0 egen minabsdprogf_m = min(absdprogf_m) gen help1 = 0 if female == 0 & new == 1 replace help1 = 1 if progf + minabsdprogf_m == meanprogf_m | progf - minabsdprogf_m == meanprogf_m gen help2 = 0 if female == 0 & new == 1 replace help2 = progf if help1 == 1 sum help2 if period == 30 & female == 0 gen nbymeanprogf_m = r(max) gen progfz_m = (progf - nbymeanprogf_m)/sdprogf_m drop absdprogf_m minabsdprogf_m help1 help2 sum progfz_m if female == 0 & period == 30 label variable progfz_m "Progesterone" sum cortf if period == 30 & female == 0 gen meancortf_m = r(mean) gen sdcortf_m = r(sd) gen dcortf_m = meancortf_m - cortf if female == 0 gen absdcortf_m = abs(dcortf_m) if female == 0 egen minabsdcortf_m = min(absdcortf_m) gen help1 = 0 if female == 0 & new == 1 replace help1 = 1 if cortf + minabsdcortf_m == meancortf_m | cortf - minabsdcortf_m == meancortf_m gen help2 = 0 if female == 0 & new == 1 replace help2 = cortf if help1 == 1 sum help2 if period == 30 & female == 0 gen nbymeancortf_m = r(max) gen cortfz_m = (cortf - nbymeancortf_m)/sdcortf_m drop absdcortf_m minabsdcortf_m help1 help2 sum cortfz_m if female == 0 & period == 30 label variable cortfz_m "Cortisol" * females sum testf if period == 30 & female == 1 gen meantestf_f = r(mean) gen sdtestf_f = r(sd) gen dtestf_f = meantestf_f - testf if female == 1 gen absdtestf_f = abs(dtestf_f) if female == 1 egen minabsdtestf_f = min(absdtestf_f) gen help1 = 0 if female == 1 & new == 1 replace help1 = 1 if testf + minabsdtestf_f == meantestf_f | testf - minabsdtestf_f == meantestf_f gen help2 = 0 if female == 1 & new == 1 replace help2 = testf if help1 == 1 sum help2 if period == 30 & female == 1 gen nbymeantestf_f = r(max) gen testfz_f = (testf - nbymeantestf_f)/sdtestf_f drop absdtestf_f minabsdtestf_f help1 help2 sum testfz_f if female == 1 & period == 30 label variable testfz_f "Testosterone" sum estf if period == 30 & female == 1 gen meanestf_f = r(mean) gen sdestf_f = r(sd) gen destf_f = meanestf_f - estf if female == 1 gen absdestf_f = abs(destf_f) if female == 1 egen minabsdestf_f = min(absdestf_f) gen help1 = 0 if female == 1 & new == 1 replace help1 = 1 if estf + minabsdestf_f == meanestf_f | estf - minabsdestf_f == meanestf_f gen help2 = 0 if female == 1 & new == 1 replace help2 = estf if help1 == 1 sum help2 if period == 30 & female == 1 gen nbymeanestf_f = r(max) gen estfz_f = (estf - nbymeanestf_f)/sdestf_f drop absdestf_f minabsdestf_f help1 help2 sum estfz_f if female == 1 & period == 30 label variable estfz_f "Estradiol" sum progf if period == 30 & female == 1 gen meanprogf_f = r(mean) gen sdprogf_f = r(sd) gen dprogf_f = meanprogf_f - progf if female == 1 gen absdprogf_f = abs(dprogf_f) if female == 1 egen minabsdprogf_f = min(absdprogf_f) gen help1 = 0 if female == 1 & new == 1 replace help1 = 1 if progf + minabsdprogf_f == meanprogf_f | progf - minabsdprogf_f == meanprogf_f gen help2 = 0 if female == 1 & new == 1 replace help2 = progf if help1 == 1 sum help2 if period == 30 & female == 1 gen nbymeanprogf_f = r(max) gen progfz_f = (progf - nbymeanprogf_f)/sdprogf_f drop absdprogf_f minabsdprogf_f help1 help2 sum progfz_f if female == 1 & period == 30 label variable progfz_f "Progesterone" sum cortf if period == 30 & female == 1 gen meancortf_f = r(mean) gen sdcortf_f = r(sd) gen dcortf_f = meancortf_f - cortf if female == 1 gen absdcortf_f = abs(dcortf_f) if female == 1 egen minabsdcortf_f = min(absdcortf_f) gen help1 = 0 if female == 1 & new == 1 replace help1 = 1 if cortf + minabsdcortf_f == meancortf_f | cortf - minabsdcortf_f == meancortf_f gen help2 = 0 if female == 1 & new == 1 replace help2 = cortf if help1 == 1 sum help2 if period == 30 & female == 1 gen nbymeancortf_f = r(max) gen cortfz_f = (cortf - nbymeancortf_f)/sdcortf_f drop absdcortf_f minabsdcortf_f help1 help2 sum cortfz_f if female == 1 & period == 30 label variable cortfz_f "Cortisol" * Computing real empirical z-scores (for risk paper) ***************************************************** * males gen ztestf_m = (testf - meantestf_m)/sdtestf_m label variable ztestf_m "Testosterone" gen zestf_m = (estf - meanestf_m)/sdestf_m label variable zestf_m "Estradiol" gen zprogf_m = (progf - meanprogf_m)/sdprogf_m label variable zprogf_m "Progesterone" gen zcortf_m = (cortf - meancortf_m)/sdcortf_m label variable zcortf_m "Cortisol" * females gen ztestf_f = (testf - meantestf_f)/sdtestf_f label variable ztestf_f "Testosterone" gen zestf_f = (estf - meanestf_f)/sdestf_f label variable zestf_f "Estradiol" gen zprogf_f = (progf - meanprogf_f)/sdprogf_f label variable zprogf_f "Progesterone" gen zcortf_f = (cortf - meancortf_f)/sdcortf_f label variable zcortf_f "Cortisol" * Interaction terms gen notpill = 1 - pill label variable notpill "Naturally Cycling" gen male = 1 - female label variable male "Male" gen femalepill = female*pill label variable femalepill "Female x Contraceptives" gen femalenotpill = female*notpill label variable femalenotpill "Female x Naturally Cycling" gen malepill = male*pill * collinear with male gen malenotpill = male*notpill gen pilltestfz_f = pill*testfz_f gen pillestfz_f = pill*estfz_f gen pillprogfz_f = pill*progfz_f gen pillcortfz_f = pill*cortfz_f label variable pilltestfz_f "Contracept. x Testost." label variable pillestfz_f "Contracept. x Estradiol" label variable pillprogfz_f "Contracept. x Progest." label variable pillcortfz_f "Contracept. x Cortisol" gen notpilltestfz_f = notpill*testfz_f gen notpillestfz_f = notpill*estfz_f gen notpillprogfz_f = notpill*progfz_f gen notpillcortfz_f = notpill*cortfz_f label variable notpilltestfz_f "Nat. Cycling x Testost." label variable notpillestfz_f "Nat. Cycling x Estradiol" label variable notpillprogfz_f "Nat. Cycling x Progest." label variable notpillcortfz_f "Nat. Cycling x Cortisol" gen maletestfz_f = male*testfz_f gen maleestfz_f = male*estfz_f gen maleprogfz_f = male*progfz_f gen malecortfz_f = male*cortfz_f label variable maletestfz_f "Male x Testost." label variable maleestfz_f "Male x Estradiol" label variable maleprogfz_f "Male x Progest." label variable malecortfz_f "Male x Cortisol" gen femaletestfz_m = female*testfz_m gen femaleestfz_m = female*estfz_m gen femaleprogfz_m = female*progfz_m gen femalecortfz_m = female*cortfz_m label variable femaletestfz_m "Female x Testost." label variable femaleestfz_m "Female x Estradiol" label variable femaleprogfz_m "Female x Progest." label variable femalecortfz_m "Female x Cortisol" gen malepilltestfz_f = male*pill*testfz_f gen malepillestfz_f = male*pill*estfz_f gen malepillprogfz_f = male*pill*progfz_f gen malepillcortfz_f = male*pill*cortfz_f * Note: collinear with male * Hormone response *************************** gen testd = tests - testf label variable testd "Testosterone Response" gen estd = ests - estf label variable testd "Estradio Response" gen progd = progs - progf label variable testd "Progesterone Response" gen cortd = corts - cortf label variable testd "Cortisol Response" gen reltestd = (tests - testf)/testf gen relestd = (ests - estf)/estf gen relprogd = (progs - progf)/progf gen relcortd = (corts - cortf)/cortf *************************** * Risk variables *************************** gen rg1 = regexm(g1, "A") replace rg1 = . if new == 0 gen rg2 = regexm(g2, "A") replace rg2 = . if new == 0 gen rg3 = regexm(g3, "A") replace rg3 = . if new == 0 gen rg4 = regexm(g4, "A") replace rg4 = . if new == 0 gen rg5 = regexm(g5, "A") replace rg5 = . if new == 0 gen rg6 = regexm(g6, "A") replace rg6 = . if new == 0 gen rg7 = regexm(g7, "A") replace rg7 = . if new == 0 gen rg8 = regexm(g8, "A") replace rg8 = . if new == 0 gen rg9 = regexm(g9, "A") replace rg9 = . if new == 0 gen rg10 = regexm(g10, "A") replace rg10 = . if new == 0 gen rgconsistent = 1 replace rgconsistent = 0 if (rg2 > rg1 | rg3 > rg2 | rg4 > rg3 | rg5 > rg4 | rg6 > rg5 | rg7 > rg6 | rg8 > rg9 | rg10 > rg9) replace rgconsistent = . if new == 0 * regconsistent is zero if choice in the gain domain are inconsistent gen riskgain = rg1 + rg2 + rg3 + rg4 + rg5 + rg6 + rg7 + rg8 + rg9 + rg10 label variable riskgain "Risk Aversion Gains" gen rgxincon = riskgain*(1 - rgconsistent) label variable rgxincon "Inaccessible x Risk Aversion Gains" * interaction variable to be used jointly with riskgain to measure the effect of * risk aversion on the gain domain if subject is not inconsistent gen rginconsistent = 1 - rgconsistent label variable rginconsistent "Inaccessible Risk Pref. Gains" gen rl1 = regexm(l1, "A") replace rl1 = . if new == 0 gen rl2 = regexm(l2, "A") replace rl2 = . if new == 0 gen rl3 = regexm(l3, "A") replace rl3 = . if new == 0 gen rl4 = regexm(l4, "A") replace rl4 = . if new == 0 gen rl5 = regexm(l5, "A") replace rl5 = . if new == 0 gen rl6 = regexm(l6, "A") replace rl6 = . if new == 0 gen rl7 = regexm(l7, "A") replace rl7 = . if new == 0 gen rl8 = regexm(l8, "A") replace rl8 = . if new == 0 gen rl9 = regexm(l9, "A") replace rl9 = . if new == 0 gen rl10 = regexm(l10, "A") replace rl10 = . if new == 0 gen rlconsistent = 1 replace rlconsistent = 0 if (rl2 < rl1 | rl3 < rl2 | rl4 < rl3 | rl5 < rl4 | rl6 < rl5 | rl7 < rl6 | rl8 < rl9 | rl10 < rl9) replace rlconsistent = . if new == 0 * rlconsistent is zero if choices in the loss domain are inconsistent gen riskloss = rl1 + rl2 + rl3 + rl4 + rl5 + rl6 + rl7 + rl8 + rl9 + rl10 label variable riskloss "Risk Aversion Losses" gen rlxincon = riskloss*(1 - rlconsistent) label variable rlxincon "Inaccessible x Risk Aversion Losses" * interaction variable to be used jointly with riskloss to measure the effect of * risk aversion on the loss domain if subject is not inconsistent gen rlinconsistent = 1 - rlconsistent label variable rlinconsistent "Inaccessible Risk Pref. Losses" gen consistent = 0 replace consistent = 1 if rgconsistent == 1 & rlconsistent == 1 replace consistent = . if new == 0 gen girrational = 0 replace girrational = 1 if rg1 == 0 gen lirrational = 0 replace lirrational = 1 if rl1 == 1 gen irrational = 0 replace irrational = 1 if rg1 == 0 | rl1 == 1 gen rgaverse = 0 replace rgaverse = 1 if riskgain > 5 & rgconsistent == 1 & girrational == 0 gen rlaverse = 0 replace rlaverse = 1 if riskloss > 5 & rlconsistent == 1 & lirrational == 0 gen rgneutral = 0 replace rgneutral = 1 if riskgain == 5 & rgconsistent == 1 & girrational == 0 gen rlneutral = 0 replace rlneutral = 1 if riskloss == 5 & rlconsistent == 1 & lirrational == 0 gen rgseeking = 0 replace rgseeking = 1 if riskgain < 5 & rgconsistent == 1 & girrational == 0 gen rlseeking = 0 replace rlseeking = 1 if riskloss < 5 & rlconsistent == 1 & lirrational == 0 gen rgmulti = 0 replace rgmulti = 1 if rgaverse == 1 replace rgmulti = 0 if rgneutral == 1 replace rgmulti = -1 if rgseeking == 1 gen rlmulti = 0 replace rlmulti = 1 if rlaverse == 1 replace rlmulti = 0 if rlneutral == 1 replace rlmulti = -1 if rlseeking == 1 gen rmulti = 0 replace rmulti = 1 if rgaverse == 1 & rlseeking == 1 replace rmulti = 2 if rgaverse == 1 & rlneutral == 1 replace rmulti = 3 if rgaverse == 1 & rlaverse == 1 replace rmulti = 4 if rgneutral == 1 & rlseeking == 1 replace rmulti = 5 if rgneutral == 1 & rlneutral == 1 replace rmulti = 6 if rgneutral == 1 & rlaverse == 1 replace rmulti = 7 if rgseeking == 1 & rlseeking == 1 replace rmulti = 8 if rgseeking == 1 & rlneutral == 1 replace rmulti = 9 if rgseeking == 1 & rlaverse == 1 replace rmulti = . if consistent == 0 | irrational == 1 gen reflection = 0 replace reflection = 1 if rgaverse == 1 & rlseeking == 1 gen greflection = 0 replace greflection = 1 if (rgseeking == 0 & rlaverse == 0) | (rgaverse == 0 & rlseeking == 0) ************************* * Other ************************* * Treatment of male sample ************************** replace mens = 0 if mens==. replace mens_adj = 0 if mens_adj==. replace mens = 0 if female==0 replace mens_adj = 0 if female==0 replace duration = 0 if female==0 *gen std_days = mens_adj*28 *replace std_days = std_days-28 *replace std_days = sqrt(std_days^2) *gen std_days2 = std_days^2 *gen std_days3 = std_days^3 *gen std_days4 = std_days^4 *gen std_days5 = std_days^5 *gen std_days6 = std_days^6 *gen std_days7 = std_days^7 *la var std_days "Adjusted Days of Menstrual Cycle" *gen mens2 = mens^2 *gen mens3 = mens^3 *gen mens4 = mens^4 *gen mens5 = mens^5 *gen mens6 = mens^6 *gen mens7 = mens^7 foreach v in v v2 v3 period age siblings asian otherethn math sciothersci econ socsci humanities pill female { replace `v' = 0 if `v' == . } **************************************************************************************************** **************************************************************************************************** **************************************************************************************************** **************************************************************************************************** **************************************************************************************************** **************************************************************************************************** **************************************************************************************************** * Save data set ********************* save "C:\bcschipper\burkhard\research\hormone\stats\hormone_data_20120906.dta", replace * Log use ********** log using "C:\bcschipper\burkhard\research\hormone\stats\hormrisk_log_20140214.smcl", replace *********************************** * Basic Demographics *********************************** sum female age white asian hispanic black other gpa math sciothersci econ socsci humanities if period == 30 & new == 1 * Hormones * Summary statistics graph twoway histogram progf if period == 30, color(*.5) || kdensity progf if period == 30, by(gender1) graph export "C:\bcschipper\burkhard\research\hormone\stats\histprogf.png", replace graph twoway histogram estf if period == 30, color(*.5) || kdensity estf if period == 30, by(gender1) graph export "C:\bcschipper\burkhard\research\hormone\stats\histestf.png", replace graph twoway histogram testf if period == 30, color(*.5) || kdensity testf if period == 30, by(gender1) graph export "C:\bcschipper\burkhard\research\hormone\stats\histtestf.png", replace graph twoway histogram cortf if period == 30, color(*.5) || kdensity cortf if period == 30, by(gender1) graph export "C:\bcschipper\burkhard\research\hormone\stats\histcortf.png", replace sum testf estf progf cortf pill if period == 30 & female == 1 & new == 1 sum testf estf progf cortf pill if period == 30 & female == 0 & new == 1 sum testf estf progf cortf if period == 30 & female == 1 & rgaverse == 1 & rlaverse == 1 sum testf estf progf cortf if period == 30 & female == 1 & rgaverse == 1 & rlneutral == 1 sum testf estf progf cortf if period == 30 & female == 1 & rgaverse == 1 & rlseeking == 1 sum testf estf progf cortf if period == 30 & female == 1 & rgneutral == 1 & rlaverse == 1 sum testf estf progf cortf if period == 30 & female == 1 & rgneutral == 1 & rlneutral == 1 sum testf estf progf cortf if period == 30 & female == 1 & rgneutral == 1 & rlseeking == 1 sum testf estf progf cortf if period == 30 & female == 1 & rgseeking == 1 & rlaverse == 1 sum testf estf progf cortf if period == 30 & female == 1 & rgseeking == 1 & rlneutral == 1 sum testf estf progf cortf if period == 30 & female == 1 & rgseeking == 1 & rlseeking == 1 sum testf estf progf cortf if period == 30 & female == 0 & rgaverse == 1 & rlaverse == 1 sum testf estf progf cortf if period == 30 & female == 0 & rgaverse == 1 & rlneutral == 1 sum testf estf progf cortf if period == 30 & female == 0 & rgaverse == 1 & rlseeking == 1 sum testf estf progf cortf if period == 30 & female == 0 & rgneutral == 1 & rlaverse == 1 sum testf estf progf cortf if period == 30 & female == 0 & rgneutral == 1 & rlneutral == 1 sum testf estf progf cortf if period == 30 & female == 0 & rgneutral == 1 & rlseeking == 1 sum testf estf progf cortf if period == 30 & female == 0 & rgseeking == 1 & rlaverse == 1 sum testf estf progf cortf if period == 30 & female == 0 & rgseeking == 1 & rlneutral == 1 sum testf estf progf cortf if period == 30 & female == 0 & rgseeking == 1 & rlseeking == 1 sum testf estf progf cortf if period == 30 & female == 1 & rgconsistent == 1 & rlconsistent == 1 sum testf estf progf cortf if period == 30 & female == 1 & rgconsistent == 1 & rlconsistent == 0 sum testf estf progf cortf if period == 30 & female == 1 & rgconsistent == 0 & rlconsistent == 1 sum testf estf progf cortf if period == 30 & female == 1 & rgconsistent == 0 & rlconsistent == 0 sum testf estf progf cortf if period == 30 & female == 0 & rgconsistent == 1 & rlconsistent == 1 sum testf estf progf cortf if period == 30 & female == 0 & rgconsistent == 1 & rlconsistent == 0 sum testf estf progf cortf if period == 30 & female == 0 & rgconsistent == 0 & rlconsistent == 1 sum testf estf progf cortf if period == 30 & female == 0 & rgconsistent == 0 & rlconsistent == 0 * quality check hormone variables reg ntestf age asian otherethn bmi dr preg pill nfolic_fol nperi_fol nlut_fol npremens_fol siblings livalone workout homobi relation luncht foodtime drinktime mouth brush blast smoking vegetarian soybean supplement dst2010 if period == 30 & female == 1, vce(robust) reg ntestf luncht foodtime drinktime mouth brush blast dst2010 if period == 30 & female == 1, vce(robust) reg ntestf age asian otherethn bmi dr preg pill nfolic_fol nperi_fol nlut_fol npremens_fol siblings livalone workout homobi relation luncht foodtime drinktime mouth brush blast smoking vegetarian soybean supplement dst2010 if period == 30 & female == 0, vce(robust) reg ntestf luncht foodtime drinktime mouth brush blast dst2010 if period == 30 & female == 0, vce(robust) reg nestf age asian otherethn bmi dr preg pill nfolic_fol nperi_fol nlut_fol npremens_fol siblings livalone workout homobi relation luncht foodtime drinktime mouth brush blast smoking vegetarian soybean supplement dst2010 if period == 30 & female == 1, vce(robust) reg nestf luncht foodtime drinktime mouth brush blast dst2010 if period == 30 & female == 1, vce(robust) reg nestf age asian otherethn bmi dr preg pill nfolic_fol nperi_fol nlut_fol npremens_fol siblings livalone workout homobi relation luncht foodtime drinktime mouth brush blast smoking vegetarian soybean supplement dst2010 if period == 30 & female == 0, vce(robust) reg nestf luncht foodtime drinktime mouth brush blast dst2010 if period == 30 & female == 0, vce(robust) reg nprogf age asian otherethn bmi dr preg pill nfolic_fol nperi_fol nlut_fol npremens_fol siblings livalone workout homobi relation luncht foodtime drinktime mouth brush blast smoking vegetarian soybean supplement dst2010 if period == 30 & female == 1, vce(robust) reg nprogf luncht foodtime drinktime mouth brush blast dst2010 if period == 30 & female == 1, vce(robust) reg nprogf age asian otherethn bmi dr preg pill nfolic_fol nperi_fol nlut_fol npremens_fol siblings livalone workout homobi relation luncht foodtime drinktime mouth brush blast smoking vegetarian soybean supplement dst2010 if period == 30 & female == 0, vce(robust) reg nprogf luncht foodtime drinktime mouth brush blast dst2010 if period == 30 & female == 0, vce(robust) reg ncortf age asian otherethn bmi dr preg pill nfolic_fol nperi_fol nlut_fol npremens_fol siblings livalone workout homobi relation luncht foodtime drinktime mouth brush blast smoking vegetarian soybean supplement dst2010 if period == 30 & female == 1, vce(robust) reg ncortf luncht foodtime drinktime mouth brush blast dst2010 if period == 30 & female == 1, vce(robust) reg ncortf age asian otherethn bmi dr preg pill nfolic_fol nperi_fol nlut_fol npremens_fol siblings livalone workout homobi relation luncht foodtime drinktime mouth brush blast smoking vegetarian soybean supplement dst2010 if period == 30 & female == 0, vce(robust) reg ncortf luncht foodtime drinktime mouth brush blast dst2010 if period == 30 & female == 0, vce(robust) * need to include session time! scatter testf tests || lfit testf tests scatter estf ests || lfit estf ests scatter progf progs || lfit progf progs scatter cortf corts || lfit cortf corts * correlation between hormones corr testf estf progf cortf if period == 30 eststo qtf: reg ntestf luncht foodtime drinktime mouth brush blast smoking dst2010 if period == 30 & female == 1, vce(robust) eststo qtm: reg ntestf luncht foodtime drinktime mouth brush blast smoking dst2010 if period == 30 & female == 0, vce(robust) eststo qef: reg nestf luncht foodtime drinktime mouth brush blast smoking dst2010 if period == 30 & female == 1, vce(robust) eststo qem: reg nestf luncht foodtime drinktime mouth brush blast smoking dst2010 if period == 30 & female == 0, vce(robust) eststo qpf: reg nprogf luncht foodtime drinktime mouth brush blast smoking dst2010 if period == 30 & female == 1, vce(robust) eststo qpm: reg nprogf luncht foodtime drinktime mouth brush blast smoking dst2010 if period == 30 & female == 0, vce(robust) eststo qcf: reg ncortf luncht foodtime drinktime mouth brush blast smoking dst2010 if period == 30 & female == 1, vce(robust) eststo qcm: reg ncortf luncht foodtime drinktime mouth brush blast smoking dst2010 if period == 30 & female == 0, vce(robust) estout qtf qtm qef qem qpf qpm qcf qcm using "C:\bcschipper\burkhard\research\hormone\quality.tex", replace /// cells(b(star fmt(%9.4f)) se(par fmt(%9.4f))) stats(N r2, fmt(%9.0g %9.4f) labels("\textit{Number of Observations}" "R\textsuperscript{2}") layout("\multicolumn{1}{c}{\quad@}")) label /// starlevels(\text{*} 0.05 \text{**} 0.01) keep(luncht foodtime drinktime mouth brush blast smoking dst2010) /// mlabels((TF) (TM) (EF) (EM) (PF) (PM) (CF) (CM), prefix(\multicolumn{1}{c}{) suffix(})) order(luncht foodtime drinktime mouth brush blast smoking dst2010) /// style(tex) collabels(,none) prehead("\begin{tabular}{l*{@span}{d}}\toprule") /// posthead(\midrule) prefoot(\bottomrule) postfoot("\addlinespace[.1cm]\multicolumn{9}{l}{Robust standard errors in parentheses; Significance levels: * 5\%; ** 1\%.} \\" "\end{tabular}") notype ******************************************* ******************************************* * Analysis for the Original 2012 Version ******************************************* ******************************************* ************************************ * Analysis Risk & Gender ************************************ eststo g0c: reg riskgain female age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & rgconsistent == 1 & girrational == 0, vce(robust) eststo g0: reg riskgain female age asian otherethn gpa math sciothersci econ socsci humanities if period == 30, vce(robust) eststo g1c: reg riskgain female age asian otherethn gpa if period == 30 & rgconsistent == 1 & girrational == 0, vce(robust) eststo l0c: reg riskloss female age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & rlconsistent == 1 & lirrational == 0, vce(robust) eststo l0: reg riskloss female age asian otherethn gpa math sciothersci econ socsci humanities if period == 30, vce(robust) eststo l1c: reg riskloss female age asian otherethn gpa if period == 30 & rlconsistent == 1 & lirrational == 0, vce(robust) estout g0c g0 g1c l0c l0 l1c using "C:\bcschipper\burkhard\research\hormone\gender1.tex", replace /// cells(b(star fmt(%9.4f)) se(par fmt(%9.4f))) stats(N r2, fmt(%9.0g %9.4f) labels("\textit{Number of Observations}" "R\textsuperscript{2}") layout("\multicolumn{1}{c}{\quad@}")) label /// starlevels(\text{*} 0.05 \text{**} 0.01) keep(female age asian otherethn gpa math sciothersci econ socsci humanities) /// mlabels((G0A) (G0) (G1A) (L0A) (L0) (L1A), prefix(\multicolumn{1}{c}{) suffix(})) order(age asian otherethn gpa math sciothersci econ socsci humanities female) /// style(tex) collabels(,none) prehead("\begin{tabular}{l*{@span}{d}}\toprule") /// posthead(\midrule) prefoot(\bottomrule) postfoot("\addlinespace[.1cm]\multicolumn{7}{l}{Robust standard errors in parentheses; Significance levels: * 5\%; ** 1\%.} \\" "\end{tabular}") notype estout g0c g0 g1c l0c l0 l1c using "C:\bcschipper\burkhard\research\hormone\s_gender1.tex", replace /// cells(b(star fmt(%9.4f)) se(par fmt(%9.4f))) stats(N r2, fmt(%9.0g %9.4f) labels("\textit{Number of Observations}" "R\textsuperscript{2}") layout("\multicolumn{1}{c}{\quad@}")) label /// starlevels(\text{*} 0.05 \text{**} 0.01) keep(female) /// mlabels((G0A) (G0) (G1A) (L0A) (L0) (L1A), prefix(\multicolumn{1}{c}{) suffix(})) indicate("Majors of Study = math sciothersci econ socsci humanities", labels(\textmd{Yes} \textmd{No})) order(female) /// style(tex) collabels(,none) prehead("\begin{tabular}{l*{@span}{d}}\toprule") /// posthead(\midrule) prefoot(\bottomrule) postfoot("\addlinespace[.1cm]\multicolumn{7}{l}{Robust standard errors in parentheses; Significance levels: * 5\%; ** 1\%.} \\" /// "\multicolumn{7}{l}{We suppress from the report the coefficients for Age, Asian, Other, and GPA.} \\""\end{tabular}") notype * Robustness checks ************************** * gender1 standard errors clustered on session reg riskgain female age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & rgconsistent == 1 & girrational == 0, robust cluster(session) reg riskgain female age asian otherethn gpa math sciothersci econ socsci humanities if period == 30, robust cluster(session) reg riskloss female age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & rlconsistent == 1 & lirrational == 0, robust cluster(session) reg riskloss female age asian otherethn gpa math sciothersci econ socsci humanities if period == 30, robust cluster(session) * gender1 with session dummies reg riskgain sessiondum* female age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & rgconsistent == 1 & girrational == 0, vce(robust) reg riskgain sessiondum* female age asian otherethn gpa math sciothersci econ socsci humanities if period == 30, vce(robust) reg riskloss sessiondum* female age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & rlconsistent == 1 & lirrational == 0, vce(robust) reg riskloss sessiondum* female age asian otherethn gpa math sciothersci econ socsci humanities if period == 30, vce(robust) * fixed effects model xtset, clear xtset datum xtreg riskgain female age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & rgconsistent == 1 & girrational == 0, fe xtreg riskgain female age asian otherethn gpa math sciothersci econ socsci humanities if period == 30, fe xtreg riskloss female age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & rlconsistent == 1 & lirrational == 0, fe xtreg riskloss female age asian otherethn gpa math sciothersci econ socsci humanities if period == 30, fe * minimum reg riskgain female if period == 30 & rgconsistent == 1 & girrational == 0 reg riskloss female if period == 30 & rlconsistent == 1 & lirrational == 0 * gender1 ordered logit ologit riskgain female age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & rgconsistent == 1 & girrational == 0, vce(robust) ologit riskgain female age asian otherethn gpa math sciothersci econ socsci humanities if period == 30, vce(robust) ologit riskloss female age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & rlconsistent == 1 & lirrational == 0, vce(robust) ologit riskloss female age asian otherethn gpa math sciothersci econ socsci humanities if period == 30, vce(robust) * gender1 ordered probit oprobit riskgain female age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & rgconsistent == 1 & girrational == 0, vce(robust) oprobit riskgain female age asian otherethn gpa math sciothersci econ socsci humanities if period == 30, vce(robust) oprobit riskloss female age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & rlconsistent == 1 & lirrational == 0, vce(robust) oprobit riskloss female age asian otherethn gpa math sciothersci econ socsci humanities if period == 30, vce(robust) * gender1 with further controls (bmi siblings) reg riskgain female age siblings bmi asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & rgconsistent == 1 & girrational == 0, vce(robust) reg riskgain female age siblings bmi asian otherethn gpa math sciothersci econ socsci humanities if period == 30, vce(robust) reg riskloss female age siblings bmi asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & rlconsistent == 1 & lirrational == 0, vce(robust) reg riskloss female age siblings bmi asian otherethn gpa math sciothersci econ socsci humanities if period == 30, vce(robust) * does not update automatically in tex file eststo r0: logit reflection female age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & consistent == 1 & irrational == 0, vce(robust) eststo r1: logit reflection female age asian otherethn gpa if period == 30 & consistent == 1 & irrational == 0, vce(robust) eststo c0: logit consistent female age asian otherethn gpa sciothersci econ socsci humanities if period == 30 & new == 1, vce(robust) estout r0 r1 c0 using "C:\bcschipper\burkhard\research\hormone\gender2.tex", replace /// cells(b(star fmt(%9.4f)) se(par fmt(%9.4f))) stats(N, fmt(%9.0g %9.4f) labels("\textit{Number of Observations}") layout("\multicolumn{1}{c}{\quad@}")) label /// starlevels(\text{*} 0.05 \text{**} 0.01) keep(female age asian otherethn gpa math sciothersci econ socsci humanities) /// mlabels((R0) (R1) (A0), prefix(\multicolumn{1}{c}{) suffix(})) order(age asian otherethn gpa math sciothersci econ socsci humanities female) /// style(tex) collabels(,none) prehead("\begin{tabular}{l*{@span}{d}}\toprule") /// posthead(\midrule) prefoot(\bottomrule) postfoot("\addlinespace[.1cm]\multicolumn{4}{l}{Robust standard errors in parentheses; Significance levels: * 5\%; ** 1\%.} \\" "\end{tabular}") notype estout r0 r1 c0 using "C:\bcschipper\burkhard\research\hormone\s_gender2.tex", replace /// cells(b(star fmt(%9.4f)) se(par fmt(%9.4f))) stats(N, fmt(%9.0g %9.4f) labels("\textit{Number of Observations}") layout("\multicolumn{1}{c}{\quad@}")) label /// starlevels(\text{*} 0.05 \text{**} 0.01) keep(female) /// mlabels((R0) (R1) (A0), prefix(\multicolumn{1}{c}{) suffix(})) indicate("Majors of Study = math sciothersci econ socsci humanities", labels(\textmd{Yes} \textmd{No})) order(female) /// style(tex) collabels(,none) prehead("\begin{tabular}{l*{@span}{d}}\toprule") /// posthead(\midrule) prefoot(\bottomrule) postfoot("\addlinespace[.1cm]\multicolumn{4}{l}{Robust standard errors in parentheses; Significance levels: * 5\%; ** 1\%.} \\" /// "\multicolumn{4}{l}{We suppress from the report the coefficients for Age, Asian, Other, and GPA.} \\""\end{tabular}") notype * Robustness checks ********************* * gender2 standard errors clustered on session logit reflection female age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & consistent == 1 & irrational == 0, robust cluster(session) logit reflection female age asian otherethn gpa if period == 30 & consistent == 1 & irrational == 0, robust cluster(session) logit consistent female age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & new == 1, robust cluster(session) * gender2 sessiondummies logit reflection sessiondum* female age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & consistent == 1 & irrational == 0, vce(robust) logit reflection sessiondum* female age asian otherethn gpa if period == 30 & consistent == 1 & irrational == 0, vce(robust) logit consistent sessiondum* female age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & new == 1, vce(robust) * session fixed effects xtset, clear xtset datum xtlogit reflection female age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & consistent == 1 & irrational == 0, fe xtlogit reflection female age asian otherethn gpa if period == 30 & consistent == 1 & irrational == 0, fe xtlogit consistent female age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & new == 1, fe * gender2 probit probit reflection female age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & consistent == 1 & irrational == 0, vce(robust) probit reflection female age asian otherethn gpa if period == 30 & consistent == 1 & irrational == 0, vce(robust) probit consistent female age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & new == 1, vce(robust) * gender2 with further controls (bmi siblings) logit reflection female age siblings bmi asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & consistent == 1 & irrational == 0, vce(robust) logit reflection female age siblings bmi asian otherethn gpa if period == 30 & consistent == 1 & irrational == 0, vce(robust) logit consistent female age siblings bmi asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & new == 1, vce(robust) * gender 2 minimum logit reflection female if period == 30 & consistent == 1 & irrational == 0, vce(robust) logit consistent female if period == 30 & new == 1, vce(robust) logit consistent female age asian otherethn gpa if period == 30 & new == 1, vce(robust) **************************************** * Analysis Risk & Sexual Behavior **************************************** reg riskgain female homobi dating livalone age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & rgconsistent == 1 & girrational == 0, vce(robust) reg riskloss female homobi dating livalone age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & rlconsistent == 1 & lirrational == 0, vce(robust) reg riskgain female homobi dating livalone if period == 30 & rgconsistent == 1 & girrational == 0, vce(robust) reg riskloss female homobi dating livalone if period == 30 & rlconsistent == 1 & lirrational == 0, vce(robust) reg riskgain homobi dating livalone age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & rgconsistent == 1 & girrational == 0 & female == 1, vce(robust) reg riskloss homobi dating livalone age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & rlconsistent == 1 & lirrational == 0 & female == 1, vce(robust) reg riskgain homobi dating livalone age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & rgconsistent == 1 & girrational == 0 & female == 0, vce(robust) reg riskloss homobi dating livalone age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & rlconsistent == 1 & lirrational == 0 & female == 0, vce(robust) reg riskgain homo dating livalone age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & rgconsistent == 1 & girrational == 0 & female == 0, vce(robust) reg riskloss homo dating livalone age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & rlconsistent == 1 & lirrational == 0 & female == 0, vce(robust) **************************************** * Analysis Risk & Contraceptives **************************************** * risk aversion and contraceptives eststo g2: reg riskgain female pill age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & rgconsistent == 1 & girrational == 0, vce(robust) eststo l2: reg riskloss female pill age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & rlconsistent == 1 & lirrational == 0, vce(robust) eststo g2f: reg riskgain pill age asian otherethn gpa sciothersci econ socsci humanities if period == 30 & rgconsistent == 1 & girrational == 0 & female == 1, vce(robust) eststo l2f: reg riskloss pill age asian otherethn gpa sciothersci econ socsci humanities if period == 30 & rlconsistent == 1 & lirrational == 0 & female == 1, vce(robust) estout g2 l2 g2f l2f using "C:\bcschipper\burkhard\research\hormone\contraceptives.tex", replace /// cells(b(star fmt(%9.4f)) se(par fmt(%9.4f))) stats(N r2, fmt(%9.0g %9.4f) labels("\textit{Number of Observations}" "R\textsuperscript{2}") layout("\multicolumn{1}{c}{\quad@}")) label /// starlevels(\text{*} 0.05 \text{**} 0.01) keep(female pill age asian otherethn gpa math sciothersci econ socsci humanities) /// mlabels((G2) (L2) (G2F) (L2F), prefix(\multicolumn{1}{c}{) suffix(})) order(age asian otherethn gpa math sciothersci econ socsci humanities female pill) /// style(tex) collabels(,none) prehead("\begin{tabular}{l*{@span}{d}}\toprule") /// posthead(\midrule) prefoot(\bottomrule) postfoot("\addlinespace[.1cm]\multicolumn{5}{l}{Robust standard errors in parentheses; Significance levels: * 5\%; ** 1\%.} \\" "\end{tabular}") notype estout g2 l2 g2f l2f using "C:\bcschipper\burkhard\research\hormone\s_contraceptives.tex", replace /// cells(b(star fmt(%9.4f)) se(par fmt(%9.4f))) stats(N r2, fmt(%9.0g %9.4f) labels("\textit{Number of Observations}" "R\textsuperscript{2}") layout("\multicolumn{1}{c}{\quad@}")) label /// starlevels(\text{*} 0.05 \text{**} 0.01) keep(female pill) /// mlabels((G2) (L2) (G2F) (L2F), prefix(\multicolumn{1}{c}{) suffix(})) order(female pill) /// style(tex) collabels(,none) prehead("\begin{tabular}{l*{@span}{d}}\toprule") /// posthead(\midrule) prefoot(\bottomrule) postfoot("\addlinespace[.1cm]\multicolumn{5}{l}{Robust standard errors in parentheses; Significance levels: * 5\%; ** 1\%.} \\" /// "\multicolumn{5}{l}{We suppress from the report the coefficients for Age, Asian, Other, GPA, and Majors of Study.} \\""\end{tabular}") notype * robustness checks * ologit ologit riskgain female pill age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & rgconsistent == 1 & girrational == 0, vce(robust) ologit riskloss female pill age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & rlconsistent == 1 & lirrational == 0, vce(robust) ologit riskgain pill age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & rgconsistent == 1 & girrational == 0 & female == 1, vce(robust) ologit riskloss pill age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & rlconsistent == 1 & lirrational == 0 & female == 1, vce(robust) * oprobit oprobit riskgain female pill age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & rgconsistent == 1 & girrational == 0, vce(robust) oprobit riskloss female pill age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & rlconsistent == 1 & lirrational == 0, vce(robust) oprobit riskgain pill age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & rgconsistent == 1 & girrational == 0 & female == 1, vce(robust) oprobit riskloss pill age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & rlconsistent == 1 & lirrational == 0 & female == 1, vce(robust) * cluster standard errors reg riskgain female pill age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & rgconsistent == 1 & girrational == 0, robust cluster(session) reg riskloss female pill age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & rlconsistent == 1 & lirrational == 0, robust cluster(session) reg riskgain pill age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & rgconsistent == 1 & girrational == 0 & female == 1, robust cluster(session) reg riskloss pill age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & rlconsistent == 1 & lirrational == 0 & female == 1, robust cluster(session) * add session dummies reg riskgain sessiondum* female pill age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & rgconsistent == 1 & girrational == 0, vce(robust) reg riskloss sessiondum* female pill age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & rlconsistent == 1 & lirrational == 0, vce(robust) reg riskgain sessiondum* pill age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & rgconsistent == 1 & girrational == 0 & female == 1, vce(robust) reg riskloss sessiondum* pill age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & rlconsistent == 1 & lirrational == 0 & female == 1, vce(robust) * fixed effects model xtset, clear xtset datum xtreg riskgain female pill age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & rgconsistent == 1 & girrational == 0, fe xtreg riskloss female pill age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & rlconsistent == 1 & lirrational == 0, fe xtreg riskgain pill age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & rgconsistent == 1 & girrational == 0 & female == 1, fe xtreg riskloss pill age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & rlconsistent == 1 & lirrational == 0 & female == 1, fe * add further controls bmi siblings workout livalone dating homobi reg riskgain female pill bmi siblings workout livalone dating homobi age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & rgconsistent == 1 & girrational == 0, vce(robust) reg riskloss female pill bmi siblings workout livalone dating homobi age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & rlconsistent == 1 & lirrational == 0, vce(robust) reg riskgain pill bmi siblings workout livalone dating homobi age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & rgconsistent == 1 & girrational == 0 & female == 1, vce(robust) reg riskloss pill bmi siblings workout livalone dating homobi age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & rlconsistent == 1 & lirrational == 0 & female == 1, vce(robust) * drop demographics controls reg riskgain female pill age asian otherethn gpa if period == 30 & rgconsistent == 1 & girrational == 0, vce(robust) reg riskloss female pill age asian otherethn gpa if period == 30 & rlconsistent == 1 & lirrational == 0, vce(robust) reg riskgain pill age asian otherethn gpa if period == 30 & rgconsistent == 1 & girrational == 0 & female == 1, vce(robust) reg riskloss pill age asian otherethn gpa if period == 30 & rlconsistent == 1 & lirrational == 0 & female == 1, vce(robust) * does not update automatically in tex file eststo r2: logit reflection female pill age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & consistent == 1 & irrational == 0, vce(robust) eststo c2: logit consistent female pill age asian otherethn gpa sciothersci econ socsci humanities if period == 30 & new == 1, vce(robust) eststo r2f: logit reflection pill age asian otherethn gpa sciothersci econ socsci if period == 30 & consistent == 1 & irrational == 0 & female == 1, vce(robust) eststo c2f: logit consistent pill age asian otherethn gpa sciothersci econ socsci humanities if period == 30 & new == 1 & female == 1, vce(robust) estout r2 r2f c2 c2f using "C:\bcschipper\burkhard\research\hormone\contraceptives2.tex", replace /// cells(b(star fmt(%9.4f)) se(par fmt(%9.4f))) stats(N, fmt(%9.0g %9.4f) labels("\textit{Number of Observations}") layout("\multicolumn{1}{c}{\quad@}")) label /// starlevels(\text{*} 0.05 \text{**} 0.01) keep(female pill age asian otherethn gpa math sciothersci econ socsci humanities) /// mlabels((R2) (R2F) (A2) (A2F), prefix(\multicolumn{1}{c}{) suffix(})) order(age asian otherethn gpa math sciothersci econ socsci humanities female pill) /// style(tex) collabels(,none) prehead("\begin{tabular}{l*{@span}{d}}\toprule") /// posthead(\midrule) prefoot(\bottomrule) postfoot("\addlinespace[.1cm]\multicolumn{5}{l}{Robust standard errors in parentheses; Significance levels: * 5\%; ** 1\%.} \\" "\end{tabular}") notype estout r2 r2f c2 c2f using "C:\bcschipper\burkhard\research\hormone\s_contraceptives2.tex", replace /// cells(b(star fmt(%9.4f)) se(par fmt(%9.4f))) stats(N, fmt(%9.0g %9.4f) labels("\textit{Number of Observations}") layout("\multicolumn{1}{c}{\quad@}")) label /// starlevels(\text{*} 0.05 \text{**} 0.01) keep(female pill) /// mlabels((R2) (R2F) (A2) (A2F), prefix(\multicolumn{1}{c}{) suffix(})) order(female pill) /// style(tex) collabels(,none) prehead("\begin{tabular}{l*{@span}{d}}\toprule") /// posthead(\midrule) prefoot(\bottomrule) postfoot("\addlinespace[.1cm]\multicolumn{5}{l}{Robust standard errors in parentheses; Significance levels: * 5\%; ** 1\%.} \\" /// "\multicolumn{5}{l}{We suppress from the report the coefficients for Age, Asian, Other, GPA, and Majors of Study.} \\""\end{tabular}") notype * Robustness checks ************************** * probit probit reflection female pill age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & consistent == 1 & irrational == 0, vce(robust) probit consistent female pill age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & new == 1, vce(robust) probit reflection pill age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & consistent == 1 & irrational == 0 & female == 1, vce(robust) probit consistent pill age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & new == 1 & female == 1, vce(robust) * cluster standard errors logit reflection female pill age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & consistent == 1 & irrational == 0, robust cluster(session) logit consistent female pill age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & new == 1, robust cluster(session) logit reflection pill age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & consistent == 1 & irrational == 0 & female == 1, robust cluster(session) logit consistent pill age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & new == 1 & female == 1, robust cluster(session) * add session dummies logit reflection sessiondum* female pill age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & consistent == 1 & irrational == 0, vce(robust) logit consistent sessiondum* female pill age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & new == 1, vce(robust) *logit reflection sessiondum* pill age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & consistent == 1 & irrational == 0 & female == 1, vce(robust) logit consistent sessiondum* pill age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & new == 1 & female == 1, vce(robust) * session fixed effects xtset, clear xtset datum xtlogit reflection female pill age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & consistent == 1 & irrational == 0, fe xtlogit consistent female pill age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & new == 1, fe *xtlogit reflection pill age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & consistent == 1 & irrational == 0 & female == 1, fe xtlogit consistent pill age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & new == 1 & female == 1, fe * bmi siblings workout livalone dating homobi reg riskgain female pill bmi siblings workout livalone dating homobi age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & rgconsistent == 1 & girrational == 0, vce(robust) reg riskloss female pill bmi siblings workout livalone dating homobi age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & rlconsistent == 1 & lirrational == 0, vce(robust) reg riskgain pill bmi siblings workout livalone dating homobi age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & rgconsistent == 1 & girrational == 0 & female == 1, vce(robust) reg riskloss pill bmi siblings workout livalone dating homobi age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & rlconsistent == 1 & lirrational == 0 & female == 1, vce(robust) logit reflection female pill bmi siblings workout livalone dating homobi age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & consistent == 1 & irrational == 0, vce(robust) logit consistent female pill bmi siblings workout livalone dating homobi age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & new == 1, vce(robust) logit reflection pill bmi siblings workout livalone dating homobi age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & consistent == 1 & irrational == 0 & female == 1, vce(robust) logit consistent pill bmi siblings workout livalone dating homobi age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & new == 1 & female == 1, vce(robust) * dropping demographics logit reflection female pill age asian otherethn gpa if period == 30 & consistent == 1 & irrational == 0, vce(robust) logit consistent female pill age asian otherethn gpa if period == 30 & new == 1, vce(robust) logit reflection pill age asian otherethn gpa if period == 30 & consistent == 1 & irrational == 0 & female == 1, vce(robust) logit consistent pill age asian otherethn gpa if period == 30 & new == 1 & female == 1, vce(robust) ******************************************* * Analysis Risk & Sex Hormones ******************************************* * Graphs ************* scatter testf riskloss || lfit testf riskloss, by(gender1) scatter estf riskloss || lfit estf riskloss, by(gender1) scatter progf riskloss || lfit progf riskloss, by(gender1) scatter cortf riskloss || lfit cortf riskloss, by(gender1) scatter testf riskgain || lfit testf riskgain, by(gender1) scatter estf riskgain || lfit estf riskgain, by(gender1) scatter progf riskgain || lfit progf riskgain, by(gender1) scatter cortf riskgain || lfit cortf riskgain, by(gender1) * we exclude now the pregnant lady (preg == 0) * (always restricted to consistent & rational) * Risk aversion and hormones: full ************************************ eststo g3: reg riskgain ntestf nestf nprogf ncortf pill female if period == 30 & preg == 0 & rgconsistent == 1 & girrational == 0, vce(robust) eststo g4: reg riskgain ntestf nestf nprogf ncortf pill female age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & preg == 0 & rgconsistent == 1 & girrational == 0, vce(robust) eststo l3: reg riskloss ntestf nestf nprogf ncortf pill female if period == 30 & preg == 0 & rlconsistent == 1 & lirrational == 0, vce(robust) eststo l4: reg riskloss ntestf nestf nprogf ncortf pill female age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & preg == 0 & rlconsistent == 1 & lirrational == 0, vce(robust) estout g3 g4 l3 l4 using "C:\bcschipper\burkhard\research\hormone\full.tex", replace /// cells(b(star fmt(%9.4f)) se(par fmt(%9.4f))) stats(N r2, fmt(%9.0g %9.4f) labels("\textit{Number of Observations}" "R\textsuperscript{2}") layout("\multicolumn{1}{c}{\quad@}")) label /// starlevels(\text{*} 0.05 \text{**} 0.01) keep(ntestf nestf nprogf ncortf pill female age asian otherethn gpa math sciothersci econ socsci humanities) /// mlabels((G3) (G4) (L3) (L4), prefix(\multicolumn{1}{c}{) suffix(})) order(age asian otherethn gpa math sciothersci econ socsci humanities female ntestf nestf nprogf ncortf pill) /// style(tex) collabels(,none) prehead("\begin{tabular}{l*{@span}{d}}\toprule") /// posthead(\midrule) prefoot(\bottomrule) postfoot("\addlinespace[.1cm]\multicolumn{5}{l}{Robust standard errors in parentheses; Significance levels: * 5\%; ** 1\%.} \\" "\end{tabular}") notype estout g3 g4 l3 l4 using "C:\bcschipper\burkhard\research\hormone\s_full.tex", replace /// cells(b(star fmt(%9.4f)) se(par fmt(%9.4f))) stats(N r2, fmt(%9.0g %9.4f) labels("\textit{Number of Observations}" "R\textsuperscript{2}") layout("\multicolumn{1}{c}{\quad@}")) label /// starlevels(\text{*} 0.05 \text{**} 0.01) keep(ntestf nestf nprogf ncortf pill female) /// mlabels((G3) (G4) (L3) (L4), prefix(\multicolumn{1}{c}{) suffix(})) indicate("Demographics = age asian otherethn" "Majors of Study \& GPA = math sciothersci econ socsci humanities", labels(\textmd{Yes} \textmd{No})) order(female ntestf nestf nprogf ncortf pill) /// style(tex) collabels(,none) prehead("\begin{tabular}{l*{@span}{d}}\toprule") /// posthead(\midrule) prefoot(\bottomrule) postfoot("\addlinespace[.1cm]\multicolumn{5}{l}{Robust standard errors in parentheses; Significance levels: * 5\%; ** 1\%.} \\" "\end{tabular}") notype * Robustness checks * full standard errors clustered on sessions reg riskgain ntestf nestf nprogf ncortf pill female if period == 30 & preg == 0 & rgconsistent == 1 & girrational == 0, robust cluster(session) reg riskgain ntestf nestf nprogf ncortf pill female age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & preg == 0 & rgconsistent == 1 & girrational == 0, robust cluster(session) reg riskloss ntestf nestf nprogf ncortf pill female if period == 30 & preg == 0 & rlconsistent == 1 & lirrational == 0, robust cluster(session) reg riskloss ntestf nestf nprogf ncortf pill female age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & preg == 0 & rlconsistent == 1 & lirrational == 0, robust cluster(session) * full session dummies reg riskgain sessiondum* ntestf nestf nprogf ncortf pill female if period == 30 & preg == 0 & rgconsistent == 1 & girrational == 0, vce(robust) reg riskgain sessiondum* ntestf nestf nprogf ncortf pill female age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & preg == 0 & rgconsistent == 1 & girrational == 0, vce(robust) reg riskloss sessiondum* ntestf nestf nprogf ncortf pill female if period == 30 & preg == 0 & rlconsistent == 1 & lirrational == 0, vce(robust) reg riskloss sessiondum* ntestf nestf nprogf ncortf pill female age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & preg == 0 & rlconsistent == 1 & lirrational == 0, vce(robust) * full session fixed effects xtset, clear xtset datum xtreg riskgain ntestf nestf nprogf ncortf pill female if period == 30 & preg == 0 & rgconsistent == 1 & girrational == 0, fe xtreg riskgain ntestf nestf nprogf ncortf pill female age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & preg == 0 & rgconsistent == 1 & girrational == 0, fe xtreg riskloss ntestf nestf nprogf ncortf pill female if period == 30 & preg == 0 & rlconsistent == 1 & lirrational == 0, fe xtreg riskloss ntestf nestf nprogf ncortf pill female age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & preg == 0 & rlconsistent == 1 & lirrational == 0, fe * without majors of study reg riskgain ntestf nestf nprogf ncortf pill female age asian otherethn gpa if period == 30 & preg == 0 & rgconsistent == 1 & girrational == 0, vce(robust) reg riskloss ntestf nestf nprogf ncortf pill female age asian otherethn gpa if period == 30 & preg == 0 & rlconsistent == 1 & lirrational == 0, vce(robust) * full ologit ologit riskgain ntestf nestf nprogf ncortf pill female if period == 30 & preg == 0 & rgconsistent == 1 & girrational == 0, vce(robust) ologit riskgain ntestf nestf nprogf ncortf pill female age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & preg == 0 & rgconsistent == 1 & girrational == 0, vce(robust) ologit riskloss ntestf nestf nprogf ncortf pill female if period == 30 & preg == 0 & rlconsistent == 1 & lirrational == 0, vce(robust) ologit riskloss ntestf nestf nprogf ncortf pill female age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & preg == 0 & rlconsistent == 1 & lirrational == 0, vce(robust) * full probit oprobit riskgain ntestf nestf nprogf ncortf pill female if period == 30 & preg == 0 & rgconsistent == 1 & girrational == 0, vce(robust) oprobit riskgain ntestf nestf nprogf ncortf pill female age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & preg == 0 & rgconsistent == 1 & girrational == 0, vce(robust) oprobit riskloss ntestf nestf nprogf ncortf pill female if period == 30 & preg == 0 & rlconsistent == 1 & lirrational == 0, vce(robust) oprobit riskloss ntestf nestf nprogf ncortf pill female age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & preg == 0 & rlconsistent == 1 & lirrational == 0, vce(robust) * full with further demographic controls (bmi siblings workout livalone dating homobi) reg riskgain ntestf nestf nprogf ncortf pill female bmi siblings workout livalone dating homobi age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & preg == 0 & rgconsistent == 1 & girrational == 0, vce(robust) reg riskloss ntestf nestf nprogf ncortf pill female bmi siblings workout livalone dating homobi age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & preg == 0 & rlconsistent == 1 & lirrational == 0, vce(robust) * full with further controls on the hormone quality (sessiondum* luncht foodtime drinktime mouth brush blast dst2010) reg riskgain ntestf nestf nprogf ncortf sessiondum* luncht foodtime drinktime mouth brush blast dst2010 pill female age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & preg == 0 & rgconsistent == 1 & girrational == 0, vce(robust) reg riskloss ntestf nestf nprogf ncortf sessiondum* luncht foodtime drinktime mouth brush blast dst2010 pill female age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & preg == 0 & rlconsistent == 1 & lirrational == 0, vce(robust) * full with further controls for dietary preferences (smoking vegetarian soybean) reg riskgain ntestf nestf nprogf ncortf smoking vegetarian soybean pill female age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & preg == 0 & rgconsistent == 1 & girrational == 0, vce(robust) reg riskloss ntestf nestf nprogf ncortf smoking vegetarian soybean pill female age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & preg == 0 & rlconsistent == 1 & lirrational == 0, vce(robust) * full with average hormone levels reg riskgain ntestav nestav nprogav ncortav pill female if period == 30 & preg == 0 & rgconsistent == 1 & girrational == 0, vce(robust) reg riskgain ntestav nestav nprogav ncortav pill female age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & preg == 0 & rgconsistent == 1 & girrational == 0, vce(robust) reg riskloss ntestav nestav nprogav ncortav pill female if period == 30 & preg == 0 & rlconsistent == 1 & lirrational == 0, vce(robust) reg riskloss ntestav nestav nprogav ncortav pill female age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & preg == 0 & rlconsistent == 1 & lirrational == 0, vce(robust) * full stepwise forward stepwise reg riskgain ntestf nestf nprogf ncortf luncht foodtime drinktime mouth brush blast pill female bmi siblings workout livalone dating homobi nconception smoking vegetarian soybean age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & preg == 0 & rgconsistent == 1 & girrational == 0, pe(0.05) stepwise reg riskloss ntestf nestf nprogf ncortf luncht foodtime drinktime mouth brush blast pill female bmi siblings workout livalone dating homobi nconception smoking vegetarian soybean age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & preg == 0 & rlconsistent == 1 & lirrational == 0, pe(0.05) * full stepwise backward stepwise reg riskgain ntestf nestf nprogf ncortf luncht foodtime drinktime mouth brush blast pill female bmi siblings workout livalone dating homobi nconception smoking vegetarian soybean age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & preg == 0 & rgconsistent == 1 & girrational == 0, pr(0.05) stepwise reg riskloss ntestf nestf nprogf ncortf luncht foodtime drinktime mouth brush blast pill female bmi siblings workout livalone dating homobi nconception smoking vegetarian soybean age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & preg == 0 & rlconsistent == 1 & lirrational == 0, pr(0.05) * reflection & hormones **************************** eststo r3: logit reflection ntestf nestf nprogf ncortf pill female age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & preg == 0 & consistent == 1 & irrational == 0, vce(robust) eststo fr3: logit reflection ntestf nestf nprogf ncortf pill age asian otherethn gpa sciothersci econ socsci humanities if period == 30 & preg == 0 & consistent == 1 & irrational == 0 & female == 1, vce(robust) eststo fr4: logit reflection ntestf nestf nprogf ncortf pill pilltestf pillestf pillprogf age asian otherethn gpa sciothersci econ socsci humanities if period == 30 & preg == 0 & consistent == 1 & irrational == 0 & female == 1, vce(robust) halfstep eststo mr3: logit reflection ntestf nestf nprogf ncortf age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & preg == 0 & consistent == 1 & irrational == 0 & female == 0, vce(robust) estout r3 fr3 fr4 mr3 using "C:\bcschipper\burkhard\research\hormone\reflection.tex", replace /// cells(b(star fmt(%9.4f)) se(par fmt(%9.4f))) stats(N r2, fmt(%9.0g %9.4f) labels("\textit{Number of Observations}" "R\textsuperscript{2}") layout("\multicolumn{1}{c}{\quad@}")) label /// starlevels(\text{*} 0.05 \text{**} 0.01) keep(ntestf nestf nprogf ncortf pill pilltestf pillestf pillprogf female age asian otherethn gpa math sciothersci econ socsci humanities) /// mlabels((R3) (R3F) (R4F) (R3M), prefix(\multicolumn{1}{c}{) suffix(})) order(age asian otherethn gpa math sciothersci econ socsci humanities female ntestf nestf nprogf ncortf pill pilltestf pillestf pillprogf) /// style(tex) collabels(,none) prehead("\begin{tabular}{l*{@span}{d}}\toprule") /// posthead(\midrule) prefoot(\bottomrule) postfoot("\addlinespace[.1cm]\multicolumn{5}{l}{Robust standard errors in parentheses; Significance levels: * 5\%; ** 1\%.} \\" "\end{tabular}") notype estout r3 fr3 fr4 mr3 using "C:\bcschipper\burkhard\research\hormone\s_reflection.tex", replace /// cells(b(star fmt(%9.4f)) se(par fmt(%9.4f))) stats(N r2, fmt(%9.0g %9.4f) labels("\textit{Number of Observations}" "R\textsuperscript{2}") layout("\multicolumn{1}{c}{\quad@}")) label /// starlevels(\text{*} 0.05 \text{**} 0.01) keep(ntestf nestf nprogf ncortf pill pilltestf pillestf pillprogf female) /// mlabels((R3) (R3F) (R4F) (R3M), prefix(\multicolumn{1}{c}{) suffix(})) order(female ntestf nestf nprogf ncortf pill pilltestf pillestf pillprogf) /// style(tex) collabels(,none) prehead("\begin{tabular}{l*{@span}{d}}\toprule") /// posthead(\midrule) prefoot(\bottomrule) postfoot("\addlinespace[.1cm]\multicolumn{5}{l}{Robust standard errors in parentheses; Significance levels: * 5\%; ** 1\%.} \\" /// "\multicolumn{5}{l}{We suppress from the report the coefficients for Age, Asian, Other, GPA, and Majors of Study.} \\""\end{tabular}") notype * Robustness checks * reflection with logit logit reflection ntestf nestf nprogf ncortf pill female age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & preg == 0 & consistent == 1 & irrational == 0, vce(robust) logit reflection ntestf nestf nprogf ncortf pill age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & preg == 0 & consistent == 1 & irrational == 0 & female == 1, vce(robust) * logit reflection ntestf nestf nprogf ncortf pill pilltestf pillestf pillprogf age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & preg == 0 & consistent == 1 & irrational == 0 & female == 1, vce(robust) logit reflection ntestf nestf nprogf ncortf age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & preg == 0 & consistent == 1 & irrational == 0 & female == 0, vce(robust) * cluster standard errors on session level logit reflection ntestf nestf nprogf ncortf pill female age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & preg == 0 & consistent == 1 & irrational == 0, robust cluster(session) logit reflection ntestf nestf nprogf ncortf pill age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & preg == 0 & consistent == 1 & irrational == 0 & female == 1, robust cluster(session) logit reflection ntestf nestf nprogf ncortf pill pilltestf pillestf pillprogf age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & preg == 0 & consistent == 1 & irrational == 0 & female == 1, robust cluster(session) halfstep logit reflection ntestf nestf nprogf ncortf age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & preg == 0 & consistent == 1 & irrational == 0 & female == 0, robust cluster(session) * add session dummies logit reflection sessiondum* ntestf nestf nprogf ncortf pill female age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & preg == 0 & consistent == 1 & irrational == 0, vce(robust) logit reflection sessiondum* ntestf nestf nprogf ncortf pill age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & preg == 0 & consistent == 1 & irrational == 0 & female == 1, vce(robust) logit reflection sessiondum* ntestf nestf nprogf ncortf pill pilltestf pillestf pillprogf age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & preg == 0 & consistent == 1 & irrational == 0 & female == 1, vce(robust) halfstep logit reflection sessiondum* ntestf nestf nprogf ncortf age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & preg == 0 & consistent == 1 & irrational == 0 & female == 0, vce(robust) * fixed effects model xtset, clear xtset datum xtlogit reflection ntestf nestf nprogf ncortf pill female age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & preg == 0 & consistent == 1 & irrational == 0, fe *xtlogit reflection ntestf nestf nprogf ncortf pill age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & preg == 0 & consistent == 1 & irrational == 0 & female == 1, fe *xtlogit reflection ntestf nestf nprogf ncortf pill pilltestf pillestf pillprogf age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & preg == 0 & consistent == 1 & irrational == 0 & female == 1, fe halfstep xtlogit reflection ntestf nestf nprogf ncortf age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & preg == 0 & consistent == 1 & irrational == 0 & female == 0, fe * drop controls logit reflection ntestf nestf nprogf ncortf pill female gpa if period == 30 & preg == 0 & consistent == 1 & irrational == 0, vce(robust) logit reflection ntestf nestf nprogf ncortf pill age gpa if period == 30 & preg == 0 & consistent == 1 & irrational == 0 & female == 1, vce(robust) logit reflection ntestf nestf nprogf ncortf pill pilltestf pillestf pillprogf gpa if period == 30 & preg == 0 & consistent == 1 & irrational == 0 & female == 1, vce(robust) halfstep logit reflection ntestf nestf nprogf ncortf gpa if period == 30 & preg == 0 & consistent == 1 & irrational == 0 & female == 0, vce(robust) * reflection with further demographic controls (dr bmi siblings workout livalone dating homobi nconception) probit reflection ntestf nestf nprogf ncortf pill female dr bmi siblings workout livalone dating homobi nconception age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & preg == 0 & consistent == 1 & irrational == 0, vce(robust) * probit reflection ntestf nestf nprogf ncortf pill dr bmi siblings workout livalone dating homobi nconception age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & preg == 0 & consistent == 1 & irrational == 0 & female == 1, vce(robust) halfstep probit reflection ntestf nestf nprogf ncortf pill pilltestf pillestf pillprogf dr bmi siblings workout livalone dating homobi nconception age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & preg == 0 & consistent == 1 & irrational == 0 & female == 1, vce(robust) probit reflection ntestf nestf nprogf ncortf dr bmi siblings workout livalone dating homobi age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & preg == 0 & consistent == 1 & irrational == 0 & female == 0, vce(robust) * reflection with further controls on the hormone quality (luncht foodtime drinktime mouth brush blast dst2010) probit reflection ntestf nestf nprogf ncortf luncht foodtime drinktime mouth brush blast dst2010 pill female age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & preg == 0 & consistent == 1 & irrational == 0, vce(robust) * probit reflection ntestf nestf nprogf ncortf luncht foodtime drinktime mouth brush blast dst2010 pill age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & preg == 0 & consistent == 1 & irrational == 0 & female == 1, vce(robust) * probit reflection ntestf nestf nprogf ncortf luncht foodtime drinktime mouth brush blast dst2010 pill pilltestf pillestf pillprogf age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & preg == 0 & consistent == 1 & irrational == 0 & female == 1, vce(robust) probit reflection ntestf nestf nprogf ncortf luncht foodtime drinktime mouth brush blast dst2010 age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & preg == 0 & consistent == 1 & irrational == 0 & female == 0, vce(robust) * reflection with further controls for dietary preferences (smoking vegetarian soybean) probit reflection ntestf nestf nprogf ncortf pill female smoking vegetarian soybean age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & preg == 0 & consistent == 1 & irrational == 0, vce(robust) probit reflection ntestf nestf nprogf ncortf pill smoking vegetarian soybean age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & preg == 0 & consistent == 1 & irrational == 0 & female == 1, vce(robust) probit reflection ntestf nestf nprogf ncortf pill pilltestf pillestf pillprogf smoking vegetarian soybean age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & preg == 0 & consistent == 1 & irrational == 0 & female == 1, vce(robust) probit reflection ntestf nestf nprogf ncortf smoking vegetarian soybean age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & preg == 0 & consistent == 1 & irrational == 0 & female == 0, vce(robust) * reflection female with further controls for menstrual cycle phases probit reflection ntestf nestf nprogf ncortf nfolic_fol nperi_fol nlut_fol npremens_fol pill age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & preg == 0 & consistent == 1 & irrational == 0 & female == 1, vce(robust) *probit reflection ntestf nestf nprogf ncortf nfolic_fol nperi_fol nlut_fol npremens_fol pill pilltestf pillestf pillprogf age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & preg == 0 & consistent == 1 & irrational == 0 & female == 1, vce(robust) * reflection with average hormones probit reflection ntestav nestav nprogav ncortav pill female smoking vegetarian soybean age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & preg == 0 & consistent == 1 & irrational == 0, vce(robust) probit reflection ntestav nestav nprogav ncortav pill smoking vegetarian soybean age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & preg == 0 & consistent == 1 & irrational == 0 & female == 1, vce(robust) probit reflection ntestav nestav nprogav ncortav pill pilltestf pillestf pillprogf smoking vegetarian soybean age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & preg == 0 & consistent == 1 & irrational == 0 & female == 1, vce(robust) probit reflection ntestav nestav nprogav ncortav smoking vegetarian soybean age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & preg == 0 & consistent == 1 & irrational == 0 & female == 0, vce(robust) * reflection stepwise forward * stepwise probit reflection ntestf nestf nprogf ncortf luncht foodtime drinktime mouth brush blast pill female smoking vegetarian soybean dr bmi siblings workout livalone dating homobi nconception age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & preg == 0 & consistent == 1 & irrational == 0, pe(0.05) * stepwise probit reflection ntestf nestf nprogf ncortf nfolic_fol nperi_fol nlut_fol npremens_fol luncht foodtime drinktime mouth brush blast pill smoking vegetarian soybean dr bmi siblings workout livalone dating homobi nconception age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & preg == 0 & consistent == 1 & irrational == 0 & female == 1, pe(0.05) * stepwise probit reflection ntestf nestf nprogf ncortf nfolic_fol nperi_fol nlut_fol npremens_fol luncht foodtime drinktime mouth brush blast pill pilltestf pillestf pillprogf smoking vegetarian soybean dr bmi siblings workout livalone dating homobi nconception age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & preg == 0 & consistent == 1 & irrational == 0 & female == 1, pe(0.05) * stepwise probit reflection ntestf nestf nprogf ncortf luncht foodtime drinktime mouth brush blast smoking vegetarian soybean dr bmi siblings workout livalone dating homobi age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & preg == 0 & consistent == 1 & irrational == 0 & female == 0, pe(0.05) * reflection stepwise backward * stepwise probit reflection ntestf nestf nprogf ncortf luncht foodtime drinktime mouth brush blast pill female smoking vegetarian soybean dr bmi siblings workout livalone dating homobi nconception age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & preg == 0 & consistent == 1 & irrational == 0, pr(0.05) * stepwise probit reflection ntestf nestf nprogf ncortf nfolic_fol nperi_fol nlut_fol npremens_fol luncht foodtime drinktime mouth brush blast pill smoking vegetarian soybean dr bmi siblings workout livalone dating homobi nconception age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & preg == 0 & consistent == 1 & irrational == 0 & female == 1, pr(0.05) * stepwise probit reflection ntestf nestf nprogf ncortf nfolic_fol nperi_fol nlut_fol npremens_fol luncht foodtime drinktime mouth brush blast pill pilltestf pillestf pillprogf smoking vegetarian soybean dr bmi siblings workout livalone dating homobi nconception age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & preg == 0 & consistent == 1 & irrational == 0 & female == 1, pr(0.05) * stepwise probit reflection ntestf nestf nprogf ncortf luncht foodtime drinktime mouth brush blast smoking vegetarian soybean dr bmi siblings workout livalone dating homobi age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & preg == 0 & consistent == 1 & irrational == 0 & female == 0, pr(0.05) * consistent & hormones ************************** eststo c3: logit consistent ntestf nestf nprogf ncortf pill female age asian otherethn gpa sciothersci econ socsci humanities if period == 30 & preg == 0, vce(robust) eststo fc3: logit consistent ntestf nestf nprogf ncortf pill age asian otherethn gpa sciothersci econ socsci humanities if period == 30 & preg == 0 & female == 1, vce(robust) eststo fc4: logit consistent ntestf nestf nprogf ncortf pill pilltestf pillestf pillprogf age asian otherethn gpa sciothersci econ socsci humanities if period == 30 & preg == 0 & female == 1, vce(robust) eststo mc3: logit consistent ntestf nestf nprogf ncortf age asian otherethn gpa sciothersci econ socsci if period == 30 & preg == 0 & female == 0, vce(robust) estout c3 fc3 fc4 mc3 using "C:\bcschipper\burkhard\research\hormone\consistent.tex", replace /// cells(b(star fmt(%9.4f)) se(par fmt(%9.4f))) stats(N r2, fmt(%9.0g %9.4f) labels("\textit{Number of Observations}" "R\textsuperscript{2}") layout("\multicolumn{1}{c}{\quad@}")) label /// starlevels(\text{*} 0.05 \text{**} 0.01) keep(ntestf nestf nprogf ncortf pill pilltestf pillestf pillprogf female age asian otherethn gpa sciothersci econ socsci humanities) /// mlabels((A3) (A3F) (A4F) (A3M), prefix(\multicolumn{1}{c}{) suffix(})) order(age asian otherethn gpa sciothersci econ socsci humanities female ntestf nestf nprogf ncortf pill pilltestf pillestf pillprogf) /// style(tex) collabels(,none) prehead("\begin{tabular}{l*{@span}{d}}\toprule") /// posthead(\midrule) prefoot(\bottomrule) postfoot("\addlinespace[.1cm]\multicolumn{5}{l}{Robust standard errors in parentheses; Significance levels: * 5\%; ** 1\%.} \\" "\end{tabular}") notype estout c3 fc3 fc4 mc3 using "C:\bcschipper\burkhard\research\hormone\s_consistent.tex", replace /// cells(b(star fmt(%9.4f)) se(par fmt(%9.4f))) stats(N r2, fmt(%9.0g %9.4f) labels("\textit{Number of Observations}" "R\textsuperscript{2}") layout("\multicolumn{1}{c}{\quad@}")) label /// starlevels(\text{*} 0.05 \text{**} 0.01) keep(ntestf nestf nprogf ncortf pill pilltestf pillestf pillprogf female) /// mlabels((A3) (A3F) (A4F) (A3M), prefix(\multicolumn{1}{c}{) suffix(})) order(female ntestf nestf nprogf ncortf pill pilltestf pillestf pillprogf) /// style(tex) collabels(,none) prehead("\begin{tabular}{l*{@span}{d}}\toprule") /// posthead(\midrule) prefoot(\bottomrule) postfoot("\addlinespace[.1cm]\multicolumn{5}{l}{Robust standard errors in parentheses; Significance levels: * 5\%; ** 1\%.} \\" /// "\multicolumn{5}{l}{We suppress from the report the coefficients for Age, Asian, Other, GPA, and Majors of Study.} \\""\end{tabular}") notype * Robustness checks * consistent with probit probit consistent ntestf nestf nprogf ncortf pill female age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & preg == 0, vce(robust) probit consistent ntestf nestf nprogf ncortf pill age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & preg == 0 & female == 1, vce(robust) probit consistent ntestf nestf nprogf ncortf pill pilltestf pillestf pillprogf age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & preg == 0 & female == 1, vce(robust) probit consistent ntestf nestf nprogf ncortf age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & preg == 0 & female == 0, vce(robust) * standard errors clusters on session level logit consistent ntestf nestf nprogf ncortf pill female age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & preg == 0, robust cluster(session) logit consistent ntestf nestf nprogf ncortf pill age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & preg == 0 & female == 1, robust cluster(session) logit consistent ntestf nestf nprogf ncortf pill pilltestf pillestf pillprogf age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & preg == 0 & female == 1, robust cluster(session) logit consistent ntestf nestf nprogf ncortf age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & preg == 0 & female == 0, robust cluster(session) * session dummies logit consistent sessiondum* ntestf nestf nprogf ncortf pill female age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & preg == 0, vce(robust) *logit consistent sessiondum* ntestf nestf nprogf ncortf pill age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & preg == 0 & female == 1, vce(robust) logit consistent sessiondum* ntestf nestf nprogf ncortf pill pilltestf pillestf pillprogf age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & preg == 0 & female == 1, vce(robust) logit consistent sessiondum* ntestf nestf nprogf ncortf age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & preg == 0 & female == 0, vce(robust) * session fixed effects xtset, clear xtset datum xtlogit consistent ntestf nestf nprogf ncortf pill female age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & preg == 0, fe *xtlogit consistent ntestf nestf nprogf ncortf pill age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & preg == 0 & female == 1, fe *xtlogit consistent ntestf nestf nprogf ncortf pill pilltestf pillestf pillprogf age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & preg == 0 & female == 1, fe *xtlogit consistent ntestf nestf nprogf ncortf age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & preg == 0 & female == 0, fe * drop controls logit consistent ntestf nestf nprogf ncortf pill female gpa if period == 30 & preg == 0, vce(robust) logit consistent ntestf nestf nprogf ncortf pill gpa if period == 30 & preg == 0 & female == 1, vce(robust) logit consistent ntestf nestf nprogf ncortf pill pilltestf pillestf pillprogf gpa if period == 30 & preg == 0 & female == 1, vce(robust) logit consistent ntestf nestf nprogf ncortf gpa if period == 30 & preg == 0 & female == 0, vce(robust) * consistent with further demographic controls (bmi siblings workout livalone dating homobi nconception) logit consistent ntestf nestf nprogf ncortf pill female bmi siblings workout livalone dating homobi nconception age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & preg == 0, vce(robust) logit consistent ntestf nestf nprogf ncortf pill bmi siblings workout livalone dating homobi nconception age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & preg == 0 & female == 1, vce(robust) logit consistent ntestf nestf nprogf ncortf pill pilltestf pillestf pillprogf bmi siblings workout livalone dating homobi nconception age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & preg == 0 & female == 1, vce(robust) logit consistent ntestf nestf nprogf ncortf bmi siblings workout livalone dating homobi age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & preg == 0 & female == 0, vce(robust) * consistent with further controls on the hormone quality (luncht foodtime drinktime mouth brush blast dst2010) logit consistent ntestf nestf nprogf ncortf luncht foodtime drinktime mouth brush blast dst2010 pill female age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & preg == 0, vce(robust) logit consistent ntestf nestf nprogf ncortf luncht foodtime drinktime mouth brush blast dst2010 pill age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & preg == 0 & female == 1, vce(robust) logit consistent ntestf nestf nprogf ncortf luncht foodtime drinktime mouth brush blast dst2010 pill pilltestf pillestf pillprogf age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & preg == 0 & female == 1, vce(robust) logit consistent ntestf nestf nprogf ncortf luncht foodtime drinktime mouth brush blast dst2010 age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & preg == 0 & female == 0, vce(robust) * consistent with further controls for dietary preferences (smoking vegetarian soybean) logit consistent ntestf nestf nprogf ncortf pill female smoking vegetarian soybean age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & preg == 0, vce(robust) logit consistent ntestf nestf nprogf ncortf pill smoking vegetarian soybean age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & preg == 0 & female == 1, vce(robust) logit consistent ntestf nestf nprogf ncortf pill pilltestf pillestf pillprogf smoking vegetarian soybean age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & preg == 0 & female == 1, vce(robust) logit consistent ntestf nestf nprogf ncortf smoking vegetarian soybean age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & preg == 0 & female == 0, vce(robust) * consistent female with further controls for menstrual cycle phases logit consistent ntestf nestf nprogf ncortf nfolic_fol nperi_fol nlut_fol npremens_fol pill age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & preg == 0 & female == 1, vce(robust) logit consistent ntestf nestf nprogf ncortf nfolic_fol nperi_fol nlut_fol npremens_fol pill pilltestf pillestf pillprogf age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & preg == 0 & female == 1, vce(robust) * consistent with average hormones logit consistent ntestav nestav nprogav ncortav pill female smoking vegetarian soybean age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & preg == 0, vce(robust) logit consistent ntestav nestav nprogav ncortav pill smoking vegetarian soybean age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & preg == 0 & female == 1, vce(robust) logit consistent ntestav nestav nprogav ncortav pill pilltestf pillestf pillprogf smoking vegetarian soybean age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & preg == 0 & female == 1, vce(robust) logit consistent ntestav nestav nprogav ncortav smoking vegetarian soybean age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & preg == 0 & female == 0, vce(robust) * consistent stepwise forward * stepwise logit consistent ntestf nestf nprogf ncortf luncht foodtime drinktime mouth brush blast pill female smoking vegetarian soybean dr bmi siblings workout livalone dating homobi nconception age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & preg == 0, pe(0.05) * stepwise logit consistent ntestf nestf nprogf ncortf nfolic_fol nperi_fol nlut_fol npremens_fol luncht foodtime drinktime mouth brush blast pill smoking vegetarian soybean dr bmi siblings workout livalone dating homobi nconception age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & preg == 0 & female == 1, pe(0.05) * stepwise logit consistent ntestf nestf nprogf ncortf nfolic_fol nperi_fol nlut_fol npremens_fol luncht foodtime drinktime mouth brush blast pill pilltestf pillestf pillprogf smoking vegetarian soybean dr bmi siblings workout livalone dating homobi nconception age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & preg == 0 & female == 1, pe(0.05) * stepwise logit consistent ntestf nestf nprogf ncortf luncht foodtime drinktime mouth brush blast smoking vegetarian soybean dr bmi siblings workout livalone dating homobi age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & preg == 0 & female == 0, pe(0.05) * consistent stepwise backward * stepwise logit consistent ntestf nestf nprogf ncortf luncht foodtime drinktime mouth brush blast pill female smoking vegetarian soybean dr bmi siblings workout livalone dating homobi nconception age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & preg == 0, pr(0.05) * stepwise logit consistent ntestf nestf nprogf ncortf nfolic_fol nperi_fol nlut_fol npremens_fol luncht foodtime drinktime mouth brush blast pill smoking vegetarian soybean dr bmi siblings workout livalone dating homobi nconception age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & preg == 0 & female == 1, pr(0.05) * stepwise logit consistent ntestf nestf nprogf ncortf nfolic_fol nperi_fol nlut_fol npremens_fol luncht foodtime drinktime mouth brush blast pill pilltestf pillestf pillprogf smoking vegetarian soybean dr bmi siblings workout livalone dating homobi nconception age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & preg == 0 & female == 1, pr(0.05) * stepwise logit consistent ntestf nestf nprogf ncortf luncht foodtime drinktime mouth brush blast smoking vegetarian soybean dr bmi siblings workout livalone dating homobi age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & preg == 0 & female == 0, pr(0.05) * Risk averison and hormones by gender ***************************************** eststo fg3: reg riskgain ntestf nestf nprogf ncortf pill age asian otherethn gpa sciothersci econ socsci humanities if period == 30 & preg == 0 & female == 1 & rgconsistent == 1 & girrational == 0, vce(robust) eststo fg4: reg riskgain ntestf nestf nprogf ncortf pill pilltestf pillestf pillprogf age asian otherethn gpa sciothersci econ socsci humanities if period == 30 & preg == 0 & female == 1 & rgconsistent == 1 & girrational == 0, vce(robust) eststo fl3: reg riskloss ntestf nestf nprogf ncortf pill age asian otherethn gpa sciothersci econ socsci humanities if period == 30 & preg == 0 & female == 1 & rlconsistent == 1 & lirrational == 0, vce(robust) eststo mg3: reg riskgain ntestf nestf nprogf ncortf age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & preg == 0 & female == 0 & rgconsistent == 1 & girrational == 0, vce(robust) eststo mg4: reg riskgain ntestf nestf nprogf ncortf age asian otherethn gpa math sciothersci econ socsci humanities lunchtestf lunchcortf if period == 30 & preg == 0 & female == 0 & rgconsistent == 1 & girrational == 0, vce(robust) eststo ml3: reg riskloss ntestf nestf nprogf ncortf age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & preg == 0 & female == 0 & rlconsistent == 1 & lirrational == 0, vce(robust) estout fg3 fg4 fl3 mg3 mg4 ml3 using "C:\bcschipper\burkhard\research\hormone\genderhorm.tex", replace /// cells(b(star fmt(%9.4f)) se(par fmt(%9.4f))) stats(N r2, fmt(%9.0g %9.4f) labels("\textit{Number of Observations}" "R\textsuperscript{2}") layout("\multicolumn{1}{c}{\quad@}")) label /// starlevels(\text{*} 0.05 \text{**} 0.01) keep(ntestf nestf nprogf ncortf pill pilltestf pillestf pillprogf age asian otherethn gpa math sciothersci econ socsci humanities lunchtestf lunchcortf) /// mlabels((G3F) (G4F) (L3F) (G3M) (G4M) (L3M), prefix(\multicolumn{1}{c}{) suffix(})) order(age asian otherethn gpa math sciothersci econ socsci humanities ntestf nestf nprogf ncortf pill pilltestf pillestf pillprogf lunchtestf lunchcortf) /// style(tex) collabels(,none) prehead("\begin{tabular}{l*{@span}{d}}\toprule") /// posthead(\midrule) prefoot(\bottomrule) postfoot("\addlinespace[.1cm]\multicolumn{7}{l}{Robust standard errors in parentheses; Significance levels: * 5\%; ** 1\%.} \\" "\end{tabular}") notype estout fg3 fg4 fl3 mg3 mg4 ml3 using "C:\bcschipper\burkhard\research\hormone\s_genderhorm.tex", replace /// cells(b(star fmt(%9.4f)) se(par fmt(%9.4f))) stats(N r2, fmt(%9.0g %9.4f) labels("\textit{Number of Observations}" "R\textsuperscript{2}") layout("\multicolumn{1}{c}{\quad@}")) label /// starlevels(\text{*} 0.05 \text{**} 0.01) keep(ntestf nestf nprogf ncortf pill pilltestf pillestf pillprogf lunchtestf lunchcortf) /// mlabels((G3F) (G4F) (L3F) (G3M) (G4M) (L3M), prefix(\multicolumn{1}{c}{) suffix(})) order(ntestf nestf nprogf ncortf pill pilltestf pillestf pillprogf lunchtestf lunchcortf) /// style(tex) collabels(,none) prehead("\begin{tabular}{l*{@span}{d}}\toprule") /// posthead(\midrule) prefoot(\bottomrule) postfoot("\addlinespace[.1cm]\multicolumn{7}{l}{Robust standard errors in parentheses; Significance levels: * 5\%; ** 1\%.} \\" /// "\multicolumn{7}{l}{We suppress from the report the coefficients for Age, Asian, Other, GPA, and Majors of Study.} \\""\end{tabular}") notype * Robustness checks * by gender with standard errors clustered on session reg riskgain ntestf nestf nprogf ncortf pill age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & preg == 0 & female == 1 & rgconsistent == 1 & girrational == 0, robust cluster(session) reg riskgain ntestf nestf nprogf ncortf pill pilltestf pillestf pillprogf age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & preg == 0 & female == 1 & rgconsistent == 1 & girrational == 0, robust cluster(session) reg riskloss ntestf nestf nprogf ncortf pill age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & preg == 0 & female == 1 & rlconsistent == 1 & lirrational == 0, robust cluster(session) reg riskloss ntestf nestf nprogf ncortf pill pilltestf pillestf pillprogf age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & preg == 0 & female == 1 & rlconsistent == 1 & lirrational == 0, robust cluster(session) reg riskgain ntestf nestf nprogf ncortf pill age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & preg == 0 & female == 0 & rgconsistent == 1 & girrational == 0, robust cluster(session) reg riskloss ntestf nestf nprogf ncortf pill age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & preg == 0 & female == 0 & rlconsistent == 1 & lirrational == 0, robust cluster(session) * sessiondummies reg riskgain sessiondum* ntestf nestf nprogf ncortf pill age asian otherethn gpa sciothersci econ socsci humanities if period == 30 & preg == 0 & female == 1 & rgconsistent == 1 & girrational == 0, vce(robust) reg riskgain sessiondum* ntestf nestf nprogf ncortf pill pilltestf pillestf pillprogf age asian otherethn gpa sciothersci econ socsci humanities if period == 30 & preg == 0 & female == 1 & rgconsistent == 1 & girrational == 0, vce(robust) reg riskloss sessiondum* ntestf nestf nprogf ncortf pill age asian otherethn gpa sciothersci econ socsci humanities if period == 30 & preg == 0 & female == 1 & rlconsistent == 1 & lirrational == 0, vce(robust) reg riskgain sessiondum* ntestf nestf nprogf ncortf age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & preg == 0 & female == 0 & rgconsistent == 1 & girrational == 0, vce(robust) reg riskgain sessiondum* ntestf nestf nprogf ncortf age asian otherethn gpa math sciothersci econ socsci humanities lunchtestf lunchcortf if period == 30 & preg == 0 & female == 0 & rgconsistent == 1 & girrational == 0, vce(robust) reg riskloss sessiondum* ntestf nestf nprogf ncortf age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & preg == 0 & female == 0 & rlconsistent == 1 & lirrational == 0, vce(robust) * session fixed effects xtset, clear xtset datum xtreg riskgain ntestf nestf nprogf ncortf pill age asian otherethn gpa sciothersci econ socsci humanities if period == 30 & preg == 0 & female == 1 & rgconsistent == 1 & girrational == 0, fe xtreg riskgain ntestf nestf nprogf ncortf pill pilltestf pillestf pillprogf age asian otherethn gpa sciothersci econ socsci humanities if period == 30 & preg == 0 & female == 1 & rgconsistent == 1 & girrational == 0, fe xtreg riskloss ntestf nestf nprogf ncortf pill age asian otherethn gpa sciothersci econ socsci humanities if period == 30 & preg == 0 & female == 1 & rlconsistent == 1 & lirrational == 0, fe xtreg riskgain ntestf nestf nprogf ncortf age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & preg == 0 & female == 0 & rgconsistent == 1 & girrational == 0, fe xtreg riskgain ntestf nestf nprogf ncortf age asian otherethn gpa math sciothersci econ socsci humanities lunchtestf lunchcortf if period == 30 & preg == 0 & female == 0 & rgconsistent == 1 & girrational == 0, fe xtreg riskloss ntestf nestf nprogf ncortf age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & preg == 0 & female == 0 & rlconsistent == 1 & lirrational == 0, fe * minimal demographics controls reg riskgain ntestf nestf nprogf ncortf pill gpa if period == 30 & preg == 0 & female == 1 & rgconsistent == 1 & girrational == 0, vce(robust) reg riskgain ntestf nestf nprogf ncortf pill pilltestf pillestf pillprogf gpa sciothersci econ socsci humanities if period == 30 & preg == 0 & female == 1 & rgconsistent == 1 & girrational == 0, vce(robust) reg riskloss ntestf nestf nprogf ncortf pill gpa if period == 30 & preg == 0 & female == 1 & rlconsistent == 1 & lirrational == 0, vce(robust) reg riskgain ntestf nestf nprogf ncortf gpa if period == 30 & preg == 0 & female == 0 & rgconsistent == 1 & girrational == 0, vce(robust) reg riskgain ntestf nestf nprogf ncortf gpa lunchtestf lunchcortf if period == 30 & preg == 0 & female == 0 & rgconsistent == 1 & girrational == 0, vce(robust) reg riskloss ntestf nestf nprogf ncortf gpa if period == 30 & preg == 0 & female == 0 & rlconsistent == 1 & lirrational == 0, vce(robust) * by gender ologit ologit riskgain ntestf nestf nprogf ncortf pill age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & preg == 0 & female == 1 & rgconsistent == 1 & girrational == 0, vce(robust) ologit riskgain ntestf nestf nprogf ncortf pill pilltestf pillestf pillprogf age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & preg == 0 & female == 1 & rgconsistent == 1 & girrational == 0, vce(robust) ologit riskloss ntestf nestf nprogf ncortf pill age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & preg == 0 & female == 1 & rlconsistent == 1 & lirrational == 0, vce(robust) ologit riskloss ntestf nestf nprogf ncortf pill pilltestf pillestf pillprogf age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & preg == 0 & female == 1 & rlconsistent == 1 & lirrational == 0, vce(robust) ologit riskgain ntestf nestf nprogf ncortf pill age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & preg == 0 & female == 0 & rgconsistent == 1 & girrational == 0, vce(robust) ologit riskloss ntestf nestf nprogf ncortf pill age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & preg == 0 & female == 0 & rlconsistent == 1 & lirrational == 0, vce(robust) * by gender oprobit oprobit riskgain ntestf nestf nprogf ncortf pill age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & preg == 0 & female == 1 & rgconsistent == 1 & girrational == 0, vce(robust) oprobit riskgain ntestf nestf nprogf ncortf pill pilltestf pillestf pillprogf age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & preg == 0 & female == 1 & rgconsistent == 1 & girrational == 0, vce(robust) oprobit riskloss ntestf nestf nprogf ncortf pill age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & preg == 0 & female == 1 & rlconsistent == 1 & lirrational == 0, vce(robust) oprobit riskloss ntestf nestf nprogf ncortf pill pilltestf pillestf pillprogf age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & preg == 0 & female == 1 & rlconsistent == 1 & lirrational == 0, vce(robust) oprobit riskgain ntestf nestf nprogf ncortf pill age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & preg == 0 & female == 0 & rgconsistent == 1 & girrational == 0, vce(robust) oprobit riskloss ntestf nestf nprogf ncortf pill age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & preg == 0 & female == 0 & rlconsistent == 1 & lirrational == 0, vce(robust) * by gender with further demographic controls (dr bmi siblings workout livalone dating homobi) reg riskgain ntestf nestf nprogf ncortf pill dr bmi siblings workout livalone dating homobi age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & preg == 0 & female == 1 & rgconsistent == 1 & girrational == 0, vce(robust) reg riskgain ntestf nestf nprogf ncortf pill pilltestf pillestf pillprogf dr bmi siblings workout livalone dating homobi age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & preg == 0 & female == 1 & rgconsistent == 1 & girrational == 0, vce(robust) reg riskloss ntestf nestf nprogf ncortf pill dr bmi siblings workout livalone dating homobi age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & preg == 0 & female == 1 & rlconsistent == 1 & lirrational == 0, vce(robust) reg riskloss ntestf nestf nprogf ncortf pill pilltestf pillestf pillprogf dr bmi siblings workout livalone dating homobi age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & preg == 0 & female == 1 & rlconsistent == 1 & lirrational == 0, vce(robust) reg riskgain ntestf nestf nprogf ncortf pill dr bmi siblings workout livalone dating homobi age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & preg == 0 & female == 0 & rgconsistent == 1 & girrational == 0, vce(robust) reg riskloss ntestf nestf nprogf ncortf pill dr bmi siblings workout livalone dating homobi age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & preg == 0 & female == 0 & rlconsistent == 1 & lirrational == 0, vce(robust) * by gender with further controls on the hormone quality (luncht foodtime drinktime mouth brush blast dst2010) reg riskgain ntestf nestf nprogf ncortf luncht foodtime drinktime mouth brush blast dst2010 pill age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & preg == 0 & female == 1 & rgconsistent == 1 & girrational == 0, vce(robust) reg riskgain ntestf nestf nprogf ncortf luncht foodtime drinktime mouth brush blast dst2010 pill pilltestf pillestf pillprogf age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & preg == 0 & female == 1 & rgconsistent == 1 & girrational == 0, vce(robust) reg riskloss ntestf nestf nprogf ncortf luncht foodtime drinktime mouth brush blast dst2010 pill age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & preg == 0 & female == 1 & rlconsistent == 1 & lirrational == 0, vce(robust) reg riskloss ntestf nestf nprogf ncortf luncht foodtime drinktime mouth brush blast dst2010 pill pilltestf pillestf pillprogf age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & preg == 0 & female == 1 & rlconsistent == 1 & lirrational == 0, vce(robust) reg riskgain ntestf nestf nprogf ncortf luncht foodtime drinktime mouth brush blast dst2010 pill age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & preg == 0 & female == 0 & rgconsistent == 1 & girrational == 0, vce(robust) reg riskloss ntestf nestf nprogf ncortf luncht foodtime drinktime mouth brush blast dst2010 pill age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & preg == 0 & female == 0 & rlconsistent == 1 & lirrational == 0, vce(robust) * by gender with further controls for dietary preferences (smoking vegetarian soybean) reg riskgain ntestf nestf nprogf ncortf smoking vegetarian soybean pill age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & preg == 0 & female == 1 & rgconsistent == 1 & girrational == 0, vce(robust) reg riskgain ntestf nestf nprogf ncortf smoking vegetarian soybean pill pilltestf pillestf pillprogf age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & preg == 0 & female == 1 & rgconsistent == 1 & girrational == 0, vce(robust) reg riskloss ntestf nestf nprogf ncortf smoking vegetarian soybean pill age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & preg == 0 & female == 1 & rlconsistent == 1 & lirrational == 0, vce(robust) reg riskloss ntestf nestf nprogf ncortf smoking vegetarian soybean pill pilltestf pillestf pillprogf age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & preg == 0 & female == 1 & rlconsistent == 1 & lirrational == 0, vce(robust) reg riskgain ntestf nestf nprogf ncortf smoking vegetarian soybean pill age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & preg == 0 & female == 0 & rgconsistent == 1 & girrational == 0, vce(robust) reg riskloss ntestf nestf nprogf ncortf smoking vegetarian soybean pill age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & preg == 0 & female == 0 & rlconsistent == 1 & lirrational == 0, vce(robust) * by gender female with further controls for menstrual cycle phases reg riskgain ntestf nestf nprogf ncortf nfolic_fol nperi_fol nlut_fol npremens_fol pill age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & preg == 0 & female == 1 & rgconsistent == 1 & girrational == 0, vce(robust) reg riskgain ntestf nestf nprogf ncortf nfolic_fol nperi_fol nlut_fol npremens_fol pill pilltestf pillestf pillprogf age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & preg == 0 & female == 1 & rgconsistent == 1 & girrational == 0, vce(robust) reg riskloss ntestf nestf nprogf ncortf nfolic_fol nperi_fol nlut_fol npremens_fol pill age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & preg == 0 & female == 1 & rlconsistent == 1 & lirrational == 0, vce(robust) reg riskloss ntestf nestf nprogf ncortf nfolic_fol nperi_fol nlut_fol npremens_fol pill pilltestf pillestf pillprogf age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & preg == 0 & female == 1 & rlconsistent == 1 & lirrational == 0, vce(robust) * by gender with average hormone levels reg riskgain ntestav nestav nprogav ncortav pill age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & preg == 0 & female == 1 & rgconsistent == 1 & girrational == 0, vce(robust) reg riskgain ntestav nestav nprogav ncortav pill pilltestf pillestf pillprogf age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & preg == 0 & female == 1 & rgconsistent == 1 & girrational == 0, vce(robust) reg riskloss ntestav nestav nprogav ncortav pill age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & preg == 0 & female == 1 & rlconsistent == 1 & lirrational == 0, vce(robust) reg riskloss ntestav nestav nprogav ncortav pill pilltestf pillestf pillprogf age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & preg == 0 & female == 1 & rlconsistent == 1 & lirrational == 0, vce(robust) reg riskgain ntestav nestav nprogav ncortav pill age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & preg == 0 & female == 0 & rgconsistent == 1 & girrational == 0, vce(robust) reg riskloss ntestav nestav nprogav ncortav pill age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & preg == 0 & female == 0 & rlconsistent == 1 & lirrational == 0, vce(robust) * by gender stepwise forward stepwise reg riskgain ntestf nestf nprogf ncortf nfolic_fol nperi_fol nlut_fol npremens_fol luncht foodtime drinktime mouth brush blast pill dr bmi siblings workout livalone dating homobi nconception smoking vegetarian soybean age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & preg == 0 & rgconsistent == 1 & girrational == 0 & female == 1, pe(0.05) stepwise reg riskloss ntestf nestf nprogf ncortf nfolic_fol nperi_fol nlut_fol npremens_fol luncht foodtime drinktime mouth brush blast pill dr bmi siblings workout livalone dating homobi nconception smoking vegetarian soybean age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & preg == 0 & rlconsistent == 1 & lirrational == 0 & female == 1, pe(0.05) stepwise reg riskgain ntestf nestf nprogf ncortf luncht foodtime drinktime mouth brush blast pill dr bmi siblings workout livalone dating homobi nconception smoking vegetarian soybean age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & preg == 0 & rgconsistent == 1 & girrational == 0 & female == 0, pe(0.05) stepwise reg riskloss ntestf nestf nprogf ncortf luncht foodtime drinktime mouth brush blast pill dr bmi siblings workout livalone dating homobi nconception smoking vegetarian soybean age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & preg == 0 & rlconsistent == 1 & lirrational == 0 & female == 0, pe(0.05) * by gender stepwise backward stepwise reg riskgain ntestf nestf nprogf ncortf nfolic_fol nperi_fol nlut_fol npremens_fol luncht foodtime drinktime mouth brush blast pill dr bmi siblings workout livalone dating homobi nconception smoking vegetarian soybean age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & preg == 0 & rgconsistent == 1 & girrational == 0 & female == 1, pr(0.05) stepwise reg riskloss ntestf nestf nprogf ncortf nfolic_fol nperi_fol nlut_fol npremens_fol luncht foodtime drinktime mouth brush blast pill dr bmi siblings workout livalone dating homobi nconception smoking vegetarian soybean age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & preg == 0 & rlconsistent == 1 & lirrational == 0 & female == 1, pr(0.05) stepwise reg riskgain ntestf nestf nprogf ncortf luncht foodtime drinktime mouth brush blast pill dr bmi siblings workout livalone dating homobi nconception smoking vegetarian soybean age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & preg == 0 & rgconsistent == 1 & girrational == 0 & female == 0, pr(0.05) stepwise reg riskloss ntestf nestf nprogf ncortf luncht foodtime drinktime mouth brush blast pill dr bmi siblings workout livalone dating homobi nconception smoking vegetarian soybean age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & preg == 0 & rlconsistent == 1 & lirrational == 0 & female == 0, pr(0.05) ************************************* * Analysis Risk \& Menstrual Cycle ************************************* sum nmenst nfolic nperi nlut npremens nduration pill pilluse pillbreak if period == 30 & preg == 0 & female == 1 & new == 1 sum nmenst_adj nfolic_adj nperi_adj nlut_adj npremens_adj nduration if period == 30 & preg == 0 & female == 1 & durationmissing == 0 & new == 1 sum nmenst_fol nfolic_fol nperi_fol nlut_fol npremens_fol nduration if period == 30 & preg == 0 & female == 1 & durationmissing == 0 & new == 1 * full eststo g5: reg riskgain age asian otherethn gpa math sciothersci econ socsci humanities nmenst_fol nfolic_fol nperi_fol nlut_fol npremens_fol pill if period == 30 & preg == 0 & rgconsistent == 1 & girrational == 0, vce(robust) eststo l5: reg riskloss age asian otherethn gpa math sciothersci econ socsci humanities nmenst_fol nfolic_fol nperi_fol nlut_fol npremens_fol pill if period == 30 & preg == 0 & rlconsistent == 1 & lirrational == 0, vce(robust) eststo r5: logit reflection age asian otherethn gpa math sciothersci econ socsci humanities nmenst_fol nfolic_fol npremens_fol pill if period == 30 & preg == 0 & consistent == 1 & irrational == 0, vce(robust) eststo c5: logit consistent age asian otherethn gpa math sciothersci econ socsci humanities nmenst_fol nfolic_fol nperi_fol nlut_fol npremens_fol pill if period == 30 & preg == 0, vce(robust) estout g5 l5 r5 c5 using "C:\bcschipper\burkhard\research\hormone\menstruation.tex", replace /// cells(b(star fmt(%9.4f)) se(par fmt(%9.4f))) stats(N r2, fmt(%9.0g %9.4f) labels("\textit{Number of Observations}" "R\textsuperscript{2}") layout("\multicolumn{1}{c}{\quad@}")) label /// starlevels(\text{*} 0.05 \text{**} 0.01) keep(nmenst_fol nfolic_fol nperi_fol nlut_fol npremens_fol pill) /// mlabels((G5) (L5) (R5) (A5), prefix(\multicolumn{1}{c}{) suffix(})) order(nmenst_fol nfolic_fol nperi_fol nlut_fol npremens_fol pill) /// style(tex) collabels(,none) prehead("\begin{tabular}{l*{@span}{d}}\toprule") /// posthead(\midrule) prefoot(\bottomrule) postfoot("\addlinespace[.1cm]\multicolumn{5}{l}{Robust standard errors in parentheses; Significance levels: * 5\%; ** 1\%.} \\" /// "\multicolumn{5}{l}{We suppress from the report coefficients for Age, Asian, Other, GPA, and Majors of Studies. } \\" "\end{tabular}") notype * Starred in 2022 run of the program. Stata ask for Version 16 now to run parmest/parmby even though it run fine back in 2015 with Stata 14 *parmby "reg riskgain age asian otherethn gpa math sciothersci econ socsci humanities nmenst_fol nfolic_fol nperi_fol nlut_fol npremens_fol pill if period == 30 & preg == 0 & rgconsistent == 1 & girrational == 0, vce(robust)", saving(junk, replace) *preserve *use junk, clear *drop if parm=="age" | parm=="asian" | parm=="otherethn" | parm=="gpa" | parm=="math" | parm=="sciothersci" | parm=="econ" | parm=="socsci" | parm=="humanities" | parm=="pill" | parm=="_cons" *gen days = . *la var days "Menstrual Phase" *replace days = 2.75 if parm== "nmenst_fol" *replace days = 9.5 if parm== "nfolic_fol" *replace days = 15 if parm== "nperi_fol" *replace days = 20 if parm== "nlut_fol" *replace days = 25.75 if parm== "npremens_fol" *twoway (scatter estimate days) (rcap min95 max95 days, sort), legend(on order(1 "Female/Male Difference" 2 "95% Confidence Interval")) xtitle(Menstrual Phase ,margin(small)) xscale(range(0 28)) xlabel(2.75 "Menstrual " 9.5 "Folicular" 15 "Peri-Ovulatory" 20 "Luteal" 25.75 "Pre-Menstrual", labsize(small)) xline(5.5 13.5 16.5 23.5, lwidth(thin) lcolor(black)) *graph export "C:\Users\bcschipper\burkhard\research\hormone\phases_nfolgain.png", replace *restore *parmby "reg riskloss age asian otherethn gpa math sciothersci econ socsci humanities nmenst_fol nfolic_fol nperi_fol nlut_fol npremens_fol pill if period == 30 & preg == 0 & rlconsistent == 1 & lirrational == 0, vce(robust)", saving(junk, replace) *preserve *use junk, clear *drop if parm=="age" | parm=="asian" | parm=="otherethn" | parm=="gpa" | parm=="math" | parm=="sciothersci" | parm=="econ" | parm=="socsci" | parm=="humanities" | parm=="pill" | parm=="_cons" *gen days = . *la var days "Menstrual Phase" *replace days = 2.75 if parm== "nmenst_fol" *replace days = 9.5 if parm== "nfolic_fol" *replace days = 15 if parm== "nperi_fol" *replace days = 20 if parm== "nlut_fol" *replace days = 25.75 if parm== "npremens_fol" *twoway (scatter estimate days) (rcap min95 max95 days, sort), legend(on order(1 "Female/Male Difference" 2 "95% Confidence Interval")) xtitle(Menstrual Phase ,margin(small)) xscale(range(0 28)) xlabel(2.75 "Menstrual " 9.5 "Folicular" 15 "Peri-Ovulatory" 20 "Luteal" 25.75 "Pre-Menstrual", labsize(small)) xline(5.5 13.5 16.5 23.5, lwidth(thin) lcolor(black)) *graph export "C:\Users\bcschipper\burkhard\research\hormone\phases_nfolloss.png", replace *restore * accuracy of self-reports scatter progf nmens_adj if female == 1 & period == 30 & preg == 0 & pill == 0 & mens < 50 || qfit progf nmens_adj if female == 1 & period == 30 & preg == 0 & pill == 0 & mens < 50 , by(calendar) * Robustness Checks ************************* * with smaller subset of demographic variables reg riskgain age asian otherethn nmenst_fol nfolic_fol nperi_fol nlut_fol npremens_fol pill if period == 30 & preg == 0 & rgconsistent == 1 & girrational == 0, vce(robust) reg riskloss age asian otherethn nmenst_fol nfolic_fol nperi_fol nlut_fol npremens_fol pill if period == 30 & preg == 0 & rlconsistent == 1 & lirrational == 0, vce(robust) logit reflection age asian otherethn nmenst_fol nfolic_fol nperi_fol nlut_fol npremens_fol pill if period == 30 & preg == 0 & consistent == 1 & irrational == 0, vce(robust) logit consistent age asian otherethn nmenst_fol nfolic_fol nperi_fol nlut_fol npremens_fol pill if period == 30 & preg == 0, vce(robust) * with unadjusted menstrual cycle phases reg riskgain age asian otherethn gpa math sciothersci econ socsci humanities nmenst nfolic nperi nlut npremens pill if period == 30 & preg == 0 & rgconsistent == 1 & girrational == 0, vce(robust) reg riskloss age asian otherethn gpa math sciothersci econ socsci humanities nmenst nfolic nperi nlut npremens pill if period == 30 & preg == 0 & rlconsistent == 1 & lirrational == 0, vce(robust) logit reflection age asian otherethn gpa math sciothersci econ socsci humanities nmenst nfolic nperi nlut npremens pill if period == 30 & preg == 0 & consistent == 1 & irrational == 0, vce(robust) logit consistent age asian otherethn gpa math sciothersci econ socsci humanities nmenst nfolic nperi nlut npremens pill if period == 30 & preg == 0, vce(robust) * with uniformly adjusted menstrual cycle phases reg riskgain age asian otherethn gpa math sciothersci econ socsci humanities nmenst_adj nfolic_adj nperi_adj nlut_adj npremens_adj pill if period == 30 & preg == 0 & rgconsistent == 1 & girrational == 0, vce(robust) reg riskloss age asian otherethn gpa math sciothersci econ socsci humanities nmenst_adj nfolic_adj nperi_adj nlut_adj npremens_adj pill if period == 30 & preg == 0 & rlconsistent == 1 & lirrational == 0, vce(robust) logit reflection age asian otherethn gpa math sciothersci econ socsci humanities nmenst_adj nfolic_adj nperi_adj nlut_adj npremens_adj pill if period == 30 & preg == 0 & consistent == 1 & irrational == 0, vce(robust) logit consistent age asian otherethn gpa math sciothersci econ socsci humanities nmenst_adj nfolic_adj nperi_adj nlut_adj npremens_adj pill if period == 30 & preg == 0, vce(robust) * ordered logit ologit riskgain age asian otherethn gpa math sciothersci econ socsci humanities nmenst_fol nfolic_fol nperi_fol nlut_fol npremens_fol pill if period == 30 & preg == 0 & rgconsistent == 1 & girrational == 0, vce(robust) ologit riskloss age asian otherethn gpa math sciothersci econ socsci humanities nmenst_fol nfolic_fol nperi_fol nlut_fol npremens_fol pill if period == 30 & preg == 0 & rlconsistent == 1 & lirrational == 0, vce(robust) * ordered probit oprobit riskgain age asian otherethn gpa math sciothersci econ socsci humanities nmenst_fol nfolic_fol nperi_fol nlut_fol npremens_fol pill if period == 30 & preg == 0 & rgconsistent == 1 & girrational == 0, vce(robust) oprobit riskloss age asian otherethn gpa math sciothersci econ socsci humanities nmenst_fol nfolic_fol nperi_fol nlut_fol npremens_fol pill if period == 30 & preg == 0 & rlconsistent == 1 & lirrational == 0, vce(robust) * probit probit reflection age asian otherethn gpa math sciothersci econ socsci humanities nmenst_fol nfolic_fol nperi_fol nlut_fol npremens_fol pill if period == 30 & preg == 0 & consistent == 1 & irrational == 0, vce(robust) probit consistent age asian otherethn gpa math sciothersci econ socsci humanities nmenst_fol nfolic_fol nperi_fol nlut_fol npremens_fol pill if period == 30 & preg == 0, vce(robust) * with demographic variables (dr bmi siblings workout livalone dating homobi smoking) reg riskgain age dr bmi siblings workout livalone dating homobi smoking asian otherethn gpa math sciothersci econ socsci humanities nmenst_fol nfolic_fol nperi_fol nlut_fol npremens_fol pill if period == 30 & preg == 0 & rgconsistent == 1 & girrational == 0, vce(robust) reg riskloss age dr bmi siblings workout livalone dating homobi smoking asian otherethn gpa math sciothersci econ socsci humanities nmenst_fol nfolic_fol nperi_fol nlut_fol npremens_fol pill if period == 30 & preg == 0 & rlconsistent == 1 & lirrational == 0, vce(robust) logit reflection age dr bmi siblings workout livalone dating homobi smoking asian otherethn gpa math sciothersci econ socsci humanities nmenst_fol nfolic_fol nperi_fol nlut_fol npremens_fol pill if period == 30 & preg == 0 & consistent == 1 & irrational == 0, vce(robust) logit consistent age dr bmi siblings workout livalone dating homobi smoking asian otherethn gpa math sciothersci econ socsci humanities nmenst_fol nfolic_fol nperi_fol nlut_fol npremens_fol pill if period == 30 & preg == 0, vce(robust) * with session dummies reg riskgain sessiondum* age asian otherethn gpa math sciothersci econ socsci humanities nmenst_fol nfolic_fol nperi_fol nlut_fol npremens_fol pill if period == 30 & preg == 0 & rgconsistent == 1 & girrational == 0, vce(robust) reg riskloss sessiondum* age asian otherethn gpa math sciothersci econ socsci humanities nmenst_fol nfolic_fol nperi_fol nlut_fol npremens_fol pill if period == 30 & preg == 0 & rlconsistent == 1 & lirrational == 0, vce(robust) logit reflection sessiondum* age asian otherethn gpa math sciothersci econ socsci humanities nmenst_fol nfolic_fol nperi_fol nlut_fol npremens_fol pill if period == 30 & preg == 0 & consistent == 1 & irrational == 0, vce(robust) logit consistent sessiondum* age asian otherethn gpa math sciothersci econ socsci humanities nmenst_fol nfolic_fol nperi_fol nlut_fol npremens_fol pill if period == 30 & preg == 0, vce(robust) * session fixed effects xtset, clear xtset datum xtreg riskgain age asian otherethn gpa math sciothersci econ socsci humanities nmenst_fol nfolic_fol nperi_fol nlut_fol npremens_fol pill if period == 30 & preg == 0 & rgconsistent == 1 & girrational == 0, fe xtreg riskloss age asian otherethn gpa math sciothersci econ socsci humanities nmenst_fol nfolic_fol nperi_fol nlut_fol npremens_fol pill if period == 30 & preg == 0 & rlconsistent == 1 & lirrational == 0, fe xtlogit reflection age asian otherethn gpa math sciothersci econ socsci humanities nmenst_fol nfolic_fol nperi_fol nlut_fol npremens_fol pill if period == 30 & preg == 0 & consistent == 1 & irrational == 0, fe xtlogit consistent age asian otherethn gpa math sciothersci econ socsci humanities nmenst_fol nfolic_fol nperi_fol nlut_fol npremens_fol pill if period == 30 & preg == 0, fe *************************************** * Analysis Risk & Digit Ratio *************************************** sum dr if new == 1 & period == 30 & female == 1 & white == 1 sum dr if new == 1 & period == 30 & female == 1 & asian == 1 sum dr if new == 1 & period == 30 & female == 1 & otherethn == 1 sum dr if new == 1 & period == 30 & female == 1 sum dr if new == 1 & period == 30 & female == 0 & white == 1 sum dr if new == 1 & period == 30 & female == 0 & asian == 1 sum dr if new == 1 & period == 30 & female == 0 & otherethn == 1 sum dr if new == 1 & period == 30 & female == 0 scatter dr riskgain if period == 30 || lfit dr riskgain if period == 30, by(gender1) scatter dr riskloss if period == 30 || lfit dr riskloss if period == 30, by(gender1) scatter dr riskgain if period == 30 & asian == 1 || lfit dr riskgain if period == 30 & asian == 1, by(gender1) scatter dr riskloss if period == 30 & asian == 1 || lfit dr riskloss if period == 30 & asian == 1, by(gender1) scatter dr riskgain if period == 30 & white == 1 || lfit dr riskgain if period == 30 & white == 1, by(gender1) scatter dr riskloss if period == 30 & white == 1 || lfit dr riskloss if period == 30 & white == 1, by(gender1) * for paper * risk aversion eststo gdr1: reg riskgain age asian otherethn female dr if period == 30 & rgconsistent == 1 & girrational == 0, vce(robust) eststo gdr2: reg riskgain age asian otherethn dr if period == 30 & rgconsistent == 1 & girrational == 0 & female == 1, vce(robust) eststo gdr3: reg riskgain age asian otherethn dr if period == 30 & rgconsistent == 1 & girrational == 0 & female == 0, vce(robust) eststo ldr1: reg riskloss age asian otherethn female dr if period == 30 & rlconsistent == 1 & lirrational == 0, vce(robust) eststo ldr2: reg riskloss age asian otherethn dr if period == 30 & rlconsistent == 1 & lirrational == 0 & female == 1, vce(robust) eststo ldr3: reg riskloss age asian otherethn dr if period == 30 & rlconsistent == 1 & lirrational == 0 & female == 0, vce(robust) estout gdr1 gdr2 gdr3 ldr1 ldr2 ldr3 using "C:\bcschipper\burkhard\research\hormone\raversedr.tex", replace /// cells(b(star fmt(%9.4f)) se(par fmt(%9.4f))) stats(N r2, fmt(%9.0g %9.4f) labels("\textit{Number of Observations}" "R\textsuperscript{2}") layout("\multicolumn{1}{c}{\quad@}")) label /// starlevels(\text{*} 0.05 \text{**} 0.01) keep(asian otherethn female dr) /// mlabels((GDR) (GDRF) (GDRM) (LDR) (LDRF) (LDRM), prefix(\multicolumn{1}{c}{) suffix(})) order(asian otherethn female dr) /// style(tex) collabels(,none) prehead("\begin{tabular}{l*{@span}{d}}\toprule") /// posthead(\midrule) prefoot(\bottomrule) postfoot("\addlinespace[.1cm]\multicolumn{7}{l}{Robust standard errors in parentheses; Significance levels: * 5\%; ** 1\%.} \\" /// "\multicolumn{7}{l}{We suppress from the report the coefficient for age.} \\" "\end{tabular}") notype * reflection, consistency eststo rdr1: logit reflection age asian otherethn female dr if period == 30 & consistent == 1 & irrational == 0, vce(robust) eststo rdr2: logit reflection age asian otherethn dr if period == 30 & consistent == 1 & irrational == 0 & female == 1, vce(robust) eststo rdr3: logit reflection age asian otherethn dr if period == 30 & consistent == 1 & irrational == 0 & female == 0, vce(robust) eststo cdr1: logit consistent age asian otherethn female dr if period == 30, vce(robust) eststo cdr2: logit consistent age asian otherethn dr if period == 30 & female == 1, vce(robust) eststo cdr3: logit consistent age asian otherethn dr if period == 30 & female == 0, vce(robust) estout rdr1 rdr2 rdr3 cdr1 cdr2 cdr3 using "C:\bcschipper\burkhard\research\hormone\refcondr.tex", replace /// cells(b(star fmt(%9.4f)) se(par fmt(%9.4f))) stats(N r2, fmt(%9.0g %9.4f) labels("\textit{Number of Observations}" "R\textsuperscript{2}") layout("\multicolumn{1}{c}{\quad@}")) label /// starlevels(\text{*} 0.05 \text{**} 0.01) keep(asian otherethn female dr) /// mlabels((RDR) (RDRF) (RDRM) (ADR) (ADRF) (ADRM), prefix(\multicolumn{1}{c}{) suffix(})) order(asian otherethn female dr) /// style(tex) collabels(,none) prehead("\begin{tabular}{l*{@span}{d}}\toprule") /// posthead(\midrule) prefoot(\bottomrule) postfoot("\addlinespace[.1cm]\multicolumn{7}{l}{Robust standard errors in parentheses; Significance levels: * 5\%; ** 1\%.} \\" /// "\multicolumn{7}{l}{We suppress from the report the coefficient for age.} \\" "\end{tabular}") notype * whites risk aversion eststo gdr4: reg riskgain age female dr if period == 30 & rgconsistent == 1 & girrational == 0 & white == 1, vce(robust) eststo gdr5: reg riskgain age dr if period == 30 & rgconsistent == 1 & girrational == 0 & white == 1 & female == 1, vce(robust) eststo gdr6: reg riskgain age dr if period == 30 & rgconsistent == 1 & girrational == 0 & white == 1 & female == 0, vce(robust) eststo ldr4: reg riskloss age female dr if period == 30 & rlconsistent == 1 & lirrational == 0 & white == 1, vce(robust) eststo ldr5: reg riskloss age dr if period == 30 & rlconsistent == 1 & lirrational == 0 & white == 1 & female == 1, vce(robust) eststo ldr6: reg riskloss age dr if period == 30 & rlconsistent == 1 & lirrational == 0 & white == 1 & female == 0, vce(robust) estout gdr4 gdr5 gdr6 ldr4 ldr5 ldr6 using "C:\bcschipper\burkhard\research\hormone\raversewdr.tex", replace /// cells(b(star fmt(%9.4f)) se(par fmt(%9.4f))) stats(N r2, fmt(%9.0g %9.4f) labels("\textit{Number of Observations}" "R\textsuperscript{2}") layout("\multicolumn{1}{c}{\quad@}")) label /// starlevels(\text{*} 0.05 \text{**} 0.01) keep(female dr) /// mlabels((GDRW) (GDRWF) (GDRWM) (LDRW) (LDRWF) (LDRWM), prefix(\multicolumn{1}{c}{) suffix(})) order(female dr) /// style(tex) collabels(,none) prehead("\begin{tabular}{l*{@span}{d}}\toprule") /// posthead(\midrule) prefoot(\bottomrule) postfoot("\addlinespace[.1cm]\multicolumn{7}{l}{Robust standard errors in parentheses; Significance levels: * 5\%; ** 1\%.} \\" /// "\multicolumn{7}{l}{We suppress from the report the coefficient for age.} \\" "\end{tabular}") notype * white reflection, consistency eststo rdr4: logit reflection age female dr if period == 30 & consistent == 1 & irrational == 0 & white == 1, vce(robust) eststo rdr5: logit reflection age dr if period == 30 & consistent == 1 & irrational == 0 & white == 1 & female == 1, vce(robust) eststo rdr6: logit reflection age dr if period == 30 & consistent == 1 & irrational == 0 & white == 1 & female == 0, vce(robust) eststo cdr4: logit consistent age female dr if period == 30 & white == 1, vce(robust) eststo cdr5: logit consistent age dr if period == 30 & white == 1 & female == 1, vce(robust) eststo cdr6: logit consistent age dr if period == 30 & white == 1 & female == 0, vce(robust) estout rdr4 rdr5 rdr6 cdr4 cdr5 cdr6 using "C:\bcschipper\burkhard\research\hormone\refconwdr.tex", replace /// cells(b(star fmt(%9.4f)) se(par fmt(%9.4f))) stats(N r2, fmt(%9.0g %9.4f) labels("\textit{Number of Observations}" "R\textsuperscript{2}") layout("\multicolumn{1}{c}{\quad@}")) label /// starlevels(\text{*} 0.05 \text{**} 0.01) keep(female dr) /// mlabels((RDRW) (RDRWF) (RDRWM) (ADRW) (ADRWF) (ADRWM), prefix(\multicolumn{1}{c}{) suffix(})) order(female dr) /// style(tex) collabels(,none) prehead("\begin{tabular}{l*{@span}{d}}\toprule") /// posthead(\midrule) prefoot(\bottomrule) postfoot("\addlinespace[.1cm]\multicolumn{7}{l}{Robust standard errors in parentheses; Significance levels: * 5\%; ** 1\%.} \\" /// "\multicolumn{7}{l}{We suppress from the report the coefficient for age.} \\" "\end{tabular}") notype * asian risk aversion eststo gdr7: reg riskgain age female dr if period == 30 & rgconsistent == 1 & girrational == 0 & asian == 1, vce(robust) eststo gdr8: reg riskgain age dr if period == 30 & rgconsistent == 1 & girrational == 0 & asian == 1 & female == 1, vce(robust) eststo gdr9: reg riskgain age dr if period == 30 & rgconsistent == 1 & girrational == 0 & asian == 1 & female == 0, vce(robust) eststo ldr7: reg riskloss age female dr if period == 30 & rlconsistent == 1 & lirrational == 0 & asian == 1, vce(robust) eststo ldr8: reg riskloss age dr if period == 30 & rlconsistent == 1 & lirrational == 0 & asian == 1 & female == 1, vce(robust) eststo ldr9: reg riskloss age dr if period == 30 & rlconsistent == 1 & lirrational == 0 & asian == 1 & female == 0, vce(robust) estout gdr7 gdr8 gdr9 ldr7 ldr8 ldr9 using "C:\bcschipper\burkhard\research\hormone\raverseadr.tex", replace /// cells(b(star fmt(%9.4f)) se(par fmt(%9.4f))) stats(N r2, fmt(%9.0g %9.4f) labels("\textit{Number of Observations}" "R\textsuperscript{2}") layout("\multicolumn{1}{c}{\quad@}")) label /// starlevels(\text{*} 0.05 \text{**} 0.01) keep(female dr) /// mlabels((GDRA) (GDRAF) (GDRAM) (LDRA) (LDRAF) (LDRAM), prefix(\multicolumn{1}{c}{) suffix(})) order(female dr) /// style(tex) collabels(,none) prehead("\begin{tabular}{l*{@span}{d}}\toprule") /// posthead(\midrule) prefoot(\bottomrule) postfoot("\addlinespace[.1cm]\multicolumn{7}{l}{Robust standard errors in parentheses; Significance levels: * 5\%; ** 1\%.} \\" /// "\multicolumn{7}{l}{We suppress from the report the coefficient for age.} \\" "\end{tabular}") notype * white reflection, consistency eststo rdr7: logit reflection age female dr if period == 30 & consistent == 1 & irrational == 0 & asian == 1, vce(robust) eststo rdr8: logit reflection age dr if period == 30 & consistent == 1 & irrational == 0 & asian == 1 & female == 1, vce(robust) eststo rdr9: logit reflection age dr if period == 30 & consistent == 1 & irrational == 0 & asian == 1 & female == 0, vce(robust) eststo cdr7: logit consistent age female dr if period == 30 & asian == 1, vce(robust) eststo cdr8: logit consistent age dr if period == 30 & asian == 1 & female == 1, vce(robust) eststo cdr9: logit consistent age dr if period == 30 & asian == 1 & female == 0, vce(robust) estout rdr7 rdr8 rdr9 cdr7 cdr8 cdr9 using "C:\bcschipper\burkhard\research\hormone\refconadr.tex", replace /// cells(b(star fmt(%9.4f)) se(par fmt(%9.4f))) stats(N r2, fmt(%9.0g %9.4f) labels("\textit{Number of Observations}" "R\textsuperscript{2}") layout("\multicolumn{1}{c}{\quad@}")) label /// starlevels(\text{*} 0.05 \text{**} 0.01) keep(female dr) /// mlabels((RDRA) (RDRAF) (RDRAM) (ADRA) (ADRAF) (ADRAM), prefix(\multicolumn{1}{c}{) suffix(})) order(female dr) /// style(tex) collabels(,none) prehead("\begin{tabular}{l*{@span}{d}}\toprule") /// posthead(\midrule) prefoot(\bottomrule) postfoot("\addlinespace[.1cm]\multicolumn{7}{l}{Robust standard errors in parentheses; Significance levels: * 5\%; ** 1\%.} \\" /// "\multicolumn{7}{l}{We suppress from the report the coefficient for age.} \\" "\end{tabular}") notype * Robustness Checks ********************* * ordered logit ologit riskgain age asian otherethn female dr if period == 30 & rgconsistent == 1 & girrational == 0, vce(robust) ologit riskgain age asian otherethn dr if period == 30 & rgconsistent == 1 & girrational == 0 & female == 1, vce(robust) ologit riskgain age asian otherethn dr if period == 30 & rgconsistent == 1 & girrational == 0 & female == 0, vce(robust) ologit riskloss age asian otherethn female dr if period == 30 & rlconsistent == 1 & lirrational == 0, vce(robust) ologit riskloss age asian otherethn dr if period == 30 & rlconsistent == 1 & lirrational == 0 & female == 1, vce(robust) ologit riskloss age asian otherethn dr if period == 30 & rlconsistent == 1 & lirrational == 0 & female == 0, vce(robust) * ordered probit oprobit riskgain age asian otherethn dr if period == 30 & rgconsistent == 1 & girrational == 0, vce(robust) oprobit riskgain age asian otherethn dr if period == 30 & rgconsistent == 1 & girrational == 0 & female == 1, vce(robust) oprobit riskgain age asian otherethn dr if period == 30 & rgconsistent == 1 & girrational == 0 & female == 0, vce(robust) oprobit riskloss age asian otherethn dr if period == 30 & rlconsistent == 1 & lirrational == 0, vce(robust) oprobit riskloss age asian otherethn dr if period == 30 & rlconsistent == 1 & lirrational == 0 & female == 1, vce(robust) oprobit riskloss age asian otherethn dr if period == 30 & rlconsistent == 1 & lirrational == 0 & female == 0, vce(robust) * reflection probit probit reflection age asian otherethn female dr if period == 30 & consistent == 1 & irrational == 0, vce(robust) probit reflection age asian otherethn dr if period == 30 & consistent == 1 & irrational == 0 & female == 1, vce(robust) probit reflection age asian otherethn dr if period == 30 & consistent == 1 & irrational == 0 & female == 0, vce(robust) * consistency probit probit consistent age asian otherethn female dr if period == 30, vce(robust) probit consistent age asian otherethn dr if period == 30 & female == 1, vce(robust) probit consistent age asian otherethn dr if period == 30 & female == 0, vce(robust) * allow inconsistent choices reg riskgain age asian otherethn female dr if period == 30, vce(robust) reg riskgain age asian otherethn dr if period == 30 & female == 1, vce(robust) reg riskgain age asian otherethn dr if period == 30 & female == 0, vce(robust) reg riskloss age asian otherethn female dr if period == 30, vce(robust) reg riskloss age asian otherethn dr if period == 30 & female == 1, vce(robust) reg riskloss age asian otherethn dr if period == 30 & female == 0, vce(robust) * drop controls reg riskgain female dr if period == 30 & rgconsistent == 1 & girrational == 0, vce(robust) reg riskgain dr if period == 30 & rgconsistent == 1 & girrational == 0 & female == 1, vce(robust) reg riskgain dr if period == 30 & rgconsistent == 1 & girrational == 0 & female == 0, vce(robust) reg riskloss female dr if period == 30 & rlconsistent == 1 & lirrational == 0, vce(robust) reg riskloss dr if period == 30 & rlconsistent == 1 & lirrational == 0 & female == 1, vce(robust) reg riskloss dr if period == 30 & rlconsistent == 1 & lirrational == 0 & female == 0, vce(robust) logit reflection female dr if period == 30 & consistent == 1 & irrational == 0, vce(robust) logit reflection dr if period == 30 & consistent == 1 & irrational == 0 & female == 1, vce(robust) logit reflection dr if period == 30 & consistent == 1 & irrational == 0 & female == 0, vce(robust) logit consistent female dr if period == 30, vce(robust) logit consistent dr if period == 30 & female == 1, vce(robust) logit consistent dr if period == 30 & female == 0, vce(robust) * add controls reg riskgain age asian otherethn gpa math sciothersci econ socsci humanities female dr if period == 30 & rgconsistent == 1 & girrational == 0, vce(robust) reg riskgain age asian otherethn gpa math sciothersci econ socsci humanities dr if period == 30 & rgconsistent == 1 & girrational == 0 & female == 1, vce(robust) reg riskgain age asian otherethn gpa math sciothersci econ socsci humanities dr if period == 30 & rgconsistent == 1 & girrational == 0 & female == 0, vce(robust) reg riskloss age asian otherethn gpa math sciothersci econ socsci humanities female dr if period == 30 & rlconsistent == 1 & lirrational == 0, vce(robust) reg riskloss age asian otherethn gpa math sciothersci econ socsci humanities dr if period == 30 & rlconsistent == 1 & lirrational == 0 & female == 1, vce(robust) reg riskloss age asian otherethn gpa math sciothersci econ socsci humanities dr if period == 30 & rlconsistent == 1 & lirrational == 0 & female == 0, vce(robust) logit reflection age asian otherethn gpa math sciothersci econ socsci humanities female dr if period == 30 & consistent == 1 & irrational == 0, vce(robust) logit reflection age asian otherethn gpa math sciothersci econ socsci humanities dr if period == 30 & consistent == 1 & irrational == 0 & female == 1, vce(robust) logit reflection age asian otherethn gpa math sciothersci econ socsci humanities dr if period == 30 & consistent == 1 & irrational == 0 & female == 0, vce(robust) logit consistent age asian otherethn gpa math sciothersci econ socsci humanities female dr if period == 30, vce(robust) logit consistent age asian otherethn gpa math sciothersci econ socsci humanities dr if period == 30 & female == 1, vce(robust) logit consistent age asian otherethn gpa math sciothersci econ socsci humanities dr if period == 30 & female == 0, vce(robust) * add even more controls reg riskgain age bmi siblings workout livalone dating homobi smoking asian otherethn gpa math sciothersci econ socsci humanities female dr if period == 30 & rgconsistent == 1 & girrational == 0, vce(robust) reg riskgain age bmi siblings workout livalone dating homobi smoking asian otherethn gpa math sciothersci econ socsci humanities dr if period == 30 & rgconsistent == 1 & girrational == 0 & female == 1, vce(robust) reg riskgain age bmi siblings workout livalone dating homobi smoking asian otherethn gpa math sciothersci econ socsci humanities dr if period == 30 & rgconsistent == 1 & girrational == 0 & female == 0, vce(robust) reg riskloss age bmi siblings workout livalone dating homobi smoking asian otherethn gpa math sciothersci econ socsci humanities female dr if period == 30 & rlconsistent == 1 & lirrational == 0, vce(robust) reg riskloss age bmi siblings workout livalone dating homobi smoking asian otherethn gpa math sciothersci econ socsci humanities dr if period == 30 & rlconsistent == 1 & lirrational == 0 & female == 1, vce(robust) reg riskloss age bmi siblings workout livalone dating homobi smoking asian otherethn gpa math sciothersci econ socsci humanities dr if period == 30 & rlconsistent == 1 & lirrational == 0 & female == 0, vce(robust) logit reflection age bmi siblings workout livalone dating homobi smoking asian otherethn gpa math sciothersci econ socsci humanities female dr if period == 30 & consistent == 1 & irrational == 0, vce(robust) logit reflection age bmi siblings workout livalone dating homobi smoking asian otherethn gpa math sciothersci econ socsci humanities dr if period == 30 & consistent == 1 & irrational == 0 & female == 1, vce(robust) logit reflection age bmi siblings workout livalone dating homobi smoking asian otherethn gpa math sciothersci econ socsci humanities dr if period == 30 & consistent == 1 & irrational == 0 & female == 0, vce(robust) logit consistent age bmi siblings workout livalone dating homobi smoking asian otherethn gpa math sciothersci econ socsci humanities female dr if period == 30, vce(robust) logit consistent age bmi siblings workout livalone dating homobi smoking asian otherethn gpa math sciothersci econ socsci humanities dr if period == 30 & female == 1, vce(robust) logit consistent age bmi siblings workout livalone dating homobi smoking asian otherethn gpa math sciothersci econ socsci humanities dr if period == 30 & female == 0, vce(robust) * white subpopulation reg riskgain age dr if period == 30 & rgconsistent == 1 & girrational == 0 & female == 1 & white == 1, vce(robust) reg riskgain age dr if period == 30 & rgconsistent == 1 & girrational == 0 & female == 0 & white == 1, vce(robust) reg riskloss age dr if period == 30 & rlconsistent == 1 & lirrational == 0 & female == 1 & white == 1, vce(robust) reg riskloss age dr if period == 30 & rlconsistent == 1 & lirrational == 0 & female == 0 & white == 1, vce(robust) logit reflection age dr if period == 30 & consistent == 1 & irrational == 0 & female == 1 & white == 1, vce(robust) logit reflection age dr if period == 30 & consistent == 1 & irrational == 0 & female == 0 & white == 1, vce(robust) logit consistent age dr if period == 30 & female == 1 & white == 1, vce(robust) logit consistent age dr if period == 30 & female == 0 & white == 1, vce(robust) * asian subpopulation reg riskgain age dr if period == 30 & rgconsistent == 1 & girrational == 0 & female == 1 & asian == 1, vce(robust) reg riskgain age dr if period == 30 & rgconsistent == 1 & girrational == 0 & female == 0 & asian == 1, vce(robust) reg riskloss age dr if period == 30 & rlconsistent == 1 & lirrational == 0 & female == 1 & asian == 1, vce(robust) reg riskloss age dr if period == 30 & rlconsistent == 1 & lirrational == 0 & female == 0 & asian == 1, vce(robust) logit reflection age dr if period == 30 & consistent == 1 & irrational == 0 & female == 1 & asian == 1, vce(robust) logit reflection age dr if period == 30 & consistent == 1 & irrational == 0 & female == 0 & asian == 1, vce(robust) logit consistent age dr if period == 30 & female == 1 & asian == 1, vce(robust) logit consistent age dr if period == 30 & female == 0 & asian == 1, vce(robust) * exclude subjects with hands that we found difficult to measure and broken2d4d reg riskgain age asian otherethn female dr if period == 30 & rgconsistent == 1 & girrational == 0 & uncleardr == 0 & broken2d4d == 0, vce(robust) reg riskgain age asian otherethn dr if period == 30 & rgconsistent == 1 & girrational == 0 & female == 1 & uncleardr == 0 & broken2d4d == 0, vce(robust) reg riskgain age asian otherethn dr if period == 30 & rgconsistent == 1 & girrational == 0 & female == 0 & uncleardr == 0 & broken2d4d == 0, vce(robust) reg riskloss age asian otherethn female dr if period == 30 & rlconsistent == 1 & lirrational == 0 & uncleardr == 0 & broken2d4d == 0, vce(robust) reg riskloss age asian otherethn dr if period == 30 & rlconsistent == 1 & lirrational == 0 & female == 1 & uncleardr == 0 & broken2d4d == 0, vce(robust) reg riskloss age asian otherethn dr if period == 30 & rlconsistent == 1 & lirrational == 0 & female == 0 & uncleardr == 0 & broken2d4d == 0, vce(robust) logit reflection age asian otherethn female dr if period == 30 & consistent == 1 & irrational == 0 & uncleardr == 0 & broken2d4d == 0, vce(robust) logit reflection age asian otherethn dr if period == 30 & consistent == 1 & irrational == 0 & female == 1 & uncleardr == 0 & broken2d4d == 0, vce(robust) logit reflection age asian otherethn dr if period == 30 & consistent == 1 & irrational == 0 & female == 0 & uncleardr == 0 & broken2d4d == 0, vce(robust) logit consistent age asian otherethn female dr if period == 30 & uncleardr == 0 & broken2d4d == 0, vce(robust) logit consistent age asian otherethn dr if period == 30 & female == 1 & uncleardr == 0 & broken2d4d == 0, vce(robust) logit consistent age asian otherethn dr if period == 30 & female == 0 & uncleardr == 0 & broken2d4d == 0, vce(robust) reg riskgain age dr if period == 30 & rgconsistent == 1 & girrational == 0 & female == 1 & white == 1 & uncleardr == 0 & broken2d4d == 0, vce(robust) reg riskgain age dr if period == 30 & rgconsistent == 1 & girrational == 0 & female == 0 & white == 1 & uncleardr == 0 & broken2d4d == 0, vce(robust) reg riskloss age dr if period == 30 & rlconsistent == 1 & lirrational == 0 & female == 1 & white == 1 & uncleardr == 0 & broken2d4d == 0, vce(robust) reg riskloss age dr if period == 30 & rlconsistent == 1 & lirrational == 0 & female == 0 & white == 1 & uncleardr == 0 & broken2d4d == 0, vce(robust) logit reflection age dr if period == 30 & consistent == 1 & irrational == 0 & female == 1 & white == 1 & uncleardr == 0 & broken2d4d == 0, vce(robust) logit reflection age dr if period == 30 & consistent == 1 & irrational == 0 & female == 0 & white == 1 & uncleardr == 0 & broken2d4d == 0, vce(robust) logit consistent age dr if period == 30 & female == 1 & white == 1 & uncleardr == 0 & broken2d4d == 0, vce(robust) logit consistent age dr if period == 30 & female == 0 & white == 1 & uncleardr == 0 & broken2d4d == 0, vce(robust) reg riskgain age dr if period == 30 & rgconsistent == 1 & girrational == 0 & female == 1 & asian == 1 & uncleardr == 0 & broken2d4d == 0, vce(robust) reg riskgain age dr if period == 30 & rgconsistent == 1 & girrational == 0 & female == 0 & asian == 1 & uncleardr == 0 & broken2d4d == 0, vce(robust) reg riskloss age dr if period == 30 & rlconsistent == 1 & lirrational == 0 & female == 1 & asian == 1 & uncleardr == 0 & broken2d4d == 0, vce(robust) reg riskloss age dr if period == 30 & rlconsistent == 1 & lirrational == 0 & female == 0 & asian == 1 & uncleardr == 0 & broken2d4d == 0, vce(robust) logit reflection age dr if period == 30 & consistent == 1 & irrational == 0 & female == 1 & asian == 1 & uncleardr == 0 & broken2d4d == 0, vce(robust) logit reflection age dr if period == 30 & consistent == 1 & irrational == 0 & female == 0 & asian == 1 & uncleardr == 0 & broken2d4d == 0, vce(robust) logit consistent age dr if period == 30 & female == 1 & asian == 1 & uncleardr == 0 & broken2d4d == 0, vce(robust) logit consistent age dr if period == 30 & female == 0 & asian == 1 & uncleardr == 0 & broken2d4d == 0, vce(robust) * drop left hand scan * session dummies * session fixed effects * exclude outlier for asian female sum dr if period == 30 & new == 1 & asian == 1 & female == 1 reg riskloss age dr if period == 30 & rlconsistent == 1 & lirrational == 0 & asian == 1 & female == 1 & dr > 0.881, vce(robust) *log using "C:\Users\bcschipper\burkhard\research\hormone\stats\aldo20110926.smcl", replace * Branas-Garzas and Rustichini type analysis * riskgain is the number of "safe" choices in Holt-Laury for gains * riskloss is the number of "safe" choices in Holt-Laury for losses * dr refers to 2D:4D * gpa refers to GPA * rgconsistent == 1 means we consider only subjects who choose consistently in the gain domain * girrational == 0 means we consider only subjects who did not choose a dominated choice in the gain domain * female == 1 refers to being female * vce(robust) controls standard errors for heteroskedasticity etc. * period == 30 is needed because our auction data is merged with the choice under risk data set * males gains reg riskgain dr if period == 30 & rgconsistent == 1 & girrational == 0 & female == 0, vce(robust) reg gpa dr if period == 30 & rgconsistent == 1 & girrational == 0 & female == 0, vce(robust) reg riskgain dr gpa if period == 30 & rgconsistent == 1 & girrational == 0 & female == 0, vce(robust) * males losses reg riskloss dr if period == 30 & rlconsistent == 1 & lirrational == 0 & female == 0, vce(robust) reg gpa dr if period == 30 & rlconsistent == 1 & lirrational == 0 & female == 0, vce(robust) reg riskloss dr gpa if period == 30 & rlconsistent == 1 & lirrational == 0 & female == 0, vce(robust) * females gains reg riskgain dr if period == 30 & rgconsistent == 1 & girrational == 0 & female == 1, vce(robust) reg gpa dr if period == 30 & rgconsistent == 1 & girrational == 0 & female == 1, vce(robust) reg riskgain dr gpa if period == 30 & rgconsistent == 1 & girrational == 0 & female == 1, vce(robust) * females losses reg riskloss dr if period == 30 & rlconsistent == 1 & lirrational == 0 & female == 1, vce(robust) reg gpa dr if period == 30 & rlconsistent == 1 & lirrational == 0 & female == 1, vce(robust) reg riskloss dr gpa if period == 30 & rlconsistent == 1 & lirrational == 0 & female == 1, vce(robust) * whites * males gains reg riskgain dr if period == 30 & rgconsistent == 1 & girrational == 0 & female == 0 & white == 1, vce(robust) reg gpa dr if period == 30 & rgconsistent == 1 & girrational == 0 & female == 0 & white == 1, vce(robust) reg riskgain dr gpa if period == 30 & rgconsistent == 1 & girrational == 0 & female == 0 & white == 1, vce(robust) * males losses reg riskloss dr if period == 30 & rlconsistent == 1 & lirrational == 0 & female == 0 & white == 1, vce(robust) reg gpa dr if period == 30 & rlconsistent == 1 & lirrational == 0 & female == 0 & white == 1, vce(robust) reg riskloss dr gpa if period == 30 & rlconsistent == 1 & lirrational == 0 & female == 0 & white == 1, vce(robust) * females gains reg riskgain dr if period == 30 & rgconsistent == 1 & girrational == 0 & female == 1 & white == 1, vce(robust) reg gpa dr if period == 30 & rgconsistent == 1 & girrational == 0 & female == 1 & white == 1, vce(robust) reg riskgain dr gpa if period == 30 & rgconsistent == 1 & girrational == 0 & female == 1 & white == 1, vce(robust) * females losses reg riskloss dr if period == 30 & rlconsistent == 1 & lirrational == 0 & female == 1 & white == 1, vce(robust) reg gpa dr if period == 30 & rlconsistent == 1 & lirrational == 0 & female == 1 & white == 1, vce(robust) reg riskloss dr gpa if period == 30 & rlconsistent == 1 & lirrational == 0 & female == 1 & white == 1, vce(robust) * asians * males gains reg riskgain dr if period == 30 & rgconsistent == 1 & girrational == 0 & female == 0 & asian == 1, vce(robust) reg gpa dr if period == 30 & rgconsistent == 1 & girrational == 0 & female == 0 & asian == 1, vce(robust) reg riskgain dr gpa if period == 30 & rgconsistent == 1 & girrational == 0 & female == 0 & asian == 1, vce(robust) * males losses reg riskloss dr if period == 30 & rlconsistent == 1 & lirrational == 0 & female == 0 & asian == 1, vce(robust) reg gpa dr if period == 30 & rlconsistent == 1 & lirrational == 0 & female == 0 & asian == 1, vce(robust) reg riskloss dr gpa if period == 30 & rlconsistent == 1 & lirrational == 0 & female == 0 & asian == 1, vce(robust) * females gains reg riskgain dr if period == 30 & rgconsistent == 1 & girrational == 0 & female == 1 & asian == 1, vce(robust) reg gpa dr if period == 30 & rgconsistent == 1 & girrational == 0 & female == 1 & asian == 1, vce(robust) reg riskgain dr gpa if period == 30 & rgconsistent == 1 & girrational == 0 & female == 1 & asian == 1, vce(robust) * females losses reg riskloss dr if period == 30 & rlconsistent == 1 & lirrational == 0 & female == 1 & asian == 1, vce(robust) reg gpa dr if period == 30 & rlconsistent == 1 & lirrational == 0 & female == 1 & asian == 1, vce(robust) reg riskloss dr gpa if period == 30 & rlconsistent == 1 & lirrational == 0 & female == 1 & asian == 1, vce(robust) *log using "C:\bcschipper\burkhard\research\hormone\stats\aldo_alldata_20110927.smcl", replace * Branas-Garzas and Rustichini type analysis * riskgain is the number of "safe" choices in Holt-Laury for gains * riskloss is the number of "safe" choices in Holt-Laury for losses * dr refers to 2D:4D * gpa refers to current Grade Point Average * female == 1 refers to being female * vce(robust) controls standard errors for heteroskedasticity etc. * period == 30 is needed because our auction data is merged with the choice under risk data set *************************** * using all data not just consistent data *************************** * males gains reg riskgain dr if period == 30 & female == 0, vce(robust) reg gpa dr if period == 30 & female == 0, vce(robust) reg riskgain dr gpa if period == 30 & female == 0, vce(robust) * males losses reg riskloss dr if period == 30 & female == 0, vce(robust) reg gpa dr if period == 30 & female == 0, vce(robust) reg riskloss dr gpa if period == 30 & female == 0, vce(robust) * females gains reg riskgain dr if period == 30 & female == 1, vce(robust) reg gpa dr if period == 30 & female == 1, vce(robust) reg riskgain dr gpa if period == 30 & female == 1, vce(robust) * females losses reg riskloss dr if period == 30 & female == 1, vce(robust) reg gpa dr if period == 30 & female == 1, vce(robust) reg riskloss dr gpa if period == 30 & female == 1, vce(robust) * whites * males gains reg riskgain dr if period == 30 & female == 0 & white == 1, vce(robust) reg gpa dr if period == 30 & female == 0 & white == 1, vce(robust) reg riskgain dr gpa if period == 30 & female == 0 & white == 1, vce(robust) * males losses reg riskloss dr if period == 30 & female == 0 & white == 1, vce(robust) reg gpa dr if period == 30 & female == 0 & white == 1, vce(robust) reg riskloss dr gpa if period == 30 & female == 0 & white == 1, vce(robust) * females gains reg riskgain dr if period == 30 & female == 1 & white == 1, vce(robust) reg gpa dr if period == 30 & female == 1 & white == 1, vce(robust) reg riskgain dr gpa if period == 30 & female == 1 & white == 1, vce(robust) * females losses reg riskloss dr if period == 30 & female == 1 & white == 1, vce(robust) reg gpa dr if period == 30 & female == 1 & white == 1, vce(robust) reg riskloss dr gpa if period == 30 & female == 1 & white == 1, vce(robust) * asians * males gains reg riskgain dr if period == 30 & female == 0 & asian == 1, vce(robust) reg gpa dr if period == 30 & female == 0 & asian == 1, vce(robust) reg riskgain dr gpa if period == 30 & female == 0 & asian == 1, vce(robust) * males losses reg riskloss dr if period == 30 & female == 0 & asian == 1, vce(robust) reg gpa dr if period == 30 & female == 0 & asian == 1, vce(robust) reg riskloss dr gpa if period == 30 & female == 0 & asian == 1, vce(robust) * females gains reg riskgain dr if period == 30 & female == 1 & asian == 1, vce(robust) reg gpa dr if period == 30 & female == 1 & asian == 1, vce(robust) reg riskgain dr gpa if period == 30 & female == 1 & asian == 1, vce(robust) * females losses reg riskloss dr if period == 30 & female == 1 & asian == 1, vce(robust) reg gpa dr if period == 30 & female == 1 & asian == 1, vce(robust) reg riskloss dr gpa if period == 30 & female == 1 & asian == 1, vce(robust) *log close ***************************** * For presentations ***************************** eststo gp1: reg riskgain female if period == 30 & preg == 0 & rgconsistent == 1 & girrational == 0, vce(robust) eststo gp2: reg riskgain ntestf nestf nprogf ncortf pill female if period == 30 & preg == 0 & rgconsistent == 1 & girrational == 0, vce(robust) eststo lp1: reg riskloss female if period == 30 & preg == 0 & rlconsistent == 1 & lirrational == 0, vce(robust) eststo lp2: reg riskloss ntestf nestf nprogf ncortf pill female if period == 30 & preg == 0 & rlconsistent == 1 & lirrational == 0, vce(robust) estout gp1 gp2 lp1 lp2 using "C:\bcschipper\burkhard\research\hormone\present_gender.tex", replace /// cells(b(star fmt(%9.4f)) se(par fmt(%9.4f))) stats(N r2, fmt(%9.0g %9.4f) labels("\textit{Number of Observations}" "R\textsuperscript{2}") layout("\multicolumn{1}{c}{\quad@}")) label /// starlevels(\text{*} 0.05 \text{**} 0.01) keep(ntestf nestf nprogf ncortf pill female) /// mlabels((Gains) (Gains) (Losses) (Losses), prefix(\multicolumn{1}{c}{) suffix(})) order(female ntestf nestf nprogf ncortf pill) /// style(tex) collabels(,none) prehead("\begin{tabular}{l*{@span}{d}}\toprule") /// posthead(\midrule) prefoot(\bottomrule) postfoot("\addlinespace[.1cm]\multicolumn{5}{l}{Robust standard errors in parentheses; Significance levels: * 5\%; ** 1\%.} \\" "\end{tabular}") notype *************************************** *************************************** * Revisions February 2014 *************************************** *************************************** ************************************ * Analysis Risk & Gender ************************************ eststo g0c: reg riskgain female age asian otherethn gpa sciothersci math econ socsci humanities if period == 30 & rgconsistent == 1 & girrational == 0, vce(robust) eststo g0: reg riskgain female age asian otherethn gpa sciothersci math econ socsci humanities if period == 30, vce(robust) eststo g1c: reg riskgain female age asian otherethn gpa if period == 30 & rgconsistent == 1 & girrational == 0, vce(robust) eststo l0c: reg riskloss female age asian otherethn gpa sciothersci math econ socsci humanities if period == 30 & rlconsistent == 1 & lirrational == 0, vce(robust) eststo l0: reg riskloss female age asian otherethn gpa sciothersci math econ socsci humanities if period == 30, vce(robust) eststo l1c: reg riskloss female age asian otherethn gpa if period == 30 & rlconsistent == 1 & lirrational == 0, vce(robust) estout g0c g0 g1c l0c l0 l1c using "C:\bcschipper\burkhard\research\hormone\gender1_collapsed.tex", replace /// cells(b(star fmt(%9.4f)) se(par fmt(%9.4f))) stats(N r2, fmt(%9.0g %9.4f) labels("\textit{Number of Observations}" "R\textsuperscript{2}") layout("\multicolumn{1}{c}{\quad@}")) label /// starlevels(\text{*} 0.05 \text{**} 0.01) keep(female) /// mlabels((G0A) (G0) (G1A) (L0A) (L0) (L1A), prefix(\multicolumn{1}{c}{) suffix(})) indicate("Demographics = age asian otherethn gpa" "Majors of Study = math sciothersci econ socsci humanities", labels(\textmd{Yes} \textmd{No})) order(female) /// style(tex) collabels(,none) prehead("\begin{tabular}{l*{@span}{d}}\toprule") /// posthead(\midrule) prefoot(\bottomrule) postfoot("\addlinespace[.1cm]\multicolumn{7}{l}{Robust standard errors in parentheses; Significance levels: * 5\%; ** 1\%.}\\""\end{tabular}") notype * does not update automatically in tex file eststo r0: logit reflection female age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & consistent == 1 & irrational == 0, vce(robust) eststo r1: logit reflection female age asian otherethn gpa if period == 30 & consistent == 1 & irrational == 0, vce(robust) eststo c0: logit consistent female age asian otherethn gpa math sciothersci econ socsci humanities if period == 30 & new == 1, vce(robust) eststo c1: logit consistent female age asian otherethn gpa if period == 30 & new == 1, vce(robust) estout r0 r1 c0 c1 using "C:\bcschipper\burkhard\research\hormone\gender2_collapsed.tex", replace /// cells(b(star fmt(%9.4f)) se(par fmt(%9.4f))) stats(N r2_p, fmt(%9.0g %9.4f) labels("\textit{Number of Observations}" "Pseudo R\textsuperscript{2}") layout("\multicolumn{1}{c}{\quad@}")) label /// starlevels(\text{*} 0.05 \text{**} 0.01) keep(female) /// eqlabels(none) mlabels((R0) (R1) (A0) (A1), prefix(\multicolumn{1}{c}{) suffix(})) indicate("Demographics = age asian otherethn gpa" "Majors of Study = math sciothersci econ socsci humanities", labels(\textmd{Yes} \textmd{No})) order(female) /// style(tex) collabels(,none) prehead("\begin{tabular}{l*{@span}{d}}\toprule") /// posthead(\midrule) prefoot(\bottomrule) postfoot("\addlinespace[.1cm]\multicolumn{5}{l}{Robust standard errors in parentheses; Significance levels: * 5\%; ** 1\%.} \\""\end{tabular}") notype **************************************** * Analysis Risk & Contraceptives **************************************** * risk aversion and contraceptives eststo g2f: reg riskgain pill age asian otherethn gpa if period == 30 & rgconsistent == 1 & girrational == 0 & female == 1, vce(robust) eststo l2f: reg riskloss pill age asian otherethn gpa if period == 30 & rlconsistent == 1 & lirrational == 0 & female == 1, vce(robust) eststo r2f: logit reflection pill age asian otherethn gpa if period == 30 & consistent == 1 & irrational == 0 & female == 1, vce(robust) eststo c2f: logit consistent pill age asian otherethn gpa if period == 30 & new == 1 & female == 1, vce(robust) * does not update automatically in tex file estout g2f l2f r2f c2f using "C:\bcschipper\burkhard\research\hormone\contraceptives_collapsed.tex", replace /// cells(b(star fmt(%9.4f)) se(par fmt(%9.4f))) stats(N r2 r2_p, fmt(%9.0g %9.4f) labels("\textit{Number of Observations}" "R\textsuperscript{2}" "Pseudo R\textsuperscript{2}") layout("\multicolumn{1}{c}{\quad@}")) label /// starlevels(\text{*} 0.05 \text{**} 0.01) keep(pill) /// eqlabels(none) mlabels((G2F) (L2F) (R2F) (A2F), prefix(\multicolumn{1}{c}{) suffix(})) indicate("Demographics = age asian otherethn gpa", labels(\textmd{Yes} \textmd{No})) order(pill) /// style(tex) collabels(,none) prehead("\begin{tabular}{l*{@span}{d}}\toprule") /// posthead(\midrule) prefoot(\bottomrule) postfoot("\addlinespace[.1cm]\multicolumn{5}{l}{Robust standard errors in parentheses; Significance levels: * 5\%; ** 1\%.} \\""\end{tabular}") notype ************************************* * Analysis Risk \& Menstrual Cycle ************************************* eststo g5: reg riskgain age asian otherethn gpa nfolic_fol nperi_fol nlut_fol npremens_fol if female == 1 & pill == 0 & period == 30 & preg == 0 & rgconsistent == 1 & girrational == 0, vce(robust) eststo l5: reg riskloss age asian otherethn gpa nfolic_fol nperi_fol nlut_fol npremens_fol if female == 1 & pill == 0 & period == 30 & preg == 0 & rlconsistent == 1 & lirrational == 0, vce(robust) eststo r5: logit reflection age asian otherethn gpa nfolic_fol nperi_fol nlut_fol npremens_fol if female == 1 & pill == 0 & period == 30 & preg == 0 & consistent == 1 & irrational == 0, vce(robust) eststo c5: logit consistent age asian otherethn gpa nfolic_fol nperi_fol nlut_fol npremens_fol if female == 1 & pill == 0 & period == 30 & preg == 0, vce(robust) estout g5 l5 r5 c5 using "C:\bcschipper\burkhard\research\hormone\menstruation_collapsed.tex", replace /// cells(b(star fmt(%9.4f)) se(par fmt(%9.4f))) stats(N r2 r2_p, fmt(%9.0g %9.4f) labels("\textit{Number of Observations}" "R\textsuperscript{2}" "Pseudo R\textsuperscript{2}") layout("\multicolumn{1}{c}{\quad@}")) label /// starlevels(\text{*} 0.05 \text{**} 0.01) keep(nfolic_fol nperi_fol nlut_fol npremens_fol) /// eqlabels(none) mlabels((G5) (L5) (R5) (A5), prefix(\multicolumn{1}{c}{) suffix(})) indicate("Demographics = age asian otherethn gpa", labels(\textmd{Yes} \textmd{No})) order(nfolic_fol nperi_fol nlut_fol npremens_fol) /// style(tex) collabels(,none) prehead("\begin{tabular}{l*{@span}{d}}\toprule") /// posthead(\midrule) prefoot(\bottomrule) postfoot("\addlinespace[.1cm]\multicolumn{5}{l}{Robust standard errors in parentheses; Significance levels: * 5\%; ** 1\%.} \\" "\end{tabular}") notype ******************************************* * Analysis Risk & Hormones ******************************************* * by gender ***************************************** * males eststo mg3: reg riskgain ztestf_m zestf_m zprogf_m zcortf_m age asian otherethn gpa if period == 30 & preg == 0 & female == 0 & rgconsistent == 1 & girrational == 0, vce(robust) eststo ml3: reg riskloss ztestf_m zestf_m zprogf_m zcortf_m age asian otherethn gpa if period == 30 & preg == 0 & female == 0 & rlconsistent == 1 & lirrational == 0, vce(robust) eststo mr3: logit reflection ztestf_m zestf_m zprogf_m zcortf_m age asian otherethn gpa if period == 30 & consistent == 1 & irrational == 0 & female == 0, vce(robust) eststo mc3: logit consistent ztestf_m zestf_m zprogf_m zcortf_m age asian otherethn gpa if period == 30 & new == 1 & female == 0, vce(robust) estout mg3 ml3 mr3 mc3 using "C:\bcschipper\burkhard\research\hormone\malehorm_collapsed.tex", replace /// cells(b(star fmt(%9.4f)) se(par fmt(%9.4f))) stats(N r2 r2_p, fmt(%9.0g %9.4f) labels("\textit{Number of Observations}" "R\textsuperscript{2}" "Pseudo R\textsuperscript{2}") layout("\multicolumn{1}{c}{\quad@}")) label /// starlevels(\text{*} 0.05 \text{**} 0.01) keep(ztestf_m zestf_m zprogf_m zcortf_m) /// eqlabels(none) mlabels((G3M) (L3M) (R3M) (A3M), prefix(\multicolumn{1}{c}{) suffix(})) indicate("Demographics = age asian otherethn gpa", labels(\textmd{Yes} \textmd{No})) order(ztestf_m zestf_m zprogf_m zcortf_m) /// style(tex) collabels(,none) prehead("\begin{tabular}{l*{@span}{d}}\toprule") /// posthead(\midrule) prefoot(\bottomrule) postfoot("\addlinespace[.1cm]\multicolumn{5}{l}{Robust standard errors in parentheses; Significance levels: * 5\%; ** 1\%.} \\""\end{tabular}") notype * Robustness checks testosterone result * full standard errors clustered on sessions reg riskgain ztestf_m zestf_m zprogf_m zcortf_m age asian otherethn gpa if period == 30 & preg == 0 & female == 0 & rgconsistent == 1 & girrational == 0, robust cluster(session) * full session dummies reg riskgain ztestf_m zestf_m zprogf_m zcortf_m sessiondum* age asian otherethn gpa if period == 30 & preg == 0 & female == 0 & rgconsistent == 1 & girrational == 0, vce(robust) * full session fixed effects xtset, clear xtset datum xtreg riskgain ztestf_m zestf_m zprogf_m zcortf_m age asian otherethn gpa if period == 30 & preg == 0 & female == 0 & rgconsistent == 1 & girrational == 0, fe * with majors of study reg riskgain ztestf_m zestf_m zprogf_m zcortf_m age asian otherethn gpa math sciothersci socsci humanities econ if period == 30 & preg == 0 & female == 0 & rgconsistent == 1 & girrational == 0, vce(robust) * add further demographics/lifestyle reg riskgain ztestf_m zestf_m zprogf_m zcortf_m age asian otherethn gpa dr bmi siblings workout livalone dating homobi if period == 30 & preg == 0 & female == 0 & rgconsistent == 1 & girrational == 0, vce(robust) * add controls on the hormone quality (luncht foodtime drinktime mouth brush blast dst2010) reg riskgain ztestf_m zestf_m zprogf_m zcortf_m age asian otherethn gpa luncht foodtime drinktime mouth brush blast dst2010 if period == 30 & preg == 0 & female == 0 & rgconsistent == 1 & girrational == 0, vce(robust) * add controls for dietary preferences reg riskgain ztestf_m zestf_m zprogf_m zcortf_m age asian otherethn gpa smoking vegetarian soybean if period == 30 & preg == 0 & female == 0 & rgconsistent == 1 & girrational == 0, vce(robust) * allow inaccessible subjects reg riskgain ztestf_m zestf_m zprogf_m zcortf_m age asian otherethn gpa if period == 30 & preg == 0 & female == 0, vce(robust) * dropping other hormones except test reg riskgain ztestf_m age asian otherethn gpa if period == 30 & preg == 0 & female == 0 & rgconsistent == 1 & girrational == 0, vce(robust) * naturally cycling females eststo fg3: reg riskgain ztestf_f zestf_f zprogf_f zcortf_f age asian otherethn gpa if pill == 0 & period == 30 & preg == 0 & female == 1 & rgconsistent == 1 & girrational == 0, vce(robust) eststo fl3: reg riskloss ztestf_f zestf_f zprogf_f zcortf_f age asian otherethn gpa if pill == 0 & period == 30 & preg == 0 & female == 1 & rlconsistent == 1 & lirrational == 0, vce(robust) eststo fr3: logit reflection ztestf_f zestf_f zprogf_f zcortf_f age asian otherethn gpa if pill == 0 & period == 30 & consistent == 1 & irrational == 0 & female == 1, vce(robust) eststo fc3: logit consistent ztestf_f zestf_f zprogf_f zcortf_f age asian otherethn gpa if pill == 0 & period == 30 & new == 1 & female == 1, vce(robust) estout fg3 fl3 fr3 fc3 using "C:\bcschipper\burkhard\research\hormone\natfemalehorm_collapsed.tex", replace /// cells(b(star fmt(%9.4f)) se(par fmt(%9.4f))) stats(N r2 r2_p, fmt(%9.0g %9.4f) labels("\textit{Number of Observations}" "R\textsuperscript{2}" "Pseudo R\textsuperscript{2}") layout("\multicolumn{1}{c}{\quad@}")) label /// starlevels(\text{*} 0.05 \text{**} 0.01) keep(ztestf_f zestf_f zprogf_f zcortf_f) /// eqlabels(none) mlabels((G3F) (L3F) (R3F) (A3F), prefix(\multicolumn{1}{c}{) suffix(})) indicate("Demographics = age asian otherethn gpa", labels(\textmd{Yes} \textmd{No})) order(ztestf_f zestf_f zprogf_f zcortf_f) /// style(tex) collabels(,none) prehead("\begin{tabular}{l*{@span}{d}}\toprule") /// posthead(\midrule) prefoot(\bottomrule) postfoot("\addlinespace[.1cm]\multicolumn{5}{l}{Robust standard errors in parentheses; Significance levels: * 5\%; ** 1\%.} \\""\end{tabular}") notype * Robustness checks cortisol result * full standard errors clustered on sessions reg riskgain ztestf_m zestf_m zprogf_m zcortf_m age asian otherethn gpa if pill == 0 & period == 30 & preg == 0 & female == 1 & rgconsistent == 1 & girrational == 0, robust cluster(session) * full session dummies reg riskgain ztestf_m zestf_m zprogf_m zcortf_m sessiondum* age asian otherethn gpa if pill == 0 & period == 30 & preg == 0 & female == 1 & rgconsistent == 1 & girrational == 0, vce(robust) * full session fixed effects xtset, clear xtset datum xtreg riskgain ztestf_m zestf_m zprogf_m zcortf_m age asian otherethn gpa if pill == 0 & period == 30 & preg == 0 & female == 1 & rgconsistent == 1 & girrational == 0, fe * with majors of study reg riskgain ztestf_m zestf_m zprogf_m zcortf_m age asian otherethn gpa math sciothersci socsci humanities econ if pill == 0 & period == 30 & preg == 0 & female == 1 & rgconsistent == 1 & girrational == 0, vce(robust) * add further demographics/lifestyle reg riskgain ztestf_m zestf_m zprogf_m zcortf_m age asian otherethn gpa dr bmi siblings workout livalone dating homobi if pill == 0 & period == 30 & preg == 0 & female == 1 & rgconsistent == 1 & girrational == 0, vce(robust) * add controls on the hormone quality (luncht foodtime drinktime mouth brush blast dst2010) reg riskgain ztestf_m zestf_m zprogf_m zcortf_m age asian otherethn gpa luncht foodtime drinktime mouth brush blast dst2010 if pill == 0 & period == 30 & preg == 0 & female == 1 & rgconsistent == 1 & girrational == 0, vce(robust) * add controls on the hormone quality (foodtime dst2010) reg riskgain ztestf_m zestf_m zprogf_m zcortf_m age asian otherethn gpa luncht foodtime drinktime mouth brush blast dst2010 if pill == 0 & period == 30 & preg == 0 & female == 1 & rgconsistent == 1 & girrational == 0, vce(robust) * add controls for dietary preferences reg riskgain ztestf_m zestf_m zprogf_m zcortf_m age asian otherethn gpa smoking vegetarian soybean if pill == 0 & period == 30 & preg == 0 & female == 1 & rgconsistent == 1 & girrational == 0, vce(robust) * allow inaccessible subjects reg riskgain ztestf_m zestf_m zprogf_m zcortf_m age asian otherethn gpa if pill == 0 & period == 30 & preg == 0 & female == 1, vce(robust) * dropping other hormones except test reg riskgain zcortf_m age asian otherethn gpa if pill == 0 & period == 30 & preg == 0 & female == 1 & rgconsistent == 1 & girrational == 0, vce(robust) * females on hormonal contraceptives eststo pg3: reg riskgain ztestf_f zestf_f zprogf_f zcortf_f if pill == 1 & period == 30 & preg == 0 & female == 1, vce(robust) eststo pl3: reg riskloss ztestf_f zestf_f zprogf_f zcortf_f if pill == 1 & period == 30 & preg == 0 & female == 1, vce(robust) * eststo pr3: logit reflection ztestf_f zestf_f zprogf_f zcortf_f if pill == 1 & period == 30 & female == 1, vce(robust) eststo pc3: logit consistent ztestf_f zestf_f zprogf_f zcortf_f if pill == 1 & period == 30 & new == 1 & female == 1, vce(robust) estout pg3 pl3 pc3 using "C:\bcschipper\burkhard\research\hormone\pillhorm_collapsed.tex", replace /// cells(b(star fmt(%9.4f)) se(par fmt(%9.4f))) stats(N r2 r2_p, fmt(%9.0g %9.4f) labels("\textit{Number of Observations}" "R\textsuperscript{2}" "Pseudo R\textsuperscript{2}") layout("\multicolumn{1}{c}{\quad@}")) label /// starlevels(\text{*} 0.05 \text{**} 0.01) keep(ztestf_f zestf_f zprogf_f zcortf_f) /// eqlabels(none) mlabels((G3P) (L3P) (A3P), prefix(\multicolumn{1}{c}{) suffix(})) order(ztestf_f zestf_f zprogf_f zcortf_f) /// style(tex) collabels(,none) prehead("\begin{tabular}{l*{@span}{d}}\toprule") /// posthead(\midrule) prefoot(\bottomrule) postfoot("\addlinespace[.1cm]\multicolumn{4}{l}{Robust standard errors in parentheses; Significance levels: * 5\%; ** 1\%.} \\""\end{tabular}") notype * Stanton et al. (2011b) hypothesis: ************************************ gen ztestf_m2 = ztestf_m^2 gen ztestf_f2 = ztestf_f^2 reg riskgain ztestf_m ztestf_m2 zestf_m zprogf_m zcortf_m age asian otherethn gpa if period == 30 & preg == 0 & female == 0 & rgconsistent == 1 & girrational == 0, vce(robust) reg riskloss ztestf_m ztestf_m2 zestf_m zprogf_m zcortf_m age asian otherethn gpa if period == 30 & preg == 0 & female == 0 & rlconsistent == 1 & lirrational == 0, vce(robust) reg riskgain ztestf_f ztestf_f2 zestf_f zprogf_f zcortf_f age asian otherethn gpa if pill == 0 & period == 30 & preg == 0 & female == 1 & rgconsistent == 1 & girrational == 0, vce(robust) reg riskloss ztestf_f ztestf_f2 zestf_f zprogf_f zcortf_f age asian otherethn gpa if pill == 0 & period == 30 & preg == 0 & female == 1 & rlconsistent == 1 & lirrational == 0, vce(robust) reg riskgain ztestf_f ztestf_f2 zestf_f zprogf_f zcortf_f age asian otherethn gpa if pill == 1 & period == 30 & preg == 0 & female == 1 & rgconsistent == 1 & girrational == 0, vce(robust) reg riskloss ztestf_f ztestf_f2 zestf_f zprogf_f zcortf_f age asian otherethn gpa if pill == 1 & period == 30 & preg == 0 & female == 1 & rlconsistent == 1 & lirrational == 0, vce(robust) * Robustness checks nonlinear testosterone result * full standard errors clustered on sessions reg riskgain ztestf_m ztestf_m2 zestf_m zprogf_m zcortf_m age asian otherethn gpa if period == 30 & preg == 0 & female == 0 & rgconsistent == 1 & girrational == 0, robust cluster(session) * full session dummies reg riskgain ztestf_m ztestf_m2 zestf_m zprogf_m zcortf_m sessiondum* age asian otherethn gpa if period == 30 & preg == 0 & female == 0 & rgconsistent == 1 & girrational == 0, vce(robust) * full session fixed effects xtset, clear xtset datum xtreg riskgain ztestf_m ztestf_m2 zestf_m zprogf_m zcortf_m age asian otherethn gpa if period == 30 & preg == 0 & female == 0 & rgconsistent == 1 & girrational == 0, fe * with majors of study reg riskgain ztestf_m ztestf_m2 zestf_m zprogf_m zcortf_m age asian otherethn gpa math sciothersci socsci humanities econ if period == 30 & preg == 0 & female == 0 & rgconsistent == 1 & girrational == 0, vce(robust) * add further demographics/lifestyle reg riskgain ztestf_m ztestf_m2 zestf_m zprogf_m zcortf_m age asian otherethn gpa dr bmi siblings workout livalone dating homobi if period == 30 & preg == 0 & female == 0 & rgconsistent == 1 & girrational == 0, vce(robust) * add controls on the hormone quality (luncht foodtime drinktime mouth brush blast dst2010) reg riskgain ztestf_m ztestf_m2 zestf_m zprogf_m zcortf_m age asian otherethn gpa luncht foodtime drinktime mouth brush blast dst2010 if period == 30 & preg == 0 & female == 0 & rgconsistent == 1 & girrational == 0, vce(robust) * add controls for dietary preferences reg riskgain ztestf_m ztestf_m2 zestf_m zprogf_m zcortf_m age asian otherethn gpa smoking vegetarian soybean if period == 30 & preg == 0 & female == 0 & rgconsistent == 1 & girrational == 0, vce(robust) * allow inaccessible subjects reg riskgain ztestf_m ztestf_m2 zestf_m zprogf_m zcortf_m age asian otherethn gpa if period == 30 & preg == 0 & female == 0, vce(robust) * dropping other hormones except test reg riskgain ztestf_m ztestf_m2 age asian otherethn gpa if period == 30 & preg == 0 & female == 0 & rgconsistent == 1 & girrational == 0, vce(robust) *************************************** * Analysis Risk & Digit Ratio *************************************** * males eststo gdrm: reg riskgain dr age asian otherethn gpa if female == 0 & period == 30 & rgconsistent == 1 & girrational == 0, vce(robust) eststo ldrm: reg riskloss dr age asian otherethn gpa if female == 0 & period == 30 & rlconsistent == 1 & lirrational == 0, vce(robust) eststo rdrm: logit reflection dr age asian otherethn gpa if female == 0 & period == 30 & consistent == 1 & irrational == 0, vce(robust) eststo cdrm: logit consistent dr age asian otherethn gpa if female == 0 & period == 30, vce(robust) estout gdrm ldrm rdrm cdrm using "C:\bcschipper\burkhard\research\hormone\dr_males.tex", replace /// cells(b(star fmt(%9.4f)) se(par fmt(%9.4f))) stats(N r2 r2_p, fmt(%9.0g %9.4f) labels("\textit{Number of Observations}" "R\textsuperscript{2}" "Pseudo R\textsuperscript{2}") layout("\multicolumn{1}{c}{\quad@}")) label /// starlevels(\text{*} 0.05 \text{**} 0.01) keep(dr) /// eqlabels(none) mlabels((GDRM) (LDRM) (RDRM) (ADRM), prefix(\multicolumn{1}{c}{) suffix(})) indicate("Demographics = age asian otherethn gpa", labels(\textmd{Yes} \textmd{No})) order(dr) /// style(tex) collabels(,none) prehead("\begin{tabular}{l*{@span}{d}}\toprule") /// posthead(\midrule) prefoot(\bottomrule) postfoot("\addlinespace[.1cm]\multicolumn{5}{l}{Robust standard errors in parentheses; Significance levels: * 5\%; ** 1\%.} \\" "\end{tabular}") notype * females eststo gdrf: reg riskgain dr age asian otherethn gpa if female == 1 & period == 30 & rgconsistent == 1 & girrational == 0, vce(robust) eststo ldrf: reg riskloss dr age asian otherethn gpa if female == 1 & period == 30 & rlconsistent == 1 & lirrational == 0, vce(robust) eststo rdrf: logit reflection dr age asian otherethn gpa if female == 1 & period == 30 & consistent == 1 & irrational == 0, vce(robust) eststo cdrf: logit consistent dr age asian otherethn gpa if female == 1 & period == 30, vce(robust) estout gdrf ldrf rdrf cdrf using "C:\bcschipper\burkhard\research\hormone\dr_females.tex", replace /// cells(b(star fmt(%9.4f)) se(par fmt(%9.4f))) stats(N r2 r2_p, fmt(%9.0g %9.4f) labels("\textit{Number of Observations}" "R\textsuperscript{2}" "Pseudo R\textsuperscript{2}") layout("\multicolumn{1}{c}{\quad@}")) label /// starlevels(\text{*} 0.05 \text{**} 0.01) keep(dr) /// eqlabels(none) mlabels((GDRF) (LDRF) (RDRF) (ADRF), prefix(\multicolumn{1}{c}{) suffix(})) indicate("Demographics = age asian otherethn gpa", labels(\textmd{Yes} \textmd{No})) order(dr) /// style(tex) collabels(,none) prehead("\begin{tabular}{l*{@span}{d}}\toprule") /// posthead(\midrule) prefoot(\bottomrule) postfoot("\addlinespace[.1cm]\multicolumn{5}{l}{Robust standard errors in parentheses; Significance levels: * 5\%; ** 1\%.} \\" "\end{tabular}") notype * white males eststo gdrwm: reg riskgain dr age gpa if white == 1 & female == 0 & period == 30 & rgconsistent == 1 & girrational == 0, vce(robust) eststo ldrwm: reg riskloss dr age gpa if white == 1 & female == 0 & period == 30 & rlconsistent == 1 & lirrational == 0, vce(robust) eststo rdrwm: logit reflection dr age gpa if white == 1 & female == 0 & period == 30 & consistent == 1 & irrational == 0, vce(robust) eststo cdrwm: logit consistent dr age gpa if white == 1 & female == 0 & period == 30, vce(robust) estout gdrwm ldrwm rdrwm cdrwm using "C:\bcschipper\burkhard\research\hormone\dr_white_males.tex", replace /// cells(b(star fmt(%9.4f)) se(par fmt(%9.4f))) stats(N r2 r2_p, fmt(%9.0g %9.4f) labels("\textit{Number of Observations}" "R\textsuperscript{2}" "Pseudo R\textsuperscript{2}") layout("\multicolumn{1}{c}{\quad@}")) label /// starlevels(\text{*} 0.05 \text{**} 0.01) keep(dr) /// eqlabels(none) mlabels((GDRWM) (LDRWM) (RDRWM) (ADRWM), prefix(\multicolumn{1}{c}{) suffix(})) indicate("Demographics = age gpa", labels(\textmd{Yes} \textmd{No})) order(dr) /// style(tex) collabels(,none) prehead("\begin{tabular}{l*{@span}{d}}\toprule") /// posthead(\midrule) prefoot(\bottomrule) postfoot("\addlinespace[.1cm]\multicolumn{5}{l}{Robust standard errors in parentheses; Significance levels: * 5\%; ** 1\%.} \\" "\end{tabular}") notype * white females eststo gdrwf: reg riskgain dr age gpa if white == 1 & female == 1 & period == 30 & rgconsistent == 1 & girrational == 0, vce(robust) eststo ldrwf: reg riskloss dr age gpa if white == 1 & female == 1 & period == 30 & rlconsistent == 1 & lirrational == 0, vce(robust) eststo rdrwf: logit reflection dr age gpa if white == 1 & female == 1 & period == 30 & consistent == 1 & irrational == 0, vce(robust) eststo cdrwf: logit consistent dr age gpa if white == 1 & female == 1 & period == 30, vce(robust) estout gdrwf ldrwf rdrwf cdrwf using "C:\bcschipper\burkhard\research\hormone\dr_white_females.tex", replace /// cells(b(star fmt(%9.4f)) se(par fmt(%9.4f))) stats(N r2 r2_p, fmt(%9.0g %9.4f) labels("\textit{Number of Observations}" "R\textsuperscript{2}" "Pseudo R\textsuperscript{2}") layout("\multicolumn{1}{c}{\quad@}")) label /// starlevels(\text{*} 0.05 \text{**} 0.01) keep(dr) /// eqlabels(none) mlabels((GDRWF) (LDRWF) (RDRWF) (ADRWF), prefix(\multicolumn{1}{c}{) suffix(})) indicate("Demographics = age gpa", labels(\textmd{Yes} \textmd{No})) order(dr) /// style(tex) collabels(,none) prehead("\begin{tabular}{l*{@span}{d}}\toprule") /// posthead(\midrule) prefoot(\bottomrule) postfoot("\addlinespace[.1cm]\multicolumn{5}{l}{Robust standard errors in parentheses; Significance levels: * 5\%; ** 1\%.} \\" "\end{tabular}") notype * asian males eststo gdram: reg riskgain dr age gpa if asian == 1 & female == 0 & period == 30 & rgconsistent == 1 & girrational == 0, vce(robust) eststo ldram: reg riskloss dr age gpa if asian == 1 & female == 0 & period == 30 & rlconsistent == 1 & lirrational == 0, vce(robust) eststo rdram: logit reflection dr age gpa if asian == 1 & female == 0 & period == 30 & consistent == 1 & irrational == 0, vce(robust) eststo cdram: logit consistent dr age gpa if asian == 1 & female == 0 & period == 30, vce(robust) estout gdram ldram rdram cdram using "C:\bcschipper\burkhard\research\hormone\dr_asian_males.tex", replace /// cells(b(star fmt(%9.4f)) se(par fmt(%9.4f))) stats(N r2 r2_p, fmt(%9.0g %9.4f) labels("\textit{Number of Observations}" "R\textsuperscript{2}" "Pseudo R\textsuperscript{2}") layout("\multicolumn{1}{c}{\quad@}")) label /// starlevels(\text{*} 0.05 \text{**} 0.01) keep(dr) /// eqlabels(none) mlabels((GDRAM) (LDRAM) (RDRAM) (ADRAM), prefix(\multicolumn{1}{c}{) suffix(})) indicate("Demographics = age gpa", labels(\textmd{Yes} \textmd{No})) order(dr) /// style(tex) collabels(,none) prehead("\begin{tabular}{l*{@span}{d}}\toprule") /// posthead(\midrule) prefoot(\bottomrule) postfoot("\addlinespace[.1cm]\multicolumn{5}{l}{Robust standard errors in parentheses; Significance levels: * 5\%; ** 1\%.} \\" "\end{tabular}") notype * Robustness check reflection effect in asian males * "quality" of fingers logit reflection dr broken broken2d4d uncleardr age gpa if asian == 1 & female == 0 & period == 30 & consistent == 1 & irrational == 0, vce(robust) * dropping demographics and adding inaccessible subjects logit reflection dr if asian == 1 & female == 0 & period == 30, vce(robust) * asian females eststo gdraf: reg riskgain dr age gpa if asian == 1 & female == 1 & period == 30 & rgconsistent == 1 & girrational == 0, vce(robust) eststo ldraf: reg riskloss dr age gpa if asian == 1 & female == 1 & period == 30 & rlconsistent == 1 & lirrational == 0, vce(robust) eststo rdraf: logit reflection dr age gpa if asian == 1 & female == 1 & period == 30 & consistent == 1 & irrational == 0, vce(robust) eststo cdraf: logit consistent dr age gpa if asian == 1 & female == 1 & period == 30, vce(robust) estout gdraf ldraf rdraf cdraf using "C:\bcschipper\burkhard\research\hormone\dr_asian_females.tex", replace /// cells(b(star fmt(%9.4f)) se(par fmt(%9.4f))) stats(N r2 r2_p, fmt(%9.0g %9.4f) labels("\textit{Number of Observations}" "R\textsuperscript{2}" "Pseudo R\textsuperscript{2}") layout("\multicolumn{1}{c}{\quad@}")) label /// starlevels(\text{*} 0.05 \text{**} 0.01) keep(dr) /// eqlabels(none) mlabels((GDRAF) (LDRAF) (RDRAF) (ADRAF), prefix(\multicolumn{1}{c}{) suffix(})) indicate("Demographics = age gpa", labels(\textmd{Yes} \textmd{No})) order(dr) /// style(tex) collabels(,none) prehead("\begin{tabular}{l*{@span}{d}}\toprule") /// posthead(\midrule) prefoot(\bottomrule) postfoot("\addlinespace[.1cm]\multicolumn{5}{l}{Robust standard errors in parentheses; Significance levels: * 5\%; ** 1\%.} \\" "\end{tabular}") notype log close log using "C:\bcschipper\burkhard\research\hormone\stats\cordelia_log_20150616.smcl", replace ************************************************************************* ************************************************************************* * Creating log file for a meta-study by Cordelia Fine, Melbourne * Email dated March 4, 2015 * Updated June 16, 2015, Email from Anna Genat dated June 14, 2015 ************************************************************************* ************************************************************************* * Correlation between risk aversion and testosterone in the gain domain ************************************************************************* * Males correlate riskgain testf if female == 0 & period == 30 reg riskgain testf if female == 0 & period == 30 * Females correlate riskgain testf if female == 1 & period == 30 reg riskgain testf if female == 1 & period == 30 * Correlation between risk aversion and testosterone in the loss domain ************************************************************************* * Males correlate riskloss testf if female == 0 & period == 30 reg riskloss testf if female == 0 & period == 30 * Females correlate riskloss testf if female == 1 & period == 30 reg riskloss testf if female == 1 & period == 30 * Correlation between risk aversion and digit ratio in the gain domain ************************************************************************* * Males correlate riskgain dr if female == 0 & period == 30 reg riskgain dr if female == 0 & period == 30 * Females correlate riskgain dr if female == 1 & period == 30 reg riskgain dr if female == 1 & period == 30 * Correlation between risk aversion and digit ratio in the loss domain ************************************************************************* * Males correlate riskloss dr if female == 0 & period == 30 reg riskloss dr if female == 0 & period == 30 * Females correlate riskloss dr if female == 1 & period == 30 reg riskloss dr if female == 1 & period == 30 log close **************************************************************************** **************************************************************************** * Revision for Journal of Economic Psychology November 2022 **************************************************************************** **************************************************************************** ************************************************************ * Referee 2 wants correlations of all biological measures ************************************************************ corrtex female dr premens menst folic peri lut testf estf progf cortf tests ests progs corts if period == 30, land file(C:\bcschipper\burkhard\research\hormone\biocorr.tex) replace ************************************************************* * Referee 2 wants test of differences of digit ratios ************************************************************* gen drwhite = dr if white==1 & new == 1 gen drasian = dr if asian==1 & new == 1 gen drothers = dr if (hispanic==1 | black == 1 | other == 1) & new == 1 sum drwhite drasian drothers if period == 30 & new == 1 sum drwhite drasian drothers if period == 30 & female==1 & new == 1 sum drwhite drasian drothers if period == 30 & female==0 & new == 1 ttest dr if period == 30 & new == 1, by(female) ttest drwhite if period == 30 & new == 1, by(female) ttest drasian if period == 30 & new == 1, by(female) ttest drothers if period == 30 & new == 1, by(female) ************************************************************************ * Dual hormone hypothesis Request by Reviewer 2022 (Post hoc analysis) ************************************************************************ gen ztestcort_m = ztestf_m*zcortf_m reg riskgain ztestf_m zcortf_m ztestcort_m age asian otherethn gpa if period == 30 & preg == 0 & female == 0 & rgconsistent == 1 & girrational == 0, vce(robust) reg riskloss ztestf_m zcortf_m ztestcort_m age asian otherethn gpa if period == 30 & preg == 0 & female == 0 & rlconsistent == 1 & lirrational == 0, vce(robust) sum zcortf_m gen avzcortf_m = r(mean) gen lowcortf_m = 0 replace lowcortf_m = 1 if zcortf_m < avzcortf_m reg riskgain ztestf_m zcortf_m lowcortf_m age asian otherethn gpa if period == 30 & preg == 0 & female == 0 & rgconsistent == 1 & girrational == 0, vce(robust) reg riskloss ztestf_m zcortf_m lowcortf_m age asian otherethn gpa if period == 30 & preg == 0 & female == 0 & rlconsistent == 1 & lirrational == 0, vce(robust) gen rattestcortf_m = ztestf_m/zcortf_m reg riskgain rattestcortf_m age asian otherethn gpa if period == 30 & preg == 0 & female == 0 & rgconsistent == 1 & girrational == 0, vce(robust) reg riskloss rattestcortf_m age asian otherethn gpa if period == 30 & preg == 0 & female == 0 & rlconsistent == 1 & lirrational == 0, vce(robust) gen ztestcort_f = ztestf_f*zcortf_f reg riskgain ztestf_f zcortf_f ztestcort_f age asian otherethn gpa if period == 30 & preg == 0 & female == 1 & pill == 0 & rgconsistent == 1 & girrational == 0, vce(robust) reg riskloss ztestf_f zcortf_f ztestcort_f age asian otherethn gpa if period == 30 & preg == 0 & female == 1 & pill == 0 & rlconsistent == 1 & lirrational == 0, vce(robust) sum zcortf_f gen avzcortf_f = r(mean) gen lowcortf_f = 0 replace lowcortf_f = 1 if zcortf_f < avzcortf_f reg riskgain ztestf_f zcortf_f lowcortf_f age asian otherethn gpa if period == 30 & preg == 0 & female == 1 & pill == 0 & rgconsistent == 1 & girrational == 0, vce(robust) reg riskloss ztestf_f zcortf_f lowcortf_f age asian otherethn gpa if period == 30 & preg == 0 & female == 1 & pill == 0 & rlconsistent == 1 & lirrational == 0, vce(robust) gen rattestcortf_f = ztestf_f/zcortf_f reg riskgain rattestcortf_f age asian otherethn gpa if period == 30 & preg == 0 & female == 1 & pill == 0 & rgconsistent == 1 & girrational == 0, vce(robust) reg riskloss rattestcortf_f age asian otherethn gpa if period == 30 & preg == 0 & female == 1 & pill == 0 & rlconsistent == 1 & lirrational == 0, vce(robust) ******************************************************************************* * T-tests for differences of accessibility in gain and loss domain (Referee 2) ******************************************************************************* sum rgconsistent rlconsistent if period == 30 & new == 1 ttest rgconsistent == rlconsistent if period == 30 & new == 1 ttest rgconsistent == rlconsistent if period == 30 & new == 1 & female == 0 ttest rgconsistent == rlconsistent if period == 30 & new == 1 & female == 1 ./. ******************************************************************************** **************best run separately ********************************************** ******************************************************************************** ******************************************************************************** ** Intraclass Correlation of Digit Measures (Requested by Referee 2) ******************************************************************************** clear clear matrix set more off set mem 50m use "C:\bcschipper\burkhard\research\hormone\stats\hormone_dr_20110728.dta" sort session subject save "C:\bcschipper\burkhard\research\hormone\stats\hormone_dr_ICC.dta", replace clear use "C:\bcschipper\burkhard\research\hormone\stats\hormone_dr_ICC.dta", replace ******************************** * GENERATE ARBITRARY ID VARIABLE ******************************** gen test = substr(session, 3, 7) gen test2 = regexr(test,"_16","1") gen test3 = regexr(test2,"_17","2") gen test4 = regexr(test3,"_15","3") destring test4, gen(test5) egen id=concat(test5 subject) destring id, replace drop test test2 test3 test4 test5 label variable id "Arbitrary ID number" ******************************************** * Formatting the data set for ICC Analysis ******************************************** stack id twod_bcs id fourd_bcs id twod_mp id fourd_mp, into(id laenge) clear wide gen rater = 1 replace rater = 2 if (twod_mp > 0 & !missing(twod_mp)) | (fourd_mp > 0 & !missing(fourd_mp)) label variable rater "Rater" gen finger = 2 replace finger = 4 if (fourd_bcs > 0 & !missing(fourd_bcs)) | (fourd_mp > 0 & !missing(fourd_mp)) label variable finger "Digit" label variable laenge "Length" gen laenge2 = laenge if (twod_mp > 0 & !missing(twod_mp)) | (twod_bcs > 0 & !missing(twod_bcs)) gen laenge4 = laenge if (fourd_mp > 0 & !missing(fourd_mp)) | (fourd_bcs > 0 & !missing(fourd_bcs)) ********************************************* * Two-way mixed ICC Analysis ********************************************* icc laenge2 id rater, mixed icc laenge4 id rater, mixed