Đang tải... (xem toàn văn)
Hướng dẫn sử dụng Stata 11, nắm vững nguyên tắc chọn dữ liệu, trình bày dữ liệu, đọc và xử lý dữ liệu, trình bày kết quả
Sử dụng chơng trình Stata 11 (iu chnh t Stata 7) néi dung CH¬NG I: GIÍI THIƯU CHUNG VỊ CH¬NG TR×NH STATA 10 2 Tặ CHỉC LU TRữ Dữ LIệU TRONG STATA (DATASET IN STATA) KHậI đẫNG V THOáT KHI STATA (OPEN AND EXIT) GIAO DIÖN STATA 11 (STATA INTERFACE) NHËP VΜ LU D÷ LIƯU (USE, INPUT AND AND SAVE) 3 CHơNG II: KHAI THáC Dữ LIệU CÊU TRĨC LƯNH TRONG STATA (STATA COMMAND SYNTAX) 2.TO¸N Tệ V HM Sẩ (OPERATORS AND FUNCTIONS) Mô Tả Dữ LIệU (DATA REPORTING) BIêN TậP V SệA CHữA D÷ LIƯU (DATA MANIPULATION) QUN SÈ TRONG VHLSS (WEIGHT) 13 14 27 41 CHơNG III: KIểM địNH GIả THIÕT VΜ PH©N TÝCH HÅI QUY 45 ƯC LẻNG V KIểM địNH GIả THIếT (ESTIMATION AND HYPOTHESIS TESTING) PHâN TíCH TơNG QUAN V HI QUY (CORRELATION AND REGRESSION) 45 52 CHơNG IV: Vẽ đ THị 59 Vẽ đ THị (GRAPH) MẫT Sẩ LOạI đ THị THấNG DẽNG LU TRữ V HIểN THị đ THị (SAVING AND GRAPH USING) 59 67 74 CHơNG V: LậP TRìNH TRONG STATA 76 GII THIệU CHUNG Về CHơNG TRìNH DO-FILE LOCAL VΜ GLOBAL MACROS TÝCH V« HÍNG VΜ MA TRËN (SCALAR AND MATRIX) LệNH đIềU KIệN V VSSNG LặP GII THIƯU VỊ FILE ADO 76 81 83 86 88 TΜI LIƯU THAM KH¶O .90 PHÔ LÔC 90 Ch¬ng I: Giíi thiƯu chung chơng trình Stata 10 Tổ chức lu trữ 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 tenchuho 101 Nguyen Van A 102 Le Thi B 103 Tran Van C Quan s¸t (observation) quymo ho 10 thunhapbq 2100 3210 1200 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 gán nhãn 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ữ Dung lợng Giá trị nhỏ Giá trị lớn dạng số (Byte) nhÊt byte -127 126 int -32,767 32,766 long -2,147,483,647 2,147,483,646 float -10^36 10^36 double -10^308 10^308 KiĨu Sè nguyªn Sè nguyªn Sè nguyªn Sè thùc Sè thùc C¸c biÕn sè cã thĨ 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 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 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 lệnh exit từ cưa sỉ lƯnh Stata Command, hc t chän exit thực đơn (menu) File Giao diện Stata 11 (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ỉ cđa Stata đợc mở việc lựa chọn tuỳ chọn thực đơn Window (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 sÏ thùc hiƯn c¸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 Lu file sè liƯu Save as Lu file sè liƯu díi tªn míi 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 Copy as picture Sao chÐp thnh hỡnh 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 Lu 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 Đa cửa sổ hép tho¹i phÝa tríc font Bring Result Window to §a cưa sỉ kÕt qu¶ phÝa tríc font Bring Graph Window to Đa cửa sổ vẽ đồ thị phía trớc font 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 hiƯn 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 Hc 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ỉ nµy Stata cho phÐp nhËp sè liƯu tõ 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 vµ paste ®Ĩ nhËp sè liƯu vào Stata Lu trữ số liệu Việc lu trữ sè liƯu cã thĨ thùc hiƯn b»ng c¸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] [qun 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 b»ng 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” vµ biÕn “ho_ten” lµ 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 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 cụ menu bar 10 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 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] 76 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ụ: "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 | Urban | 4269 1730 71.16 28.84 71.16 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 77 "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 + Rural | 4269 Urban | 1730 71.16 28.84 71.16 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 hc Run thực đơn Tool 1.3 Một số lu ý soạn thảo do-file version # 78 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] 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 nhớ lớn nhớ 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 79 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) 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 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 cã thể 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ù /* */ chØ dïng đợc do-file không dùng đợc cửa sổ lƯnh command 80 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 dãy 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 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’ tơng đơng với gõ: 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 81 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 ngoặc 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 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ơ: 82 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 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) 83 Stata định nghĩa ma trận A[r, c] mảng hình chữ nhật gồm r hàng (row) c cột (column) Ví dụ: Nếu ma trận A đợc tạo th× 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 b»ng 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 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 84 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 Tạo véctơ hàng 3) 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 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 85 Xoá ma trận Ma trận tích vô hớng xoá khỏi bé 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), 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 đằng sau else đợc thực hiện, trờng hợp else không đợc Stata thực hiƯn 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) { commands } else { comands 86 } - C¸c lƯnh ifelse đợ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’