3 HLM High School and Beyond HLM6 HLM6 C: Program Files HLM6S 2 C: Program MATHACH Files HLM6S Examples AppendxA school SECTOR Socio-Economic
|
|
|
- たつぞう しのしま
- 5 years ago
- Views:
Transcription
1 S. W. Raudenbush HLM6 student edition SAS/STAT MIXED R 2 HLM6 HLM HLM Hierarchical Linear A. S. Bryk S. W. Raudenbush Models HLM SSI *1 HLM6 student edition *2 student edition HLM6 (1) GUI (2) (3) okumurin@p.u-tokyo.ac.jp *1 *2
2 3 HLM High School and Beyond HLM6 HLM6 C: Program Files HLM6S 2 C: Program MATHACH Files HLM6S Examples AppendxA school SECTOR Socio-Economic Status: SES MEANSES MATHACH SES SECTOR MEANSES 2 Raudenbush and Bryk 1 (2002) MATHACH SECTOR SES MEANSES SES 1: 7185 i j SES SES.j j MATHACH SES MAT HACH ij = β 0j + β 1j (SES ij SES.j ) + r ij (1) 2: { SECTOR 1 (Catholic) SECT OR = 0 (public) MEANSES β 0j = γ 00 + γ 01 (SECT OR j ) + γ 02 (MEANSES j ) + u 0j (2) β 1j = γ 10 + γ 11 (SECT OR j ) + γ 12 (MEANSES j ) + u 1j (3) 2 2 u 0j, u 1j 1 2 SECTOR, MEANSES
3 HLM MDMT.mdmt MDM MDM.mdm STS.sts 3. HLM.hlm MLM.mlm HLM HLM ASCII 2. SPSS, SAS, SYSTAT, STATA ASCII HLM6 ID ID ID ID ID 3
4 3.3 ID ID 1. 1 (HSB1.DAT) 8 2 (HSB2.DAT) SPSS SPSS.sav 1 ASCII MDM ASCII ASCII 4
5 3.4 MDM 1 (HSB1.SAV) 2 (HSB2.SAV) ASCII 3.4 MDM HLM6 MDM Multivariate Data Matrix.mdm MDM MDM MDM HLM5 SSM Sufficient MDM Statistics Matrix 6 ASCII.dat SPSS.sav MDM 5
![3.4 MDM 3.4.1 HLM6 HLM6 HLM6 HLM6S WHLMS.exe WHLMS.exe [ ] [ ] HLM6](/docs-images/91/105372447/images/6-0.jpg "3.4.2 ASCII MDM ASCII MDM MDM 1 1 [File]-[Make new MDM file]-[ascii")
6 3.4 MDM HLM6 HLM6 HLM6 HLM6S WHLMS.exe WHLMS.exe [ ] [ ] HLM ASCII MDM ASCII MDM MDM 1 1 [File]-[Make new MDM file]-[ascii input] 6
![3.4 MDM [OK] 2 HLM2 Make MDM - HLM2 MDM MDM File Name MDM hsb-gd.](/docs-images/91/105372447/images/7-1.jpg "mdm Nesting of input data persons within groups Level-1 Specification 1 Level-1 File Name: 1 C: Program")
7 3.4 MDM [OK] 2 HLM2 Make MDM - HLM2 MDM MDM File Name MDM hsb-gd.mdm Nesting of input data persons within groups Level-1 Specification 1 Level-1 File Name: 1 C: Program Files HLM6S Examples Browse AppendxA hsb1.dat 7
![3.4 MDM Number of Variables: ID 4 MDM / Missing Data: Data Format: (A4,4F12.3) ( ) A ID A ID Labels: F 4F12.3 3 4 ID F [ ].[ ] hsb1.](/docs-images/91/105372447/images/8-0.jpg "dat X / 2 8X 8, ACH Level-2 Specification MINORITY, FEMALE, SES, MATH- 2 1 C: Program Files HLM6S Examples 2 AppendxA hsb2.dat 6 Data Format (A4,6F12.")
8 3.4 MDM Number of Variables: ID 4 MDM / Missing Data: Data Format: (A4,4F12.3) ( ) A ID A ID Labels: F 4F ID F [ ].[ ] hsb1.dat X / 2 8X 8, ACH Level-2 Specification MINORITY, FEMALE, SES, MATH- 2 1 C: Program Files HLM6S Examples 2 AppendxA hsb2.dat 6 Data Format (A4,6F12.3) Save mdmt file [Save mdmt file] MDMT MDM Template.mdmt SIZE SECTOR PRA- CAD DISCLIM HIM- MDMT MDM *3 INTY MEANSES *4 *3 MDM *4 MDMT MDM [File]-[Make new MDM from old MDM template (.mdmt) file] MDMT MDM 8
9 3.4 MDM Make MDM [Make MDM] MDM MDM Check Stats [Check Stats] STS hlm2mdm.sts STS Done [Done] MDM MDM [File]-[Display MDM stats] MDM 9
![file]-[stat package input]](/docs-images/91/105372447/images/10-2.jpg "[OK] ASCII 2 HLM2 Make MDM -")
10 3.4 MDM SPSS MDM SPSS.sav MDM [File]-[Make new MDM file]-[stat package input] [OK] ASCII 2 HLM2 Make MDM - HLM2 MDM 10
11 3.4 MDM MDM File Name MDM.mdm hsb-gd.mdm Input File Type SPSS/Windows Nesting of input data persons within groups Level-1 Specification 1 Browse: [Browse] Level-1 File C: Program Files HLM6S Examples Name: AppendxA hsb1.sav Choose Variabes: [Choose Variables] 1 SPSS ID MDM : Missing Data? Yes No Yes hsb1.sav MDM / Level-2 Specification 2 Browse: [Browse] 2 C: Program Files HLM6S Examples AppendxA hsb2.sav Choose Variables: [Choose Variables] 2 11
![3.4 MDM Save mdmt file [Save mdmt file] MDM hsb-gd.](/docs-images/91/105372447/images/12-0.jpg "mdmt Make MDM [Make MDM] MDM MDM Check Stats [Check Stats] STS")
12 3.4 MDM Save mdmt file [Save mdmt file] MDM hsb-gd.mdmt Make MDM [Make MDM] MDM MDM Check Stats [Check Stats] STS hlm2mdm.sts STS 12
![3.5 Done [Done] MDM 3.](/docs-images/91/105372447/images/13-0.jpg "5 MDM MDM [File]-[Create HLM6 (1) (3) a new")
![model using an existing GUI HLM6 MDM file] MDM](/docs-images/91/105372447/images/13-1.jpg "3.5.1 1 (1) MATHACH SES [Level-1] 1 1 MATHACH")
13 3.5 Done [Done] MDM 3.5 MDM MDM [File]-[Create HLM6 (1) (3) a new model using an existing GUI HLM6 MDM file] MDM (1) MATHACH SES [Level-1] 1 1 MATHACH Outcome variable MATHACH 13
14 3.5 SES add variable group centered β 0 2 SES MEANSES HLM6 1 SES Delete variable from model 2 [Level-2] 2 1 β 0 (2) 14
15 3.6 SECTOR MEANSES uncenterd β 1 (3) SECTOR MEANSES uncenterd 2 u 1 u Basic Settings [Basic Settings] [Basic Setting] [Level-1] [Outcome] 15
16 3.6 [Level-1 Residual File] 16
![3.6 Variables in residual file Residual File Type SPSS resfil1.sav Free Format.txt [OK] CSV.](/docs-images/91/105372447/images/17-0.jpg "csv 2 [Level-2 Residual File] SPSS Excel Variables in residual file 1 Residual File Type SPSS [OK] resfil1.sav resfil2.")
17 3.6 Variables in residual file Residual File Type SPSS resfil1.sav Free Format.txt [OK] CSV.csv 2 [Level-2 Residual File] SPSS Excel Variables in residual file 1 Residual File Type SPSS [OK] resfil1.sav resfil2.sav 2 Title no title Output file name HSB-GD Graph file name [OK] 17
![8 hsb-gd.hlm [Run Analysis] HLM 3.](/docs-images/91/105372447/images/18-1.jpg "9 [File]-[View Output] [Basing Settings]")
18 3.7 HLM Other Settings 3.7 HLM HLM.hlm HMLM.mlm [File]-[Save As] HLM 3.8 hsb-gd.hlm [Run Analysis] HLM 3.9 [File]-[View Output] [Basing Settings] hlm2.txt 18
19 3.9 SPECIFICATIONS FOR THIS HLM2 RUN Problem Title: HSB-GD MDM The data source for this run = hsb-gd.mdm.hlm The command file for this run = C:\Program Files\HLM6S\Examples\AppendxA\hsb-gd.hlm Output file name = C:\Program Files\HLM6S\Examples\AppendxA\hlm2.txt The maximum number of level-1 units = 7185 The maximum number of level-2 units = 160 The maximum number of iterations = 100 Method of estimation: restricted maximum likelihood Weighting Specification Weight Variable Weighting? Name Normalized? Level 1 no Level 2 no Precision no MATHACH The outcome variable is MATHACH The model specified for the fixed effects was: Level-1 Level-2 Coefficients Predictors INTRCPT1, B0 INTRCPT2, G00 SECTOR, G01 MEANSES, G02 * SES slope, B1 INTRCPT2, G10 SECTOR, G11 MEANSES, G12 19
20 3.9 SES * - This level-1 predictor has been centered around its group mean. The model specified for the covariance components was: Sigma squared (constant across level-2 units) Tau dimensions INTRCPT1 SES slope Summary of the model specified (in equation format) Level-1 Model Y = B0 + B1*(SES) + R Level-2 Model B0 = G00 + G01*(SECTOR) + G02*(MEANSES) + U0 B1 = G10 + G11*(SECTOR) + G12*(MEANSES) + U1 Iterations stopped due to small change in likelihood function ******* ITERATION 61 ******* 1 Sigma_squared = Tau INTRCPT1,B SES,B Tau (as correlations) INTRCPT1,B SES,B (Reliability) Random level-1 coefficient Reliability estimate INTRCPT1, B SES, B
21 3.9 2 The value of the likelihood function at iteration 61 = E The outcome variable is MATHACH Final estimation of fixed effects: Standard Approx. Fixed Effect Coefficient Error T-ratio d.f. P-value For INTRCPT1, B0 INTRCPT2, G SECTOR, G MEANSES, G For SES slope, B1 INTRCPT2, G SECTOR, G MEANSES, G The outcome variable is MATHACH Final estimation of fixed effects (with robust standard errors) Standard Approx. Fixed Effect Coefficient Error T-ratio d.f. P-value For INTRCPT1, B0 INTRCPT2, G SECTOR, G MEANSES, G For SES slope, B1 INTRCPT2, G SECTOR, G MEANSES, G Final estimation of variance components: Random Effect Standard Variance df Chi-square P-value Deviation Component INTRCPT1, U SES slope, U level-1, R
22 3.10 Deviance Deviance Statistics for current covariance components model Deviance = Number of estimated parameters = HLM6 (1) 1 SECTOR [File]-[Graph Equations]-[Level-1 equation graphing] 1 22
![3.10 X-focus Level-1: Number of groups: 1 SES groups All Z-focus SECTOR [OK]](/docs-images/91/105372447/images/23-0.jpg "Cathoric SECTOR=1 SECTOR 2 21 2 SECTOR INTRCPT1, B0 1.226384 SES slope, B1 1.")
23 3.10 X-focus Level-1: Number of groups: 1 SES groups All Z-focus SECTOR [OK] Cathoric SECTOR=1 SECTOR SECTOR INTRCPT1, B SES slope, B
24 HLM6 1 MDM l1resid 1 fitval 2 mathach 24
25 eb : ebintrcp 1 β 0j 1 empirical Bayes u 0j ebses 2 1 β 1j u 1j ol 2 ol : olintrcp olses 1 β 0j β 1j 2 ordinary least-squares u 0j u 1j ebintrcp olintrcp eb- eb ses olses fv : 1 2 ec : 1 ( fv ) ( eb 4 SAS PROC MIXED β 0j β 1j ) 23 HLM6 HLM6 1 SAS MIXED 29 R SAS PROC MIXED HLM6 25
26 4.1 SAS (2) (3) (1) MAT HACH ij = [γ 00 + γ 01 SECT OR j + γ 02 MEANSES j + γ 10 (SES ij SES.j )+ γ 11 SECT OR j (SES ij SES.j )+ γ 12 MEANSES j (SES ij SES.j )]+ [u 0j + u 1j (SES ij SES.j ) + r ij ] (4) [ ] 2 [ ] HLM6 SAS R (4) SAS PROC MIXED Singer (1998) HSB HLM6 4.1 SAS ***************************; * SAS HSB ; * By Taichi OKUMURA ; ***************************; OPTIONS ls = 80; /* */ DATA hsb; INFILE C:\Program Files\SAS Institute\SAS\V8\hsb-sas.dat ; INPUT id school ses mathach meanses cses sector; RUN; /* proc mixed */ PROC MIXED DATA=hsb NOCLPRINT COVTEST; CLASS school; MODEL mathach = sector meanses cses sector*cses meanses*cses / SOLUTION DDFM=bw NOTEST; RANDOM intercept cses / TYPE = un SUB=school; RUN; QUIT; MIXED PROC MIXED NOCLPRINT: COVTEST: school 26
27 4.2 SAS CLASS school MODEL SOLUTION: DDFM=bw: NOTEST: RANDOM TYPE = un: sub=school: school 4.2 SAS / / SAS The Mixed Procedure Model Information Data Set Dependent Variable Covariance Structure Subject Effect Estimation Method Residual Variance Method Fixed Effects SE Method Degrees of Freedom Method WORK.HSB mathach Unstructured school REML /* REML */ Profile Model-Based Between-Within Dimensions Covariance Parameters 4 Columns in X 6 Columns in Z Per Subject 2 Subjects 160 Max Obs Per Subject 67 Observations Used 7185 Observations Not Used 0 Total Observations
28 4.2 SAS Iteration History Iteration Evaluations -2 Res Log Like Criterion Convergence criteria met. Covariance Parameter Estimates Standard Z Cov Parm Subject Estimate Error Value Pr Z /* 2 */ UN(1,1) school <.0001 UN(2,1) school UN(2,2) school /* 1 */ Residual <.0001 Fit Statistics -2 Res Log Likelihood /* deviance */ AIC (smaller is better) AICC (smaller is better) BIC (smaller is better) Null Model Likelihood Ratio Test DF Chi-Square Pr > ChiSq <.0001 /* 2 */ Solution for Fixed Effects Standard Effect Estimate Error DF t Value Pr > t Intercept <.0001 /* G00 */ sector <.0001 /* G01 */ meanses <.0001 /* G02 */ cses <.0001 /* G10 */ sector*cses <.0001 /* G11 */ meanses*cses /* G12 */ 28
29 5 R R (R Development Core Team, 2005) nmle lme *5 nlme HSB 5.1 SAS HSB 10 school ses mathach meanses cses sector Public Public Public Public Public Public Public Public Public Public 5.2 R R HSB By Taichi OKUMURA nlme library(nlme) 1 data(mathachieve) MathAchieve[1:10,] 2 data(mathachschool) MathAchSchool[1:10,] SES mses attach(mathachieve) *5 29
30 5.3 R mses <- tapply(ses, School, mean) detach(mathachieve) HSB HSB <- as.data.frame(mathachieve[, c("school", "SES", "MathAch")]) names(hsb) <- c("school", "ses", "mathach") SES mses HSB$meanses <- mses[as.character(hsb$school)] HSB$cses <- HSB$ses - HSB$meanses sector HSB sector <- MathAchSchool$Sector names(sector) <- row.names(mathachschool) HSB$sector <- sector[as.character(hsb$school)] HSB$sector <- factor(hsb$sector, levels=c( Public, Catholic )) result.hsb result.hsb <- lme(mathach ~ sector + meanses + cses + sector*cses + meanses*cses, random = ~ cses school, data = HSB) summary(result.hsb) 5.3 R / / Linear mixed-effects model fit by REML /* REML */ Data: HSB AIC BIC loglik /* LogLik -2 deviance */ Random effects: Formula: ~cses school Structure: General positive-definite, Log-Cholesky parametrization StdDev Corr (Intercept) (Intr) /* 2 */ cses /* 2 */ Residual /* 1 */ /* 2 */ Fixed effects: mathach ~ sector + meanses + cses + sector * cses + meanses * cses Value Std.Error DF t-value p-value (Intercept) e+00 /* G00 */ sectorcatholic e-04 /* G01 */ meanses e+00 /* G02 */ cses e+00 /* G10 */ sectorcatholic:cses e+00 /* G11 */ meanses:cses e-04 /* G12 */ Correlation: (Intr) sctrct meanss cses sctrc: sectorcatholic meanses cses
31 REFERENCES sectorcatholic:cses meanses:cses Standardized Within-Group Residuals: Min Q1 Med Q3 Max Number of Observations: 7185 /* */ Number of Groups: 160 /* */ References R Development Core Team. (2005). R: A language and environment for statistical computing. Vienna, Austria. (ISBN ) Raudenbush, S. W., & Bryk, A. S. (2002). Hierarchical linear models: Applications and data analysis methods (2nd ed.). Sage. Singer, J. D. (1998). Using SAS PROC MIXED to fit multilevel models, hierarchical models, and individual growth models. Journal of Educational and Behavioral Statistics, 24,
: (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.
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
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
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.
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,
4.9 Hausman Test Time Fixed Effects Model vs Time Random Effects Model Two-way Fixed Effects Model
1 EViews 5 2007 7 11 2010 5 17 1 ( ) 3 1.1........................................... 4 1.2................................... 9 2 11 3 14 3.1 Pooled OLS.............................................. 14
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](/thumbs/84/89087584.jpg "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
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)
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......................
1 I EViews View Proc Freeze
EViews 2017 9 6 1 I EViews 4 1 5 2 10 3 13 4 16 4.1 View.......................................... 17 4.2 Proc.......................................... 22 4.3 Freeze & Name....................................
講義のーと : データ解析のための統計モデリング. 第５回
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
第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
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
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....................
講義のーと : データ解析のための統計モデリング. 第３回
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
% 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.
10
z c j = N 1 N t= j1 [ ( z t z ) ( )] z t j z q 2 1 2 r j /N j=1 1/ N J Q = N(N 2) 1 N j j=1 r j 2 2 χ J B d z t = z t d (1 B) 2 z t = (z t z t 1 ) (z t 1 z t 2 ) (1 B s )z t = z t z t s _ARIMA CONSUME
最小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) ( )
Microsoft PowerPoint - GLMMexample_ver pptx
Linear Mixed Model ( 以下 混合モデル ) の短い解説 この解説のPDFは http://www.lowtem.hokudai.ac.jp/plantecol/akihiro/sumida-index.html の お勉強 のページにあります. ver 20121121 と との間に次のような関係が見つかったとしよう 全体的な傾向に対する回帰直線を点線で示した ところが これらのデータは実は異なる
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
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:
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
5.2 White
1 EViews 1 : 2007/5/15 2007/5/25 1 EViews 4 2 ( 6 2.1............................................ 6 2.2 Workfile............................................ 7 2.3 Workfile............................................
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
1 EViews 2 2007/5/17 2007/5/21 4 OLS 2 4.1.............................................. 2 4.2................................................ 9 4.3.............................................. 11 4.4
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
80 X 1, X 2,, X n ( λ ) λ P(X = x) = f (x; λ) = λx e λ, x = 0, 1, 2, x! l(λ) = n f (x i ; λ) = n λ x i e λ x i! = λ n x i e nλ n x i! n n log l(λ) = log(λ) x i nλ log( x i!) log l(λ) λ = 1 λ n x i n =
<4D F736F F D20939D8C7689F090CD985F93C18EEA8D758B E646F63>
Gretl OLS omitted variable omitted variable AIC,BIC a) gretl gretl sample file Greene greene8_3 Add Define new variable l_g_percapita=log(g/pop) Pg,Y,Pnc,Puc,Ppt,Pd,Pn,Ps Add logs of selected variables
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)
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 (
![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 (](/thumbs/82/85163103.jpg "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 (")
7 2 2008 7 10 1 2 2 1.1 2............................................. 2 1.2 2.......................................... 2 1.3 2........................................ 3 1.4................................................
一般化線形 (混合) モデル (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!
こんにちは由美子です
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
Studies of Foot Form for Footwear Design (Part 9) : Characteristics of the Foot Form of Young and Elder Women Based on their Sizes of Ball Joint Girth
Studies of Foot Form for Footwear Design (Part 9) : Characteristics of the Foot Form of Young and Elder Women Based on their Sizes of Ball Joint Girth and Foot Breadth Akiko Yamamoto Fukuoka Women's University,
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
87 6.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, V(y t y t 1, y t 2, ) = σ 2 3. Thus, y t y t 1,
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
第13回：交差項を含む回帰・弾力性の推定
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 = β
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
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................................
kubostat2017b p.1 agenda I 2017 (b) probability distribution and maximum likelihood estimation :
kubostat2017b p.1 agenda I 2017 (b) probabilit distribution and maimum likelihood estimation kubo@ees.hokudai.ac.jp http://goo.gl/76c4i 2017 11 14 : 2017 11 07 15:43 1 : 2 3? 4 kubostat2017b (http://goo.gl/76c4i)
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,
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
!!! 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
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
ECCS. ECCS,. ( 2. Mac Do-file Editor. Mac Do-file Editor Windows Do-file Editor Top Do-file e
1 1 2015 4 6 1. ECCS. ECCS,. (https://ras.ecc.u-tokyo.ac.jp/guacamole/) 2. Mac Do-file Editor. Mac Do-file Editor Windows Do-file Editor Top Do-file editor, Do View Do-file Editor Execute(do). 3. Mac System
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
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
(missing data analysis) - - 1/16/2011 (missing data, missing value) (list-wise deletion) (pair-wise deletion) (full information maximum likelihood method, FIML) (multiple imputation method) 1 missing completely
σ 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)
,, Poisson 3 3. t t y,, y n Nµ, σ 2 y i µ + ɛ i ɛ i N0, σ 2 E[y i ] µ * i y i x i y i α + βx i + ɛ i ɛ i N0, σ 2, α, β *3 y i E[y i ] α + βx i
![,, Poisson 3 3. t t y,, y n Nµ, σ 2 y i µ + ɛ i ɛ i N0, σ 2 E[y i ] µ * i y i x i y i α + βx i + ɛ i ɛ i N0, σ 2, α, β *3 y i E[y i ] α + βx i](/thumbs/87/95421301.jpg ",, Poisson 3 3. t t y,, y n Nµ, σ 2 y i µ + ɛ i ɛ i N0, σ 2 E[y i ] µ * i y i x i y i α + βx i + ɛ i ɛ i N0, σ 2, α, β *3 y i E[y i ] α + βx i")
Armitage.? SAS.2 µ, µ 2, µ 3 a, a 2, a 3 a µ + a 2 µ 2 + a 3 µ 3 µ, µ 2, µ 3 µ, µ 2, µ 3 log a, a 2, a 3 a µ + a 2 µ 2 + a 3 µ 3 µ, µ 2, µ 3 * 2 2. y t y y y Poisson y * ,, Poisson 3 3. t t y,, y n Nµ,
untitled
IT (1, horiike@ml.me.titech.ac.jp) (1, jun-jun@ms.kagu.tus.ac.jp) 1. 1-1 19802000 2000ITIT IT IT TOPIX (%) 1TOPIX 2 1-2. 80 80 ( ) 2004/11/26 S-PLUS 2 1-3. IT IT IT IT 2. 2-1. a. b. (Size) c. B/M(Book
第2回：データの加工・整理
2 2018 4 13 1 / 24 1. 2. Excel 3. Stata 4. Stata 5. Stata 2 / 24 1 cross section data e.g., 47 2009 time series data e.g., 1999 2014 5 panel data e.g., 47 1999 2014 5 3 / 24 micro data aggregate data 4
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
TS001 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 postestimation 49 mwp-055 corrgram/ac/pac 56 mwp-009 dfgls
(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
/22 R MCMC R R MCMC? 3. Gibbs sampler : kubo/
2006-12-09 1/22 R MCMC R 1. 2. R MCMC? 3. Gibbs sampler : kubo@ees.hokudai.ac.jp http://hosho.ees.hokudai.ac.jp/ kubo/ 2006-12-09 2/22 : ( ) : : ( ) : (?) community ( ) 2006-12-09 3/22 :? 1. ( ) 2. ( )
Introduction Purpose This training course demonstrates the use of the High-performance Embedded Workshop (HEW), a key tool for developing software for
Introduction Purpose This training course demonstrates the use of the High-performance Embedded Workshop (HEW), a key tool for developing software for embedded systems that use microcontrollers (MCUs)
.. est table TwoSLS1 TwoSLS2 GMM het,b(%9.5f) se Variable TwoSLS1 TwoSLS2 GMM_het hi_empunion totchr
3,. Cameron and Trivedi (2010) Microeconometrics Using Stata, Revised Edition, Stata Press 6 Linear instrumentalvariables regression 9 Linear panel-data models: Extensions.. GMM xtabond., GMM(Generalized
kubostat7f p GLM! logistic regression as usual? N? GLM GLM doesn t work! GLM!! probabilit distribution binomial distribution : : β + β x i link functi
kubostat7f p statistaical models appeared in the class 7 (f) kubo@eeshokudaiacjp https://googl/z9cjy 7 : 7 : The development of linear models Hierarchical Baesian Model Be more flexible Generalized Linear
BR001
BR001 Stata 11 Stata Stata11 whitepaper mwp 3 mwp-027 22 mwp-028 / 40 mwp-001 logistic/logit 50 mwp-039 logistic/logit postestimation 60 mwp-040 margins 74 mwp-029 regress 90 mwp-037 regress postestimation
現代日本論演習／比較現代日本論研究演習I「統計分析の基礎」
URL: http://tsigeto.info/statg/ I ( ) 3 2017 2 ( 7F) 1 : (1) ; (2) 1998 (70 20% 6 8 ) (30%) ( 2) ( 2) 2 1. (4/13) 2. SPSS (4/20) 3. (4/27) [ ] 4. (5/11 6/1) [1, 4 ] 5. (6/8) 6. (6/15 6/29) [2, 5 ] 7. (7/6
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
Autumn 2005 1 9 13 14 16 16 DATA _null_; SET sashelp.class END=eof; FILE 'C: MyFiles class.txt'; /* */ PUT name sex age; IF eof THEN DO; FILE LOG; /* */ PUT '*** ' _n_ ' ***'; END; DATA _null_;
Microsoft Word - 計量研修テキスト_第5版）.doc
Q3-1-1 テキスト P59 10.8.3.2.1.0 -.1 -.2 10.4 10.0 9.6 9.2 8.8 -.3 76 78 80 82 84 86 88 90 92 94 96 98 R e s i d u al A c tual Fi tte d Dependent Variable: LOG(TAXH) Date: 10/26/05 Time: 15:42 Sample: 1975
Introduction Purpose This training course describes the configuration and session features of the High-performance Embedded Workshop (HEW), a key tool
Introduction Purpose This training course describes the configuration and session features of the High-performance Embedded Workshop (HEW), a key tool for developing software for embedded systems that
ブック
ARMA Estimation on Process of ARMA Time Series Model Sanno University Bulletin Vol.26 No. 2 February 2006 ARMA Estimation on Process of ARMA Time Series Model Many papers and books have been published
t sex N y y y Diff (1-2)
Armitage 1 1.1 2 t 1.2 SAS Proc GLM 2 2.1 1 1 2.1.1 50 1 1 t sex N y 50 116.45 119.6 122.75 11.071 1.5657 93.906 154.32 y 50 127.27 130.7 134.13 12.072 1.7073 102.68 163.37 y Diff (1-2) -15.7-11.1-6.504
Presentation Title Goes Here
SAS 9: (reprise) SAS Institute Japan Copyright 2004, SAS Institute Inc. All rights reserved. Greetings, SAS 9 SAS 9.1.3 Copyright 2004, SAS Institute Inc. All rights reserved. 2 Informations of SAS 9 SAS
H22 BioS t (i) treat1 treat2 data d1; input patno treat1 treat2; cards; ; run; 1 (i) treat = 1 treat =
H BioS t (i) treat treat data d; input patno treat treat; cards; 3 8 7 4 8 8 5 5 6 3 ; run; (i) treat treat data d; input group patno period treat y; label group patno period ; cards; 3 8 3 7 4 8 4 8 5
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
4 2 p = p(t, g) (1) r = r(t, g) (2) p r t g p r dp dt = p dg t + p g (3) dt dr dt = r dg t + r g dt 3 p t p g dt p t 3 2 4 r t = 3 4 2 Benefit view dp
( ) 62 1 1 47 2 3 47 2 e-mail:miyazaki@ngu.ac.jp 1 2000 2005 1 4 2 p = p(t, g) (1) r = r(t, g) (2) p r t g p r dp dt = p dg t + p g (3) dt dr dt = r dg t + r g dt 3 p t p g dt p t 3 2 4 r t = 3 4 2 Benefit
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
¥¤¥ó¥¿¡¼¥Í¥Ã¥È·×¬¤È¥Ç¡¼¥¿²òÀÏ Âè2²ó
2 2015 4 20 1 (4/13) : ruby 2 / 49 2 ( ) : gnuplot 3 / 49 1 1 2014 6 IIJ / 4 / 49 1 ( ) / 5 / 49 ( ) 6 / 49 (summary statistics) : (mean) (median) (mode) : (range) (variance) (standard deviation) 7 / 49
ｽﾗｲﾄﾞ ﾀｲﾄﾙなし
LightCycler Software Ver.3.5 : 200206 1/30 Windows NT Windows NT Ctrl + Alt + Delete LightCycler 3 Front Screen 2/30 LightCycler3 Front RUN Data Analysis LightCycler Data Analysis Edit Graphics Defaults
9 8 7 (x-1.0)*(x-1.0) *(x-1.0) (a) f(a) (b) f(a) Figure 1: f(a) a =1.0 (1) a 1.0 f(1.0)
E-mail: takio-kurita@aist.go.jp 1 ( ) CPU ( ) 2 1. a f(a) =(a 1.0) 2 (1) a ( ) 1(a) f(a) a (1) a f(a) a =2(a 1.0) (2) 2 0 a f(a) a =2(a 1.0) = 0 (3) 1 9 8 7 (x-1.0)*(x-1.0) 6 4 2.0*(x-1.0) 6 2 5 4 0 3-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)
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
Microsoft Word - Win-Outlook.docx
Microsoft Office Outlook での設定方法 (IMAP および POP 編 ) How to set up with Microsoft Office Outlook (IMAP and POP) 0. 事前に https://office365.iii.kyushu-u.ac.jp/login からサインインし 以下の手順で自分の基本アドレスをメモしておいてください Sign
Attendance Demand for J-League õ Shinsuke KAWAI* and Takeo HIRATA* Abstract The purpose of this study was to clarify the variables determining the attendance in J-league matches, using the 2,699 J-league
1 環境統計学ぷらす 第 5 回 一般 ( 化 ) 線形混合モデル 高木俊 2013/11/21
1 環境統計学ぷらす 第 5 回 一般 ( 化 ) 線形混合モデル 高木俊 shun.takagi@sci.toho-u.ac.jp 2013/11/21 2 予定 第 1 回 : Rの基礎と仮説検定 第 2 回 : 分散分析と回帰 第 3 回 : 一般線形モデル 交互作用 第 4.1 回 : 一般化線形モデル 第 4.2 回 : モデル選択 (11/29?) 第 5 回 : 一般化線形混合モデル
A Nutritional Study of Anemia in Pregnancy Hematologic Characteristics in Pregnancy (Part 1) Keizo Shiraki, Fumiko Hisaoka Department of Nutrition, Sc
A Nutritional Study of Anemia in Pregnancy Hematologic Characteristics in Pregnancy (Part 1) Keizo Shiraki, Fumiko Hisaoka Department of Nutrition, School of Medicine, Tokushima University, Tokushima Fetal
(X) (Y ) Y = intercept + c X + e (1) e c c M = intercept + ax + e (2) a Y = intercept + cx + bm + e (3) (1) X c c c (3) b X M Y (indirect effect) a b
21 12 23 (mediation analysis) Figure 1 X Y M (mediator) mediation model Baron and Kenny (1986) 1 1) mediated moderation ( moderated mediation) 2) (multilevel mediation model) a M b X c (c ) Y 1: 1 1.1
Solution Report
CGE 3 GAMS * Date: 2018/07/24, Version 1.1 1 2 2 GAMSIDE 3 2.1 GAMS................................. 3 2.2 GAMSIDE................................ 3 2.3 GAMSIDE............................. 7 3 GAMS 11
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.)
Validation of a Food Frequency Questionnaire Based on Food Groups for Estimating Individual Nutrient Intake Keiko Takahashi *', Yukio Yoshimura *', Ta
Validation of a Food Frequency Questionnaire Based on Food Groups for Estimating Individual Nutrient Intake Keiko Takahashi *', Yukio Yoshimura *', Tae Kaimoto *', Daisuke Kunii *', Tatsushi Komatsu *2
udc-2.dvi
13 0.5 2 0.5 2 1 15 2001 16 2009 12 18 14 No.39, 2010 8 2009b 2009a Web Web Q&A 2006 2007a20082009 2007b200720082009 20072008 2009 2009 15 1 2 2 2.1 18 21 1 4 2 3 1(a) 1(b) 1(c) 1(d) 1) 18 16 17 21 10
Mantel-Haenszelの方法
Mantel-Haenszel 2008 6 12 ) 2008 6 12 1 / 39 Mantel & Haenzel 1959) Mantel N, Haenszel W. Statistical aspects of the analysis of data from retrospective studies of disease. J. Nat. Cancer Inst. 1959; 224):
Title 生活年令による学級の等質化に関する研究 (1) - 生活年令と学業成績について - Author(s) 与那嶺, 松助 ; 東江, 康治 Citation 研究集録 (5): 33-47 Issue Date 1961-12 URL http://hdl.handle.net/20.500.12000/ Rights 46 STUDIES ON HOMOGENEOUS
13 Student Software TI-Nspire CX CAS TI Web TI-Nspire CX CAS Student Software ( ) 1 Student Software 37 Student Software Nspire Nspire Nspir
13 Student Software TI-Nspire CX CAS TI Web TI-Nspire CX CAS Student Software ( ) 1 Student Software 37 Student Software 37.1 37.1 Nspire Nspire Nspire 37.1: Student Software 13 2 13 Student Software esc
ＪＭＰ V4 による生存時間分析
V4 1 SAS 2000.11.18 4 ( ) (Survival Time) 1 (Event) Start of Study Start of Observation Died Died Died Lost End Time Censor Died Died Censor Died Time Start of Study End Start of Observation Censor
計量経済分析 2011 年度夏学期期末試験 担当 : 別所俊一郎 以下のすべてに答えなさい. 回答は日本語か英語でおこなうこと. 1. 次のそれぞれの記述が正しいかどうか判定し, 誤りである場合には理由, あるいはより適切な 記述はどのようなものかを述べなさい. (1) You have to wo
計量経済分析 2011 年度夏学期期末試験 担当 : 別所俊一郎 以下のすべてに答えなさい. 回答は日本語か英語でおこなうこと. 1. 次のそれぞれの記述が正しいかどうか判定し, 誤りである場合には理由, あるいはより適切な 記述はどのようなものかを述べなさい. (1) You have to worry about perfect multicollinearity in the multiple
7 ( 7 ( Workfile Excel hatuden 1000kWh kion_average kion_max kion_min date holiday *1 obon 7.1 Workfile 1. Workfile File - New -
1 EViews 4 2007 7 4 7 ( 2 7.1 Workfile............................................ 2 7.2........................................... 4 8 6 8.1................................................. 6 8.2................................................
浜松医科大学紀要
On the Statistical Bias Found in the Horse Racing Data (1) Akio NODA Mathematics Abstract: The purpose of the present paper is to report what type of statistical bias the author has found in the horse
15 2004.03 194
The Statistical Processing using EXCEL MIYOSHI Yoshihiko In this paper, I summarize the method of performing statistical processing using only the basic function of EXCEL without the VBA macro, add-in
Cleaner XL 1.5 クイックインストールガイド
Autodesk Cleaner XL 1.5 Contents Cleaner XL 1.5 2 1. Cleaner XL 3 2. Cleaner XL 9 3. Cleaner XL 12 4. Cleaner XL 16 5. 32 2 1. Cleaner XL 1. Cleaner XL Cleaner XL Administrators Cleaner XL Windows Media
現代日本論演習／比較現代日本論研究演習I「統計分析の基礎」
URL: http://tsigeto.info/statg/ I () 3 2016 2 ( 7F) 1 : (1); (2) 1998 (70 20% 6 9 ) (30%) ( 2) ( 2) 2 1. (4/14) 2. SPSS (4/21) 3. (4/28) [] 4. (5/126/2) [1, 4] 5. (6/9) 6. (6/166/30) [2, 5] 7. (7/78/4)
わが国企業による資金調達方法の選択問題
* 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
1 Stata SEM LightStone 1 5 SEM Stata Alan C. Acock, Discovering Structural Equation Modeling Using Stata, Revised Edition, Stata Press. Introduc
1 Stata SEM LightStone 1 5 SEM Stata Alan C. Acock, 2013. Discovering Structural Equation Modeling Using Stata, Revised Edition, Stata Press. Introduction to confirmatory factor analysis 9 Stata14 2 1
dvi
2017 65 2 185 200 2017 1 2 2016 12 28 2017 5 17 5 24 PITCHf/x PITCHf/x PITCHf/x MLB 2014 PITCHf/x 1. 1 223 8522 3 14 1 2 223 8522 3 14 1 186 65 2 2017 PITCHf/x 1.1 PITCHf/x PITCHf/x SPORTVISION MLB 30
_念3）医療2009_夏.indd
Evaluation of the Social Benefits of the Regional Medical System Based on Land Price Information -A Hedonic Valuation of the Sense of Relief Provided by Health Care Facilities- Takuma Sugahara Ph.D. Abstract
(3) (5) (6) 162 青 山 充 て ん層 厚 さ と流 量 計 水 柱 高 さの 関係 を 示 した 付 属 チ ャ ー トに よ り直 ち に平 均 粒子 径 を知 る よ う にな って い る. しか し,試 料 重 量 を 変 え て実 験 す る場 合,チ ρp/w=1を ャー
161 Tadaaki Aoyama: Measurement of Mean Particle Size of Tungsten and Tungsten Carbide Powders by Fisher Sub-Sieve Sizer not Dependent on Porosity. An equation for the specific surface area of powder,
LM35 高精度・摂氏直読温度センサIC
Precision Centigrade Temperature Sensors Literature Number: JAJSB56 IC A IC D IC IC ( ) IC ( K) 1/4 55 150 3/4 60 A 0.1 55 150 C 40 110 ( 10 ) TO-46 C CA D TO-92 C IC CA IC 19831026 24120 11800 ds005516
Microsoft Word - 計量研修テキスト_第5版）.doc
Q8-1 テキスト P131 Engle-Granger 検定 Dependent Variable: RM2 Date: 11/04/05 Time: 15:15 Sample: 1967Q1 1999Q1 Included observations: 129 RGDP 0.012792 0.000194 65.92203 0.0000 R -95.45715 11.33648-8.420349
Microsoft Word - Meta70_Preferences.doc
Image Windows Preferences Edit, Preferences MetaMorph, MetaVue Image Windows Preferences Edit, Preferences Image Windows Preferences 1. Windows Image Placement: Acquire Overlay at Top Left Corner: 1 Acquire
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
(GSYS2.4) GSYS2.4 Manual SUZUKI Ryusuke Hokkaido University Hospital Abstract GSYS2.4 is an update version of GSYS version 2. Main features added in this version are Magnifying glass function, Automatically