倒U型
最近在学习倒U型,虽然自己没有用上,但是留着作为以后备用,放假了可以时常更新之前学习的stata代码。
这里放一下检验倒U型的代码,主要是从连老师和B站视频中学到的,后期会放上相关链接。
sysuse "nlsw88.dta", clear
gen ttl_exp2 = ttl_exp * ttl_exp
reg wage ttl_exp ttl_exp2
ssc install utest
utest ttl_exp ttl_exp2, fieller
sum ttl_exp
///这是看utest结论
global tp: dis %3.1f -_b[ttl_exp]/(2*_b[ttl_exp2])
dis "Turn point = " $tp
//找到极值点
local b0 = _b[_cons]
local b1 = _b[ttl_exp]
local b2 = _b[ttl_exp2]
sum ttl_exp
local min: dis %3.1f r(min)
local max: dis %3.1f r(max)
#delimit
twoway( function y = `b2'*x^2 + `b1'*x + `b0',
range(`min' `max') )
( function y = `b2'*x^2 + `b1'*x + `b0',
range(`max' 40) lp(dash) )
,
ytitle("wage") xtitle("ttl_exp")
xline(`min',lc(red) lw(*1.5))
xline($tp, lp(dash) lc(green))
xline(`max',lc(red) lw(*1.5))
xlabel(`min' "Min (`min')" 5(5)25
`max' "Max (`max')", angle(50))
legend(off)
;
#delimit cr
///画出倒U型形状
///有一个疑问是,连老师专栏中写道倒U型的slope显示是负向,但是实际上,slope显示为正画出的图像是倒U型。
///附上三点信息:
1.一般一次项为正,二次项为负为倒U型可以画图出来;
2.在学习时候安装好了utest后可以help一下,查看代码和使用方法;
3.两个变量多为连续变量。
由于在处理时候,每个人疑惑不一样,附上学习的链接。
##链接如下
链接:link.https://www.bilibili.com/video/BV1j7411u7HT
https://www.lianxh.cn/news/c36ef8e23a704.html