楼主: Alfred_G
3299 2

[Stata] 自用Stata绘图（1） [推广有奖]

已卖：1287份资源

教授

38%

还不是VIP/贵宾

-

0%
威望
0 级
论坛币
21223 个
通用积分
362.9526
学术水平
256 点
热心指数
263 点
信用等级
167 点
经验
94930 点
帖子
905
精华
2
在线时间
1403 小时
注册时间
2013-12-19
最后登录
2025-11-9

初级学术勋章 初级热心勋章 中级学术勋章 中级热心勋章

楼主
Alfred_G 学生认证  发表于 2015-10-17 10:11:30 |AI写论文

+2 论坛币
k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群
赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

自用Stata绘图（1）

1.带标注的散点图：
  1. sysuse lifeexp, clear
  2. keep if region==2 | region==3
  3. replace gnppc = gnppc / 1000
  4. label var gnppc "GNP per capita (thousands of dollars)"
  5. gen lgnp = log(gnp)
  6. qui reg lexp lgnp
  7. predict hat
  8. label var hat "Linear prediction"
  9. replace country = "Trinidad" if country=="Trinidad and Tobago"
  10. replace country = "Para" if country == "Paraguay"
  11. gen pos = 3
  12. replace pos = 9 if lexp > hat
  13. replace pos = 3 if country == "Colombia"
  14. replace pos = 3 if country == "Para"
  15. replace pos = 3 if country == "Trinidad"
  16. replace pos = 9 if country == "United States"
  17. #delimit ;
  18. twoway (scatter lexp gnppc, mlabel(country) mlabv(pos))
  19.         (line hat gnppc, sort)
  20.         , xsca(log) xlabel(.5 5 10 15 20 25 30, grid)
  21.         legend(off)
  22.         title("Life expectancy vs. GNP per capita")
  23.         subtitle("North, Central, and South America")
  24.         note("Data source: World bank, 1998")
  25.         ytitle("Life expectancy at birth (years)") ;
  26.         #delimit cr
复制代码

lifeexp.png

2.针式图

  1. sysuse lifeexp, clear
  2. keep if region==3
  3. gen lngnp = ln(gnppc)
  4. quietly regress le lngnp
  5. predict r, resid
  6. twoway dropline r gnp, ///
  7.   yline(0, lstyle(foreground)) mlabel(country) mlabpos(9) ///
  8.   ylab(-6(1)6) ///
  9.   subtitle("Regression of life expectancy on ln(gnp)" "Residuals:" " ", pos(11)) ///
  10.   note("Residuals in years; positive values indicate" "longer than predicted life expectancy")
复制代码

Graph.png

3.直方图

  1. sysuse sp500, clear
  2. #delimit ;

  4. histogram volume, freq
  5.   xaxis(1 2)
  6.   ylabel(0(10)60, grid)
  7.   xlabel(12321 "mean"
  8.         9735 "-1 s.d."
  9.        14907 "+1 s.d."
  10.         7149 "-2 s.d."
  11.        17493 "+2 s.d."
  12.        20078 "+3 s.d."
  13.        22664 "+4 s.d.", axis(2) grid gmax)
  14.   xtitle("", axis(2))
  15.   subtitle("S&P 500, January 2001 - December 2001")
  16.   note("Source:  Yahoo!Finance and Commodity Systems, Inc.")
  17.   ;
  18. #delimit cr
复制代码

Graph.png

4.带拟合分布曲线的直方图

  1. sysuse sp500, clear
  2. #delimit ;

  4. histogram volume, freq normal
  5. xaxis(1 2)
  6. ylabel(0(10)60, grid)
  7. xlabel(12321 "mean"
  8.       9735 "-1 s.d."
  9.      14907 "+1 s.d."
  10.       7149 "-2 s.d."
  11.      17493 "+2 s.d."
  12.      20078 "+3 s.d."
  13.      22664 "+4 s.d."
  14. , axis(2) grid gmax)
  15. xtitle("", axis(2))
  16. subtitle("S&P 500, January 2001 - December 2001")
  17. note("Source:  Yahoo!Finance and Commodity Systems, Inc.")
  18. ;
  19. #delimit cr
复制代码

Graph.png

5.区间图

  1. sysuse sp500, clear
  2. gen month = month(date)
  3. sort month
  4. by month: egen lo = min(volume)
  5. by month: egen hi = max(volume)
  6. format lo hi %10.0gc
  7. by month: keep if _n==_N
  8. #delimit ;

  10. twoway rcap lo hi month,
  11.   xlabel(1 "J"  2 "F"  3 "M"  4 "A"  5 "M"  6 "J"
  12.                  7 "J"  8 "A"  9 "S" 10 "O" 11 "N" 12 "D")
  13.   xtitle("Month of 2001")
  14.   ytitle("High and Low Volume")
  15.   yaxis(1 2) ylabel(12321 "12,321 (mean)", axis(2) angle(0))
  16.   ytitle("", axis(2))
  17.   yline(12321, lstyle(foreground))
  18.   msize(*2)
  19.   title("Volume of the S&P 500", margin(b+2.5))
  20.   note("Source:  Yahoo!Finance and Commodity Systems Inc.")
  21. ;
  22. #delimit cr
复制代码

区间图.png

6.区域图：
  1. sysuse gnp96, clear
  2. #delimit ;

  4. twoway area d.gnp96 date, xlabel(36(8)164, angle(90))
  5.   ylabel(-100(50)200, angle(0))
  6.   ytitle("Billions of 1996 Dollars")
  7.   xtitle("")
  8.   subtitle("Change in U.S. GNP", position(11))
  9.   note("Source: U.S. Department of Commerce, Bureau of Economic Analysis")
  10.   ;
  11. #delimit cr
复制代码

Graph.png

5.根据函数制作的曲线图
  1. #d;
  2. twoway function y=exp(-x/6)*sin(x), range(0 12.57)
  3.   yline(0, lstyle(foreground))
  4.   xlabel(0 3.14 "pi" 6.28 "2 pi" 9.42 "3 pi" 12.57 "4 pi")
  5.   plotregion(style(none))
  6.   xscale(noline)
  7. ;
  8. #d cr
复制代码

Graph.png



6.正态曲线下的面积图

  1. #delimit ;

  3. twoway
  4.         function y=normalden(x), range(-4 -1.96) bcolor(gs12) recast(area) ||
  5.         function y=normalden(x), range(1.96 4) bcolor(gs12) recast(area) ||
  6.         function y=normalden(x), range(-4 4) clstyle(foreground) ||,
  7.         plotregion(style(none))
  8.         yscale(off) xscale(noline)
  9.         legend(off)
  10.         xlabel(-4 "-4 sd" -3 "-3 sd" -2 "-2 sd" -1 "-1 sd" 0 "mean" 1 "1 sd"   2 "2 sd"   3 "3 sd"   4 "4 sd", grid gmin gmax)
  11.         xtitle("")
  12.         ;
  13. #delimit cr
复制代码

Graph.png




二维码

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0
关键词：stata绘图 Stata tata EXPECTANCY World Bank prediction thousands replace country dollars

相关帖子

已有 2 人评分经验 论坛币 学术水平 热心指数 信用等级 收起 理由
niuniuyiwan + 100 + 100 + 5 + 5 + 5 感谢支持代码库建设！
xddlovejiao1314 + 100 + 100 + 3 + 3 + 3 精彩帖子

总评分: 经验 + 200  论坛币 + 200  学术水平 + 8  热心指数 + 8  信用等级 + 8   查看全部评分

本帖被以下文库推荐

知识和能力是在交流中增加的，平淡做人，认真做事，不功利，不急躁~

沙发
xddlovejiao1314 学生认证  发表于 2015-10-17 11:22:02
谢谢分享。
已有 2 人评分论坛币 学术水平 热心指数 信用等级 收起 理由
于役 + 1 + 1 + 1 精彩帖子
niuniuyiwan + 10 + 1 热心帮助其他会员

总评分: 论坛币 + 10  学术水平 + 1  热心指数 + 2  信用等级 + 1   查看全部评分

藤椅
于役 发表于 2019-1-20 14:26:46
复式散点图要怎么做啊 求指教
返回列表
高级模式
B Color Image Link Quote Code Smilies
您需要登录后才可以回帖 登录 | 我要注册

本版微信群
加好友,备注jltj
拉您入交流群

京ICP备16021002号-2 京B2-20170662号 京公网安备 11010802022788号 论坛法律顾问：王进律师 知识产权保护声明   免责及隐私声明

GMT+8, 2025-12-6 01:32