第 60 章 模型的边际效应 |
您所在的位置:网站首页 › 边际效应怎么求 › 第 60 章 模型的边际效应 |
第 60 章 模型的边际效应
本章介绍模型的边际效应,主要围绕marginaleffects宏包,本章的内容也是来源该宏包的说明文档。 60.1 边际效应边际效应,测量的是某一个预测因子变化一个单位与伴随的响应变量的变化之间的关联。用数学语言表示,就是回归方程对x求偏导。 假定我们建立的回归方程是一个二次函数, \[ y = -x^2 \] 那么对x的偏导数 \[ \frac{\partial y}{\partial x} = -2x \] 可以看到,此时的边际效应就是曲线的斜率 当\(x0\),斜率为负,x增加y也减少,在这个位置上边际效应为负 60.2 marginaleffects function最简单的线性模型,每个因子的边际效应就是预测因子的系数,与因子的取值无关。但是复杂点模型,因子边际效应不仅仅与因子的取值有关,而且还与其它因子的值也有关。 我们下面用企鹅数据来说明。 library(tidyverse) library(palmerpenguins) library(marginaleffects)我们先构建一个二元变量fat_penguin(是否为胖子), 1表示是,0表示不是。并建立logitisc回归模型 dat % drop_na() %>% mutate( fat_penguin = if_else(body_mass_g > median(body_mass_g), 1, 0) ) mod % count(term, contrast) ## ## Term Contrast ## bill_length_mm dY/dX ## flipper_length_mm dY/dX ## species Chinstrap - Adelie ## species Gentoo - Adelie ## ## Columns: term, contrast, n边际效应对连续变量非常适合。但离散变量的边际效应,不太好理解,因此采用对照方法,具体为,以某一层级为基线,那么从基线切换到其它层级,伴随着响应变量的变化,就为离散变量的边际效应。 60.3 平均边际效应 summary(mfx) ## ## Term Contrast Estimate Std. Error z ## bill_length_mm mean(dY/dX) 0.0268 0.00591 4.538 ## flipper_length_mm mean(dY/dX) 0.0108 0.00241 4.506 ## species mean(Chinstrap) - mean(Adelie) -0.4076 0.05786 -7.044 ## species mean(Gentoo) - mean(Adelie) 0.0658 0.10956 0.601 ## Pr(>|z|) 2.5 % 97.5 % ## |
CopyRight 2018-2019 办公设备维修网 版权所有 豫ICP备15022753号-3 |