如何画（倒）U 型图？

黄老师，我想问一下，为什么我安装marginscontplot，安装不上去？就一直显示
ssc install: "marginscontplot" not found at SSC, type search marginscontplot
(To find all packages at SSC that start with m, type ssc describe m)

小丸子ddx 发表于 2018-11-27 16:14
黄老师，我想问一下，为什么我安装marginscontplot，安装不上去？就一直显示
ssc install: "marginscontpl ...

请 search marginscontplot 后安装。

黃河泉 发表于 2018-11-27 17:12
请 search marginscontplot 后安装。

谢谢黄老师，我已经安装上了

放一个我在 Stata 寒假现场班中的 Stata 范例，供参考 ：https://gitee.com/arlionn/stata_training/blob/master/README.md

1. *---------------------------
   
2. *-4.5.2  特殊交乘项: 平方项  (Wooldridge, Sec6.3, p.189-194)
   
3.   
   
4.     *-二次函数回顾
   
5.         *
   
6.         *  y = a*x^2 + b*x + c
   
7.         *
   
8.         *  转折点: tpoint = -b/2a
   
9.         *
   
10.         *  边际效果: dy/dx = b + 2a*x
    
11.         
    
12.         
    
13.     *-转折点(The turnning point)
    
14.       
    
15.           sysuse "nlsw88.dta", clear
    
16.           global xx "hours age tenure married south i.race"
    
17.           gen ttl_exp2 = ttl_exp^2
    
18.           reg wage ttl_exp ttl_exp2 $xx
    
19.       *------------------------------------------------
    
20.       *     wage |   Coef.   Std. Err.      t    P>|t|
    
21.       *----------+-------------------------------------
    
22.       *  ttl_exp |   0.492      0.110     4.49   0.000
    
23.       * ttl_exp2 |  -0.009      0.005    -1.99   0.046
    
24.       *                  (ommited)                     
    
25.       *------------------------------------------------
    
26.          
    
27.           keep if e(sample)==1  //删除缺漏值
    
28.         
    
29.           predict wage_hat if e(sample)
    
30.         
    
31.           twoway (scatter wage ttl_exp) (qfit wage_hat ttl_exp )
    
32.   
    
33.     *-转折点的确定
    
34.           global tp: dis %3.1f -_b[ttl_exp]/(2*_b[ttl_exp2])
    
35.           dis "Turn point = " $tp
    
36.         
    
37.         *---------------------------图示-----------------------------
    
38.           local b0 = _b[_cons]
    
39.           local b1 = _b[ttl_exp]
    
40.           local b2 = _b[ttl_exp2]
    
41.           sum ttl_exp
    
42.           local min:  dis %3.1f r(min)
    
43.           local max:  dis %3.1f r(max)
    
44.       #delimit        
    
45.           twoway ( function y = `b2'*x^2 + `b1'*x + `b0',
    
46.                    range(`min' `max') )
    
47.                  ( function y = `b2'*x^2 + `b1'*x + `b0',
    
48.                    range(`max' 40) lp(dash) )
    
49.                          ,
    
50.                          ytitle("wage") xtitle("ttl_exp")
    
51.                          xline(`min',lc(red) lw(*1.5))
    
52.                          xline($tp,  lp(dash) lc(green))
    
53.                          xline(`max',lc(red) lw(*1.5))
    
54.                          xlabel(`min' "Min (`min')" 5(5)25
    
55.                                  $tp "TP ($tp)" `max' "Max (`max')", angle(50))
    
56.                          legend(off)
    
57.                          caption("The turn point is: $tp")
    
58.                          ;
    
59.           #delimit cr
    
60.           graph export "$Out\fig_Ushape.png", replace
    
61.     *---------------------------图示-----------------------------
    
62.         
    
63.         *-转折点右侧的观察值占比：要注意其实际意义
    
64.           keep if e(sample)
    
65.           count if ttl_exp>$tp
    
66.           sum ttl_exp
    
67.           *-结论: 此时的转折点不具有实际意义                                                  
    
68. 
    
69. 
    
70.         *-使用 margin 命令分析边际效应 (推荐使用)
    
71.          
    
72.           reg wage c.ttl_exp##c.ttl_exp $xx
    
73.                
    
74.           sum ttl_exp
    
75.           global min:  dis %4.1f r(min)
    
76.           global max:  dis %4.1f r(max)
    
77.           margins, dydx(ttl_exp) at(ttl_exp=($min 1(3)28 $max))
    
78.       *-----------------------------------------------------------
    
79.       *       |         Delta-method
    
80.       *       |  dy/dx     Std. Err.    t    P>|t|      [95% CI]
    
81.       *-------+---------------------------------------------------
    
82.       *ttl_exp|      
    
83.       *   _at |
    
84.       *    1  |  0.490      0.109     4.51   0.000    0.277  0.703
    
85.       *    2  |  0.474      0.101     4.69   0.000    0.276  0.672
    
86.       *    3  |  0.419      0.076     5.56   0.000    0.271  0.568
    
87.       *    4  |  0.365      0.052     7.03   0.000    0.263  0.467
    
88.       *    5  |  0.311      0.034     9.04   0.000    0.243  0.378
    
89.       *    6  |  0.257      0.034     7.56   0.000    0.190  0.323
    
90.       *    7  |  0.202      0.051     3.97   0.000    0.102  0.302
    
91.       *    8  |  0.148      0.074     1.99   0.047    0.002  0.294
    
92.       *    9  |  0.094      0.100     0.94   0.347   -0.102  0.289
    
93.       *   10  |  0.040      0.126     0.31   0.754   -0.207  0.286
    
94.       *   11  | -0.015      0.152    -0.10   0.923   -0.314  0.284
    
95.       *   12  | -0.031      0.160    -0.19   0.847   -0.346  0.284
    
96.       *-----------------------------------------------------------
    
97.          
    
98.           marginsplot, xlabel(, format(%3.1f) angle(60)) //绘图,附带 95% 置信区间
    
99.           graph export "$Out\fig_Wage_marginplot_A.png", replace
    
100.          
     
101.           marginsplot, addplot(function y=0, range($min $max)  ///
     
102.                            lcolor(red) lpattern(dash) legend(off)) ///
     
103.                                    xlabel(, format(%3.1f) angle(60)) ///
     
104.                                    ylabel(, format(%2.1f) angle(0))
     
105.           graph export "$Out\fig_Wage_marginplot_B.png", replace  

复制代码

fig_Ushape.png






Stata 寒假现场班    https://gitee.com/arlionn/stata_training/blob/master/README.md

arlionn 发表于 2018-12-2 09:51
放一个我在 Stata 寒假现场班中的 Stata 范例，供参考 ：https://gitee.com/arlionn/stata_training/blob/m ...

1. 感谢连老师提供的精彩例子。2. 我已经将其设定为精彩回复，这样一来，其他论坛朋友就可以更清楚并容易注意到此精彩例子。3. 不知您有无注意到，好像"大部分"的情况与例子，都是倒 U 而非正 U 型关系？

其实，我那个例子是想说明：多数情况下我们只看到了 U 型或 倒 U 型的半边。也就是说，在转折点左侧或右侧的样本数很少，此时下结论时要谨慎。

arlionn 发表于 2018-12-5 11:15
其实，我那个例子是想说明：多数情况下我们只看到了 U 型或 倒 U 型的半边。也就是说，在转折点左侧或右侧的 ...

也没错！

黄老师，我用了你给的画U型关系图的命令，
xi:reg ca3 sub c.sub#c.sub down bm lev size roa indr growth ci lnky top1 lnage mng com ex gdprate ppi i.year i.indu,robust
marginscontplot sub, ci
为什么stata一直提示invalid numlist has too many elements？

zyj592991367 发表于 2019-7-5 16:26
黄老师，我用了你给的画U型关系图的命令，
xi:reg ca3 sub c.sub#c.sub down bm lev size roa indr growth ...

看不出问题！

感谢分享，想请问一下这个图的含义是什么？是pcgdp的边际作用么？如果是的话，跟后面您给出的marginsplot命令中的结果图有什么区别吖？谢谢您！