stata中对数命令
倾向值匹配、IPW等
Stata的因果推断功能允许您从观察数据中估计实验类型的因果效应(治疗效应)。是否对连续、二元、计数或生存结果感兴趣? 是在模拟结果过程还是治疗过程?Stata可以估算您的治疗效果。有了这样一套全面的估计量,您就会找到适合您的那一个。
估计量
• 双重差分(DID)
• 异质性DID(新增)
• 逆概率权重(IPW)
• 倾向得分匹配
• 协变量匹配
• 回归调整
• 加权回归
• 双重稳健方法
– 增强型IPW(AIPW)
– 带回归调整的IPW
– 具有lasso选择控件的AIPW
• 因果中介(新增)
统计
• 平均处理效应(ATEs)
• 处理后的ATEs(ATETs)
• 潜在结果手段(POM)
• 直接影响、间接影响
结果
• 连续—线性
• 二进制—logistic、概率、异方差
• 计数—Poisson
• 分数
• 非负,包括指数平均值
• 生存率—指数、Weibull、伽玛、对数正态
治疗
• 二进制—logistic、概率、异方差
• 多值——多项logistic
诊断
• 重叠图
• 协变量平衡
内生性处理效应
• 连续、截尾、二进制、有序和计数结果
• ATEs、ATETs和POMs
• 结合内生协变量、样本选择和面板数据
• 内生性检验
也许你是一名医学研究人员,知道医生在决定是否开药治疗高血压时考虑的变量,但你对应该用来模拟血压的变量却远没有那么自信。
您可以使用IPW估计器,
• teffects ipw (bp) (drug x1 x2)
或倾向性得分匹配
• teffects psmatch (bp) (drug x1 x2)
了解更多影响血压的变量,但不了解哪些决定是否开药的变量?
使用回归调整,
• teffects ra (bp x1 x3) (drug)
或最近邻匹配,
• teffects nnmatch (bp x1 x3) (drug)
如果你同时知道血压建模和是否开具药物处方,你可以使用双重稳健估计量中的一个。
使用增强的IPW,
• teffects aipw (bp x1 x3) (drug x1 x2)
或采用回归调整的IPW,
• teffects ipwra (bp x1 x3) (drug x1 x2)
令人惊讶的是,使用这些双重稳健的方法,我们只需要正确对待两个模型规范中的一个。
让我们看一个使用带回归调整的IPW例子。我们将收缩压建模为年龄、高血压家族史、吸烟和体重的函数。是否开药被模拟为年龄、家族史和体重的函数。
• teffects ipwra
(bp age fam_history smokes weight)
(drug age fam_history weight)
ATE为-21.4。若所有患者均按医嘱服药,平均血压将下降21.4mm Hg。在没有人给药的情况下,167.5的POM为我们提供了估计的收缩压。
让我们开看一些诊断。协变量是否平衡?我们可以查看年龄的核密度图,例如,比较处理过的和未处理过的。
• tebalance density age
加权后,年龄的密度非常相似。
我们可以对所有协变量进行平衡性检验。
• tebalance overid
我们并不反对所有协变量都是平衡的这一原假设。
我们只谈到了可用的治疗效果估计量的范围。
你有二元结果吗?
• teffects ipwra (y x1 x2, probit) (treat x1 x3)
或是生存结果?
• stteffects ipwra (x1 x2) (treat x1 x3)
是否存在内生性问题?
• eteffects (y x1 x2) (treat x1 x3)
你有很多潜在的控制因素吗?
• telasso (y x1 x2 x3 x4 ...)
(treat x1 x2 x3 x4 ...)
你有中介变量吗?
• mediate (y x1 x2) (m x1 x3) (treat)
等等……
北京天演融智软件有限公司(科学软件网)是STATA软件在中国的授权经销商,为中国的软件用户提供优质的软件销售和培训服务。
","gnid":"90368c0f7b5da87ea","img_data":[{"flag":2,"img":[{"desc":"","height":"383","title":"","url":"https://p0.ssl.img.360kuai.com/t0173f3da782f36af1c.png","width":"900"},{"desc":"","height":"360","title":"","url":"https://p0.ssl.img.360kuai.com/t015bd43ed3d5c90515.png","width":"606"},{"desc":"","height":"698","title":"","url":"https://p0.ssl.img.360kuai.com/t011010d5360808171c.png","width":"680"},{"desc":"","height":"374","title":"","url":"https://p0.ssl.img.360kuai.com/t01bae8dac627cd765c.png","width":"706"},{"desc":"","height":"323","title":"","url":"https://p0.ssl.img.360kuai.com/t01bce6f80d24249957.png","width":"529"},{"desc":"","height":"154","title":"","url":"https://p0.ssl.img.360kuai.com/t01b1d9efaff5e2cc46.png","width":"531"}]}],"original":0,"pat":"art_src_0,fts0,sts0","powerby":"pika","pub_time":1712479719000,"pure":"","rawurl":"http://zm.news.so.com/c4f62740a72fdcf19f2007480070a569","redirect":0,"rptid":"09a640408d231305","rss_ext":[],"s":"t","src":"天演融智","tag":[],"title":"Stata软件之因果推断
凌往宗4671如何在STATA中添加新的命令 -
麻皆须15912418229 ______ 首先,打开stata软件.我们需要确定自己安装的命令是否已经在stata系统之中,简单的一个测试方法便是在输入框中输入“help XXX”,以“fsum”名利为例,输入“h fsum”.如果没有该命令,则会出现图片的提示.这里也可以选择“是”来软件自动搜索,但是我们提供另外一个下载方式.我们可以直接在stata命令窗口中输入下载命令:“ssc install fsum”命令.4 等待下载,当出现以下图样的时候,命令便下载完成了.5 重复上面第二个步骤,便可以查看新命令的功能了.
凌往宗4671stata怎么用命令计算有几个变量 -
麻皆须15912418229 ______ stata怎么用命令计算有几个变量 dis `c(k)'可以看到数据中含有的变量的个数.不过这个命令会把“股票代码”和“日期”也包含在内.
凌往宗4671如何在stata中做GMM -
麻皆须15912418229 ______ 输入命令: gmm moment_prog [if] [in] [weight], /// {equations(namelist)|nequations(#)} /// {parameters(namelist)|nparameters(#)} [options] /// [program_options]
凌往宗4671stata 怎么建command log -
麻皆须15912418229 ______ 你要先log close或cap clear后才能新建,因为内存中在某一时段只能打开一个,你要把内存中的log文件clear后,才能新建.
凌往宗4671stata 中ppml命令中变量能去对数吗 -
麻皆须15912418229 ______ foreach var of varlist x1-x100{ gen log`var'=log(`var') }
凌往宗4671如何在STATA中实现nested Logit回归 -
麻皆须15912418229 ______ 在stata中多元的logit命令是:mlogit y x,base(1) y是你的因变量 x你的自变量 base(1)的意思是你选择第一选项为参照项
凌往宗4671stata什么命令可以列 -
麻皆须15912418229 ______ 问题没有问明白,以下是常用命令 更改目录命令 ·cd c:/statadata stata的工作目录将由默认工作目录切换到c盘的statadata文件夹下.在stata中可以用正斜线(/)也可以用反斜线(\),但建议用正斜线.在实际中由于文件名路径太长所以通常...
凌往宗4671在matlab中取模和取对数分别应该用什么函数? -
麻皆须15912418229 ______ 取模命令:abs 对数命令:log或log10(分别是自然对数和以10为低的对数)
凌往宗4671如何用stata做协整 -
麻皆须15912418229 ______ 命令格式: vecrank depvarlist [if] [in] [, options] 例: webuse balance2 vecrank y i c vecrank y i c, lags(5) vecrank y i c, lags(5) level99
凌往宗4671如何用stata 做一个相关性分析的矩阵? -
麻皆须15912418229 ______ 在stata里help cor. stata的命令名是correlate [varlist] [if] [in] [weight] [, correlate_options] stata 里面分析相关性的命令是 pwcorr a b c d e , sig 结果就有了包括了显著性的判断标准,stata里面没有星星,直接根据sig,也就是p的值来判断是否显著...