untitled

Size: px
Start display at page:

Download "untitled"

Transcription

1 IT (1, (1, ITIT IT IT TOPIX (%) 1TOPIX 2

2 ( ) 2004/11/26 S-PLUS 2

3 1-3. IT IT IT IT a. b. (Size) c. B/M(Book value to Market value) d. (VOL) 2004/11/26 S-PLUS 3

4 2-1-a. (Mkt) Mkt = (TOPIX ) () IT IT IT 2-1-b. (Size) Blanz and Reinganum, % Reinganum, 1992 () () 2004/11/26 S-PLUS 4

5 2-1-c. B/M() B/M B/M () () 2-1-d. (VOL) (=1100) 2004/11/26 S-PLUS 5

6 R 2-2. Fama and FrenchB/M Fama and French Fama and French rf Size B / t R rf Mkt M = t β Mkt + β Size + β B / M + β Vol + ε 0 t 1 B/M () 2004/11/26 S-PLUS 6 2 Vol β i ε 3 (1)

7 3. 12JASDAQ(1) () (2004.3) IT IT IT 3-a,b IT / // 3-a 2004/11/26 S-PLUS 3-b 7

8 3-2. ITIT 1999IT IT a2000ITIT 4ITIT1 5-b5-b2000IT IT1 IT /11/26 S-PLUS 8 5-a 5-b1

9 (Ri-rf) (Mkt) B/M(VOL)B/M B/M(<) () 3 Mkt(size)B/M (VOL) B/M(<) IT(5) () 2004/11/26 S-PLUS 9

10 4-2. *** Linear Model *** Call: lm(formula = Ri-rf ~ -1 + Mkt + Size + B/M + VOL, data = all, na.action = na.exclude) Residuals: Min 1Q Median 3Q Max Coefficients: Value Std. Error t value Pr(> t ) Mkt *** Size *** B/M * VOL *** Residual standard error: on degrees of freedom Multiple R-Squared: F-statistic: on 4 and degrees of freedom, the p-value is observations deleted due to missing values (***, **, * 1%, 5%, 10%) 1% B/M B/M10%1% /11/26 S-PLUS 10

11 4-3. : u u = ( R rf ) E[ β Mkt + β Size + β B / M + β VOL ] t 6 0 t 1 (2) 6IT IT IT IT1 3-aIT 2004/11/26 S-PLUS

12 (***, **, * 1%, 5%, 10%) ~ IT ~ T-test Sign test 4t( ) IT total ITtotal IT 2004/11/26 S-PLUS 12

13 4-5. F test for variance equality data: IT$Residuals and IT$Residuals F = , num df = 38312, denom df = , p-value = 0 alternative hypothesis: true ratio of variances is not equal to 1 data: IT$Residuals and IT$Residuals F = , num df = 38312, denom df = 34422, p-value = 0 alternative hypothesis: true ratio of variances is not equal to a IT ITITIT 4IT 7-a,b,c a IT 1998IT 3 IT 2004/11/26 S-PLUS 13

14 b c b7-a IT 1999 IT IT IT2000 IT ITIT1999 IT 2004/11/26 S-PLUS 14

15 IT IT (Eli Bartov, Partha Mohanram, and Chandrakanth Seethamraju) (SGR, Sales Growth Rate) (GPOS, Gross Profit On Sales)(CEOS, Capital Expenditure On Sales) (ROIC, Return On Invested Capital)(3) u = α + α CEOS + α GPOS + α ROIC + α SGR (3) 2004/11/26 S-PLUS 15

16 (***, **, * 1%, 5%, 10%) 5IT1999 (CEOS)(GPOS) 1999 IT 2004/11/26 S-PLUS 16

17 B/M4 () IT 2000(35) (4) (7) (5) IT (6) IT (64) 2004/11/26 S-PLUS 17

18 IT ITIT IT ITIT 5-3. IT () 1998ITIT IT 2004/11/26 S-PLUS 18

19 [1] Eugene F.Fama and Kenneth R.French Multifactor Explanations of Pricing Anomalies, The Journal of Finance(19963) [2] Eli Bartov, Partha Mohanram, and Chandrakanth Seethamraju, Valuation of Internet Stocks An IPO Perspective Journal of Accounting Research, Vol. 40 No.2 (20023) [3] James L.Davis and Eugene F.Fama and Kenneth R.French Characteristic, Covariance,and Average Returns The Center for Research insecurity Prices Working Paper No.471(1999) [4] Luc Keuleneer and Willem VerhoogRecent Trends in Valuation FromStrategy to ValueRoyal NIVRA (2003) [5] (2003) [6] (1997) [7] 2004/11/26 S-PLUS 19

20 (S-Plus) menudescribe(data = all, variables = "RiRf,Mkt,LNB,BtoM,VOL") menucor(data = all, variables = "RiRf,Mkt,LNB,BtoM,VOL",na.method = "Available") menulm(data = all, formula = "RiRf ~ -1 + Mkt + LNB + BtoM + VOL") 4 menudescribe(data = IT96.R, variables = "Residuals") ( XX96.RXX03.RIT, subit, nonit24total ) menuttest1(data.x = IT96.R, x = Residuals) 4-5. Fvar.test(IT$Residuals, subit$residuals) ( subitnonit2 ) 5 menulm(data = IT98.R, formula = Residuals ~ CEOS98 + GPOS98 + ROIC98 + SGR98) ( XX98.RXX00.RIT, subit, nonit9 ) 2004/11/26 S-PLUS 20

21 2 guiplot( PlotType = "", DataSet = "growth", Columns = "year,growth") guimodify( "LinePlot", Name = "GS1$1$1", LineStyle = "Solid", LineColor = "Lt Cyan", 6 guiplot( PlotType = "", DataSet = "Average.weekly.Residuals", Columns = "week,it, IT,IT") guimodify( "Graph2D", Name = "GS1$1", PanelType = "By Plot") LineWeight = "8", ConnectLineType = "To Y = 0") guimodify( "LinePlot", Name = "GS1$1$1", UseTextAsSymbol = T, VarySymbolText = "y Column") guimodify( "LinePlot", Name = "GS1$1$1", SymbolHeight = "0.2") guimodify( "Axis2dY", Name = "GS1$1$Axis2dY1", (ex. 1981/09/ ) AxisMin = "-2", AxisMax = "Auto") 2004/11/26 S-PLUS 21

J1順位と得点者数の関係分析

J1順位と得点者数の関係分析 2015 年度 S-PLUS & Visual R Platform 学生研究奨励賞応募 J1 順位と得点者数の関係分析 -J リーグの得点数の現状 - 目次 1. はじめに 2. 研究目的 データについて 3.J1 リーグの得点数の現状 4. 分析 5. まとめ 6. 今後の課題 - 参考文献 - 東海大学情報通信学部 経営システム工学科 山田貴久 1. はじめに 1993 年 5 月 15 日に

More information

(lm) lm AIC 2 / 1

(lm) lm AIC 2 / 1 W707 s-taiji@is.titech.ac.jp 1 / 1 (lm) lm AIC 2 / 1 : y = β 1 x 1 + β 2 x 2 + + β d x d + β d+1 + ϵ (ϵ N(0, σ 2 )) y R: x R d : β i (i = 1,..., d):, β d+1 : ( ) (d = 1) y = β 1 x 1 + β 2 + ϵ (d > 1) y

More information

DAA09

DAA09 > summary(dat.lm1) Call: lm(formula = sales ~ price, data = dat) Residuals: Min 1Q Median 3Q Max -55.719-19.270 4.212 16.143 73.454 Coefficients: Estimate Std. Error t value Pr(> t ) (Intercept) 237.1326

More information

Use R

Use R Use R! 2008/05/23( ) Index Introduction (GLM) ( ) R. Introduction R,, PLS,,, etc. 2. Correlation coefficient (Pearson s product moment correlation) r = Sxy Sxx Syy :, Sxy, Sxx= X, Syy Y 1.96 95% R cor(x,

More information

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

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 R John Fox 2006 8 26 2008 8 28 1 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 GUI R R R Console > ˆ 2 ˆ Fox(2005) jfox@mcmaster.ca

More information

201711grade2.pdf

201711grade2.pdf 2017 11 26 1 2 28 3 90 4 5 A 1 2 3 4 Web Web 6 B 10 3 10 3 7 34 8 23 9 10 1 2 3 1 (A) 3 32.14 0.65 2.82 0.93 7.48 (B) 4 6 61.30 54.68 34.86 5.25 19.07 (C) 7 13 5.89 42.18 56.51 35.80 50.28 (D) 14 20 0.35

More information

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

k3 ( :07 ) 2 (A) k = 1 (B) k = 7 y x x 1 (k2)?? x y (A) GLM (k 2012 11 01 k3 (2012-10-24 14:07 ) 1 6 3 (2012 11 01 k3) kubo@ees.hokudai.ac.jp web http://goo.gl/wijx2 web http://goo.gl/ufq2 1 3 2 : 4 3 AIC 6 4 7 5 8 6 : 9 7 11 8 12 8.1 (1)........ 13 8.2 (2) χ 2....................

More information

統計研修R分散分析(追加).indd

統計研修R分散分析(追加).indd http://cse.niaes.affrc.go.jp/minaka/r/r-top.html > mm mm TRT DATA 1 DM1 2537 2 DM1 2069 3 DM1 2104 4 DM1 1797 5 DM2 3366 6 DM2 2591 7 DM2 2211 8

More information

28

28 y i = Z i δ i +ε i ε i δ X y i = X Z i δ i + X ε i [ ] 1 δ ˆ i = Z i X( X X) 1 X Z i [ ] 1 σ ˆ 2 Z i X( X X) 1 X Z i Z i X( X X) 1 X y i σ ˆ 2 ˆ σ 2 = [ ] y i Z ˆ [ i δ i ] 1 y N p i Z i δ ˆ i i RSTAT

More information

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

k2 ( :35 ) ( k2) (GLM) web   web   1 : 2012 11 01 k2 (2012-10-26 16:35 ) 1 6 2 (2012 11 01 k2) (GLM) kubo@ees.hokudai.ac.jp web http://goo.gl/wijx2 web http://goo.gl/ufq2 1 : 2 2 4 3 7 4 9 5 : 11 5.1................... 13 6 14 6.1......................

More information

1 15 R Part : website:

1 15 R Part : website: 1 15 R Part 4 2017 7 24 4 : website: email: http://www3.u-toyama.ac.jp/kkarato/ kkarato@eco.u-toyama.ac.jp 1 2 2 3 2.1............................... 3 2.2 2................................. 4 2.3................................

More information

untitled

untitled 2011/6/22 M2 1*1+2*2 79 2F Y YY 0.0 0.2 0.4 0.6 0.8 0.000 0.002 0.004 0.006 0.008 0.010 0.012 1.0 1.5 2.0 2.5 3.0 3.5 4.0 Y 0 50 100 150 200 250 YY A (Y = X + e A ) B (YY = X + e B ) X 0.00 0.05 0.10

More information

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

1 Stata SEM LightStone 3 2 SEM. 2., 2,. Alan C. Acock, Discovering Structural Equation Modeling Using Stata, Revised Edition, Stata Press. 1 Stata SEM LightStone 3 2 SEM. 2., 2,. Alan C. Acock, 2013. Discovering Structural Equation Modeling Using Stata, Revised Edition, Stata Press. 2 3 2 Conservative Depress. 3.1 2. SEM. 1. x SEM. Depress.

More information

こんにちは由美子です

こんにちは由美子です Analysis of Variance 2 two sample t test analysis of variance (ANOVA) CO 3 3 1 EFV1 µ 1 µ 2 µ 3 H 0 H 0 : µ 1 = µ 2 = µ 3 H A : Group 1 Group 2.. Group k population mean µ 1 µ µ κ SD σ 1 σ σ κ sample mean

More information

Ł\”ƒ1PDFŠp

Ł\”ƒ1PDFŠp 48 2006 89 103 Keiji ABE This article investigates the first-day return and underwriter s ability on initial public offerings (IPOs). First, we find that, on average, offer prices of 405 IPOs listed on

More information

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

Stata11 whitepapers mwp-037 regress - regress regress. regress mpg weight foreign Source SS df MS Number of obs = 74 F( mwp-037 regress - regress 1. 1.1 1.2 1.3 2. 3. 4. 5. 1. regress. regress mpg weight foreign Source SS df MS Number of obs = 74 F( 2, 71) = 69.75 Model 1619.2877 2 809.643849 Prob > F = 0.0000 Residual

More information

最小2乗法

最小2乗法 2 2012 4 ( ) 2 2012 4 1 / 42 X Y Y = f (X ; Z) linear regression model X Y slope X 1 Y (X, Y ) 1 (X, Y ) ( ) 2 2012 4 2 / 42 1 β = β = β (4.2) = β 0 + β (4.3) ( ) 2 2012 4 3 / 42 = β 0 + β + (4.4) ( )

More information

% 10%, 35%( 1029 ) p (a) 1 p 95% (b) 1 Std. Err. (c) p 40% 5% (d) p 1: STATA (1). prtesti One-sample test of pr

% 10%, 35%( 1029 ) p (a) 1 p 95% (b) 1 Std. Err. (c) p 40% 5% (d) p 1: STATA (1). prtesti One-sample test of pr 1 1. 2014 6 2014 6 10 10% 10%, 35%( 1029 ) p (a) 1 p 95% (b) 1 Std. Err. (c) p 40% 5% (d) p 1: STATA (1). prtesti 1029 0.35 0.40 One-sample test of proportion x: Number of obs = 1029 Variable Mean Std.

More information

インターネットを活用した経済分析 - フリーソフト Rを使おう

インターネットを活用した経済分析 - フリーソフト Rを使おう R 1 1 1 2017 2 15 2017 2 15 1/64 2 R 3 R R RESAS 2017 2 15 2/64 2 R 3 R R RESAS 2017 2 15 3/64 2-4 ( ) ( (80%) (20%) 2017 2 15 4/64 PC LAN R 2017 2 15 5/64 R R 2017 2 15 6/64 3-4 R 15 + 2017 2 15 7/64

More information

> usdata01 と打ち込んでエンター キーを押すと V1 V2 V : : : : のように表示され 読み込まれていることがわかる ここで V1, V2, V3 は R が列のデータに自 動的につけた変数名である ( variable

> usdata01 と打ち込んでエンター キーを押すと V1 V2 V : : : : のように表示され 読み込まれていることがわかる ここで V1, V2, V3 は R が列のデータに自 動的につけた変数名である ( variable R による回帰分析 ( 最小二乗法 ) この資料では 1. データを読み込む 2. 最小二乗法によってパラメーターを推定する 3. データをプロットし 回帰直線を書き込む 4. いろいろなデータの読み込み方について簡単に説明する 1. データを読み込む 以下では read.table( ) 関数を使ってテキストファイル ( 拡張子が.txt のファイル ) のデー タの読み込み方を説明する 1.1

More information

こんにちは由美子です

こんにちは由美子です 1 2 . sum Variable Obs Mean Std. Dev. Min Max ---------+----------------------------------------------------- var1 13.4923077.3545926.05 1.1 3 3 3 0.71 3 x 3 C 3 = 0.3579 2 1 0.71 2 x 0.29 x 3 C 2 = 0.4386

More information

R Console >R ˆ 2 ˆ 2 ˆ Graphics Device 1 Rcmdr R Console R R Rcmdr Rcmdr Fox, 2007 Fox and Carvalho, 2012 R R 2

R Console >R ˆ 2 ˆ 2 ˆ Graphics Device 1 Rcmdr R Console R R Rcmdr Rcmdr Fox, 2007 Fox and Carvalho, 2012 R R 2 R John Fox Version 1.9-1 2012 9 4 2012 10 9 1 R R Windows R Rcmdr Mac OS X Linux R OS R R , R R Console library(rcmdr)

More information

36巻目次.indd

36巻目次.indd 36 12 論 文 知的財産評価と企業価値 岩城康史 ( 関西学院大学大学院 ) 岡田克彦 ( 関西学院大学 ) 要 1970 1977 2009 33 旨 キーワード 1 はじめに 400 120 1 2013 9 6 1,000 Griliches et al. 1986 1 2014 6 22 2 36 12 Trajtenberg 1990 Weighted Patent Count, WPC

More information

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

kubostat2018d p.2 :? bod size x and fertilization f change seed number? : a statistical model for this example? i response variable seed number : { i kubostat2018d p.1 I 2018 (d) model selection and kubo@ees.hokudai.ac.jp http://goo.gl/76c4i 2018 06 25 : 2018 06 21 17:45 1 2 3 4 :? AIC : deviance model selection misunderstanding kubostat2018d (http://goo.gl/76c4i)

More information

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

講義のーと :  データ解析のための統計モデリング. 第3回 Title 講義のーと : データ解析のための統計モデリング Author(s) 久保, 拓弥 Issue Date 2008 Doc URL http://hdl.handle.net/2115/49477 Type learningobject Note この講義資料は, 著者のホームページ http://hosho.ees.hokudai.ac.jp/~kub ードできます Note(URL)http://hosho.ees.hokudai.ac.jp/~kubo/ce/EesLecture20

More information

BMIdata.txt DT DT <- read.table("bmidata.txt") DT head(dt) names(dt) str(dt)

BMIdata.txt DT DT <- read.table(bmidata.txt) DT head(dt) names(dt) str(dt) ?read.table read.table(file, header = FALSE, sep = "", quote = "\" ", dec = ".", numerals = c("allow.loss", "warn.loss", "no.loss"), row.names, col.names, as.is =!stringsasfactors, na.strings = "NA", colclasses

More information

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

講義のーと :  データ解析のための統計モデリング. 第5回 Title 講義のーと : データ解析のための統計モデリング Author(s) 久保, 拓弥 Issue Date 2008 Doc URL http://hdl.handle.net/2115/49477 Type learningobject Note この講義資料は, 著者のホームページ http://hosho.ees.hokudai.ac.jp/~kub ードできます Note(URL)http://hosho.ees.hokudai.ac.jp/~kubo/ce/EesLecture20

More information

人工知能学会研究会資料 SIG-FPAI-B Predicting stock returns based on the time lag in information diffusion through supply chain networks 1 1 Yukinobu HA

人工知能学会研究会資料 SIG-FPAI-B Predicting stock returns based on the time lag in information diffusion through supply chain networks 1 1 Yukinobu HA 人工知能学会研究会資料 SIG-FPAI-B508-08 - - Predicting stock returns based on the time lag in information diffusion through supply chain networks 1 1 Yukinobu HAMURO 1 Katsuhiko OKADA 1 1 1 Kwansei Gakuin University

More information

1 CAPM: I-,,, I- ( ) 1 I- I- I- ( CAPM) I- CAPM I- 1 I- Jensen Fama-French 3 I- Fama-French 3 I- Fama-MacBeth I- SMB-FL, HML-FL Fama-MacBeth 1 Fama-Fr

1 CAPM: I-,,, I- ( ) 1 I- I- I- ( CAPM) I- CAPM I- 1 I- Jensen Fama-French 3 I- Fama-French 3 I- Fama-MacBeth I- SMB-FL, HML-FL Fama-MacBeth 1 Fama-Fr 1 CAPM: I-,,, I- ( ) 1 I- I- I- ( CAPM) I- CAPM I- 1 I- Jensen Fama-French 3 I- Fama-French 3 I- Fama-MacBeth I- SMB-FL, HML-FL Fama-MacBeth 1 Fama-French (FF) 3 [5] (Capital Asset Pricing Model; CAPM

More information

autocorrelataion cross-autocorrelataion Lo/MacKinlay [1988, 1990] (A)

autocorrelataion cross-autocorrelataion Lo/MacKinlay [1988, 1990] (A) Discussion Paper Series A No.425 2002 2 186-8603 iwaisako@ier.hit-u.ac.jp 14 1 24 autocorrelataion cross-autocorrelataion Lo/MacKinlay [1988, 1990] 1990 12 13 (A) 12370027 13 1 1980 Lo/MacKinlay [1988]

More information

R による共和分分析 1. 共和分分析を行う 1.1 パッケージ urca インスツールする 共和分分析をするために R のパッケージ urca をインスツールする パッケージとは通常の R には含まれていない 追加的な R のコマンドの集まりのようなものである R には追加的に 600 以上のパッ

R による共和分分析 1. 共和分分析を行う 1.1 パッケージ urca インスツールする 共和分分析をするために R のパッケージ urca をインスツールする パッケージとは通常の R には含まれていない 追加的な R のコマンドの集まりのようなものである R には追加的に 600 以上のパッ R による共和分分析 1. 共和分分析を行う 1.1 パッケージ urca インスツールする 共和分分析をするために R のパッケージ urca をインスツールする パッケージとは通常の R には含まれていない 追加的な R のコマンドの集まりのようなものである R には追加的に 600 以上のパッケージが用意されており それぞれ分析の目的に応じて標準の R にパッケージを追加していくことになる インターネットに接続してあるパソコンで

More information

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

一般化線形 (混合) モデル (2) - ロジスティック回帰と GLMM .. ( ) (2) GLMM kubo@ees.hokudai.ac.jp I http://goo.gl/rrhzey 2013 08 27 : 2013 08 27 08:29 kubostat2013ou2 (http://goo.gl/rrhzey) ( ) (2) 2013 08 27 1 / 74 I.1 N k.2 binomial distribution logit link function.3.4!

More information

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

kubostat2015e p.2 how to specify Poisson regression model, a GLM GLM how to specify model, a GLM GLM logistic probability distribution Poisson distrib kubostat2015e p.1 I 2015 (e) GLM kubo@ees.hokudai.ac.jp http://goo.gl/76c4i 2015 07 22 2015 07 21 16:26 kubostat2015e (http://goo.gl/76c4i) 2015 (e) 2015 07 22 1 / 42 1 N k 2 binomial distribution logit

More information

⑤香取 徹.indd

⑤香取 徹.indd 77 割引現在価値モデルの理論と有効性 Theory and Efficiency of Discounted Present Value Model Toru Katori FCF EVA FCF EVA RI 1 はじめに 1980 Free Cash Flow: FCF Copeland et. Al 19901990 Economic Value Added: EVA Stewart 1991

More information

上場変更と株価:株主分散と流動性変化のインパクト

上場変更と株価:株主分散と流動性変化のインパクト Merton Amihud and Mendelson NASDAQ JASDAQ JASDAQ JASDAQ QUICK E-mailjuno@waseda.jp E-mailshibata-mai@c.metro-u.ac.jp E-mailtakeshi.shimatani@boj.or.jp E-mailtokiko.shimizu@boj.or.jp JASDAQ JASDAQ Merton

More information

1 911 9001030 9:00 A B C D E F G H I J K L M 1A0900 1B0900 1C0900 1D0900 1E0900 1F0900 1G0900 1H0900 1I0900 1J0900 1K0900 1L0900 1M0900 9:15 1A0915 1B0915 1C0915 1D0915 1E0915 1F0915 1G0915 1H0915 1I0915

More information

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

1 Stata SEM LightStone 4 SEM 4.. Alan C. Acock, Discovering Structural Equation Modeling Using Stata, Revised Edition, Stata Press 3. 1 Stata SEM LightStone 4 SEM 4.. Alan C. Acock, 2013. Discovering Structural Equation Modeling Using Stata, Revised Edition, Stata Press 3. 2 4, 2. 1 2 2 Depress Conservative. 3., 3,. SES66 Alien67 Alien71,

More information

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

Microsoft Word - 計量研修テキスト_第5版).doc Q10-2 テキスト P191 1. 記述統計量 ( 変数 :YY95) 表示変数として 平均 中央値 最大値 最小値 標準偏差 観測値 を選択 A. 都道府県別 Descriptive Statistics for YY95 Categorized by values of PREFNUM Date: 05/11/06 Time: 14:36 Sample: 1990 2002 Included

More information

R分散分析06.indd

R分散分析06.indd http://cse.niaes.affrc.go.jp/minaka/r/r-top.html > mm mm TRT DATA 1 DM1 2537 2 DM1 2069 3 DM1 2104 4 DM1 1797 5 DM2 3366 6 DM2 2591 7 DM2 2211 8 DM2

More information

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

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 BR003 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 sampsi 47 mwp-044 sdtest 54 mwp-043 signrank/signtest

More information

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

kubostat2017c p (c) Poisson regression, a generalized linear model (GLM) : : kubostat2017c p.1 2017 (c), a generalized linear model (GLM) : kubo@ees.hokudai.ac.jp http://goo.gl/76c4i 2017 11 14 : 2017 11 07 15:43 kubostat2017c (http://goo.gl/76c4i) 2017 (c) 2017 11 14 1 / 47 agenda

More information

ID JJF00296 Earnings management and market timing: The case of Japanese management buyouts Yasutomo Tsukioka 2-16 Japan Journal of Finance Vol.

ID JJF00296 Earnings management and market timing: The case of Japanese management buyouts Yasutomo Tsukioka 2-16 Japan Journal of Finance Vol. ID JJF00296 Earnings management and market timing: The case of Japanese management buyouts Yasutomo Tsukioka 2-16 Japan Journal of Finance 331.2 Vol.33 / No. 1.2. 201312 Dec. 2013 Japan Finance Association

More information

yamadaiR(cEFA).pdf

yamadaiR(cEFA).pdf R 2012/10/05 Kosugi,E.Koji (Yamadai.R) Categorical Factor Analysis by using R 2012/10/05 1 / 9 Why we use... 3 5 Kosugi,E.Koji (Yamadai.R) Categorical Factor Analysis by using R 2012/10/05 2 / 9 FA vs

More information

σ t σ t σt nikkei HP nikkei4csv H R nikkei4<-readcsv("h:=y=ynikkei4csv",header=t) (1) nikkei header=t nikkei4csv 4 4 nikkei nikkei4<-dataframe(n

σ t σ t σt nikkei HP nikkei4csv H R nikkei4<-readcsv(h:=y=ynikkei4csv,header=t) (1) nikkei header=t nikkei4csv 4 4 nikkei nikkei4<-dataframe(n R 1 R R R tseries fseries 1 tseries fseries R Japan(Tokyo) R library(tseries) library(fseries) 2 t r t t 1 Ω t 1 E[r t Ω t 1 ] ɛ t r t = E[r t Ω t 1 ] + ɛ t ɛ t 2 iid (independently, identically distributed)

More information

回帰分析 単回帰

回帰分析 単回帰 回帰分析 単回帰 麻生良文 単回帰モデル simple regression model = α + β + u 従属変数 (dependent variable) 被説明変数 (eplained variable) 独立変数 (independent variable) 説明変数 (eplanator variable) u 誤差項 (error term) 撹乱項 (disturbance term)

More information

L Y L( ) Y0.15Y 0.03L 0.01L 6% L=(10.15)Y 108.5Y 6%1 Y y p L ( 19 ) [1990] [1988] 1

L Y L( ) Y0.15Y 0.03L 0.01L 6% L=(10.15)Y 108.5Y 6%1 Y y p L ( 19 ) [1990] [1988] 1 1. 1-1 00 001 9 J-REIT 1- MM CAPM 1-3 [001] [1997] [003] [001] [1999] [003] 1-4 0 . -1 18 1-1873 6 1896 L Y L( ) Y0.15Y 0.03L 0.01L 6% L=(10.15)Y 108.5Y 6%1 Y y p L 6 1986 ( 19 ) -3 17 3 18 44 1 [1990]

More information

textream ( )

textream ( ) 27 EM151010 2016 1 18 textream ( ) 1 4 2 6 3 8 4 10 5 12 5.1.......................................... 12 5.2 Comtemporary regression................................... 13 5.3........................................

More information

卒業論文

卒業論文 Y = ax 1 b1 X 2 b2...x k bk e u InY = Ina + b 1 InX 1 + b 2 InX 2 +...+ b k InX k + u X 1 Y b = ab 1 X 1 1 b 1 X 2 2...X bk k e u = b 1 (ax b1 1 X b2 2...X bk k e u ) / X 1 = b 1 Y / X 1 X 1 X 1 q YX1

More information

I- Fama-French 3, Idiosyncratic (I- ) I- ( ) 1 I- I- I- 1 I- I- Jensen Fama-French 3 SMB-FL, HML-FL I- Fama-French 3 I- Fama-MacBeth Fama-MacBeth I- S

I- Fama-French 3, Idiosyncratic (I- ) I- ( ) 1 I- I- I- 1 I- I- Jensen Fama-French 3 SMB-FL, HML-FL I- Fama-French 3 I- Fama-MacBeth Fama-MacBeth I- S I- Fama-French 3, Idiosyncratic (I- ) I- ( ) 1 I- I- I- 1 I- I- Jensen Fama-French 3 SMB-FL, HML-FL I- Fama-French 3 I- Fama-MacBeth Fama-MacBeth I- SMB-FL, HML-FL Fama-MacBeth 2, 3, 5 I- HML-FL 1 Fama-French

More information

早稲田大学現代政治経済研究所 ダブルトラック オークションの実験研究 宇都伸之早稲田大学上條良夫高知工科大学船木由喜彦早稲田大学 No.J1401 Working Paper Series Institute for Research in Contemporary Political and Ec

早稲田大学現代政治経済研究所 ダブルトラック オークションの実験研究 宇都伸之早稲田大学上條良夫高知工科大学船木由喜彦早稲田大学 No.J1401 Working Paper Series Institute for Research in Contemporary Political and Ec 早稲田大学現代政治経済研究所 ダブルトラック オークションの実験研究 宇都伸之早稲田大学上條良夫高知工科大学船木由喜彦早稲田大学 No.J1401 Working Paper Series Institute for Research in Contemporary Political and Economic Affairs Waseda University 169-8050 Tokyo,Japan

More information

11夏特集号初校.indd

11夏特集号初校.indd 1 2 3 5 50 40 7 6 3 ABC 3 5 A 5% B C 100 3 1 2 3 A 5% 5% 5% B 10% 5% 0% C 20% 10% 15% A 15.8% 15.0% 0.8% B 15.5% 15.0% 0.5% C 12.2% 15.0% 2.8% 2,000 1,500 1,000 500 0 10% 5% 3% 1% 01 5 10 15 20 25 30

More information

表紙4_1/山道 小川内 小川内 芦塚

表紙4_1/山道 小川内 小川内 芦塚 1 2008.1Vol.23 2008.1Vol.23 2 2008.1Vol.23 3 2008.1Vol.23 4 5 2008.1Vol.23 2008.1Vol.23 6 7 2008.1Vol.23 2008.1Vol.23 8 9 2008.1Vol.23 10 2008.1Vol.23 11 2008.1Vol.23 12 2008.1Vol.23 Center 13 2008.1Vol.23

More information

Page

Page Page Page 1 3 4 M&A DCF NPV 1-1 BIS 1-3 1-2 5 6 1-3 2 2 0 200 010 2-1 A B B A B 7 8 Coffee Break CAPM Capital Asset Pricing model 3 IRR NPV EVA BS PL 9 2 ROA 3-1 EBIT Earnings Before Interest,Taxes) 3-2

More information

と入力する すると最初の 25 行が表示される 1 行目は変数の名前であり 2 列目は企業番号 (1,,10),3 列目は西暦 (1935,,1954) を表している ( 他のパネルデータを分析する際もデ ータをこのように並べておかなくてはならない つまりまず i=1 を固定し i=1 の t に関

と入力する すると最初の 25 行が表示される 1 行目は変数の名前であり 2 列目は企業番号 (1,,10),3 列目は西暦 (1935,,1954) を表している ( 他のパネルデータを分析する際もデ ータをこのように並べておかなくてはならない つまりまず i=1 を固定し i=1 の t に関 R によるパネルデータモデルの推定 R を用いて 静学的パネルデータモデルに対して Pooled OLS, LSDV (Least Squares Dummy Variable) 推定 F 検定 ( 個別効果なしの F 検定 ) GLS(Generalized Least Square : 一般化最小二乗 ) 法による推定 およびハウスマン検定を行うやり方を 動学的パネルデータモデルに対して 1 階階差

More information

目次 第 1 章序論 第 2 章データの概要 第 3 章 J リーグの現状 第 4 章分析 第 5 章まとめ 第 6 章今後の課題

目次 第 1 章序論 第 2 章データの概要 第 3 章 J リーグの現状 第 4 章分析 第 5 章まとめ 第 6 章今後の課題 2016 年度 S-PLUS & Visual R Platform 学生研究奨励賞応募 J1 順位と観客動員数の関係分析 - ゴールの重要性 - 目次 第 1 章序論 第 2 章データの概要 第 3 章 J リーグの現状第 4 章分析 第 5 章まとめ 第 6 章今後の課題 - 参考文献 - 東海大学情報通信学部 経営システム工学科 山田貴久 目次 第 1 章序論 第 2 章データの概要 第 3

More information

fiúŁÄ”s‘ê‡ÌŁª”U…−…X…N…v…„…~…A…•‡Ì ”s‘ê™´›ß…−…^†[…fiŠ‚ª›Âfl’«

fiúŁÄ”s‘ê‡ÌŁª”U…−…X…N…v…„…~…A…•‡Ì ”s‘ê™´›ß…−…^†[…fiŠ‚ª›Âfl’« 2016/3/11 Realized Volatility RV 1 RV 1 Implied Volatility IV Volatility Risk Premium VRP 1 (Fama and French(1988) Campbell and Shiller(1988)) (Hodrick(1992)) (Lettau and Ludvigson (2001)) VRP (Bollerslev

More information

2 I- I- (1) 2 I- (2) 2 I- 1 [18] I- I-. 1 I- I- Jensen [11] I- FF 3 I- FF 3 2 2.1 CAPM n ( i = 1,..., n) M t R i,t, i = 1,..., n R M,t ( ) R i,t = r i

2 I- I- (1) 2 I- (2) 2 I- 1 [18] I- I-. 1 I- I- Jensen [11] I- FF 3 I- FF 3 2 2.1 CAPM n ( i = 1,..., n) M t R i,t, i = 1,..., n R M,t ( ) R i,t = r i 1 Idiosyncratic,, Idiosyncratic (I- ) I- 1 I- I- Jensen I- Fama-French 3 I- Fama-French 3 1 Fama-French (FF) 3 [6] (Capital Asset Pricing Model; CAPM [12, 15]) CAPM ( [2, 10, 14, 16]) [18] Idiosyncratic

More information

1 kawaguchi p.1/81

1 kawaguchi p.1/81 1 kawaguchi atsushi@kurume-u.ac.jp 2005 7 2 p.1/81 2.1 2.2 2.2.3 2.3 AUC 4.4 p.2/81 X Z X = α + βz + e α : Z = 0 X ( ) β : Z X ( ) e : 0 σ 2 p.3/81 2.1 Z X 1 0.045 2 0.114 4 0.215 6 0.346 7 0.41 8 0.52

More information

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

わが国企業による資金調達方法の選択問題 * takeshi.shimatani@boj.or.jp ** kawai@ml.me.titech.ac.jp *** naohiko.baba@boj.or.jp No.05-J-3 2005 3 103-8660 30 No.05-J-3 2005 3 1990 * E-mailtakeshi.shimatani@boj.or.jp ** E-mailkawai@ml.me.titech.ac.jp

More information

untitled

untitled Cross [1973]French [1980] Rogalsk [1984]Arel [1990]Arel [1987]Rozeff and Knney [1974]seasonaltycalendar structure [2004] 12 Half-Year Effect [2004] [1983] [1990]ChanHamao and Lakonshok [1991] Fama and

More information

アダストリア売り上げデータによる 現状把握と今後の方針 東海大学情報通信学部経営システム工学科佐藤健太

アダストリア売り上げデータによる 現状把握と今後の方針 東海大学情報通信学部経営システム工学科佐藤健太 アダストリア売り上げデータによる 現状把握と今後の方針 東海大学情報通信学部経営システム工学科佐藤健太 目次 1. 研究背景 2. 研究目的 3. データ概要 4. 分析手順 5. 分析結果 6. 戦略予想 7. まとめ 8. 今後の課題 参考文献 2016/10/27 2016 年 S-PLUS &Visual R Platform 学生研究奨励賞 1 1. 研究背景 Ⅰ アダストリア (¹) とは,

More information

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

kubostat2017e p.1 I 2017 (e) GLM logistic regression : : :02 1 N y count data or kubostat207e p. I 207 (e) GLM kubo@ees.hokudai.ac.jp https://goo.gl/z9ycjy 207 4 207 6:02 N y 2 binomial distribution logit link function 3 4! offset kubostat207e (https://goo.gl/z9ycjy) 207 (e) 207 4

More information

オーストラリア研究紀要 36号(P)☆/3.橋本

オーストラリア研究紀要 36号(P)☆/3.橋本 36 p.9 202010 Tourism Demand and the per capita GDP : Evidence from Australia Keiji Hashimoto Otemon Gakuin University Abstract Using Australian quarterly data1981: 2 2009: 4some time-series econometrics

More information

<4D F736F F F696E74202D BD95CF97CA89F090CD F6489F18B4195AA90CD816A>

<4D F736F F F696E74202D BD95CF97CA89F090CD F6489F18B4195AA90CD816A> 主な多変量解析 9. 多変量解析 1 ( 重回帰分析 ) 目的変数 量的 説明変数 質的 あり量的 重回帰分析 数量化 Ⅰ 類 質的 判別分析 数量化 Ⅱ 類 なし 主成分分析因子分析多次元尺度構成法 数量化 Ⅲ 類数量化 Ⅳ 類 その他 クラスタ分析共分散構造分析 説明変数 : 独立変数 予測変数 目的変数 : 従属変数 基準変数 3 1. 単回帰分析各データの構造 y b ax a α: 1,,,

More information

高齢化とマクロ投資比率―国際パネルデータを用いた分析―

高齢化とマクロ投資比率―国際パネルデータを用いた分析― 196 2017 * ** ** ** ** 160 2 2 JEL Classification Codes E21, E22, J11 Keywords * ESRI 28 ESRI 29 3 17 ESRI ** 115 196 Population Aging and Domestic Investment An Analysis Using International Panel Data

More information

!!! 2!

!!! 2! 2016/5/17 (Tue) SPSS (mugiyama@l.u-tokyo.ac.jp)! !!! 2! 3! 4! !!! 5! (Population)! (Sample) 6! case, observation, individual! variable!!! 1 1 4 2 5 2 1 5 3 4 3 2 3 3 1 4 2 1 4 8 7! (1) (2) (3) (4) categorical

More information

研修コーナー

研修コーナー l l l l l l l l l l l α α β l µ l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l l

More information

Vol. 42 No pp Headcount ratio p p A B pp.29

Vol. 42 No pp Headcount ratio p p A B pp.29 1990 2003 2005 2000 1998 2004 2001 2 2000 2001 2000 1 Vol. 42 No. 2 2005 pp.21-22 25 25-29 30-34 1999 1 Headcount ratio 2 1995 20-25 25-30 2005 p.25 2005 2000 2 15 34 2003 p.3 15 34 A B 3 4 3 3 2003 pp.29-332001

More information

Y-...W

Y-...W 259 21 1 1 2003 2 1 2 1 15 22 260 56 3 2 2. 1 2. 2 2. 3 2. 4 3 3. 1 3. 2 4 2 2. 1 1965 1999 2 Brick, Frierman, and Kim 1998 Choate 1997Fischer, Heinkel, and Zechner 1989 Gilson 1997 Mauer and Triantis

More information

130 Oct Radial Basis Function RBF Efficient Market Hypothesis Fama ) 4) 1 Fig. 1 Utility function. 2 Fig. 2 Value function. (1) (2)

130 Oct Radial Basis Function RBF Efficient Market Hypothesis Fama ) 4) 1 Fig. 1 Utility function. 2 Fig. 2 Value function. (1) (2) Vol. 47 No. SIG 14(TOM 15) Oct. 2006 RBF 2 Effect of Stock Investor Agent According to Framing Effect to Stock Exchange in Artificial Stock Market Zhai Fei, Shen Kan, Yusuke Namikawa and Eisuke Kita Several

More information

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

Microsoft Word - 計量研修テキスト_第5版).doc Q9-1 テキスト P166 2)VAR の推定 注 ) 各変数について ADF 検定を行った結果 和文の次数はすべて 1 である 作業手順 4 情報量基準 (AIC) によるラグ次数の選択 VAR Lag Order Selection Criteria Endogenous variables: D(IG9S) D(IP9S) D(CP9S) Exogenous variables: C Date:

More information

以下の内容について説明する 1. VAR モデル推定する 2. VAR モデルを用いて予測する 3. グレンジャーの因果性を検定する 4. インパルス応答関数を描く 1. VAR モデルを推定する ここでは VAR(p) モデル : R による時系列分析の方法 2 y t = c + Φ 1 y t

以下の内容について説明する 1. VAR モデル推定する 2. VAR モデルを用いて予測する 3. グレンジャーの因果性を検定する 4. インパルス応答関数を描く 1. VAR モデルを推定する ここでは VAR(p) モデル : R による時系列分析の方法 2 y t = c + Φ 1 y t 以下の内容について説明する 1. VAR モデル推定する 2. VAR モデルを用いて予測する 3. グレンジャーの因果性を検定する 4. インパルス応答関数を描く 1. VAR モデルを推定する ここでは VAR(p) モデル : R による時系列分析の方法 2 y t = c + Φ 1 y t 1 + + Φ p y t p + ε t, ε t ~ W.N(Ω), を推定することを考える (

More information

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

Rによる計量分析:データ解析と可視化 - 第3回  Rの基礎とデータ操作・管理 R 3 R 2017 Email: gito@eco.u-toyama.ac.jp October 23, 2017 (Toyama/NIHU) R ( 3 ) October 23, 2017 1 / 34 Agenda 1 2 3 4 R 5 RStudio (Toyama/NIHU) R ( 3 ) October 23, 2017 2 / 34 10/30 (Mon.) 12/11 (Mon.)

More information

(2/24) : 1. R R R

(2/24) : 1. R R R R? http://hosho.ees.hokudai.ac.jp/ kubo/ce/2004/ : kubo@ees.hokudai.ac.jp (2/24) : 1. R 2. 3. R R (3/24)? 1. ( ) 2. ( I ) : (p ) : cf. (power) p? (4/24) p ( ) I p ( ) I? ( ) (5/24)? 0 2 4 6 8 A B A B (control)

More information

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

第11回:線形回帰モデルのOLS推定 11 OLS 2018 7 13 1 / 45 1. 2. 3. 2 / 45 n 2 ((y 1, x 1 ), (y 2, x 2 ),, (y n, x n )) linear regression model y i = β 0 + β 1 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

More information

16soukatsu_p1_40.ai

16soukatsu_p1_40.ai 2 2016 DATA. 01 3 DATA. 02 4 DATA. 03 5 DATA. 04 6 DATA. 05 7 DATA. 06 8 DATA. 07 9 DATA. 08 DATA. 09 DATA. 10 DATA. 11 DATA. 12 DATA. 13 DATA. 14 10 11 12 13 COLUMN 1416 17 18 19 DATA. 15 20 DATA. 16

More information

Microsoft Word - StatsDirectMA Web ver. 2.0.doc

Microsoft Word - StatsDirectMA Web ver. 2.0.doc Web version. 2.0 15 May 2006 StatsDirect ver. 2.0 15 May 2006 2 2 2 Meta-Analysis for Beginners by using the StatsDirect ver. 2.0 15 May 2006 Yukari KAMIJIMA 1), Ataru IGARASHI 2), Kiichiro TSUTANI 2)

More information

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

H22 BioS (i) I treat1 II treat2 data d1; input group patno treat1 treat2; cards; ; run; I H BioS (i) I treat II treat data d; input group patno treat treat; cards; 8 7 4 8 8 5 5 6 ; run; I II sum data d; set d; sum treat + treat; run; sum proc gplot data d; plot sum * group ; symbol c black

More information

要約

要約 No.7, 563-574 (2006) A Consideration of the Goodyear Tire & Rubber Company Global Strategies and Valuation of These Strategies USHIYAMA Yukio Nihon University, Graduate School of Social and Cultural Studies

More information

表紙_目次.PDF

表紙_目次.PDF JIL 2 1980 vol.3 15 3 1 3 15 18 IT 26 28 33 20 OECD 80 90 (1998) (2000) 2 Aghion et al. (1999) OECD 1970 1 2 70 1980 90 1.1 Kuznets U inverted U-shaped hypothesis GNP U 1 18901940 6 3 Kuznets(1955) (1963)

More information

151021slide.dvi

151021slide.dvi : Mac I 1 ( 5 Windows (Mac Excel : Excel 2007 9 10 1 4 http://asakura.co.jp/ books/isbn/978-4-254-12172-8/ (1 1 9 1/29 (,,... (,,,... (,,, (3 3/29 (, (F7, Ctrl + i, (Shift +, Shift + Ctrl (, a i (, Enter,

More information

Ver. 1.0 (discretion) (revenue management) (fictitious revenue recognition) (window dressing) (fraudulent accounting practice) SEC (General

Ver. 1.0 (discretion) (revenue management) (fictitious revenue recognition) (window dressing) (fraudulent accounting practice) SEC (General 2004.11.25 Ver. 1.0 (discretion) (revenue management) (fictitious revenue recognition) (window dressing) (fraudulent accounting practice) SEC (Generally Accepted Accounting Principles: GAAP) GAAP (realization

More information

: (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

: (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 218 6 219 6.11: (EQS) /EQUATIONS V1 = 30*V999 + 1F1 + E1; V2 = 25*V999 +.54*F1 + E2; V3 = 16*V999 + 1.46*F1 + E3; V4 = 10*V999 + 1F2 + E4; V5 = 19*V999 + 1.29*F2 + E5; V6 = 17*V999 + 2.22*F2 + E6; CALIS.

More information

Isogai, T., Building a dynamic correlation network for fat-tailed financial asset returns, Applied Network Science (7):-24, 206,

Isogai, T., Building a dynamic correlation network for fat-tailed financial asset returns, Applied Network Science (7):-24, 206, H28. (TMU) 206 8 29 / 34 2 3 4 5 6 Isogai, T., Building a dynamic correlation network for fat-tailed financial asset returns, Applied Network Science (7):-24, 206, http://link.springer.com/article/0.007/s409-06-0008-x

More information

GLM PROC GLM y = Xβ + ε y X β ε ε σ 2 E[ε] = 0 var[ε] = σ 2 I σ 2 0 σ 2 =... 0 σ 2 σ 2 I ε σ 2 y E[y] =Xβ var[y] =σ 2 I PROC GLM

GLM PROC GLM y = Xβ + ε y X β ε ε σ 2 E[ε] = 0 var[ε] = σ 2 I σ 2 0 σ 2 =... 0 σ 2 σ 2 I ε σ 2 y E[y] =Xβ var[y] =σ 2 I PROC GLM PROC MIXED ( ) An Introdunction to PROC MIXED Junji Kishimoto SAS Institute Japan / Keio Univ. SFC / Univ. of Tokyo e-mail address: jpnjak@jpn.sas.com PROC MIXED PROC GLM PROC MIXED,,,, 1 1.1 PROC MIXED

More information

情報管理学科で学ぶ

情報管理学科で学ぶ 1/17 ` http://www.biwako.shiga-u.ac.jp/sensei/kumazawa/ 6............................................ 5 1............................... 1 1.1 I II III 1 1.2 2 1.3 2 2......................................

More information

読めば必ずわかる 分散分析の基礎 第2版

読めば必ずわかる 分散分析の基礎 第2版 2 2003 12 5 ( ) ( ) 2 I 3 1 3 2 2? 6 3 11 4? 12 II 14 5 15 6 16 7 17 8 19 9 21 10 22 11 F 25 12 : 1 26 3 I 1 17 11 x 1, x 2,, x n x( ) x = 1 n n i=1 x i 12 (SD ) x 1, x 2,, x n s 2 s 2 = 1 n n (x i x)

More information