密度函数曲线的三维图做法

6717

收藏 2013-07-12

webuse pennxrate
kdensity ppp if year==1970
这是画1970年所有个体ppp指标的密度函数曲线，我想画1970-2003年所有个体ppp指标的密度函数曲线，最好有年份这一坐标轴，估计需要做三维图，这样一张图可以看出随时间的变化不同个体的分布状态，请问该怎么做呢？

扫码加我拉你入群

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

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

全部回复

h3327156

2013-7-12 17:31:21

***Stata的三维没有比较好看，年份太细，变化不大的，年度应当要跨比较多年才有变化
***建议参考手册绘多条的kdensity
***原本ppp图没有比较好看，这边改用ln(ppp)
webuse pennxrate
g lnppp=ln(ppp)
kdensity lnppp, nogr g(x s)
forv i=1970(10)2000{
  kdensity lnppp if year==`i', nogr g(s`i') at(x)
  la var s`i' "`i'Year"
  loc m `m' s`i'
}
line `m' x, sort ytitle(Density)

扫码加我拉你入群

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

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

h3327156

2013-7-12 21:19:45

***利用scat3绘画看看***如果楼主还是坚持三维
webuse pennxrate
g lnppp=ln(ppp)
kdensity lnppp, nogr g(x s)
forv i=1970(10)2000{
  kdensity lnppp if year==`i', nogr g(s`i') at(x)
  la var s`i' "`i'Year"
  loc m `m' s`i'
}
*line `m' x, sort ytitle(Density)

g g=1970 in 1/50
replace g=1980 in 51/100
replace g=1990 in 101/150
replace g=2000 in 151/200

forv a=1/150{
replace x=x[`a'] in `=50+`a''
}

forv b=1/50{
replace s1970=s1980[`b'] in `=50+`b''
replace s1970=s1990[`b'] in `=100+`b''
replace s1970=s1990[`b'] in `=150+`b''
}
g ss=s1970

scat3 x g ss
*g是时间，图形很难看，没办法啦，kdensity的default是50点
*如果您可以改成1000点，就会很漂亮的图。