Sử dụng chơng trình Stata nội dung Chơng I: Giới thiệu chung chơng trình Stata Tæ chøc l−u tr÷ d÷ liƯu Stata (Dataset in Stata) 2 Khởi động tho¸t khái Stata (Open and exit) 3 Giao diÖn Stata (Stata interface) Biên làm việc (log file) Error! Bookmark not defined Nhập lu liệu (Use, input and and save) Chơng II: Khai thác liệu CÊu tróc lƯnh Stata (Stata command syntax) 2.To¸n tư vµ hµm sè (Operators and functions) 11 Mô tả liệu (Data reporting) 12 Biên tập sửa ch÷a d÷ liƯu (Data manipulation) 25 QuyÒn sè VHLSS (Weight) 38 Chơng III: Kiểm định giả thiết phân tích håi quy 43 Ước lợng kiểm định giả thiết (Estimation and hypothesis testing) 43 Ph©n tích tơng quan hồi quy (Correlation and regression) 49 Chơng IV: Vẽ đồ thị 55 Vẽ đồ thị (graph) 55 Một số loại đồ thÞ th−êng dïng 64 Lu trữ hiển thị đồ thị (Saving and graph using) 71 Chơng V: Lập trình Stata 73 Giíi thiƯu chung vỊ ch−¬ng tr×nh do-file 73 Local vµ global macros 77 Tích vô hớng ma trận (scalar and matrix) 80 Lệnh điều kiện vòng lặp 82 Giíi thiƯu vỊ file ado 84 Tài liệu tham khảo 85 Phô lôc 85 Chơng I: Giới thiệu chung chơng trình Stata Tỉ chøc l−u tr÷ d÷ liƯu Stata (Dataset in Stata) Stata phần mềm thống kê sử dụng để quản lý, phân tích số liệu vẽ đồ thị Stata cho phép lu trữ thông tin đặc điểm đối tợng nghiên cứu Số liệu lu trữ Stata đợc hiển thị dới dạng bảng nh ví dụ sau: maho 101 102 103 tenchuho Nguyen Van A Le Thi B Tran Van C quymoho 10 thunhapbq 2100 3210 1200 Quan sát (observation) Mỗi hàng ngang bảng số liệu đợc gọi quan sát, hay ghi (record) lu trữ số liệu đối tợng nghiên cứu ví dụ có quan sát lu trữ số liệu M hộ (maho); Tên chủ hộ (tenchuho); Quy mô hộ (quymoho); Thu nhập bình quân (thunhapbq) hộ gia đình Biến (variable) Thông tin đối tợng nghiên cứu đợc thu thập lu trữ theo đặc điểm chúng Các đặc điểm đợc gọi biến Biến đợc xem cột bảng số liệu ví dụ có biến, với tên maho, tenchuho, quymoho, thunhapbq Tên biến dài từ đến 32 ký tự, đợc bắt đầu chữ dấu gạch dới (_) Tên biến bao gồm chữ, sỗ dấu gạch dới Các ký tự đặc biệt khác dùng để đặt tên cho biến Biến xác định (identifying variables) Thông thờng biến có biến dùng để nhận dạng quan sát, đợc gọi biến xác định Nhờ có biến xác định mà quan sát phân biệt đợc với Mỗi quan sát có giá trị biến ví dụ trên, biến xác định maho, quan sát biến maho nhận giá trị Các đặc điểm biến Các biến đợc g¸n nh‰n (chó thÝch) VÝ dơ biÕn maho cã thĨ đợc gán nhn M hộ Biến đợc định dạng (format) biến số biến ký tự với loại lu trữ khác Biến số lu trữ dới loại byte; int; long; float; double Còn biến ký tự lu trữ dới dạng str1 đến str80 cho độ dài khác Kiểu lu trữ dạng số byte int long float double Dung lợng (Byte) 4 Giá trị nhỏ Giá trị lớn Kiểu -127 -32,767 -2,147,483,647 -10^36 -10^308 126 32,766 2,147,483,646 10^36 10^308 Sè nguyªn Sè nguyên Số nguyên Số thực Số thực Các biến số bao gồm biến rời rạc liên tục Các biến nh quy mô hộ gia đình, giới tính chủ hộ, vùng địa lý, trình độ giáo dục biến rời rạc (discrete) (hay gọi biến phân loại (categorical)) Các biến đợc lu trữ dới dạng byte, int, long Các biến liên tục (continuous) nh thu nhập, chi tiêu hộ lu trữ dới dạng float double Biến ký tự (string) dùng để lu trữ loại ký tù VÝ dơ biÕn tenchuho lµ biÕn kiĨu ký tự dùng để lu trữ tên chủ hộ Kiểu lu trữ dạng chữ str1 str2 str80 Byte Độ dµi lín nhÊt 2 80 80 Khởi động thoát khỏi Stata (Open and exit) Stata đợc khởi động tơng tự nh chơng trình tin học ứng dụng khác, cách kích vào biểu tợng tệp wstata.exe Windows explorer, chọn cách chọn Start -> Program -> Stata Chơng trình đợc thoát b»ng lƯnh exit tõ cưa sỉ lƯnh Stata Command, tuỳ chọn exit thực đơn (menu) File Giao diện Stata (Stata interface)1 Sau Stata đợc khởi động, giao diện Stata đợc lên, bao gồm thực đơn (menu bar) cùng, dới công cụ (tool bar) cưa sỉ (windows) C¸c cưa sỉ cđa Stata C¸c cửa sổ Stata đợc mở việc lựa chọn tuỳ chọn thực đơn Windows (menu bar) Các cửa sổ bao gồm: Results Hiển thị lệnh kết Graph Hiển thị đồ thị Viewer Hiển thị cửa sổ trợ giúp (help) hiển thị nội dung file văn (text) Command Dùng để gõ câu lệnh Review Hiển thị lệnh đ thực Variables Hiển thị danh sách biến tệp số liệu Data editor Hiển thị sửa chữa số liệu dới dạng bảng Do-file editor Hiển thị cửa sổ để soạn thảo chơng trình Thanh thực đơn (Menu bar) Bằng cách kích vào thực đơn tuỳ chọn đó, Stata thực lệnh khác Thanh thực đơn bao gồm nhóm lệnh sau đây: File Open Mở file số liệu View Xem c¸c file cđa Stata cưa sỉ Viewer Save L−u file sè liÖu Save as L−u file sè liÖu dới tên File name Chọn tên file để đa vào cửa sổ lệnh Log Đóng, mở, xem lại log file Save graph Lu giữ file đồ thị Print graph In đồ thị Print results In kết Exit Thoát khỏi Stata Edit Copy text Sao chép văn (text) Copy tables Sao chép bảng biểu Paste Dán Table copy options Lùa chän chÐp b¶ng sè liƯu Graph copy options Lựa chọn chép đồ thị (không có Stata 7) Prefs Các tuỳ chọn màu sắc, phông chữ, kích cỡ Windows Results Mở cửa sổ kết Graph Mở cửa sổ đồ thị Log Mở cưa sỉ log file Viewer Më cưa sỉ trỵ gióp (help) vµ xem néi dung file Command Më cưa sỉ câu lệnh Review Mở cửa sổ lệnh đ thực hiƯn Variables Më cưa sỉ danh s¸ch c¸c biÕn cđa tƯp sè liƯu Help/Search Më cưa sỉ trỵ gióp (help) Data editor Mở cửa để xem số liệu lu trữ dới dạng bảng Do-file editor Mở cửa sổ viết chơng trình Help Các trợ giúp liên quan đến việc sử dụng Stata Thanh công cụ (tool bar) Các tuỳ chọn công cụ đợc thiết kế để thực lệnh thông dụng Stata Nếu di chuyển trỏ đến nút lên câu huớng dẫn, bao gồm: Open (use) Mở file sè liƯu Stata Save L−u tr÷ file sè liƯu ®Üa Print results In néi dung cđa cưa sỉ kết Begin log Mở, đóng xem nội dung cđa file log Start viewer Më cưa sỉ trỵ trỵ (help) Bring Dialog Window to font §−a cưa sỉ hép thoại phía trớc Bring Result Window to font Đa cưa sỉ kÕt qu¶ phÝa tr−íc Bring Graph Window to font Đa cửa sổ vẽ đồ thị phía trớc Do-file editor Mở cửa sổ soạn thảo chơng trình Data editor Më cưa sỉ sưa ch÷a sè liƯu Data browser Më cưa sỉ xem sè liƯu Clear –more- condition Tắt lệnh more Break Dừng việc thực lệnh chơng trình Nhập lu liệu (Use, input and and save) Më tƯp sè liƯu ®ang cã: Tệp số liệu Stata đợc mở lựa chọn Open thực đơn File; nút Open (use) công cụ tool bar Nếu file số liệu có dung lợng lớn phải thiết lập bé nhí cÇn dïng cho Stata b»ng lƯnh: set memory #[k|m] VÝ dô: set mem 32m set mem 32000k NhËp số liệu Có số cách để nhập số liệu tõ bµn phÝm vµo bé nhí cđa Stata - Sư dụng cửa sổ Stata editor để nhập số liệu Hoặc tõ cưa sỉ command, gâ lƯnh edit Sau ®ã nhËp số liệu theo kiểu biểu bảng cửa sổ Stata cho phÐp nhËp sè liƯu tõ c¸c file sở liệu khác Trớc hết file số liệu cần đợc lu trữ dới dạng text (có thể chơng trình Excel), quan sát đợc dòng giá trị cách phẩy (commas) dấu cách (tab) Sau dùng lệnh copy paste để nhập số liệu vào Stata L−u tr÷ sè liƯu ViƯc l−u tr÷ sè liƯu thực tùy chọn Save Save as thực đơn (menu bar); nút Save công cụ (tool bar) Chú ý: Xem thêm lệnh infile outfile Chơng II: Khai thác liệu Cấu trúc lệnh Stata (Stata command syntax) Cấu trúc lệnh Stata nh− sau: [by danh s¸ch biÕn:] Có ph¸p lƯnh [danh sách biến] [biểu thức] [điều kiện] [phạm vi] [quyền sè] [, t chän] Trong phÇn H−íng dÉn s− dơng (Help) Stata, cú pháp lệnh trình bày tiếng Anh nh− sau: [by varlist:] command [varlist] [=exp] [if exp] [in range] [weight] [, options] Trong dấu ngoặc vuông ký hiệu tuỳ chọn Chú ý: - Các câu lệnh Stata đợc viết chữ thờng - Đối với tên biến, Stata phân biệt chữ viết thờng với chữ viÕt hoa VÝ dơ, cïng mét tƯp sè liƯu, biến Ho_ten biến ho_ten biến khác - Các tuỳ chọn đợc ký hiệu dấu ngoặc vuông [ ] Các tuỳ chọn có không câu lệnh Các tham số bắt buộc (tên biến) đợc đặt dấu ngoặc < > Các câu lệnh không thực đợc tham số bắt buộc không đợc khai báo - Một sè lƯnh Stata cho phÐp viÕt t¾t VÝ dơ lƯnh summarize viết tắt sum Trong tài liệu phần gạch chân dới cú pháp câu lệnh cú pháp viết tắt câu lệnh - Các ví dụ tài liệu sử dụng số liệu Điều tra Mức sống dân c năm 1998 Tổng cục Thống kê tiến hành Trong Tệp chi tiêu tổng hợp Hhexp98n.dta thờng xuyên đợc sử dơng by danh s¸ch biÕn (by varlist): Stata sÏ thùc câu lệnh với theo giá trị đợc danh sách biến Biến đợc danh sách biến đợc yêu cầu sẵp xếp trớc thùc hiƯn lƯnh VÝ dơ: sort sex by sex: sum rlpcex1 -> sex = Variable | Obs Mean Std Dev Min Max -+ rlpcex1 | 4375 2980.906 2430.648 357.318 45801.71 -> sex = Variable | Obs Mean Std Dev Min Max -+ rlpcex1 | 1624 3748.368 3231.241 376.9805 30624.77 Dïng c«ng menu bar Danh s¸ch biÕn (varlist) ChØ danh s¸ch biến chịu tác động câu lệnh Nếu nh biến đợc lệnh Stata có tác dụng lên tất biến (all variables) VÝ dô: sum hhsize sex reg7 Variable | Obs Mean Std Dev Min Max -+ hhsize | 5999 4.752292 1.954292 19 sex | 5999 1.270712 4443645 reg7 | 5999 4.01917 2.145305 sum Variable | Obs Mean Std Dev Min Max -+ househol | 5999 19617.86 11201.92 101 38820 year | 5999 97.94666 2247337 97 98 month | 5999 6.340723 3.011082 12 Lệnh sum hiển thị thống kê tất biến tệp số liệu Điều kiện (if exp) Stata thực câu lệnh quan sát mà giá trị cho kết biểu thức Ví dô: sum poor if reg7==1 Variable | Obs Mean Std Dev Min Max -+ poor | 859 4982538 5002882 Lệnh có tác dụng quan sát mà biến reg7 có giá trị b»ng Ph¹m vi (in range) ChØ ph¹m vi quan sát chịu tác động câu lệnh Range (phạm vi) có dạng sau: sum poor in 10 Tính giá trị trung bình biến poor cho quan sát 10 (chính giá trị biến poor quan sát thứ 10) sum poor in 10/100 Tính giá trị trung bình biến poor cho quan sát từ 10 đến 100 sum poor in f/100 Tính giá trị trung bình biến poor cho quan sát từ đến 100 sum poor in 100/l Tính giá trị trung bình biến poor cho quan sát từ thứ 100 đến quan sát cuối Quyền số (weight) Cho phÐp tÝnh to¸n sđ dơng qun sè T chọn quyền số đợc trình bày kỹ mục chơng Các tuỳ chọn (Options) Nhiều câu lệnh Stata cho phép tuỳ chọn riêng Các tuỳ chọn đợc sau dấu phẩy Ví dơ: LƯnh sum cã t chän lµ detail, cho phÐp tính toán thêm số thống kê khác giá trị trung bình độ lệnh chuẩn sum rlpcex1, detail comp.M&Reg price adj.pc tot exp Percentiles Smallest 1% 682.9575 357.318 5% 1012.433 366.2792 10% 1238.088 376.9805 Obs 5999 25% 1671.054 381.3502 Sum of Wgt 5999 50% 2397.042 75% 90% 3711.917 5940.803 Largest 26944.64 30624.77 Mean Std Dev Variance 3188.667 2692.567 7249918 10 graph using "c:\do thi 1" "c:\do thi 2" "c:\do thi 3", margin(10) title("Mot so dac diem cua ho gia dinh") region1 region2 region3 persons 1:urban 98; 0:rural 98 12% foodpoor 18% poor but still above food povert 1500 70% nonpoor region4 region7 region5 region6 1000 500 Total region1 region3 region2 region5 region4 region7 region6 24% poor1 16% poor2 16% poor3 12% poor4 10% poor5 4% poor6 18% poor7 Mot so dac diem cua ho gia dinh Chó ý: Chóng ta co thĨ kÕt hỵp lƯnh saving với using để lu trữ đồ thị VÝ dô: graph using "c:\do thi 1" "c:\do thi 2" "c:\do thi 3", margin(10) title("Mot so dac die m cua ho gia dinh") saving("c:\do thi tong hop") graph using "c:\do thi tong hop" 72 Chơng V: Lập trình Stata Giới thiệu chung chơng trình do-file 1.1 Mở lu trữ do-file Stata cho phép viết tệp đợc gọi do-file bao gồm lệnh cđa Stata Thay v× thùc hiƯn tõng lƯnh mét tõ cửa sổ lệnh command, tệp do-file lần lợt thực lệnh Chơng trình Stata đợc soạn thảo cửa sổ do-file editor Cửa sổ đợc mở cách kích vào thực đơn Windows chọn tuỳ chọn do-file editor Một cách khác để mở cửa sổ gõ lệnh doedit cửa sổ lệnh command Ví dụ: Một chơng trình đợc soạn th¶o cưa sỉ do-file editor nh− sau: -clear set mem 32m use "C:\VLSS98\Hhexp98n.dta", clear tab urban98 sum hhsize gen new=hhsizet gen new=hhsize Sau soạn thảo, do-file đợc lu trữ tuỳ chọn Save as thực đơn File cửa sổ do-file editor Tên do-file đợc lệnh doedit nh sau: doedit (tên do-file) Tệp do-file có phần mở rộng ví dụ lu trữ đoạn chơng trình dới tên chơng trình th mục Vlss98 ổ đĩa C 1.2 Thực tệp do-file Để chạy do-file cửa sổ lệnh gâ mét hai lÖnh sau: filename [, nostop] run filename [, nostop] LƯnh run thùc hiƯn c¸c lƯnh do-file nhng không hiển thị kết hình Trong trình thực do-file, có câu lệnh sai Stata báo lỗi ngừng việc thực câu lệnh sau Tuy nhiên tuỳ chọn nostop đợc Stata bỏ qua câu lệnh bị lỗi tiếp tục thực lệnh sau câu lệnh lỗi Ví dụ: 73 "c:\vlss98\chuong trinh 1" clear set mem 32m (32768k) use "C:\VLSS98\Hhexp98n.dta", clear tab urban98 1:urban 98; | 0:rural 98 | Freq Percent Cum + Rural | 4269 71.16 71.16 Urban | 1730 28.84 100.00 + Total | 5999 100.00 sum hhsize Variable | Obs Mean Std Dev Min Max -+ hhsize | 5999 4.752292 1.954292 19 gen new=hhsizet hhsizet not found r(111); end of do-file r(111); Víi tuú chän nostop "c:\vlss98\chuong trinh 1", nostop clear set mem 32m (32768k) use "C:\VLSS98\Hhexp98n.dta", clear tab urban98 1:urban 98; | 0:rural 98 | Freq Percent Cum 74 + Rural | 4269 71.16 71.16 Urban | 1730 28.84 100.00 + Total | 5999 100.00 sum hhsize Variable | Obs Mean Std Dev Min Max -+ hhsize | 5999 4.752292 1.954292 19 gen new=hhsizet hhsizet not found r(111); gen new=hhsize end of do-file Thùc hiƯn (ch¹y) b»ng lÖnh run run "c:\vlss98\chuong trinh 1", nostop hhsizet not found C¸c do-file cã thĨ thùc hiƯn b»ng t chọn Do thực đơn File, thực trực tiÕp cưa sỉ Do-file editor b»ng t chän Do Run thực đơn Tool 1.3 Một số lu ý soạn thảo do-file version # Khi soạn thảo tệp do-file nên đa dòng lệnh vào đầu chơng trình để thông báo phiên Stata đợc dùng để soạn thảo do-file Ví dụ nh dùng Stata 7.0 để soạn thảo do-file câu lệnh đợc đa vào đầu chơng trình nh− sau: version 7.0 clear use Hhexp98n.dta tab reg7 … Các phiên Stata khác có khác cú pháp ý nghĩa câu lệnh Lệnh version cho phép chơng trình Stata chạy hiểu đợc nội dung tệp do-file đợc viết phiên khác set memory #[k|m] 75 Nếu nh file số liệu đòi hỏi nhớ lớn nhớ mà Stata sử dụng phải thiết lập nhớ lớn cho Stata lệnh Chú ý không nên thiÕt lËp bé nhí lín h¬n bé nhí cđa RAM m¸y tÝnh VÝ dơ: use "C:\Hhexp98n.dta", clear no room to add more observations r(901); set mem 32m (32768k) use "C:\Hhexp98n.dta", clear set more off/on Theo chÕ ®é ngầm định, thực lệnh nh kết việc xử lý lệnh dài cửa sổ kết (Stata Results), hình dừng lại phải ấn phím (chẳng hạn Enter Space bar) để kết tiếp tục đợc hiển thị Lệnh set more off cho phép kết không bị dừng lại mà đợc hiển thị liên tục thực xong câu lệnh do-file Lệnh set more on khôi phục lại chế độ ngầm định Ký tự * /* */ Stata không thực câu lệnh đợc bắt đầu ký tự * nằm hai nhóm ký tự /* */ Các ký tự dùng để viết thích do-file VÝ dô: -version 7.0 set mem 32m use "C:\Hhexp98n.dta", clear * Tao bien thu nhap cua ho gia dinh /* Bien bang Thu nhap binh quan nhan voi Quy mo ho*/ gen hhexp = rlpcex1 * hhsize #delimit ; Khi câu lệnh do-file editor dài dùng lệnh để thông báo câu lệnh đợc kết thúc ký tự (;) Theo chế độ ngầm định câu lệnh đợc kÕt thóc xng dßng b»ng viƯc gâ phÝm Enter Để khôi phục lại chế độ ngầm định dùng lệnh #delimit cr Ví dụ: lệnh vẽ đồ thị ë ch−¬ng tr−íc: graph meanexp meanedu ahhsize, title (Do thi chi tieu va hoc van chu ho) l1title(Chi tieu binh quan (tr dong)) l2title(So nam hoc cua chu ho) b2title (Quy mo ho gia dinh) xlabel ylabel xline (5 10 to 20) yline(2 to 8) connect(ll) 76 tuơng đơng víi: #delimit ; graph meanexp meanedu ahhsize, title (Do thi chi tieu va hoc van chu ho) l1title(Chi tieu binh quan (tr dong)) l2title(So nam hoc cua chu ho) b2title (Quy mo ho gia dinh) xlabel ylabel xline (5 10 to 20) yline(2 to 8) connect(ll) ; gen hhexp = rlpcex1 * hhsize ; …… Sau ®ã chóng ta nên khôi phục lại chế độ ngầm định nh câu lệnh sau viết dßng b»ng lƯnh: #delimit cr Chó ý: - Chóng ta dùng ký tự /* */ để viết câu lƯnh dµi nh− sau: graph meanexp meanedu ahhsize, title (Do thi chi tieu va hoc van chu ho) /* */ l1title(Chi tieu binh quan (tr dong)) l2title(So nam hoc cua chu ho) /* */ b2title (Quy mo ho gia dinh) xlabel ylabel xline (5 10 to 20) yline(2 to 8) connect(ll); - Các lệnh # delimit cách viết câu lệnh dài sử dụng ký tự /* */ dùng đợc do-file không dùng đợc cửa sỉ lƯnh command Local vµ global macros Macros lµ biến đợc dùng chơng trình Stata Biến macros đợc xem nh đoạn ký tự gọi macroname (tên macros) - tơng ứng với dy ký tự khác - đợc gọi macro contents (nội dung macro) Có hai loại macros local macros (macros néi bé) vµ global macros (macros toµn bé) 2.1 Local macros NÕu chóng ta gâ: local hogd “age hhsize rlpcex1” (DÊu nh¸y kÐp co thĨ bá qua, tøc lµ cã thĨ gâ: local hogd age hhsize rlpcex1) Khi `hogd đợc hiểu tơng đơng với: age hhsize rlpcex1 hogd đợc gọi tên macros, age hhsize rlpcex1 nội dung macros Để sử dơng néi dung cđa macros, chóng ta gâ tªn cđa macros dấu trích dẫn bên trái ( ) nằm phía bên trái bàn phím - dấu trích dẫn bên phải ( ) nằm phía phải bên dới bàn phím Nh− vËy nÕu chóng ta gâ: summarize `hogd’ th× tơng đơng với gõ: 77 summarize age hhsize rlpcex1 NÕu chóng ta gâ: local tb “summarize” th× chóng ta cã thĨ thùc hiƯn lƯnh summarize age hhsize rlpcex1 b»ng c¸ch gâ: `tb' `hogd' Variable | Obs Mean Std Dev Min Max -+ -age | 5999 48.01284 13.7702 16 95 hhsize | 5999 4.752292 1.954292 19 rlpcex1 | 5999 3188.667 2692.567 357.318 45801.71 Để hiển thị nội dung local macros gâ lƯnh macros list _(tªn local macros) VÝ dơ: macro list _hogd _hogd: age hhsize rlpcex1 Để xoá local macros dung lênh macros drop _(tªn local macros) VÝ dơ: macro drop _hogd macro list _hogd local macro `hogd' not found r(111); 2.2 Global macros NÕu chóng ta gâ: global diaban “reg7 province commune” (hc cã thĨ bá qua dÊu ngc kép: global diaban reg7 province commune) Khi $diaban tơng đơng với: reg7 province commune diaban đợc gọi tên macros, reg7 province commune nội dung macros Để sử dụng đợc nội dung global macros chóng ta gâ ký hiƯu $ liỊn tr−íc tªn cđa macros Nh− vËy nÕu chóng ta gâ: describe $diaban tơng đơng với gõ: describe : reg7 province commune describe $diaban 78 storage display value variable name type format label variable label reg7 int %8.0g Code by regions province float %9.0g Province code commune float %9.0g commune code PSU-SVY commands global mota "describe" $mota $diaban storage display value variable name type format label variable label reg7 int %8.0g Code by regions province float %9.0g Province code commune float %9.0g commune code PSU-SVY commands §Ĩ hiển thị nội dung global macros gâ lƯnh macros list (tªn global macros) VÝ dơ: global diaban "reg7 province commune" macro list diaban diaban: reg7 province commune Để xoá global macros cã thĨ dïng lƯnh macros drop (tªn local macros) VÝ dô: macro drop diaban macro list diaban global macro $diaban not found r(111); 2.3 Sự khác local macros global macros Local macros tồn chơng trình Một chơng trình không hiểu đợc local macros đợc sử dụng chơng trình khác Trong đó, đ đợc khai báo, global macros đợc hiểu tất chơng trình tồn nhớ Stata suốt trình hoạt động Ví dụ: Thực đoạn chơng trình khai báo local macros a Sau thực lệnh hiển thị nội dung local macros này, nhng macros không tồn đoạn chơng trinh khác hay ë bé nhí cđa Stata "C:\WINDOWS\TEMP\STD010000.tmp" local a "chuong trinh thong ke Stata" end of do-file 79 macro list _a local macro `a' not found r(111); Trong ®ã ®èi víi global macros "C:\WINDOWS\TEMP\STD010000.tmp" global b "chuong trinh thong ke Stata" end of do-file macro list b b: chuong trinh thong ke Stata Tích vô hớng ma trận (scalar and matrix) 3.1 Ma trận (matrix) Stata định nghĩa ma trận A[r, c] mảng hình chữ nhËt gåm r hµng (row) vµ c cét (column) VÝ dụ: Nếu ma trận A đ đợc tạo chóng ta cã thĨ xem néi dung cđa ma trËn ®ã nh− sau: matrix list A A[3,3] c1 c2 c3 r1 r2 r3 10 11 14 ma trận A bao gồm phÇn tư (element): 1, 2, 4, 3, 4, 7, 10, 11, 14 Các cột đợc đặt tên c1, c2, c3, hàng r1, r2, r3 Phần tử giao điểm dòng cột đợc ký hiệu A[1, 2] Trong ví dụ A[1, 2] chứa giá trị 3.2 Tích vô hớng (scalar) Tích vô hớng chứa phần tử số Tích vô hớng đợc định nghĩa lƯnh sau: scalar scalar_name = expression VÝ dơ: scalar a = 10 scalar list a a = 10 scalar b = a* scalar list b b= 20 80 Trong chừng mực đó, tích vô hớng xem nh trờng hợp đặc biệt ma trận có phân tử (một hàng vµ mét cét) 3.3 Mét sè lƯnh lµm viƯc víi ma trËn ThiĨt lËp kÝch th−íc ma trËn Gia trÞ ngầm định kích thớc ma trận tối đa 40 hµng vµ 40 cét Chóng ta cã thĨ thay ®ỉi kÝch th−íc tèi ®a nµy b»ng lƯnh: set matsize 500 Lệnh cho phép ma trận đợc tạo bao gồm 500 hàng 500 cét T¹o ma trËn Ma trËn cã thĨ t¹o câu lệnh trực tiếp Ví dụ: matrix mymat = (1,2\3,4) Các phần tử đợc phân biệt dấu phẩy, hàng đợc phân biệt dấu gạch chÐo matrix myvec = (1 3) T¹o véctơ hàng matrix mycol = (1/5/3/1/3) Tạo véctơ cột Ma trận đợc tạo từ số liệu lệnh: mkmat [điều kiện] [phạm vi] [, matrix(tên ma trận) ] Ví dụ: input maho quymo thunhap maho quymo thunhap 101 1200 103 1400 105 3200 107 1000 109 2500 end mkmat maho quymo thunhap, matrix(A) matrix list A A[5,3] maho quymo thunhap r1 101 1200 r2 103 1400 r3 105 3200 r4 107 1000 81 r5 109 2500 TÝnh to¸n ma trËn matrix D =B T¹o ma trËn D b»ng ma trËn B matrix C = (C+C)/2 TÝnh l¹i ma trËn C dùa giá trị của matrix D = A*A T¹o ma trËn D b»ng tÝch ma trËn A ma trận chuyển vị A Xoá ma trận Ma trận tích vô hớng xoá khỏi nhí b»ng lƯnh: matrix drop scalar drop VÝ dơ: matrix drop A scalar drop B Lệnh điều kiện vòng lặp 4.1 Lệnh ifelse Cú pháp: điều kiện (điểu kiện logic) { Nhãm c©u lƯnh } else C©u lƯnh Stata sÏ kiểm tra điều kiện logic (expression), điều kiện lệnh Nhóm câu lệnh đợc thực hiện, điều kiện sai lệnh ®»ng sau else sÏ ®−ỵc thùc hiƯn, tr−êng hỵp else không đợc Stata thực c¸c lƯnh sau lƯnh if {…} VÝ dơ: local a=invnorm(uniform()) if `a'>=0 { display "So ngau nhien tao lon hon hoac bang 0" } else di "So ngau nhien tao nho hon 0" macro list _a Chó ý: - Sư ký hiƯu { } cho phÐp viÕt nhiỊu c©u lƯnh sau else ®iỊu kiƯn (®iĨu kiƯn) { 82 commands } else { comands } - C¸c lƯnh if…else đợc sử dụng lồng với điều kiện (điểu kiện) { Nhóm câu lệnh } else điều kiện (điều kiện) { 4.2 Lệnh while Cú pháp: while { Nhóm câu lệnh } Stata sÏ kiĨm tra ®iỊu kiƯn logic (expression), nÕu ®iỊu kiƯn lệnh Nhóm câu lệnh đợc thực hiện, điều kiện sai lệnh không đợc thực Ví dụ: local i=1 while `i’