Dear Statalist:
I am learning a bit more about the user-written -sensemakr- package. If you're not familiar with it, it's an implemention of some ideas in Cinelli & Hazlett (2020, in the Journal of the Royal Statistical Society: Statistical Methodology, Series B). You can run an OLS regression within it, and it will report how big the effect size of a potential confounder would have to be in order to reduce your coefficient of interest to zero.
The problem is that it can only run basic OLS commands within sensemakr itself. I am wondering if there is a way to run more complex OLS models and then pass them to -sensemakr-. I think it may be possible because that functionality is available for its sister package in R, but I'm pretty sure that if it's a possibility in Stata that it is undocumented.
This obviously doesn't work, but I'm wondering if there's a way to use some of the options in regress that can then be used within -sensemakr- in some way. I also tried:
And that also didn't work, so we can rule that one out.
Thank you for your time,
Jonathan
(Apologies in advance: I will probably post a combined R / Stata question to stackexchange soon, but I'm going to focus mostly on R there and I wanted to see what this community had to say first)
I am learning a bit more about the user-written -sensemakr- package. If you're not familiar with it, it's an implemention of some ideas in Cinelli & Hazlett (2020, in the Journal of the Royal Statistical Society: Statistical Methodology, Series B). You can run an OLS regression within it, and it will report how big the effect size of a potential confounder would have to be in order to reduce your coefficient of interest to zero.
Code:
sysuse auto2.dta regress price mpg weight rep78 sensemakr price mpg weight rep78, treat(rep78) bench(weight)
Code:
sysuse auto2.dta regress price mpg weight rep78, vce(bootstrap, reps(100)) cluster(foreign) est sto m1 sensemakr m1, treat(rep78) bench(weight)
Code:
sysuse auto2.dta bootstrap, cluster(foreign) reps(100): sensemakr price mpg weight rep78, treat(rep78) bench(weight)
Thank you for your time,
Jonathan
(Apologies in advance: I will probably post a combined R / Stata question to stackexchange soon, but I'm going to focus mostly on R there and I wanted to see what this community had to say first)
Comment