第13回:交差項を含む回帰・弾力性の推定

Similar documents
第11回:線形回帰モデルのOLS推定

第6回:データセットの結合

第2回:データの加工・整理

Stata11 whitepapers mwp-037 regress - regress regress. regress mpg weight foreign Source SS df MS Number of obs = 74 F(

最小2乗法

4 OLS 4 OLS 4.1 nurseries dual c dual i = c + βnurseries i + ε i (1) 1. OLS Workfile Quick - Estimate Equation OK Equation specification dual c nurser

y i OLS [0, 1] OLS x i = (1, x 1,i,, x k,i ) β = (β 0, β 1,, β k ) G ( x i β) 1 G i 1 π i π i P {y i = 1 x i } = G (

BR001

kubostat2017c p (c) Poisson regression, a generalized linear model (GLM) : :

Rによる計量分析:データ解析と可視化 - 第3回 Rの基礎とデータ操作・管理

!!! 2!

k2 ( :35 ) ( k2) (GLM) web web 1 :

<4D F736F F D20939D8C7689F090CD985F93C18EEA8D758B E646F63>

1. Stata( ステータ ) Stata は,StataCorp 社の販売している統計ソフトウェアで, 計量経済学においてもっともよく使われています 最新の計量経済学的手法の論文を執筆する際に,Stata による推定方法 ( コマンド ) も同時に発表されることがよくあり, 基本的な分析からより

2 Tobin (1958) 2 limited dependent variables: LDV 2 corner solution 2 truncated censored x top coding censor from above censor from below 2 Heck

28

Stata 11 Stata ts (ARMA) ARCH/GARCH whitepaper mwp 3 mwp-083 arch ARCH 11 mwp-051 arch postestimation 27 mwp-056 arima ARMA 35 mwp-003 arima postestim

一般化線形 (混合) モデル (2) - ロジスティック回帰と GLMM

4.9 Hausman Test Time Fixed Effects Model vs Time Random Effects Model Two-way Fixed Effects Model

ECCS. ECCS,. ( 2. Mac Do-file Editor. Mac Do-file Editor Windows Do-file Editor Top Do-file e

自由集会時系列part2web.key

第9回 日経STOCKリーグレポート 審査委員特別賞<地域の元気がでるで賞>

k3 ( :07 ) 2 (A) k = 1 (B) k = 7 y x x 1 (k2)?? x y (A) GLM (k

Stata 11 Stata ROC whitepaper mwp anova/oneway 3 mwp-042 kwallis Kruskal Wallis 28 mwp-045 ranksum/median / 31 mwp-047 roctab/roccomp ROC 34 mwp-050 s

講義のーと : データ解析のための統計モデリング. 第5回

kubostat2017e p.1 I 2017 (e) GLM logistic regression : : :02 1 N y count data or

講義のーと : データ解析のための統計モデリング. 第3回

1 Stata SEM LightStone 4 SEM 4.. Alan C. Acock, Discovering Structural Equation Modeling Using Stata, Revised Edition, Stata Press 3.

kubostat2015e p.2 how to specify Poisson regression model, a GLM GLM how to specify model, a GLM GLM logistic probability distribution Poisson distrib

: (EQS) /EQUATIONS V1 = 30*V F1 + E1; V2 = 25*V *F1 + E2; V3 = 16*V *F1 + E3; V4 = 10*V F2 + E4; V5 = 19*V99

Microsoft Word - Text_5_STATA_1_Jp_ doc

「スウェーデン企業におけるワーク・ライフ・バランス調査 」報告書

PDF


2010 II / y = e x y = log x = log e x 2. ( e x ) = e x 3. ( ) log x = 1 x 1.2 Warming Up 1 u = log a M a u = M a 0

Excelにおける回帰分析(最小二乗法)の手順と出力

80 X 1, X 2,, X n ( λ ) λ P(X = x) = f (x; λ) = λx e λ, x = 0, 1, 2, x! l(λ) = n f (x i ; λ) = i=1 i=1 n λ x i e λ i=1 x i! = λ n i=1 x i e nλ n i=1 x

計量経済分析 2011 年度夏学期期末試験 担当 : 別所俊一郎 以下のすべてに答えなさい. 回答は日本語か英語でおこなうこと. 1. 次のそれぞれの記述が正しいかどうか判定し, 誤りである場合には理由, あるいはより適切な 記述はどのようなものかを述べなさい. (1) You have to wo

H22 BioS t (i) treat1 treat2 data d1; input patno treat1 treat2; cards; ; run; 1 (i) treat = 1 treat =

JSP58-program


.. est table TwoSLS1 TwoSLS2 GMM het,b(%9.5f) se Variable TwoSLS1 TwoSLS2 GMM_het hi_empunion totchr

151021slide.dvi

kubostat2018d p.2 :? bod size x and fertilization f change seed number? : a statistical model for this example? i response variable seed number : { i

2 G(k) e ikx = (ik) n x n n! n=0 (k ) ( ) X n = ( i) n n k n G(k) k=0 F (k) ln G(k) = ln e ikx n κ n F (k) = F (k) (ik) n n= n! κ n κ n = ( i) n n k n

untitled

応力とひずみ.ppt

TS002

Stata 11 Stata VAR VEC whitepaper mwp 4 mwp-084 var VAR 14 mwp-004 varbasic VAR 25 mwp-005 svar VAR 31 mwp-007 vec intro VEC 47 mwp-008 vec VEC 75 mwp

Stata 11 whitepaper mwp 4 mwp mwp-028 / 41 mwp mwp mwp-079 functions 72 mwp-076 insheet 89 mwp-030 recode 94 mwp-033 reshape wide

renshumondai-kaito.dvi

2

1 Tokyo Daily Rainfall (mm) Days (mm)

Word 2000 Standard

68 A mm 1/10 A. (a) (b) A.: (a) A.3 A.4 1 1


R John Fox R R R Console library(rcmdr) Rcmdr R GUI Windows R R SDI *1 R Console R 1 2 Windows XP Windows * 2 R R Console R ˆ R

1 15 R Part : website:

Microsoft Word - 計量研修テキスト_第5版).doc

(pdf) (cdf) Matlab χ ( ) F t

1 Stata SEM LightStone 3 2 SEM. 2., 2,. Alan C. Acock, Discovering Structural Equation Modeling Using Stata, Revised Edition, Stata Press.

Solution Report

60 (W30)? 1. ( ) 2. ( ) web site URL ( :41 ) 1/ 77

7. フィリップス曲線 経済統計分析 (2014 年度秋学期 ) フィリップス曲線の推定 ( 経済理論との関連 ) フィリップス曲線とは何か? 物価と失業の関係 トレード オフ 政策運営 ( 財政 金融政策 ) への含意 ( 計量分析の手法 ) 関数形の選択 ( 関係が直線的でない場合の推定 ) 推

(lm) lm AIC 2 / 1

AR(1) y t = φy t 1 + ɛ t, ɛ t N(0, σ 2 ) 1. Mean of y t given y t 1, y t 2, E(y t y t 1, y t 2, ) = φy t 1 2. Variance of y t given y t 1, y t

2

2


2


9.6.3 SPSS (optional) SPSS jpn2003_jpn.sav CSV GRETL pisa2003stu_jpn.csv, pisa2003sch_jpn.csv gretl label SAS SPSS script OECD PISA SAS SPSS csv SPSS

東海道新幹線でDS



JA2008

tokei01.dvi

kubostat2017b p.1 agenda I 2017 (b) probability distribution and maximum likelihood estimation :

untitled

1 kawaguchi p.1/81

2 1,2, , 2 ( ) (1) (2) (3) (4) Cameron and Trivedi(1998) , (1987) (1982) Agresti(2003)

9 1 (1) (2) (3) (4) (5) (1)-(5) (i) (i + 1) 4 (1) (2) (3) (4) (5) (1)-(2) (1)-(5) (5) 1


1 I EViews View Proc Freeze

1 1.1 (JCPRG) 30 Nuclear Reaction Data File (NRDF) PC GSYS2.4 JCPRG GSYS2.4 Java Windows, Linux, Max OS X, FreeBSD GUI PNG, GIF, JPEG X Y GSYS2

X X X Y R Y R Y R MCAR MAR MNAR Figure 1: MCAR, MAR, MNAR Y R X 1.2 Missing At Random (MAR) MAR MCAR MCAR Y X X Y MCAR 2 1 R X Y Table 1 3 IQ MCAR Y I

I L01( Wed) : Time-stamp: Wed 07:38 JST hig e, ( ) L01 I(2017) 1 / 19

2017_Eishin_Style_H01

81

Introduction Purpose This training course demonstrates the use of the High-performance Embedded Workshop (HEW), a key tool for developing software for

The effect of smoking habit on the labor productivities

10:30 12:00 P.G. vs vs vs 2

Microsoft Word - 計量研修テキスト_第5版).doc

わが国企業による資金調達方法の選択問題

H22 BioS (i) I treat1 II treat2 data d1; input group patno treat1 treat2; cards; ; run; I


こんにちは由美子です

Microsoft Word - 計量研修テキスト_第5版).doc

グラフ数値読み取りシステム (GSYS2.4) 利用の手引

nakata/nakata.html p.1/20

Transcription:

13 2018 7 27 1 / 31

1. 2. 2 / 31

y i = β 0 + β X x i + β Z z i + β XZ x i z i + u i, E(u i x i, z i ) = 0, E(u i u j x i, z i ) = 0 (i j), V(u i x i, z i ) = σ 2, i = 1, 2,, n x i z i 1 3 / 31

y i = β 0 + β X x i + β Z z i + β XZ x i z i + u i x i y i x i = β X + β XZ z i. x i x i 1 y i z i y i = β 0 + β X x i + u i y i = β 1 x i x i z i 4 / 31

β X β XZ 0 z i 0 β X > 0, β XZ > 0 z i = 0 z i β X > 0, β XZ < 0 zi = 0 z i β X < 0, β XZ > 0 zi = 0 z i β X < 0, β XZ < 0 zi = 0 x i 5 / 31

β X β XZ 0 x i z i y i 6 / 31

Stata Stata generate c. #c. c. # 2 3 generate 7 / 31

c i = β 0 + β Y y i + β M m i + β D d i + β YM y i m i + u i (1) c i : y i : m i : d i : 0 8 / 31

1 1. Stata 2. File Log Begin... 2018microdata1 lecture20180727.smcl 3. File Open 4. consumption2009.dta 9 / 31

5. Statistics Linear models and related Linear regression 6. Model Dependent variable: expenditure_th 7. Independent variables: income_th, deposit_th, female 8. Model Independent variable: c.income_th#c.deposit_th c. # 10 / 31

9. SE/Robust Standard error type: Robust White 10. Reporting Set table formats 11. Coef/SE/CI Decimal format 4 decimals 4 12. p-value Decimal format 4 decimals p 4 13. Test statistic Decimal format 2 decimals t 2 11 / 31

14. Format settings for coefficient tables OK 15. regress - Linear regression OK 12 / 31

1% H 0 H0 1 H0 H0 13 / 31

ln y i = β 0 + β 1 ln x i + u i, E(u i x i ) = 0, E(u i u j x i ) = 0 (i j), V(u i x i ) = σ 2, i = 1, 2,, n ln y i ln x i ln 1 14 / 31

β 1 ln y i ln x i β 1 = d ln y i d ln x i. ln xi ln x i 1 15 / 31

d ln y i d ln x i = d ln yi d ln x i = d ln y i dy i = d ln y i dy i dy i y i. dx i x i dy i dx i dy i dx i dx i d ln x i 1. d ln x i dx i d ln y i dy i = 1 y i, d ln x i dx i = 1 x i 16 / 31

d ln y i dy i dy i dx i 1 d ln x i dx i = 1 dy i 1 y i dx 1 i x i = 1 dy i y dx i i x i 1 = dy i y i = dy i y i. dx i x i dx i x i d ln y i d ln x i = dy i y i dx i x i 17 / 31

β 1 = d ln y i d ln x i = dy i y i. dx i x i dx i : x i x i dy i : y i y i dx i x i : x i x i dy i y i : y i y i 18 / 31

β 1 = dy i y i dx i x i β 1, = y i x i. dy i y i dx i x i x i 1% y i % y i x i elasticity y i x i e.g., β 1 x i 1% y i β 1 % 19 / 31

ln y i = β 0 + β 1 ln x i + u i, i = 1, 2,..., n ln x i β 1 y i x i β 1 y i x i e.g., β 1 OLS ˆβ 1 y i x i H 0 : β 1 = 0 vs H 1 : β 1 0 0 y i x i 20 / 31

level Level-Level Model: y i = β 0 + β 1 x i + u i β1 x i 1 y i β 1 Log-Log Model: ln y i = β 0 + β 1 ln x i + u i β 1 x i 1% y i β 1 % 100 β 1 % 21 / 31

Log-Level Model: ln y i = β 0 + β 1 x i + u i β1 x i 1 y i 100 β 1 % β 1 y i x i semi-elasticity Level-Log Model: y i = β 0 + β 1 ln x i + u i β1 x i 1% y i 1 100 β 1 22 / 31

Log-Log Model Log-Level Model Log-Log Model Level-Level Model R 2 R 2 23 / 31

Stata generate Command generate ( )=ln(( )) generate ln ( 0 0 0 0 24 / 31

ln c i = β 0 + β Y ln y i + β M ln m i + β D d i + u i (2) % d i 0 25 / 31

2 1. Command generate lnexpenditure=ln(expenditure) Enter 2. Command generate lnincome=ln(income) Enter 3. Command generate lndeposit=ln(deposit) Enter 4. 26 / 31

5. Command list prefecture male female expenditure lnexpenditure income lnincome deposit lndeposit Enter list more 6 2 27 / 31

6. Statistics Linear models and related Linear regression 7. Model Dependent variable: lnexpenditure 8. Independent variables: Independent variables: lnincome, lndeposit, female 9. OK 28 / 31

0.7197 1% H0 0.7197 1% 0.7197% 1 0.039 5% H0 0.7197 1% 0.039% 29 / 31

0.1246 1% H0 12.46% 30 / 31

3 1. File Log Close 2. Stata lecture20180727.smcl 2018microdata1 3. lecture20180727.smcl 31 / 31