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

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

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

扫码加我 拉你入群

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

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

分享0 收藏2 回帖

关键词：密度函数 三维图 kdensity Density Rate 绑架

***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)

***利用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点，就会很漂亮的图。

h3327156 发表于 2013-7-12 21:19
***利用scat3绘画看看***如果楼主还是坚持三维
webuse pennxrate
g lnppp=ln(ppp)

谢谢！谢谢！虽然结果不怎么好看，不过还是很感激！可能要求助R了，看是否能画出更好的三维图...

学习了

学习了，呵呵