yamadaiR(cEFA).pdf
|
|
- さやな くぬぎ
- 5 years ago
- Views:
Transcription
1 R 2012/10/05 Kosugi,E.Koji (Yamadai.R) Categorical Factor Analysis by using R 2012/10/05 1 / 9
2 Why we use Kosugi,E.Koji (Yamadai.R) Categorical Factor Analysis by using R 2012/10/05 2 / 9
3 FA vs categorical FA, Kosugi,E.Koji (Yamadai.R) Categorical Factor Analysis by using R 2012/10/05 3 / 9
4 One of the reasons 3 7 Kosugi,E.Koji (Yamadai.R) Categorical Factor Analysis by using R 2012/10/05 4 / 9
5 (2002) 1 2 polychoric correlation coefficient (polyserial correlation coefficient) 3 Kosugi,E.Koji (Yamadai.R) Categorical Factor Analysis by using R 2012/10/05 5 / 9
6 Polychoric Correlation Polyserial Correlation Tetrachoric Correlation Kosugi,E.Koji (Yamadai.R) Categorical Factor Analysis by using R 2012/10/05 6 / 9
7 images of latent continuity Figure : image of latent continuity and expression x ξ x ξ x = 1 ξ < a 1 x = 2 a 1 ξ < a 2 x = 3 a 2 ξ < a 3. x = s. a s 1 ξ (1) Kosugi,E.Koji (Yamadai.R) Categorical Factor Analysis by using R 2012/10/05 7 / 9
8 ρ X,Y (2step-ML) Kosugi,E.Koji (Yamadai.R) Categorical Factor Analysis by using R 2012/10/05 8 / 9
9 Follow me with R code... Kosugi,E.Koji (Yamadai.R) Categorical Factor Analysis by using R 2012/10/05 9 / 9
10 > library(psych) > library(polycor) > # sample statistics > sample <- read.csv("cefasample.csv",head=f,na.strings="*") > head(sample) V1 V2 V3 V4 V5 V6 V7 V > summary(sample) V1 V2 V3 V4 Min. :1.000 Min. :1.000 Min. :1.000 Min. : st Qu.: st Qu.: st Qu.: st Qu.:3.000 Median :4.000 Median :4.000 Median :4.000 Median :4.000 Mean :3.913 Mean :4.127 Mean :3.901 Mean : rd Qu.: rd Qu.: rd Qu.: rd Qu.:4.000 Max. :5.000 Max. :5.000 Max. :5.000 Max. :5.000 NA's :2 NA's :1 V5 V6 V7 V8 Min. :1.000 Min. :1.000 Min. :1.00 Min. : st Qu.: st Qu.: st Qu.:2.00 1st Qu.:3.000 Median :4.000 Median :3.000 Median :3.00 Median :4.000 Mean :3.955 Mean :3.138 Mean :2.78 Mean : rd Qu.: rd Qu.: rd Qu.:4.00 3rd Qu.:4.000 Max. :5.000 Max. :5.000 Max. :5.00 Max. :5.000 NA's :1 > table(sample$v1) > describe(sample) var n mean sd median trimmed mad min max range skew kurtosis se V V V V V V V V
11 > # peason cor > peason.cor <- cor(sample,use="complete.obs") > print(peason.cor,digit=2) V1 V2 V3 V4 V5 V6 V7 V8 V V V V V V V V > # polychoric cor > polychoric.cor <- polychoric(sample) > print(polychoric.cor$rho) V1 V2 V3 V4 V5 V6 V7 V V V V V V V V V8 V V V V V V V V > # > # compare, peason vs polycor > # > > # FA > fa.parallel(peason.cor,n.obs=355) Parallel analysis suggests that the number of factors = 3 and the number of components = 2
12 > fa.parallel(polychoric.cor$rho,n.obs=355) Parallel analysis suggests that the number of factors = 3 and the number of components = > fa.result.peason <- fa(peason.cor,n.obs=355,fm="gls",nfactors=3,rotate="promax") > fa.result.polych <- fa(polychoric.cor$rho,n.obs=355,fm="gls",nfactors=3,rotate="promax") > print(fa.result.peason,digit=3,sort=t) Factor Analysis using method = gls Call: fa(r = peason.cor, nfactors = 3, n.obs = 355, rotate = "promax", fm = "gls") Standardized loadings (pattern matrix) based upon correlation matrix item GLS2 GLS1 GLS3 h2 u2 V V V V V V V V GLS2 GLS1 GLS3 SS loadings Proportion Var Cumulative Var Proportion Explained Cumulative Proportion With factor correlations of GLS2 GLS1 GLS3 GLS GLS GLS Test of the hypothesis that 3 factors are sufficient. The degrees of freedom for the null model are 28 and the objective function was wit The degrees of freedom for the model are 7 and the objective function was 0.03 The root mean square of the residuals (RMSR) is The df corrected root mean square of the residuals is 0.04 The number of observations was 355 with Chi Square = with prob < Tucker Lewis Index of factoring reliability = RMSEA index = and the 90 % confidence intervals are NA BIC =
13 Fit based upon off diagonal values = Measures of factor score adequacy GLS2 GLS1 GLS3 Correlation of scores with factors Multiple R square of scores with factors Minimum correlation of possible factor scores > print(fa.result.polych,digit=3,sort=t) Factor Analysis using method = gls Call: fa(r = polychoric.cor$rho, nfactors = 3, n.obs = 355, rotate = "promax", fm = "gls") Standardized loadings (pattern matrix) based upon correlation matrix item GLS3 GLS1 GLS2 h2 u2 V V V V V V V V GLS3 GLS1 GLS2 SS loadings Proportion Var Cumulative Var Proportion Explained Cumulative Proportion With factor correlations of GLS3 GLS1 GLS2 GLS GLS GLS Test of the hypothesis that 3 factors are sufficient. The degrees of freedom for the null model are 28 and the objective function was wit The degrees of freedom for the model are 7 and the objective function was The root mean square of the residuals (RMSR) is The df corrected root mean square of the residuals is The number of observations was 355 with Chi Square = with prob < Tucker Lewis Index of factoring reliability =
14 RMSEA index = and the 90 % confidence intervals are BIC = Fit based upon off diagonal values = Measures of factor score adequacy GLS3 GLS1 GLS2 Correlation of scores with factors Multiple R square of scores with factors Minimum correlation of possible factor scores > # > # sample <- subset(sample,select=c("v11","v13","v20","v5","v4","v17","v12","v15")) > # write.table(sample,"cefasample.csv",sep=",",row.name=f,col.name=f,na="*") > > > # mixed pattern > sample.cat <- data.frame(lapply(sample[1:3],factor),sample[4:8]) > summary(sample.cat) V1 V2 V3 V4 V5 V6 1: 2 1: 3 1 : 2 Min. :1.000 Min. :1.000 Min. : : 26 2: 13 2 : 18 1st Qu.: st Qu.: st Qu.: : 61 3: 31 3 : 60 Median :4.000 Median :4.000 Median : :178 4:197 4 :206 Mean :3.853 Mean :3.955 Mean : : 88 5:111 5 : 67 3rd Qu.: rd Qu.: rd Qu.:4.000 NA's: 2 Max. :5.000 Max. :5.000 Max. :5.000 NA's :1 V7 V8 Min. :1.00 Min. : st Qu.:2.00 1st Qu.:3.000 Median :3.00 Median :4.000 Mean :2.78 Mean : rd Qu.:4.00 3rd Qu.:4.000 Max. :5.00 Max. :5.000 NA's :1 > hetcor.cor <- hetcor(sample.cat) > hetcor.cor$correlations V1 V2 V3 V4 V5 V6 V7 V V V V V V V V
15 V8 V V V V V V V V > hetcor.cor$type [,1] [,2] [,3] [,4] [,5] [1,] "" "Polychoric" "Polychoric" "Polyserial" "Polyserial" [2,] "Polychoric" "" "Polychoric" "Polyserial" "Polyserial" [3,] "Polychoric" "Polychoric" "" "Polyserial" "Polyserial" [4,] "Polyserial" "Polyserial" "Polyserial" "" "Pearson" [5,] "Polyserial" "Polyserial" "Polyserial" "Pearson" "" [6,] "Polyserial" "Polyserial" "Polyserial" "Pearson" "Pearson" [7,] "Polyserial" "Polyserial" "Polyserial" "Pearson" "Pearson" [8,] "Polyserial" "Polyserial" "Polyserial" "Pearson" "Pearson" [,6] [,7] [,8] [1,] "Polyserial" "Polyserial" "Polyserial" [2,] "Polyserial" "Polyserial" "Polyserial" [3,] "Polyserial" "Polyserial" "Polyserial" [4,] "Pearson" "Pearson" "Pearson" [5,] "Pearson" "Pearson" "Pearson" [6,] "" "Pearson" "Pearson" [7,] "Pearson" "" "Pearson" [8,] "Pearson" "Pearson" "" > fa.parallel(hetcor.cor$correlations,n.obs=355) Parallel analysis suggests that the number of factors = 3 and the number of components = > fa.result.hetcor <- fa(hetcor.cor$correlations,n.obs=355,fm="gls",nfactors=3,rotate="proma > print(fa.result.hetcor,digit=3,sort=t) Factor Analysis using method = gls Call: fa(r = hetcor.cor$correlations, nfactors = 3, n.obs = 355, rotate = "promax", fm = "gls") Standardized loadings (pattern matrix) based upon correlation matrix item GLS1 GLS2 GLS3 h2 u2 V V V V
16 V V V V GLS1 GLS2 GLS3 SS loadings Proportion Var Cumulative Var Proportion Explained Cumulative Proportion With factor correlations of GLS1 GLS2 GLS3 GLS GLS GLS Test of the hypothesis that 3 factors are sufficient. The degrees of freedom for the null model are 28 and the objective function was wit The degrees of freedom for the model are 7 and the objective function was The root mean square of the residuals (RMSR) is The df corrected root mean square of the residuals is The number of observations was 355 with Chi Square = with prob < Tucker Lewis Index of factoring reliability = RMSEA index = and the 90 % confidence intervals are BIC = Fit based upon off diagonal values = Measures of factor score adequacy GLS1 GLS2 GLS3 Correlation of scores with factors Multiple R square of scores with factors Minimum correlation of possible factor scores > 7
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 informationDAA09
> 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 information1 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( ) 1.1 Polychoric Correlation Polyserial Correlation Graded Response Model Partial Credit Model Tetrachoric Correlation ( ) 2 x y x y s r 1 x 2
1 (,2007) SPSSver8 1997 (2002) 1. 2. polychoric correlation coefficient (polyserial correlation coefficient) 3. (1999) M-plus R 1 ( ) 1.1 Polychoric Correlation Polyserial Correlation Graded Response Model
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(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 information1 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
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
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 informationk3 ( :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講義のーと : データ解析のための統計モデリング. 第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統計研修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 informationUse 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 informationk2 ( :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 informationkubostat2017c 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講義のーと : データ解析のための統計モデリング. 第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σ 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一般化線形 (混合) モデル (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 informationMicrosoft 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 informationuntitled
1 Hitomi s English Tests 1 2 3 4 5 6 7 8 9 10 11 12 13 14 1 1 0 1 1 0 1 0 0 0 1 0 0 1 0 2 0 0 1 1 0 0 0 0 0 1 1 1 1 0 3 1 1 0 0 0 0 1 0 1 0 1 0 1 1 4 1 1 0 1 0 1 1 1 1 0 0 0 1 1 5 1 1 0 1 1 1 1 0 0 1 0
More information浜松医科大学紀要
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
More informationIsogai, 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¥¤¥ó¥¿¡¼¥Í¥Ã¥È·×¬¤È¥Ç¡¼¥¿²òÀÏ Âè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
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 informationkubostat2015e 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塗装深み感の要因解析
17 Analysis of Factors for Paint Depth Feeling Takashi Wada, Mikiko Kawasumi, Taka-aki Suzuki ( ) ( ) ( ) The appearance and quality of objects are controlled by paint coatings on the surfaces of the objects.
More informationMales are Susceptible to Unrequited Love: Development of the Benefits and Tsukasa KATO MY ALL Jungle Smile SEVENTH HEAVENPerfumeOh! YeahOR- ANGE RANGE
Males are Susceptible to Unrequited Love: Development of the Benefits and Tsukasa KATO MY ALL Jungle Smile SEVENTH HEAVENPerfumeOh! YeahOR- ANGE RANGEw-inds Google 115 Table BULS: the Benefits of Unrequited
More informationTable 1 Experimental conditions Fig. 1 Belt sanded surface model Table 2 Factor loadings of final varimax criterion 5 6
JSPE-54-04 Factor Analysis of Relationhsip between One's Visual Estimation and Three Dimensional Surface Roughness Properties on Belt Sanded Surface Motoyoshi HASEGAWA and Masatoshi SHIRAYAMA This paper
More information最小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 informationAR(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,
More information<30315F985F95B65F90B490852E696E6464>
Modeling for Change by Latent Difference Score Model: Adapting Process of the Student of Freshman at Half Year Intervals Kazuaki SHIMIZU and Norihiro MIHO Abstract The purpose of this paper is to present
More information0801391,繊維学会ファイバ12月号/報文-01-西川
Pattern Making Method and Evaluation by Dots of Monochrome Shigekazu Nishikawa 1,MarikoYoshizumi 1,andHajime Miyake 2 1 Miyagi University of Education, 149, Aramaki-aza-Aoba, Aoba-ku, Sendai-shi, Miyagi
More informationkubostat2018d 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 informationStata11 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こんにちは由美子です
Sample size power calculation Sample Size Estimation AZTPIAIDS AIDSAZT AIDSPI AIDSRNA AZTPr (S A ) = π A, PIPr (S B ) = π B AIDS (sampling)(inference) π A, π B π A - π B = 0.20 PI 20 20AZT, PI 10 6 8 HIV-RNA
More informationuntitled
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 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 informationStata User Group Meeting in Kyoto / ( / ) Stata User Group Meeting in Kyoto / 21
Stata User Group Meeting in Kyoto / 2017 9 16 ( / ) Stata User Group Meeting in Kyoto 2017 9 16 1 / 21 Rosenbaum and Rubin (1983) logit/probit, ATE = E [Y 1 Y 0 ] ( / ) Stata User Group Meeting in Kyoto
More informationStata 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 information80 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 =
More information¥¤¥ó¥¿¡¼¥Í¥Ã¥È·×¬¤È¥Ç¡¼¥¿²òÀÏ Âè2²ó
2 212 4 13 1 (4/6) : ruby 2 / 35 ( ) : gnuplot 3 / 35 ( ) 4 / 35 (summary statistics) : (mean) (median) (mode) : (range) (variance) (standard deviation) 5 / 35 (mean): x = 1 n (median): { xr+1 m, m = 2r
More informationACLI-EBC-CLHIA Interim Proposal _J_ June Final.PDF
/ ACLI-EBC-CLHIA 2 5800 2004 3 3 1 2004 3 55 7 ACLI EBC / CLHIA 3 3 20047 CTE 3 ACLI-EBC-CLHIA 3 CTE CTE CTE60 CTE 1) CTE(60) CTE(80) CTE(90) 2) 100 3) 8) 3 Mercer Oliver Wyman Actuarial Consulting ACLI-EBC-CLHIA
More informationuntitled
005 n X i i 1 5 i 1 5 i 1 X i 3 X i 40 n i1 i i n X i 40 1 005 95 005 testing statistical hypothesis - A B A B 5 ()()() ()()() 3 005 ( ) null hypothesis 5 pp0.01p1 4 005 (1) 1 4 6 () N i 1 ( X i X ) N
More information28
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 informationNewgarten, BL., Havighrst, RJ., & Tobin, S.Life Satisfaction Index-A LSIDiener. E.,Emmons,R.A.,Larsen,R.J.,&Griffin,S. The Satisfaction With Life Scal
青年期における人生に対する 積極的態度に関する研究 KJ 法による検討と尺度の構成を中心として 海老根 理 絵 臨床心理学コース The research of the positive attitude toward life in adolescence Rie EBINE The purpose of this study is to assess qualitatively the structure
More informationMicrosoft PowerPoint - 譫礼峩荵
2011 3 11 323 25 HP 2013 JOC 2014 JOC HP 2014 Go for it 2011 J A 1 2 3 2013bj 2013 Human Welfare 7 1 2015 CiNii 1980 3 1990 17 2000 63 2010 46 2014 7 29 2004 2012 20132007 2003 2003 2002 3 1 20 70 1 2013
More informationMicrosoft Word - 計量研修テキスト_第5版).doc
Q4-1 テキスト P83 多重共線性が発生する回帰 320000 280000 240000 200000 6000 4000 160000 120000 2000 0-2000 -4000 74 76 78 80 82 84 86 88 90 92 94 96 98 R e s i dual A c tual Fi tted Dependent Variable: C90 Date: 10/27/05
More informationwin版8日目
8 日目 : 項目のチェック (2) 1 日 30 分くらい,30 日で何とか R をそこそこ使えるようになるための練習帳 :Win 版 昨日は, 平均値などの基礎統計量を計算する試行錯誤へご招待しましたが (?), 今日は簡 単にやってみます そのためには,psych というパッケージが必要となりますが, パッケー ジのインストール & 読み込みの詳しい方法は, 後で説明します 以下の説明は,psych
More information15 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
More informationMicrosoft Word - News 18 本文.doc
Argonauta 18: 17 30 (2010) 2007 1971 Cohen et al. 2003 A, B, C X, Y, Z 17 error error 1 0 1 1 X Y 1981 i) ii) 0 18 iii) iv) Yi = 0 + 1 X1i + 2 X2i + + p Xpi + i Xp pxpi + i i i Xp i Xp i 2 X Y 0 Xi Xi
More informationA 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
More informationJ. Jpn. Acad. Nurs. Sci. 35: (2015)
J. Jpn. Acad. Nurs. Sci., Vol. 35, pp. 257 266, 2015 DOI: 10.5630/jans.35.257 Development of Natural Disaster Preparedness Scale for Nursing Department of Hospital: Reliability and Validity as Scale Ayumi
More information現代日本論演習/比較現代日本論研究演習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
More informationudc-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
More informationuntitled
WinLD R (16) WinLD https://www.biostat.wisc.edu/content/lan-demets-method-statistical-programs-clinical-trials WinLD.zip 2 2 1 α = 5% Type I error rate 1 5.0 % 2 9.8 % 3 14.3 % 5 22.6 % 10 40.1 % 3 Type
More information9 2 1 f(x, y) = xy sin x cos y x y cos y y x sin x d (x, y) = y cos y (x sin x) = y cos y(sin x + x cos x) x dx d (x, y) = x sin x (y cos y) = x sin x
2009 9 6 16 7 1 7.1 1 1 1 9 2 1 f(x, y) = xy sin x cos y x y cos y y x sin x d (x, y) = y cos y (x sin x) = y cos y(sin x + x cos x) x dx d (x, y) = x sin x (y cos y) = x sin x(cos y y sin y) y dy 1 sin
More informationVisual Evaluation of Polka-dot Patterns Yoojin LEE and Nobuko NARUSE * Granduate School of Bunka Women's University, and * Faculty of Fashion Science,
Visual Evaluation of Polka-dot Patterns Yoojin LEE and Nobuko NARUSE * Granduate School of Bunka Women's University, and * Faculty of Fashion Science, Bunka Women's University, Shibuya-ku, Tokyo 151-8523
More informationMicrosoft 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日立金属技報 Vol.34
Influence of Misorientation Angle between Adjacent Grains on Magnetization Reversal in Nd-Fe-B Sintered Magnet Tomohito Maki Rintaro Ishii Mitsutoshi Natsumeda Takeshi Nishiuchi Ryo Uchikoshi Masaaki Takezawa
More information10
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
More informationオーストラリア研究紀要 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 informationMicrosoft 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
More information10:30 12:00 P.G. vs vs vs 2
1 10:30 12:00 P.G. vs vs vs 2 LOGIT PROBIT TOBIT mean median mode CV 3 4 5 0.5 1000 6 45 7 P(A B) = P(A) + P(B) - P(A B) P(B A)=P(A B)/P(A) P(A B)=P(B A) P(A) P(A B) P(A) P(B A) P(B) P(A B) P(A) P(B) P(B
More informationRによる計量分析:データ解析と可視化 - 第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 informationDO(Sep2014)
Mathematics Diagnostic Examination Guidance Examination Overview The Mathematics Examination will contain two sections. Section A contains fifteen 2-point questions and Section B contains five 4-point
More informationIPSJ SIG Technical Report Vol.2009-BIO-17 No /5/26 DNA 1 1 DNA DNA DNA DNA Correcting read errors on DNA sequences determined by Pyrosequencing
DNA 1 1 DNA DNA DNA DNA Correcting read errors on DNA sequences determined by Pyrosequencing Youhei Namiki 1 and Yutaka Akiyama 1 Pyrosequencing, one of the DNA sequencing technologies, allows us to determine
More informationWeb Stamps 96 KJ Stamps Web Vol 8, No 1, 2004
The Journal of the Japan Academy of Nursing Administration and Policies Vol 8, No 1, pp 43 _ 57, 2004 The Literature Review of the Japanese Nurses Job Satisfaction Research Which the Stamps-Ozaki Scale
More informationThe Japanese Journal of Health Psychology, 29(S): (2017)
Journal of Health Psychology Research 2017, Vol. 29, Special issue, 139 149Journal of Health Psychology Research 2016, J-STAGE Vol. Advance 29, Special publication issue, 139 149 date : 5 December, 2016
More information分布
(normal distribution) 30 2 Skewed graph 1 2 (variance) s 2 = 1/(n-1) (xi x) 2 x = mean, s = variance (variance) (standard deviation) SD = SQR (var) or 8 8 0.3 0.2 0.1 0.0 0 1 2 3 4 5 6 7 8 8 0 1 8 (probability
More informationkubostat2017e 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パーソナリティ研究2006 第14巻 第2号 214–226
2006 14 2 214 226 2006 (NPI-35) 1) Raskin & Hall (1979) Emmons (1984) (Narcissistic Personality Inventory; NPI) Raskin & Terry (1988) NPI (Narcissistic Personality Inventory-35; NPI-35) 1 35 5 2 3 4 NPI-35
More information02[021-046]小山・池田(責)岩.indd
Developing a Japanese Enryo-Sasshi Communication Scale: Revising a Trial Version of a Scale Based on Results of a Pilot Survey KOYAMA Shinji and IKEDA Yutaka Toward exploring Japanese Enryo-Sasshi communication
More informationA comparative study of the team strengths calculated by mathematical and statistical methods and points and winning rate of the Tokyo Big6 Baseball Le
Powered by TCPDF (www.tcpdf.org) Title 東京六大学野球リーグ戦において勝敗結果から計算する優勝チームと勝点 勝率との比較研究 Sub Title A comparative study of the team strengths calculated by mathematical and statistical methods and points and winning
More information1 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 information149 (Newell [5]) Newell [5], [1], [1], [11] Li,Ryu, and Song [2], [11] Li,Ryu, and Song [2], [1] 1) 2) ( ) ( ) 3) T : 2 a : 3 a 1 :
Transactions of the Operations Research Society of Japan Vol. 58, 215, pp. 148 165 c ( 215 1 2 ; 215 9 3 ) 1) 2) :,,,,, 1. [9] 3 12 Darroch,Newell, and Morris [1] Mcneil [3] Miller [4] Newell [5, 6], [1]
More information.3 ˆβ1 = S, S ˆβ0 = ȳ ˆβ1 S = (β0 + β1i i) β0 β1 S = (i β0 β1i) = 0 β0 S = (i β0 β1i)i = 0 β1 β0, β1 ȳ β0 β1 = 0, (i ȳ β1(i ))i = 0 {(i ȳ)(i ) β1(i ))
Copright (c) 004,005 Hidetoshi Shimodaira 1.. 3. 4. 004-10-01 16:15:07 shimo cat(" 1: "); c(mea(), mea()) cat(" : "); mmea
More information!!! 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(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,
More informationR分散分析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日本統計学会誌, 第44巻, 第2号, 251頁-270頁
44, 2, 205 3 25 270 Multiple Comparison Procedures for Checking Differences among Sequence of Normal Means with Ordered Restriction Tsunehisa Imada Lee and Spurrier (995) Lee and Spurrier (995) (204) (2006)
More information紀要1444_大扉&目次_初.indd
1. Hi, How are you? / What s up? / How s it going? A / Nice talking to you. 2. Oh really? / That s great! / That s A, B interesting! / Are you serious? / Sounds good. / You too! / That s too bad. 3. Sorry?
More informationuntitled
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
More informationStudies 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,
More information151021slide.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 information66-1 田中健吾・松浦紗織.pwd
Abstract The aim of this study was to investigate the characteristics of a psychological stress reaction scale for home caregivers, using Item Response Theory IRT. Participants consisted of 337 home caregivers
More informationAkhtar, S., & Thomson, J. A. 1982 Overview: Elson, M (Ed.) 1987 The Kohut seminars on self psychology and psychotherapy with adolescents and young adults. New York: W.W.Norton & Company. Gabbard, G.O.
More information8日目
8 日 目 : 項 目 のチェック(2) 1 日 30 分 くらい,30 日 で 何 とかRをそこそこ 使 えるようになるための 練 習 帳 :Mac 版 昨 日 は, 平 均 値 などの 基 礎 統 計 量 を 計 算 する 試 行 錯 誤 へご 招 待 しましたが(?), 今 日 は 簡 単 にやってみます そのためには,psych というパッケージが 必 要 となりますので,R を 起 動 したら,まずこれをとってきてください
More information: : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : :
1 1 1 1.1 : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : 1 1.2 : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : 1 1.2.1 t :
More information(1) (2) (3) (1) (2) (3) (1) (2) (3) (4) (1) (2)
(1) (2) (3) (1) (2) (3) (1) (2) (3) (4) (1) (2) (1) (2) (3) (1) (2) (3) (1) (2) (3) vs. (2) (1) 16 22 143,977 24,030 55,769 600 (2) (3) (4) (1) (2) (1) (2) 2 (3) (abuse) (1) (2) (3) (1) (2) (3)
More information12日目
12 日目 : 因子分析 (2) 1 日 30 分くらい,30 日で何とか R をそこそこ使えるようになるための練習帳 :Mac 版 さて, 昨日はいろいろと分析を試していただけたでしょうか? 因子分析は, かなりの試行錯誤を必要とするので, じっくりと構えて取り組む必要があります 投入する変数を削除することにも挑戦した人はわかるでしょうが, 以下のような変数をまとめておく作業は, ここで結構役立ちます
More information講義のーと : データ解析のための統計モデリング. 第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
More informationR 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療養病床に勤務する看護職の職務関与の構造分析
原著 :. JDS Job Diagnostic SurveyHackman & OldhamStamps, Herzberg Ⅰ. 諸言,, 10.,, 11 Ⅱ. 方法 1. 概念枠組みと質問紙の測定尺度 Hackman & Oldham Hackman & Oldham JDS 内発的動機づけ職務特性 技能多様性 タスク明確性 タスク重要性 自律性 職務からのフィードバック 他者からのフィードバック
More informationMicrosoft 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
More informationMeas- urement Angoff, W. H. 19654 Equating non-parallel tests. Journal of Educational Measurement, 1, 11-14. Angoff, W. H. 1971a Scales, norms and equivalent scores. In R. L. Thorndike (Ed.) Educational
More information- 1-128 - 2 -
127 - 1-128 - 2 - - 3-129 - 4 - 2-5 - 130-6 - - 7-131 - 8 - - 9-132 - 10 - 6041 3 () 1 ( ) () 6041 (1010) 1041 (192) 1941 () 2 (1) (2) (3) () 3 1 1 () 4 2 () 5 1 2 3 4 () 6 () 7-11 - 133-12 - 134 135 136
More information撮 影
DC cathode ray tube, 2.2 log log log + log log / / / A method determining tone conversion characteristics of digital still camera from two pictorial images Tone conversion characteristic Luminance
More information10_08.dvi
476 67 10 2011 pp. 476 481 * 43.72.+q 1. MOS Mean Opinion Score ITU-T P.835 [1] [2] [3] Subjective and objective quality evaluation of noisereduced speech. Takeshi Yamada, Shoji Makino and Nobuhiko Kitawaki
More information