楼主: 匿名
20369 212

[Stata] [Stata] 如何添加直角坐标系及极坐标系及SDAS_绘图实例 [推广有奖]

181楼
niuniuyiwan 在职认证  发表于 2015-10-15 20:49:55
  1.      
  2.                
  3.                                 //命令安装
  4.                                 ssc install spgrid
  5.                                 ssc install spkde
  6.                                 ssc install spmap
  7.                                 
  8.                                 //数据集
  9.                                 net get spkde, from(http://fmwww.bc.edu/RePEc/bocode/s)
  10.                                 net get spmap, from(http://fmwww.bc.edu/RePEc/bocode/s)
  11.                                 net get spgrid, from(http://fmwww.bc.edu/RePEc/bocode/s)
  12.                                 
  13.                                 //修改变量值单位
  14.                                 use "italy-outlinecoordinates",clear
  15.                                 replace _X=_X/1000
  16.                                 replace _Y=_Y/1000
  17.                                 save,replace
  18.                                 
  19.                                 //略去(kilometers)
  20.                                 spgrid using "italy-outlinecoordinates",   ///
  21.                                 resolution(w10)  ///
  22.                                 cells("GridCells.dta")      ///
  23.                                 points("GridPoints.dta")   ///
  24.                                 replace compress dots
  25.                                 
  26.                                 use "Italy-DataPoints.dta", clear
  27.                                 spkde dcvd95 pop95 using "GridPoints.dta",   ///
  28.                                 xcoord(xcoord) ycoord(ycoord)              ///
  29.                                 bandwidth(fbw) fbw(100) dots               ///
  30.                                 saving("Kde.dta", replace)
复制代码

Graph.png

data.JPG


182楼
niuniuyiwan 在职认证  发表于 2015-10-21 19:31:51
  1.                
  2. //命令安装
  3. ssc install spgrid
  4. ssc install spkde
  5. ssc install spmap

  7. //数据集(安装在当前路径下)
  8. net get spkde, from(http://fmwww.bc.edu/RePEc/bocode/s)
  9. net get spmap, from(http://fmwww.bc.edu/RePEc/bocode/s)
  10. net get spgrid, from(http://fmwww.bc.edu/RePEc/bocode/s)

  12. //(对原例数据集中变量值单位大小重新设定，以加快运行效率)
  13. use "italy-outlinecoordinates",clear
  14. replace _X=_X/1000
  15. replace _Y=_Y/1000
  16. save,replace
复制代码
  1. //例1
  2. spgrid, shape(square) resolution(w10) xrange(0 500) yrange(0 200) ///
  3. verbose replace cells("Rectangle-GridCells(Square).dta") ///
  4. points("Rectangle-GridPoints(Square).dta")
  5. use "Rectangle-GridPoints(Square).dta", clear
  6. spmap using "Rectangle-GridCells(Square).dta", id(spgrid_id)
复制代码
Graph.png

183楼
niuniuyiwan 在职认证  发表于 2015-10-21 19:33:13
  1. //例2
  2. spgrid, shape(hexagonal) resolution(w10) xrange(0 500) yrange(0 200) ///
  3. verbose replace cells("Rectangle-GridCells(Hexagonal).dta") ///
  4. points("Rectangle-GridPoints(Hexagonal).dta")
  5. use "Rectangle-GridPoints(Hexagonal).dta", clear
  6. spmap using "Rectangle-GridCells(Hexagonal).dta", id(spgrid_id)
复制代码
Graph.png

184楼
niuniuyiwan 在职认证  发表于 2015-10-21 19:36:39
  1. spgrid using "Italy-OutlineCoordinates.dta", ///
  2. shape(hexagonal) resolution(w10) ///
  3. verbose replace ///
  4. cells("Italy-GridCells(Hexagonal).dta") ///
  5. points("Italy-GridPoints(Hexagonal).dta")

  7. use "Italy-GridPoints(Hexagonal).dta", clear

  9. spmap using "Italy-GridCells(Hexagonal).dta", ///
  10. id(spgrid_id) ///
  11. poly(data("Italy-OutlineCoordinates.dta") ///
  12. ocolor(red) osize(thick))
复制代码

Graph.png



185楼
niuniuyiwan 在职认证  发表于 2015-10-21 19:38:30
  1. spgrid using "Italy-OutlineCoordinates.dta", ///
  2. shape(hexagonal) resolution(w10) ///
  3. verbose replace compress ///
  4. cells("Italy-GridCells(HexValid).dta") ///
  5. points("Italy-GridPoints(HexValid).dta")
  6. use "Italy-GridPoints(HexValid).dta", clear
  7. spmap using "Italy-GridCells(HexValid).dta", ///
  8. id(spgrid_id) ///
  9. poly(data("Italy-OutlineCoordinates.dta") ///
  10. ocolor(red) osize(medium))
复制代码

Graph.png



186楼
niuniuyiwan 在职认证  发表于 2015-10-21 19:39:56
  1. spgrid using "Italy-OutlineCoordinates.dta", ///
  2. shape(hexagonal) resolution(w10) ///
  3. mapexclude("Italy-Exclude.dta") ///
  4. verbose replace compress ///
  5. cells("Italy2-GridCells(HexValid).dta") ///
  6. points("Italy2-GridPoints(HexValid).dta")

  8. use "Italy2-GridPoints(HexValid).dta", clear

  10. spmap using "Italy2-GridCells(HexValid).dta", ///
  11. id(spgrid_id) ///
  12. poly(data("Italy-OutlineCoordinates.dta") ///
  13. ocolor(red) osize(medium))
复制代码

Graph.png



187楼
niuniuyiwan 在职认证  发表于 2015-10-21 19:41:22
  1.    
  2.      use "Italy-DataPoints.dta", clear
  3.      spkde using "GridPoints.dta",     ///
  4.         xcoord(xcoord) ycoord(ycoord)   ///
  5.         bandwidth(fbw) fbw(100) dots    ///
  6.         edgecorrect                     ///
  7.         saving("Kde.dta", replace)
  8.      use "Kde.dta", clear
  9.      spmap lambda using "GridCells.dta",     ///
  10.          id(spgrid_id) clnum(20)              ///
  11.          fcolor(Rainbow) ocolor(none ..)      ///
  12.          legend(off)                          ///
  13.          point(data("Italy-DataPoints.dta")   ///
  14.          x(xcoord) y(ycoord))
复制代码

Graph.png



188楼
niuniuyiwan 在职认证  发表于 2015-10-21 19:43:48
  1.      
  2.       use "Italy-DataPoints.dta", clear
  3.      spkde using "GridPoints.dta",     ///
  4.         xcoord(xcoord) ycoord(ycoord)   ///
  5.         kernel(normal)                  ///
  6.         bandwidth(fbw) fbw(ad5) dots    ///
  7.         saving("Kde.dta", replace)
  8.      use "Kde.dta", clear
  9.      spmap lambda using "GridCells.dta",     ///
  10.          id(spgrid_id) clnum(20)              ///
  11.          fcolor(Rainbow) ocolor(none ..)      ///
  12.          legend(off)                          ///
  13.          point(data("Italy-DataPoints.dta")   ///
  14.          x(xcoord) y(ycoord))
复制代码

Graph.png


189楼
niuniuyiwan 在职认证  发表于 2015-10-22 07:11:46
  1. set scheme s1color
  2. u cumbrian_cirques, clear
  3. scatter width length , name(g1_1, replace)
  4. scatter width length , ms(Oh) name(g1_2, replace)
  5. graph combine g1_1 g1_2, saving(g1, asis replace)
复制代码

Graph.png



190楼
niuniuyiwan 在职认证  发表于 2015-10-22 08:12:37
  1. set scheme s1color
  2. u cumbrian_cirques, clear


  5. scatter maxgrad grade, ms(Oh) yla(40(10)80, format(%4.0g)) xla(, valuelabel) xsc(r(0.8 5.2)) name(g2_1, replace)
  6. scatter maxgrad grade, ms(Oh) yla(40(10)80, format(%4.0g)) xla(, valuelabel) xsc(r(0.8 5.2)) jitter(1) name(g2_2, replace)
  7. graph combine g2_1 g2_2, saving(g2, replace asis)
复制代码

Graph.png



返回列表
1 ...141516171819202122下一页
高级模式
B Color Image Link Quote Code Smilies
您需要登录后才可以回帖 登录 | 我要注册

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

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

GMT+8, 2025-12-31 22:19