***** STATA/SE 16.0 DO-FILE CODE FOR: "Reflection increases belief in God through self-questioning among non-believers" (Yilmaz & Isler, 2019) ***** * Contents: Analysis code for Experiment 1 and Experiment 2 * ** EXPERIMENT 1 ** * Data preparation * Treatment numeric gen TreatmentN = 0 replace TreatmentN = 1 if Treatment=="TP" replace TreatmentN = 2 if Treatment=="TD" * Religious believer (1) vs atheists and agnostics (2) [excluding those who chose "other" in religious affliation] gen Believer = 1 if Affliation!="Other" replace Believer = 0 if Affliation=="Agnostic" | Affliation=="Atheist" * Composite intuition manipulation score gen AvMan_Check = (Man_Check_Q1+ Man_Check_Q2)/2 * Descriptive statistics tab Gender tabstat Age, stat(mean sd) tab Treatment tab Affliation * Manipulation check * Questions tabstat AvMan_Check, by(Treatment) stat(mean sd) ci means AvMan_Check if TreatmentN==0 ci means AvMan_Check if TreatmentN==1 ci means AvMan_Check if TreatmentN==2 anova AvMan_Check TreatmentN estat esize pwmean AvMan_Check , over(TreatmentN) mcompare(tukey) effects * Response times tabstat Belief_T, by(Treatment) stat(mean sd) ci means Belief_T if TreatmentN==0 ci means Belief_T if TreatmentN==1 ci means Belief_T if TreatmentN==2 anova Belief_T TreatmentN estat esize pwmean Belief_T , over(TreatmentN) mcompare(tukey) effects * Hypothesis tests * Belief in God tabstat Belief_In_God, by(Treatment) stat(mean sd) ci means Belief_In_God if TreatmentN==0 ci means Belief_In_God if TreatmentN==1 ci means Belief_In_God if TreatmentN==2 anova Belief_In_God TreatmentN estat esize pwmean Belief_In_God , over(TreatmentN) mcompare(tukey) effects * Exploratory analysis * Controlling for religiosity anova Belief_In_God TreatmentN c.Religiosity estat esize * Controlling for religiosity and comparing time-pressure vs. time-delay anova Belief_In_God TreatmentN c.Religiosity if Treatment!="Control" estat esize * Believers only tabstat Belief_In_God if Believer==1, by(Treatment) stat(mean sd) ci means Belief_In_God if TreatmentN==0 & Believer==1 ci means Belief_In_God if TreatmentN==1 & Believer==1 ci means Belief_In_God if TreatmentN==2 & Believer==1 anova Belief_In_God TreatmentN if Believer==1 estat esize * Atheists and agnostics only tabstat Belief_In_God if Believer==0, by(Treatment) stat(mean sd) ci means Belief_In_God if TreatmentN==0 & Believer==0 ci means Belief_In_God if TreatmentN==1 & Believer==0 ci means Belief_In_God if TreatmentN==2 & Believer==0 anova Belief_In_God TreatmentN if Believer==0 estat esize * Error terms anova Belief_In_God TreatmentN if Believer==0 predict w if Believer==0, resid kdensity w, normal swilk w ** EXPERIMENT 2 ** ** Analysis in "wide" data format * Data preparation * CRT 2 Score gen CRT_1C = 0 replace CRT_1C = 1 if CRT2_1=="2" gen CRT_2C = 0 replace CRT_2C = 1 if CRT2_2==8 gen CRT_3C = 0 replace CRT_3C = 1 if CRT2_3=="Emily" gen CRT_4C = 0 replace CRT_4C = 1 if CRT2_4=="0" gen CRT_Score = CRT_1C+CRT_2C+CRT_3C+CRT_4C * Change in Belief gen Change_in_Belief = Belief_Reflection-Belief_Intuition * Religious believer (1) vs atheists and agnostics (2) [excluding those who chose "other" in religious affliation] gen Believer = 1 if Affliation!="Other" replace Believer = 0 if Affliation=="Agnostic" | Affliation=="Atheist" * Belief revision types * Upward revision gen Upward = 0 replace Upward = 1 if Belief_Reflection>Belief_Intuition * Downward revision gen Downward = 0 replace Downward = 1 if Belief_Reflection