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":"cache","pub_time":1712479719000,"pure":"","rawurl":"http://zm.news.so.com/c4f62740a72fdcf19f2007480070a569","redirect":0,"rptid":"09a640408d231305","rss_ext":[],"s":"t","src":"天演融智","tag":[],"title":"Stata软件之因果推断
文美榕3505stata数据匹配求和 求助 -
郎庄琰18388004314 ______ 你要问什么方面的赋值? 一些简单的 (1)gen x=1 (创造X变量,使X变量为1) (2)gen x=. (创造X变量,不赋值) gen y=1 replace x=y-1 (使x变量变为y变量减一,也就是0)
文美榕3505如何求对数平均值 -
郎庄琰18388004314 ______[答案] 已知:x1,x2,.,xn 均大于零;取对数:lgx1,lgx2,.,lgxn,那么对数平均值为:(lgx1+lgx2+.+lgxn)/n (1) //:把n个数据取对数加起来再除以数据个数n由于:lgx1+lgx2+...+lgxn=lg(x1*x2*...*xn) 因此对数平均值也可写成...
文美榕3505怎样对一个数求对数 -
郎庄琰18388004314 ______ 对数函数其实就是指数函数的反函数,所以如果不用计算器的话,只能求比较特殊的,比如log以2为底,真数是4,其实求的就是2的多少次幂是4,所以log2 4 = 2,同理,求log以2为底,真数是16,就是求2的多少次幂是16,log2 16 = 4,当然,如果是log4 2,那就是说2是4的二分之一次幂,所以log4 2 =1/2,log16 2=1/4.
文美榕3505如何用STATA求相关系数并且还有显著性检验 -
郎庄琰18388004314 ______ regress 这个命令就可以啊,然后执行完后用test检验就可以了. 本文来自: 人大经济论坛 Stata专版
文美榕3505用stata做logit回归,Log likelihood 值为多少比较正常?说明什么问题?对300多个数据做回归,十几个变量,log likelihood是负的一百多 -
郎庄琰18388004314 ______[答案] -100多可以的,如果significant的话问题不大
文美榕3505如何用stata进行面板数据的回归分析 -
郎庄琰18388004314 ______ 面板数据用xi即可 熟练掌握 我可以代分析
文美榕3505stata怎么弄线性回归 当应变量为log时 -
郎庄琰18388004314 ______ 所谓回归分析法,是在掌握大量观察数据的基础上,利用数理统计方法建立因变量与自变量之间的回归关系函数表达式(称回归方程式).回归分析中,当研究的因果关系只涉及因变量和一个自变量时,叫做一元回归分析;当研究的因果关系涉及因变量和两个或两个以上自变量时,叫做多元回归分析.此外,回归分析中,又依据描述自变量与因变量之间因果关系的函数表达式是线性的还是非线性的,分为线性回归分析和非线性回归分析.通常线性回归分析法是最基本的分析方法,遇到非线性回归问题可以借助数学手段化为线性回归问题处理.
文美榕3505Stata问题求助:如何按类型填充相同的数值 -
郎庄琰18388004314 ______ 先提取六位数的前两位数字,可以用subinstr函数,并命名为class 接着对class变量对频率分析,并生成新的分类变量type 命令为tab class,gen(type)
文美榕3505在stata中如何实现系数约束检验 -
郎庄琰18388004314 ______ t值和p>t都是测你的系数是否significant用的,stata会自动对所有系数做t检验,比方说你的回归结果里除了xexper以外都是有效的,xexper是无效的,因为p>t这栏里大于0.1了.95%那个就是95%的置信区间
文美榕3505求对数运算的一些公式 -
郎庄琰18388004314 ______[答案] logaM+logaN=logMN logaN^n=nlogaN logaM-logaN=logaM/N