楼主: 匿名
20361 212

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

31楼
niuniuyiwan 在职认证  发表于 2015-8-27 07:53:48
  1. //make up some data
  2. set seed 1
  3. clear
  4. set obs 50
  5. generate a=runiform()*1000
  6. generate b=runiform()*10
  7. generate s=runiform()<.5

  9. label define sexlab 1 `"{fontface "Female and male symbols":M }"' 0 ///
  10.       `"{fontface "Female and male symbols":F }"'
  11. label values s sexlab

  13. scatter a b , ms(i)  mlab(s) mlabpos(c) ///
  14. mlabsize(*2) title( "Example" `"{fontface "Female and male symbols":M F}"')
复制代码

Graph.png


以上为未加装字体时的图像 FEMALE AND MALE SYMBOLS FONT.zip (4.56 KB) 本附件包括：

  • FMSYMB.TTF
  • FMSYMB.WRI
  • LICENSE.TXT

解压后双击>-----安装-----即可

字体下载地址：http://www.fonts4free.net/female-and-male-symbols-font.html#

代码同上


Graph.png



32楼
niuniuyiwan 在职认证  发表于 2015-8-27 07:55:09
  1. // Fourier plot

  3. clear
  4. set obs 8
  5. gen x=(_n-1)*(2*_pi)/8
  6. gen y=sin(x)
  7. mata: h=st_data(.,"y")
  8. mata: fft(h)

  10. # delimit ;
  11. twoway scatter y x ||       
  12. function y=1/8*(4*sin(1*x)-4*sin( 7*x)), range(x) ||       
  13. function y=1/8*(4*sin(1*x)-4*sin(-1*x)), range(x) ||,
  14. legend(cols(1) label(1 "sampled values")       
  15.   label(2 "sum of high frequency components")       
  16.   label(3 "sum of low frequency components"))       
  17. title("Reconstructed signal for y=sin(x) sampled at 8 points");       
  18. # delimit cr       
  19.                        
复制代码

Graph.png



33楼
niuniuyiwan 在职认证  发表于 2015-8-27 07:58:00
  1. *修改后，添加#delimit执行无误
  2. // Circle using function
  3. #delimit ;
  4. twoway (function y = sqrt(1 -(x)^2),
  5. range(-1 1) lwidth(thick) lcolor(red))
  6. (function y = -sqrt(1 - (x)^2),
  7. range(-1 1) lwidth(thick) lcolor(red)), aspect(1);
复制代码

Graph.png


34楼
niuniuyiwan 在职认证  发表于 2015-8-27 08:39:36
  1. // normalden using function
  2. clear
  3. set obs 100

  5. generate x1 = rnormal(0,1)
  6. generate x2 = rnormal(0.5,1.3)

  8. summarize x1
  9. local M1 = r(mean)
  10. local SD1 = r(sd)

  12. summarize x2
  13. local M2 = r(mean)
  14. local SD2 = r(sd)

  16. twoway function x1 = normalden(x, `M1', `SD1') , range(-4 4) || ///
  17. function x2 = normalden(x, `M2', `SD2'), range(-4 5)
复制代码

Graph.png



35楼
niuniuyiwan 在职认证  发表于 2015-8-27 12:18:46
  1. //Mandelbrot Set
  2. *Example taken from Stata News
  3. *http://www.stata.com/stata-news/statanews.27.4.pdf
  4. *Please note that on my computer this took 3 hours to run.
  5. //Stata——14 MP,大约半小时
  6. clear all

  8. mata:
  9. real scalar escape(complex scalar Z, C, real scalar R, maxiter)
  10. {
  11. real scalar i
  12. for (i=0; i < maxiter; i++) {
  13.   if (norm(Z) <= R) {
  14.    Z = Z*Z+C
  15. }
  16.   else return(i)
  17.   }
  18.   return(i)
  19. }
  20. end

  22. mata:
  23. mandelbrot_set = J(201*201, 3, .)
  24. cnt = 1
  25. C = -0.8+0.156i
  26. for(i=0; i<=200; i++) {
  27. for(j=0; j<=200; j++) {
  28.   Z = C(-2+i*4/200, -2+j*4/200)
  29.   n = escape(Z, C, 100, 400)
  30.   mandelbrot_set[cnt,.]= (i, j, n)
  31.   cnt++
  32. }
  33. }
  34. end

  36. getmata (x y escape)= mandelbrot_set
  37. quiet summarize escape
  38. local min = r(min)-1
  39. local max = r(max)+1
  40. twoway contour escape y x, ccuts(`min'(1)`max')  ///
  41.   clegend(off) graphr(m(zero)) xscale(off)        ///
  42.   xlab(,nogrid)  yscale(off) ylab(,nogrid) scheme(s2color)
复制代码

Graph.png



36楼
niuniuyiwan 在职认证  发表于 2015-8-27 13:56:20
  1. //2nd Row of axis labels

  3. webuse turksales, clear
  4. line sales t,                              ///
  5.     xtick(119.5(4)159.5, tlength(*1.5))    ///
  6.     xmlabel(120/159,                       ///
  7.         noticks format(%tqq))              ///
  8.     xlabel(121.5(4)157.5,                  ///
  9.         noticks format(%tqCY) labgap(2.5)) ///
  10.     xtitle("")

  12. *(Maarten Buis - Stata list,  Thu, 7 Nov 2013)
复制代码


Graph.png



37楼
niuniuyiwan 在职认证  发表于 2015-8-27 14:14:28
  1. // Gap between bars in bar graph

  3. sysuse auto, clear

  5. graph bar (mean) price turn mpg weight, over(rep78)

  7. set obs 75
  8. replace rep78=2.5 in 75

  10. graph bar (mean) price turn mpg weight, over(rep78,  relabel( 3 " "  ))

  12. *(Nick Cox - Stata list,  Wed, 26 May 2010)
复制代码

Graph.png


38楼
niuniuyiwan 在职认证  发表于 2015-8-27 16:50:58
  1. // Candle Stick plot

  3. clear
  4. input ///
  5. str3 share         str20   time     price
  6.       abc          "1/5/2011:10:34" 10.34
  7.       abc          "1/5/2011:10:35" 10.25
  8.       abc          "1/5/2011:10:36" 10.15
  9.       abc          "1/5/2011:10:37" 10.73
  10.       abc          "1/5/2011:10:39" 10.95
  11.       abc          "1/5/2011:10:41" 10.53
  12.       abc          "1/5/2011:10:45" 10.99
  13.       abc          "1/5/2011:10:49" 11.05
  14.       abc          "2/5/2011:10:34" 9.34
  15.       abc          "2/5/2011:10:35" 9.25
  16.       abc          "2/5/2011:10:36" 9.15
  17.       abc          "2/5/2011:10:37" 9.73
  18.       abc          "2/5/2011:10:39" 9.95
  19.       abc          "2/5/2011:10:41" 9.53
  20.       abc          "2/5/2011:10:45" 9.99
  21.       abc          "2/5/2011:10:49" 9.05
  22.       abc          "3/5/2011:10:34" 8.34
  23.       abc          "3/5/2011:10:35" 8.25
  24.       abc          "3/5/2011:10:36" 8.15
  25.       abc          "3/5/2011:10:37" 8.73
  26.       abc          "3/5/2011:10:39" 8.95
  27.       abc          "3/5/2011:10:41" 8.53
  28.       abc          "3/5/2011:10:45" 8.99
  29.       abc          "3/5/2011:10:49" 12.05
  30.       abc          "4/5/2011:10:34" 12.34
  31.       abc          "4/5/2011:10:35" 12.25
  32.       abc          "4/5/2011:10:36" 12.15
  33.       abc          "4/5/2011:10:37" 12.73
  34.       abc          "4/5/2011:10:39" 12.95
  35.       abc          "4/5/2011:10:41" 12.53
  36.       abc          "4/5/2011:10:45" 12.99
  37.       abc          "4/5/2011:10:49" 17.05
  38. end
  39. l

  41. gen double date=dofc(clock(time,"DMY hm"))
  42. format date %td

  44. collapse  (p1) p1=price (p25) p25=price ///
  45. (p50) p50=price (p75) p75=price         ///
  46. (p99) p99=price, by(date)

  48. gen uav=p75- 1.5*(p25-p75)
  49. gen lav=p75+ 1.5*(p25-p75)
  50.                
  51. gen col=2               
  52. replace col=(p50-p50[_n-1])>0 if _n!=1       
  53.                
  54. twoway (rbar p75 p25 date if col==0, barwidth(.6) color(red))  ///
  55. (rbar p75 p25 date if col==1, barwidth(.6) color(green))       ///
  56. (rbar p75 p25 date if col==2, barwidth(.6) color(black))       ///
  57. (rspike uav lav date if col==2, lwidth(thick) color( black))   ///
  58. (rspike uav lav date  if col==1 , lwidth(thick) color( green)) ///
  59. (rspike uav lav date if col==0 , lwidth(thick) color(red)  ) , ///
  60. legend(off)
复制代码

Graph.png



39楼
niuniuyiwan 在职认证  发表于 2015-8-27 21:45:44
  1. // Adding average line to bar graph

  3. webuse grunfeld, clear   //getting data online
  4. keep if company == 1
  5. keep if year < 1945

  7. generate year2 = year + 0.3
  8. generate year3 = year - 0.3

  10. generate average = (invest + mvalue + kstock) / 3

  12. twoway bar invest year3, barw(0.3) || ///
  13. bar mvalue year, barw(0.3)         || ///
  14. bar kstock year2 , barw(0.3)       || ///
  15. line average year , name(a2, replace)

  17. *(Nick Cox - Stata list,  13/4/11)
复制代码

a2.png



40楼
niuniuyiwan 在职认证  发表于 2015-8-27 23:09:57
  1. // Graph Title - Includes years
  2. clear

  4. inp str10 date ricepr
  5. "01/01/2006"   700
  6. "01/02/2006"   700
  7. "01/03/2006"   900
  8. "01/04/2006"   900
  9. "01/05/2006"   900
  10. "01/06/2006"   900
  11. "01/07/2006"   900
  12. "01/08/2006"   933.33
  13. "01/09/2006"   1000
  14. "01/10/2006"   1000
  15. "01/11/2006"   1000
  16. "01/12/2006"   1000
  17. "01/01/2007"   1050
  18. "01/02/2007"   1087.5
  19. "01/03/2007"   1100
  20. "01/04/2007"   1100
  21. "01/05/2007"   1100
  22. "01/06/2007"   1100
  23. "01/07/2007"   1175
  24. "01/08/2007"   1200
  25. "01/09/2007"   1200
  26. "01/10/2007"   1300
  27. "01/11/2007"   1400
  28. "01/12/2007"   1400
  29. end

  31. gen date2=date(date,"DMY")
  32. gen monthlabel = substr("JFMAMJJASOND", month(date2), 1)

  34. replace monthlabel=monthlabel+ " " + string(year(date2),"%5.0f") ///
  35. if month(date2)==1

  37. labmask date2, values(monthlabel) //user written program, must be installed first
  38. levelsof date2, local(levels)
  39. line ricepr date2, xla(`levels', valuelabel angle(45))

  41. *(Nick Cox - Stata list,  Tue, 04 May 2010)
复制代码

Graph.png



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

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

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

GMT+8, 2025-12-31 20:45