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

// Graph - Slide twoway                                                                       ///   (scatteri 6 0 "Kinds of graphs Stata can make", mlabsize(huge) pstyle(p1)) ///   (scatteri 5 0 "Schemes for controlling look of graphs", mlabsize(huge)     ///       pstyle(p1))                                                            ///   (scatteri 4 0 "Structure of graph commands", mlabsize(huge) pstyle(p1))    ///   (scatteri 3 0 "graph options", mlabsize(huge) pstyle(p1))                  ///   (scatteri 2 0 "Building complex graphs", mlabsize(huge) pstyle(p1))        ///   (scatteri 1 0 "Additional examples", mlabsize(huge) pstyle(p1)),           ///       title("Overview")                                                      ///    xlabel(-.05 1) ylabel(-1 7)                                               ///    xscale(off) yscale(off)                                                   ///    legend(off)                                                               ///    note("Stata graphics are fun and easy" "Stata graphics are powerful",     ///          ring(0)) 复制代码

1. twoway function y=-3*(1-(abs(x)/2)^0.5)^0.5, range(-4 4) || ///
   
2.         function y= (1-(abs(x)-1)^2)^0.5, range(-4 4) yscale(range(-3 1)) ///
   
3.         graphregion(fcolor(pink)) title("七夕" , size(vhuge)  ///
   
4.         margin(bottom)) legend(off) yscale(off) xscale(off)

复制代码

1. // histogram
   
2. 
   
3. sysuse auto, clear
   
4. 
   
5. twoway histogram mpg if !foreign,                              ///
   
6. start(10) width(2) freq bfcolor(none) blcolor(pink)            ///
   
7. || histogram mpg if foreign, freq start(10) width(2) barw(1.8) ///
   
8. bfcolor(none) blcolor(blue)                                    ///
   
9. legend(order(1 "Domestic" 2 "Foreign") pos(2) ring(0) col(1))

复制代码

1. // Text that varies in size due to a calculated factor
   
2.   
   
3. clear programs
   
4.        
   
5.         *--------------- begin example ---------------------
   
6.         program s
   
7.        
   
8.         sysuse auto, clear
   
9.         tempname hdle
   
10.         tempfile info
    
11.         postfile `hdle' str20  variable missings using `info'
    
12.        
    
13.         qui ds, has(type numeric)
    
14.        
    
15.         qui foreach var in `r(varlist)'{
    
16.             /*want more missings
    
17.             to make this interesting */
    
18.             replace `var'=. if runiform()<0.1
    
19.             qui count if `var'==.
    
20.             post `hdle' ("`var'") (r(N))
    
21.           }
    
22.        
    
23.         postclose `hdle'
    
24.        
    
25.         use `info', clear
    
26.         gen number=_n
    
27.         list, noobs
    
28.         gen x = 1
    
29.        
    
30.         sum missings, meanonly
    
31.         gen factor = missings/r(mean)
    
32.         levelsof variable, local(aa)
    
33.        
    
34.         local a=1        
    
35.        
    
36.         foreach i of local aa {
    
37.        
    
38.         local f = factor[`a'
    
39.         local txt `"`txt' text(`a' 1 "`i'", size(*`f'))"'
    
40.         capture `++a'
    
41.         }
    
42.        
    
43.         scatter number x, msymbol(i)  ///
    
44.              xtitle("") xscale(off)   ///
    
45.              ytitle("") yscale(off)   ///
    
46.              plotregion(lstyle(none)) ///
    
47.              `txt'        
    
48.         *----------------- end example -----------------------------
    
49.         end
    
50.        
    
51.         s                 //calls program "s" above

复制代码

1	sysuse uslifeexp.dta, clear
2	graph twoway line le year
3	/* OR */
4	twoway line le year
5	/* OR */
6	line le year

1	sysuse uslifeexp.dta, clear
2	graph twoway line le year
3	/* OR */
4	twoway line le year
5	/* OR */
6	line le year, scheme(s1mono)

1	sysuse uslifeexp.dta, clear
2	line le year, scheme(economist)
3	/* to see list of
4	scheme names:
5	graph query, schemes
6	to change default scheme:
7	set scheme schemename
8	*/

1	*Multiple Dependent Variables
2
3	sysuse uslifeexp.dta, clear
4	line le_wmale le_wfemale le_bmale le_bfemale year

1	*Adding Text
2	sysuse uslifeexp.dta, clear
3	line le_wmale le_wfemale le_bmale le_bfemale year  ///
4	, text(32 1920 "{bf:1918} {it:Influenza} Pandemic", place(3))

1	// Shade Part of a Graph
2
3	sysuse auto, clear
4
5	summarize price, mean
6	local mean = r(mean)
7	kdensity price, gen(x h)
8
9	line h x, || area h x if  x < `mean'