こんにちは由美子です
|
|
- あいり たもん
- 5 years ago
- Views:
Transcription
1 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 HIV-RNA 20 n A = AZT S A = AZTn B = PI S B = PI π A = S A /n A, π B = S B /n B, AIDS =δ = π A - π B 1
2 H 0 : = 0 H A : 0 δh 0 reject H A accept δ π π π π δ Sample α/2 = α/2 = You Sample H 0 God Truth H A H 0 A=B (OK) Type II error (β error) H A Type I error A is not B (α error) truth A B A B Type I error(α error) A B Type II error(β error)fda α error β error 2
3 Truth Sample β -error Power = 1 - β -error power? If difference of means increase, If sample number increased β-error one side 2 AIDS AZT 20PI 30 HIV-RNA α = 0.05 (two sided), β = 0.20 ( power = 80%) α = 0.05 (two sided) δ = ( )/0.3 (1 0.3) (1 0.2) = n =[(Zα + Zβ)/d] 2 = [( )/0.1644] 2 = 290 AZT 290 PI
4 Power 90two sided test α = 0.05, 35% % 45% ( ) / (0.45 x x 0.65) = [( ) / 0.145] = 500 / 50% 60% ( ) / (0.60 x x 0.50) = [( ) / 0.143] = 513 / 4
5 α 5
6 6
7 7
8 success / failure AIDSAZT 24 HIV-RNA µ A PIµ B H 0 : µ A - µ B =0 H A :µ A - µ B 0 AIDS X A, X B µ A - µ B X A - X B X A - X B sample variance (Zα + Zβ) 2 δ 2 δ = µ 1 µ 2 8
9 SD (Zα + Zβ) 2 δ 2 δ = µ 1 µ SD [ 2 x S 2 (Zα + Zβ) 2 / 2 ] = [2 x 20 2 x ( ) 2 ]/10 2 = 63/arm [ 2 x S 2 (Zα + Zβ) 2 / 2 ] = [2 x 20 2 x ( ) 2 ]/6 2 = 175/arm 95(confidence interval) 95(confidence interval) SE Sample size 95CI (π B π A ) π π π π standard error (SE) phase II trial (1 arm) 95CI
10 one arm ππ When π = 0.8, n = 62 When π = 0.7, n = 81 When π = 0.6, n = 93 When π = 0.5, n = 96 When π = 0.4, n = 93 When π = 0.3, n = 81 When π = 0.2, n = 62 When π = 0.1, n = Clinical Equivalence Trials Bio-equivalence under the curve Cmax Clinical Equivalence Trials 10 Clinical Equivalence Trials one side AZT AIDS ddi AIDS D < D < 0.1 ddi AZT AZT OK (standard) 0.7 (new)
11 = πs πn 95% CI upper %CI 0.2 OK 95% CI upper 95CI 0.2 (πs πn)+ 1.65[πS (1 - πs)/n + πs (1 - πs)/n] πs = πn=0.7 sample size 1.65[πS (1 - πs)/n + πs (1 - πs)/n] = [2 x 0.7(1 0.7)/n]=0.2 n = 29 / arm 0.6 πs πn= [πS (1 - πs)/n + πs (1 - πs)/n] = [2 x 0.7(1 0.7)/n]=0.1 n = 115 / arm sample size 11
12 (πs πn)+ 1.65[πS (1 - πs)/n + πs (1 - πs)/n]=( ) [0.80 (1 0.80)/ (1 0.75)/100] = % p 1 p 2 95CI 10simple mastectomynon-equivalent α β δ δ α β δ equivalent 12
13 Hazard function Prob (T>t) = e λt H 0 : λ 1(t) = λ 2(t) H A : λ 1(t) = constant λ 2(t) 2 *hazard rate λ 1 /λ 2 * = ln(λ 1 /λ 2 )/2 d =[(Zα + Zβ)/*] 2 d: X 1 (λ 1 ) (λ 2 )α = 0.05, power = 80% * = ln(λ 1 /λ 2 )/2 = ln(1.5)/2 = d = [( ) / 0.287] 2 = Years of additional follow-up Years of accrual Accrual follow-up accrual 1 follow-up 2 3 Prob (T>t) = e λt
14 t = 11 Prob (T>1) = e λ1 = 0.5 ln (0.5) = λ 1 = 0.69 λ 1 /λ 2 = 1.5 = 0.69/λ 2 λ 2 = 0.46 Hazard function λ Accrual years = A, Follow-up years = F F A/2 A/2 + F accrual 2 2 2/2 + 2 = Prob (T>t) = e λt 3 T failure time Prob (T>3) = e λt = e 0.69 x 3 = Prob (T<3) = = Prob (T>3) = e λt = e 0.46 x 3 = Prob (T<3) = = sample size event %, 74.8% 14
15 (sample size) accrual AIDS randomized clinical trial AIDS AIDS (accrual)1 Type I error 5%, type II error 20% Prob (T>t) = e λt H 0 : λ 1(t) = λ 2(t) H A : λ 1(t) = constant λ 2(t) * = ln(λ 1 /λ 2 )/2 = ln(2.0/1.5)/2 = d = [( ) / 0.203] 2 = Prob (T>t) = e λt 1.5 Prob (T>1.5) = e λ x 1.5 = 0.5 ln (0.5) = λ 1 = 0.46 λ 1 /λ 2 = 2.0/1.5 = 1.33= 0.46/λ 2 λ 2 =
16 accrual 3 1 3/2 + 1 = Prob (T>t) = e λt 2.5 T failure time Prob (T>2.5) = e λt = e 0.46 x 2.5 = Prob (T<2.5) = = Prob (T>2.5) = e λt = e 0.35 x 2.5 = Prob (T<2.5) = = sample size event %, 58.3% (sample size)
17 STATA sample size randomized placebo controlled clinical trial randomization 4 68 pilot study placebo 498 ± 20.2 sec, 485 ± 19.5 sec 0.7 change method α = 0.05 (two sided), 90% power STATA command. sampsi , sd1(20.2) sd2(19.5) method(change) pre(1) post(3) r1(.7) Estimated sample size for two samples with repeated measures Assumptions: alpha = (two-sided) power = m1 = 498 m2 = 485 sd1 = 20.2 sd2 = 19.5 n2/n1 = 1.00 number of follow-up measurements = 3 correlation between follow-up measurements = number of baseline measurements = 1 correlation between baseline & follow-up = Method: CHANGE relative efficiency = adjustment to sd = adjusted sd1 = adjusted sd2 = Estimated required sample sizes: n1 = 20 n2 = 20 placebo 20 17
18 Clinical trials with repeated measures sampsi , sd1(20.2) sd2(19.5) method(change) pre(1) post(3) r1(.7) n1(1 5) n2(15) Estimated power for two samples with repeated measures Assumptions: alpha = (two-sided) m1 = 498 m2 = 485 sd1 = 20.2 sd2 = 19.5 sample size n1 = 15 n2 = 15 n2/n1 = 1.00 number of follow-up measurements = 3 correlation between follow-up measurements = number of baseline measurements = 1 correlation between baseline & follow-up = Method: CHANGE relative efficiency = adjustment to sd = adjusted sd1 = adjusted sd2 = Estimated power: power = placebo 20. sampsi , sd1(20.2) sd2(19.5) method(change) pre(1) post(3) r1(.7) n1(20) n2(15) Estimated power for two samples with repeated measures Assumptions: alpha = (two-sided) m1 = 498 m2 = 485 sd1 = 20.2 sd2 = 19.5 sample size n1 = 20 n2 = 15 18
19 n2/n1 = 0.75 number of follow-up measurements = 3 correlation between follow-up measurements = number of baseline measurements = 1 correlation between baseline & follow-up = Method: CHANGE relative efficiency = adjustment to sd = adjusted sd1 = adjusted sd2 = Estimated power: power =
20 Two-sample test of equality of proportions Yes/no type sample size H 0 reject α = 0.05, power 0.80 sample. sampsi , power(0.8) Estimated sample size for two-sample comparison of proportions Test Ho: p1 = p2, where p1 is the proportion in population 1 and p2 is the proportion in population 2 Assumptions: alpha = (two-sided) power = p1 = p2 = n2/n1 = 1.00 Estimated required sample sizes: n1 = 222 n2 = phase I trial 90. sampsi , power(0.9) Estimated sample size for two-sample comparison of proportions Test Ho: p1 = p2, where p1 is the proportion in population 1 and p2 is the proportion in population 2 Assumptions: alpha = (two-sided) power = p1 = p2 = n2/n1 = 1.00 Estimated required sample sizes: n1 = 287 n2 =
21 300 placebo 150. sampsi , n1(300) r(0.5) Estimated power for two-sample comparison of proportions Test Ho: p1 = p2, where p1 is the proportion in population 1 and p2 is the proportion in population 2 Assumptions: alpha = (two-sided) p1 = p2 = sample size n1 = 300 n2 = 150 n2/n1 = 0.50 Estimated power: power = placebo 2:1 80. sampsi , power(0.8) r(0.5) Estimated sample size for two-sample comparison of proportions Test Ho: p1 = p2, where p1 is the proportion in population 1 and p2 is the proportion in population 2 Assumptions: alpha = (two-sided) power = p1 = p2 = n2/n1 = 0.50 Estimated required sample sizes: n1 = 349 n2 =
22 One sample test of proportion golden standard 75 a = 0.05, 80. sampsi , power(0.8) onesample Estimated sample size for one-sample comparison of proportion to hypothesized value Test Ho: p = , where p is the proportion in the population Assumptions: alpha = (two-sided) power = alternative p = Estimated required sample size: n = historical comparison randomized clinical trial. sampsi , power(0.8) onesample Estimated sample size for one-sample comparison of proportion to hypothesized value Test Ho: p = , where p is the proportion in the population Assumptions: alpha = (two-sided) power = alternative p = Estimated required sample size: n = 29 22
23 Two sample test of equality of means endpoint 105 mmhg SD 10 mmhg 98 SD 10 placebo 2:1 80α = sampsi , p(0.8) r(2) sd1(10) sd2(10) Estimated sample size for two-sample comparison of means Test Ho: m1 = m2, where m1 is the mean in population 1 and m2 is the mean in population 2 Assumptions: alpha = (two-sided) power = m1 = 105 m2 = 98 sd1 = 10 sd2 = 10 n2/n1 = 2.00 Estimated required sample sizes: n1 = 25 n2 = 50 23
こんにちは由美子です
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分布
(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 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 informationPhase II clinical trials for patients with cancer
II II I 1 I III 14 20 14 20 14 0.05 14 20 1 10 25 20-80 100-200 95 CI (π B π A ) 1.96 [π A (1 π A )/n A + π B (1 π B )/n B ] π A = probability of response rate by treatment A π B = probability of response
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% 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 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 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 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 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 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 information1 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 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 informationJMP 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
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 informationMicrosoft 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 informationuntitled
2 : n =1, 2,, 10000 0.5125 0.51 0.5075 0.505 0.5025 0.5 0.4975 0.495 0 2000 4000 6000 8000 10000 2 weak law of large numbers 1. X 1,X 2,,X n 2. µ = E(X i ),i=1, 2,,n 3. σi 2 = V (X i ) σ 2,i=1, 2,,n ɛ>0
More informationECCS. 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
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 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 informationI L01( Wed) : Time-stamp: Wed 07:38 JST hig e, ( ) L01 I(2017) 1 / 19
I L01(2017-09-20 Wed) : Time-stamp: 2017-09-20 Wed 07:38 JST hig e, http://hig3.net ( ) L01 I(2017) 1 / 19 ? 1? 2? ( ) L01 I(2017) 2 / 19 ?,,.,., 1..,. 1,2,.,.,. ( ) L01 I(2017) 3 / 19 ? I. M (3 ) II,
More informationカテゴリ変数と独立性の検定
II L04(2015-05-01 Fri) : Time-stamp: 2015-05-01 Fri 22:28 JST hig 2, Excel 2, χ 2,. http://hig3.net () L04 II(2015) 1 / 20 : L03-S1 Quiz : 1 2 7 3 12 (x = 2) 12 (y = 3) P (X = x) = 5 12 (x = 3), P (Y =
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症例数設定? What is sample size estimation? 医療機器臨床試験のコンサルティングで最も相談件数が多いのは 症例数の設定 Many a need of consulting for device clinical trial is sample size estimat
医療機器臨床試験における 例数設計のいろは Sample size estimation for clinical trials of devices -First education- 株式会社バイオスタティスティカルリサ - チ古川敏仁 Furukawa Toshihito, Biostatistical Research 2005 年 9 月 3 日第 2 回医療機器臨床試験研究会 Copyright(C)
More informationα β *2 α α β β α = α 1 β = 1 β 2.2 α 0 β *3 2.3 * *2 *3 *4 (µ A ) (µ P ) (µ A > µ P ) 10 (µ A = µ P + 10) 15 (µ A = µ P +
Armitage 1 1.1 2 t *1 α β 1.2 µ x µ 2 2 2 α β 2.1 1 α β α ( ) β *1 t t 1 α β *2 α α β β α = α 1 β = 1 β 2.2 α 0 β 1 0 0 1 1 5 2.5 *3 2.3 *4 3 3.1 1 1 1 *2 *3 *4 (µ A ) (µ P ) (µ A > µ P ) 10 (µ A = µ P
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 information(pdf) (cdf) Matlab χ ( ) F t
(, ) (univariate) (bivariate) (multi-variate) Matlab Octave Matlab Matlab/Octave --...............3. (pdf) (cdf)...3.4....4.5....4.6....7.7. Matlab...8.7.....9.7.. χ ( )...0.7.3.....7.4. F....7.5. t-...3.8....4.8.....4.8.....5.8.3....6.8.4....8.8.5....8.8.6....8.9....9.9.....9.9.....0.9.3....0.9.4.....9.5.....0....3
More informationy 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................................................
More informationVol. 36, Special Issue, S 3 S 18 (2015) PK Phase I Introduction to Pharmacokinetic Analysis Focus on Phase I Study 1 2 Kazuro Ikawa 1 and Jun Tanaka 2
Vol. 36, Special Issue, S 3 S 18 (2015) PK Phase I Introduction to Pharmacokinetic Analysis Focus on Phase I Study 1 2 Kazuro Ikawa 1 and Jun Tanaka 2 1 2 1 Department of Clinical Pharmacotherapy, Hiroshima
More information2007-Kanai-paper.dvi
19 Estimation of Sound Source Zone using The Arrival Time Interval 1080351 2008 3 7 S/N 2 2 2 i Abstract Estimation of Sound Source Zone using The Arrival Time Interval Koichiro Kanai The microphone array
More information第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 information64 3 g=9.85 m/s 2 g=9.791 m/s 2 36, km ( ) 1 () 2 () m/s : : a) b) kg/m kg/m k
63 3 Section 3.1 g 3.1 3.1: : 64 3 g=9.85 m/s 2 g=9.791 m/s 2 36, km ( ) 1 () 2 () 3 9.8 m/s 2 3.2 3.2: : a) b) 5 15 4 1 1. 1 3 14. 1 3 kg/m 3 2 3.3 1 3 5.8 1 3 kg/m 3 3 2.65 1 3 kg/m 3 4 6 m 3.1. 65 5
More information自由集会時系列part2web.key
spurious correlation spurious regression xt=xt-1+n(0,σ^2) yt=yt-1+n(0,σ^2) n=20 type1error(5%)=0.4703 no trend 0 1000 2000 3000 4000 p for r xt=xt-1+n(0,σ^2) random walk random walk variable -5 0 5 variable
More information1 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 informationyamadaiR(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 informationkiyo5_1-masuzawa.indd
.pp. A Study on Wind Forecast using Self-Organizing Map FUJIMATSU Seiichiro, SUMI Yasuaki, UETA Takuya, KOBAYASHI Asuka, TSUKUTANI Takao, FUKUI Yutaka SOM SOM Elman SOM SOM Elman SOM Abstract : Now a small
More informationpositron 1930 Dirac 1933 Anderson m 22Na(hl=2.6years), 58Co(hl=71days), 64Cu(hl=12hour) 68Ge(hl=288days) MeV : thermalization m psec 100
positron 1930 Dirac 1933 Anderson m 22Na(hl=2.6years), 58Co(hl=71days), 64Cu(hl=12hour) 68Ge(hl=288days) 0.5 1.5MeV : thermalization 10 100 m psec 100psec nsec E total = 2mc 2 + E e + + E e Ee+ Ee-c mc
More informationスライド 1
SAS による二項比率の差の非劣性検定の正確な方法について 武藤彬正宮島育哉榊原伊織株式会社タクミインフォメーションテクノロジー Eact method of non-inferiority test for two binomial proportions using SAS Akimasa Muto Ikuya Miyajima Iori Sakakibara Takumi Information
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 information臨床試験
3000km 5000km Biomedical Research 1964 The purpose of biomedical research involving human subjects must be to improve diagnostic, therapeutic, and prophylactic procedures and the understanding of the aetiology
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 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こんにちは由美子です
Prevalence (proportion) 1. 1991 23.8 2. 1960 52 85 2477 310 cross sectional study prevalence Time referent: prevalence 1985 40 45 prevalence 0.5 80 85 43 Time referent 1985 time referent Risk and Cumulative
More information鉄鋼協会プレゼン
NN :~:, 8 Nov., Adaptive H Control for Linear Slider with Friction Compensation positioning mechanism moving table stand manipulator Point to Point Control [G] Continuous Path Control ground Fig. Positoining
More informationtokei01.dvi
2. :,,,. :.... Apr. - Jul., 26FY Dept. of Mechanical Engineering, Saga Univ., JAPAN 4 3. (probability),, 1. : : n, α A, A a/n. :, p, p Apr. - Jul., 26FY Dept. of Mechanical Engineering, Saga Univ., JAPAN
More information2 ( ) i
25 Study on Rating System in Multi-player Games with Imperfect Information 1165069 2014 2 28 2 ( ) i ii Abstract Study on Rating System in Multi-player Games with Imperfect Information Shigehiko MORITA
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 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 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²¾ÁÛ¾õ¶·É¾²ÁË¡¤Î¤¿¤á¤Î¥Ñ¥Ã¥±¡¼¥¸DCchoice ¡Ê»ÃÄêÈÇ¡Ë
DCchoice ( ) R 2013 2013 11 30 DCchoice package R 2013/11/30 1 / 19 1 (CV) CV 2 DCchoice WTP 3 DCchoice package R 2013/11/30 2 / 19 (Contingent Valuation; CV) WTP CV WTP WTP 1 1989 2 DCchoice package R
More information層化
(confounding) Confounding Exposure? Outcome Confounders adjustment (confounders) Not Confounding? 2 exposure 3rd factor 3rd factor NEJM 2000; 342: 1930-6. (adjustment) univariate or crude relative risk
More informationICU ( ) ICU ( ) ICU SpO 2 ICU ( ) -297 S.E DEX: N=30 DEX: N= 6 (J-DEX Placebo: N=30 Placebo: N= 3 65 DEX: N=16 DEX: N= 3 Placebo: N=
ICU ( ICU ( ICU SpO 2 ICU ( -297 S.E. -75-94 65 DEX: N=30 DEX: N= 6 ( Placebo: N=30 Placebo: N= 3 65 DEX: N=16 DEX: N= 3 Placebo: N=19 Placebo: N= 4 65 DEX: N=53 DEX: N=42 DEX: N=21 DEX: N=10 ( Placebo:
More informationSize Effect of Biomass on Carbonization Rate Treated in Superheated Steam Combined with Far Infrared Heating Akiko ISA Yoshio HAGURA and Kanichi Kit Graduate School of Biosphere Science, Hiroshima University,
More informationnsg04-28/ky208684356100043077
δ!!! μ μ μ γ UBE3A Ube3a Ube3a δ !!!! α α α α α α α α α α μ μ α β α β β !!!!!!!! μ! Suncus murinus μ Ω! π μ Ω in vivo! μ μ μ!!! ! in situ! in vivo δ δ !!!!!!!!!! ! in vivo Orexin-Arch Orexin-Arch !!
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 information日本製薬工業協会シンポジウム 生存時間解析の評価指標に関する最近の展開ー RMST (restricted mean survival time) を理解するー 2. RMST の定義と統計的推測 2018 年 6 月 13 日医薬品評価委員会データサイエンス部会タスクフォース 4 生存時間解析チー
日本製薬工業協会シンポジウム 生存時間解析の評価指標に関する最近の展開ー RMST (restricted mean survival time) を理解するー 2. RMST の定義と統計的推測 2018 年 6 月 13 日医薬品評価委員会データサイエンス部会タスクフォース 4 生存時間解析チーム 日本新薬 ( 株 ) 田中慎一 留意点 本発表は, 先日公開された 生存時間型応答の評価指標 -RMST(restricted
More information, (GPS: Global Positioning Systemg),.,, (LBS: Local Based Services).. GPS,.,. RFID LAN,.,.,.,,,.,..,.,.,,, i
25 Estimation scheme of indoor positioning using difference of times which chirp signals arrive 114348 214 3 6 , (GPS: Global Positioning Systemg),.,, (LBS: Local Based Services).. GPS,.,. RFID LAN,.,.,.,,,.,..,.,.,,,
More informationbosai-2002.dvi
45 B-2 14 4 Annuals of Disas. Prev. Res. Inst., Kyoto Univ., No. 45 B-2, 22 5 m 5 m :,,, 1. 2. 2.1 27 km 2 187 km 2 14 % 77 % 47 7, 9 2, 54 6 7, 9 16, 57 8 1, 9 47 2 1 57 5 2.2 45 2 Fig. 1 2 2.3 Fig. 2
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 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 information2005 2006.2.22-1 - 1 Fig. 1 2005 2006.2.22-2 - Element-Free Galerkin Method (EFGM) Meshless Local Petrov-Galerkin Method (MLPGM) 2005 2006.2.22-3 - 2 MLS u h (x) 1 p T (x) = [1, x, y]. (1) φ(x) 0.5 φ(x)
More informationII (No.2) 2 4,.. (1) (cm) (2) (cm) , (
II (No.1) 1 x 1, x 2,..., x µ = 1 V = 1 k=1 x k (x k µ) 2 k=1 σ = V. V = σ 2 = 1 x 2 k µ 2 k=1 1 µ, V σ. (1) 4, 7, 3, 1, 9, 6 (2) 14, 17, 13, 11, 19, 16 (3) 12, 21, 9, 3, 27, 18 (4) 27.2, 29.3, 29.1, 26.0,
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 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 Economic Affairs Waseda University 169-8050 Tokyo,Japan
More informationy = x 4 y = x 8 3 y = x 4 y = x 3. 4 f(x) = x y = f(x) 4 x =,, 3, 4, 5 5 f(x) f() = f() = 3 f(3) = 3 4 f(4) = 4 *3 S S = f() + f() + f(3) + f(4) () *4
Simpson H4 BioS. Simpson 3 3 0 x. β α (β α)3 (x α)(x β)dx = () * * x * * ɛ δ y = x 4 y = x 8 3 y = x 4 y = x 3. 4 f(x) = x y = f(x) 4 x =,, 3, 4, 5 5 f(x) f() = f() = 3 f(3) = 3 4 f(4) = 4 *3 S S = f()
More information漸化式のすべてのパターンを解説しましたー高校数学の達人・河見賢司のサイト
https://www.hmg-gen.com/tuusin.html https://www.hmg-gen.com/tuusin1.html 1 2 OK 3 4 {a n } (1) a 1 = 1, a n+1 a n = 2 (2) a 1 = 3, a n+1 a n = 2n a n a n+1 a n = ( ) a n+1 a n = ( ) a n+1 a n {a n } 1,
More informationA comparison of abdominal versus vaginal hysterectomy for leiomyoma and adenomyosis Kenji ARAHORI, Hisasi KATAYAMA, Suminori NIOKA Department of Obstetrics and Gnecology, National Maizuru Hospital,Kyoto,
More informationnews
ETL NEWS 1999.9 ETL NEWS 1999.11 Establishment of an Evaluation Technique for Laser Pulse Timing Fluctuations Optoelectronics Division Hidemi Tsuchida e-mail:tsuchida@etl.go.jp A new technique has been
More informationuntitled
January 2009 Rotor-Gene Q Sample & Assay Technologies 1 1-1 1.1 1-1 1.2 Rotor-Gene Q 1-1 1.3 Rotor-Gene Q 1-1 1.4 1-1 2 2-1 2.1 Run File 2-1 2.2 2-2 2.3 2-2 2.4 PCR 2-3 2.5 2-3 2.6 2-4 2.7 2-4 2.8 2-5
More information¼§À�ÍýÏÀ – Ê×ÎòÅŻҼ§À�¤È¥¹¥Ô¥ó¤æ¤é¤® - No.7, No.8, No.9
No.7, No.8, No.9 email: takahash@sci.u-hyogo.ac.jp Spring semester, 2012 Introduction (Critical Behavior) SCR ( b > 0) Arrott 2 Total Amplitude Conservation (TAC) Global Consistency (GC) TAC 2 / 25 Experimental
More informationVol. 29, No. 2, (2008) FDR Introduction of FDR and Comparisons of Multiple Testing Procedures that Control It Shin-ichi Matsuda Department of
Vol. 29, No. 2, 125 139 (2008) FDR Introduction of FDR and Comparisons of Multiple Testing Procedures that Control It Shin-ichi Matsuda Department of Information Systems and Mathematical Sciences, Faculty
More information2.2 h h l L h L = l cot h (1) (1) L l L l l = L tan h (2) (2) L l 2 l 3 h 2.3 a h a h (a, h)
1 16 10 5 1 2 2.1 a a a 1 1 1 2.2 h h l L h L = l cot h (1) (1) L l L l l = L tan h (2) (2) L l 2 l 3 h 2.3 a h a h (a, h) 4 2 3 4 2 5 2.4 x y (x,y) l a x = l cot h cos a, (3) y = l cot h sin a (4) h a
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 information数理統計学Iノート
I ver. 0/Apr/208 * (inferential statistics) *2 A, B *3 5.9 *4 *5 [6] [],.., 7 2004. [2].., 973. [3]. R (Wonderful R )., 9 206. [4]. ( )., 7 99. [5]. ( )., 8 992. [6],.., 989. [7]. - 30., 0 996. [4] [5]
More informationFig. 1.1 Annual commercial landings of masu salmon in Hokkaido during 1970-2002. Fig. 1.2 Number of hatchery-reared masu salmon stocked in Hokkaido, 1980-2001. Fig. 2.2a Masu salmon landed
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 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 informationMicrosoft Word - Power_Analysis_Jp_ docx
Power Analysis using G*Power Version 1.0 013 年 3 月 03 日 評価学博士 佐々木亮 サンプルサイズの検討方法 1. 最低のサンプルサイズサンプルサイズに関する考え方 統計分析を用いた調査報告書では サンプルサイズとして 30 あるいは 5 を用いている場合が頻繁に見られる 事前 事後比較のための 1 群の t 検定では まさに 30 あるいは 5 が必要ということになり
More informationA Study of Adaptive Array Implimentation for mobile comunication in cellular system GD133
A Study of Adaptive Array Implimentation for mobile comunication in cellular system 15 1 31 01GD133 LSI DSP CMA 10km/s i 1 1 2 LS-CMA 5 2.1 CMA... 5 2.1.1... 5 2.1.2... 7 2.1.3... 10 2.2 LS-CMA... 13 2.2.1...
More informationJ. Jpn. Soc. Soil Phys. No. 126, p (2014) ECH 2 O 1 2 Calibration of the capacitance type of ECH 2 O soil moisture sensors Shoichi MITSUISHI 1 a
J. Jpn. Soc. Soil Phys. No. 126, p.63 70 (2014) ECH 2 O 1 2 Calibration of the capacitance type of ECH 2 O soil moisture sensors Shoichi MITSUISHI 1 and Masaru MIZOGUCHI 1 Abstract: Volumetric water contents,
More information2 2 1?? 2 1 1, 2 1, 2 1, 2, 3,... 1, 2 1, 3? , 2 2, 3? k, l m, n k, l m, n kn > ml...? 2 m, n n m
2009 IA I 22, 23, 24, 25, 26, 27 4 21 1 1 2 1! 4, 5 1? 50 1 2 1 1 2 1 4 2 2 2 1?? 2 1 1, 2 1, 2 1, 2, 3,... 1, 2 1, 3? 2 1 3 1 2 1 1, 2 2, 3? 2 1 3 2 3 2 k, l m, n k, l m, n kn > ml...? 2 m, n n m 3 2
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 informationN cos s s cos ψ e e e e 3 3 e e 3 e 3 e
3 3 5 5 5 3 3 7 5 33 5 33 9 5 8 > e > f U f U u u > u ue u e u ue u ue u e u e u u e u u e u N cos s s cos ψ e e e e 3 3 e e 3 e 3 e 3 > A A > A E A f A A f A [ ] f A A e > > A e[ ] > f A E A < < f ; >
More information2 1 1 α = a + bi(a, b R) α (conjugate) α = a bi α (absolute value) α = a 2 + b 2 α (norm) N(α) = a 2 + b 2 = αα = α 2 α (spure) (trace) 1 1. a R aα =
1 1 α = a + bi(a, b R) α (conjugate) α = a bi α (absolute value) α = a + b α (norm) N(α) = a + b = αα = α α (spure) (trace) 1 1. a R aα = aα. α = α 3. α + β = α + β 4. αβ = αβ 5. β 0 6. α = α ( ) α = α
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.. 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
More informationADM-Hamiltonian Cheeger-Gromov 3. Penrose
ADM-Hamiltonian 1. 2. Cheeger-Gromov 3. Penrose 0. ADM-Hamiltonian (M 4, h) Einstein-Hilbert M 4 R h hdx L h = R h h δl h = 0 (Ric h ) αβ 1 2 R hg αβ = 0 (Σ 3, g ij ) (M 4, h ij ) g ij, k ij Σ π ij = g(k
More information2 The Bulletin of Meiji University of Integrative Medicine 3, Yamashita 10 11
1-122013 1 2 1 2 20 2,000 2009 12 1 2 1,362 68.1 2009 1 1 9.5 1 2.2 3.6 0.82.9 1.0 0.2 2 4 3 1 2 4 3 Key words acupuncture and moxibustion Treatment with acupuncture, moxibustion and Anma-Massage-Shiatsu
More information