untitled

Size: px
Start display at page:

Download "untitled"

Transcription

1 R Network Meta-Analysis

2 R Network Meta-Analysis homogeneity consistency Ranking Network Meta-Analysis(homogeneity) (consistency) Ranking 1

3 Network Meta-Analysis Network Meta-Analysis P Ranking R Bayesian Network Meta-Analysis 2

4 Network Meta-Analysis Meta-Analysis 2 Network Meta-Analysis Multiple Treatments Comparison Mixed Treatments Comparison 3 3

5 Network Meta-Analysis Network Meta-Analysis Meta-Analysis Network Meta-Analysis I. II. homogeneity consistency III. IV I II R 4

6 Network Meta-Analysis Network Meta-Analysis Network Meta-Analysis 3 multi-arm study Network A vs. B A vs. C B vs. C + A vs. D A B C D 5

7 Network Meta-Analysis Incomplete Block Design Sequence Incomplete Block Design IBD Period A B Patient 11,...,1n 1 2 B A Patient 21,...,2n 2 3 C B Patient 31,...,3n 3 4 B C Patient 41,...,4n 4 Network Meta-Analysis NMA 1: A vs. B Study 11,...,1n 1 IBD A B B C A C B A C GLM MIXED A 2: B vs. C Study 21,...,2n 2 B C 6

8 Network Meta-Analysis Incomplete Block Design Network Meta-Analysis NMA Incomplete Block Design IBD Senn et al IBD Sequence Sequence NMA IBD 1 NMA IBD NMA transitivity consistency 7

9 Network Meta-Analysis Network Meta-Analysis homogeneity consistency Ranking Bayesian Network Meta-Analysis B Treatment Mean Diff. 95% CI A D A B C D [-1.94; -1.57] [-0.61; -0.19] [-1.15; -0.74] C

10 Network Meta-Analysis homogeneity heterogeneity transitivity similarity Network Meta-Analysis Network Meta-Analysis effect modifier consistency Network Meta-Analysis consistency inconsistency 9 Caldwell (2014) Krahn (2013)

11 Network Meta-Analysis Network Meta-Analysis P Ranking R Bayesian Network Meta-Analysis 10

12 Senn2013 Senn et al. (2013) Senn et al. (2013) multi-arm study 1 10 acarbose benfluorex metformin miglitol pioglitazone rosiglitazone sitagliptin sulfonylurea alone vildagliptin placebo HbA1c % Senn2013 studlab treat1 1 treat2 2 TE HbA1c % TREAT1 TREAT2 sete TE 11

13 Senn2013 metf 1 benf migl 1 1 acar piog vild plac 2 1 sulf 6 rosi 1 sita multi-arm Willms

14 Senn2013 TE sete treat1 treat2 studlab metf sulf Alex rosi plac Baksi acar plac Costa rosi plac Davidson metf plac DeFronzo piog rosi Derosa vild plac Garber metf plac Gonzalez-Ortiz piog metf Hanefeld sita plac Hermansen migl plac Johnston migl plac Johnston1998a migl plac Johnston1998b rosi plac Kerenyi rosi metf Kim piog plac Kipnes metf plac Lewin benf plac Moulin acar sulf Oyama rosi plac Rosenstock benf plac Stucci rosi sulf Vongthavaravat acar metf Willms metf plac Willms acar plac Willms rosi plac Wolffenbuttel rosi metf Yang rosi plac Zhu

15 Senn2013 ID TE sete treat1 treat2 studlab s d metf sulf Alex acar plac Costa metf plac DeFronzo metf plac Gonzalez-Ortiz acar metf Willms metf plac Willms acar plac Willms acar metf sulf Senn metf vs. plac 3 multi-arm 1 3 ID=5 plac 14

16 Network Meta-Analysis 1 2 metf acar sulf plac metf % [-1.10, -0.43] % [-0.60, -0.14] % [-1.88, -1.40] acar 0.76 % [ 0.42, 1.10] 0.39 % [-0.02, 0.80] % [-1.14, -0.62] sulf 0.37 % [ 0.14, 0.60] % [-0.80, 0.02] % [-1.61, -0.94] plac 1.64 % [ 1.40, 1.88] 0.88 % [ 0.62, 1.14] 1.27 % [ 0.94, 1.61] Treatment metf acar sulf plac Mean Diff % -0.88% -1.27% 95% CI [-1.88, -1.40] [-1.14, -0.62] [-1.61, -0.94] [95%, 95% ]

17 Network Meta-Analysis, Cov, diag,,,,, plac 1,,5 1,,4 metf acar sulf V 1 V 4 V 5 Y 1 Y 4 Y 5 Willms metf acar plac M A P ID=6 SETE Var(M-P) = Var(M) + Var(P) - 2Cov(M,P) ID=5 ID=7 Cov(M-P, A-P) = Cov(M,A) Cov(M,P) Cov(P,A) +Var(P) 16

18 Network Meta-Analysis 1-step Network Estimate: 1.64, 0.88, step metf acar sulf, Cov Direct Estimate:,,,, Cov, diag,,, Network Estimate: 1.64, 0.88, 1.27 X 1 X 3 X 4 17

19 Network Meta-Analysis,Cov,,Cov 1,,5 1,,4 DerSimonian and Laird Cov 18 multi-arm τ /2 ID=5(acar vs. metf)

20 homogeneity consistency heterogeneity Within-designs Q statistic, inconsistency Between-designs Q statistic, / heterogeneity

21 P Ranking Cov( ) Ranking metf metformin P 1. metf p,,,, 2. P : 1, metf metformin P Ranking P P Bayesian Network Meta-Analysis SUCRA Surface Under the Cumulative RAnking 20 Rucker G and Schwarzer G (2015)

22 Network Meta-Analysis > # 1-step > y <- matrix(c(-0.37, -0.8, -1.9, -0.4, -1.2, -1), ncol=1) > X <- matrix(c(1, 0, -1, 0, 1, 0, 1, 0, 0, 1, 0, 0, 1, 0, 0, 0, 1, 0), ncol=3, by=t) > V <- matrix(c( , 0, 0, 0, 0, 0, + 0, , 0, 0, 0, 0, + 0, 0, , 0, 0, 0, + 0, 0, 0, , 0, 0, + 0, 0, 0, 0, , , + 0, 0, 0, 0, , ), ncol=6, by=t) > ( theta <- solve(t(x) %*% solve(v) %*% X) %*% t(x) %*% solve(v) %*% y ) [,1] [1,] [2,] [3,] > ( Cov <- solve(t(x) %*% solve(v) %*% X) ) [,1] [,2] [,3] [1,] [2,] [3,] > ( theta_ci <- cbind(theta-qnorm(0.975)*sqrt(diag(cov)), + theta+qnorm(0.975)*sqrt(diag(cov))) ) [,1] [,2] [1,] [2,] [3,]

23 Network Meta-Analysis > # 2-step > theta_dir <- matrix(rep(0, 5), ncol=1) > Xa <- X[c(1, 2, 3, 5, 6),] > Va <- matrix(rep(0,25), ncol=5) > theta_dir[1,1] <- sum(solve(v[1,1]) %*% y[1]) * V[1,1]; Va[1,1] <- V[1,1] > theta_dir[2,1] <- sum(solve(v[2,2]) %*% y[2]) * V[2,2]; Va[2,2] <- V[2,2] > theta_dir[3,1] <- sum(solve(v[3:4,3:4]) %*% y[3:4]) / sum(diag(solve(v[3:4,3:4]))) > Va[3,3] <- 1/sum(diag(solve(V[3:4,3:4]))) > theta_dir[4:5,1] <- y[5:6]; Va[4:5,4:5] <- V[5:6,5:6] > ( theta_a <- solve(t(xa) %*% solve(va) %*% Xa) %*% t(xa) %*% solve(va) %*% theta_dir ) [,1] [1,] [2,] [3,] > ( Cov_a <- solve(t(xa) %*% solve(va) %*% Xa) ) [,1] [,2] [,3] [1,] [2,] [3,] > ( theta_ci_a <- cbind(theta-qnorm(0.975)*sqrt(diag(cov_a)), + theta+qnorm(0.975)*sqrt(diag(cov_a))) ) [,1] [,2] [1,] [2,] [3,]

24 Network Meta-Analysis > # Qnet, Qhet, Qinc > ( Q_net <- t(y - X %*% theta) %*% solve(v) %*% (y-x %*% theta) ) [,1] [1,] > ( Q_inc <- t(theta_dir - Xa %*% theta_a) %*% solve(va) %*% (theta_dir - Xa %*% theta_a) ) [,1] [1,] > ( Q_het <- Q_net - Q_inc ) [,1] [1,] > # P score for metf > CC <- -1*diag(3); CC[,1] = 1; E <- CC %*% theta; z <- rep(0,3) > for (i in 1:3) { + c <- CC[i,,drop=F] + z[i] <- sqrt(t(c %*% theta) %*% solve(c %*% Cov %*% t(c)) %*% (c %*% theta)) + if (E[i] < 0) z[i] <- -1*z[i] + } > p_metf <- 1-(pnorm(z)) > ( pscore_metf <- mean(p_metf) ) [1]

25 inconsistency Detaching a single design inconsistency Krahn 2013 Detaching a single design, 1, 0, Cov 1 0 detach, 24

26 inconsistency Detaching a single design,,,,,,,,,,,, 0,,,,, 0 inconsistency, 0 inconsistency 25

27 inconsistency network estimate plac:metf_ plac:acar:metf plac:acar plac:metf plac:acar_ plac:acar:metf plac:metf_ plac:acar:metf direct estimate plac:acar plac:metf plac:acar_ plac:acar:metf (inconsistency) direct estimate network estimate Network estimate 26

28 inconsistency plac:metf_ plac:acar:metf Detach plac:acar plac:metf plac:acar_ plac:acar:metf Detach plac:metf_ plac:acar:metf plac:acar plac:metf plac:acar_ plac:acar:metf = inconsistent evidence = supportive evidence 27

29 inconsistency Detach plac:metf_ plac:acar:metf plac:acar plac:metf plac:acar_ plac:acar:metf plac:metf_ plac:acar:metf Detach plac:acar plac:metf plac:acar_ plac:acar:metf multi-arm metf vs. plac inconsistency metf vs. plac network estimate direct estimate metf vs. plac 3 2-arm 2 metf vs. plac 3 inconsistency 28

30 Network Meta-Analysis Network Meta-Analysis P Ranking R Bayesian Network Meta-Analysis 29

31 Network Meta-Analysis netmeta() netmeta(te, sete, treat1, treat2, studlab, data=senn2013, sm, level=0.95, level.comb=0.95, comb.fixed=true, comb.random=true, reference="", seq=null, title="",...) data TE sete treat1 treat2 studlab sm TE RD Risk Difference RR Risk Ratio OR Odds Ratio MD Mean Difference SMD Standardized Mean Difference IRR Incidence Rate Ratio IRD Incidence Rate Difference level level.comb pooled estimate comb.fixed comb.random reference reference="plac" seq title 30

32 Network Meta-Analysis # install.packages("netmeta", dep=t) library(netmeta) data(senn2013) options(width=1000) x <- netmeta(te, sete, treat1, treat2, studlab, data=senn2013, sm="md", comb.fixed=true, comb.random=true) d1 <- decomp.design(x) # Q d2 <- netmeta:::decomp.tau(x, tau=x$tau) # Q r2 <- netrank(x, small="good") # P forest(x, pooled='fixed') # Forest Plot() forest(x, pooled='random') # Forest Plot() netgraph(x, thickness='se.fixed'); title('fixed Effect Model') # () netgraph(x, thickness='se.random'); title('random Effect Model') # () netheat(x, random=f); mtext('fixed Effect Model', 1, 1, cex=1.8) # Heatmap ( ) netheat(x,random=t,tau=x$tau); mtext('random Effect Model',1,1,cex=1.8) #( ) x$te.random <- x$te.fixed; x$pval.random <- x$pval.fixed r1 <- netrank(x, small="good") # P ( small="good" or "bad" 31

33 Network Meta-Analysis x$te.fixed x$lower.fixed x$upper.fixed x$te.random x$lower.random x$upper.random d1$q.decomp Q homogeneity / consistency d1$q.het.design Design-specific decomposition of within-designs Q statistic d1$q.inc.detach Between-designs Q statistic after detaching of single designs d1$xxx 3 d2$q.decomp d2$q.het.design d2$q.inc.detach Q r1$pscore r2$pscore Network Ranking / 32

34 Network Meta-Analysis x <- netmeta(te, sete, treat1, treat2, studlab, data=senn2013, sm="md", comb.fixed=true, comb.random=true, reference="plac", seq=c('plac', 'acar', 'benf', 'metf', 'migl', 'piog','rosi', 'sita', 'sulf', 'vild')) # reference reference="plac" seq 33

35 Network Meta-Analysis > x$te.fixed # acar benf metf migl piog plac rosi sita sulf vild acar benf metf migl piog plac rosi sita sulf vild > x$lower.fixed # acar benf metf migl piog plac rosi sita sulf vild acar benf metf migl piog plac rosi sita sulf vild > r2$pscore # P acar benf metf migl piog plac rosi sita sulf vild

36 Network Meta-Analysis metf:sulf rosi:sulf metf:piog plac:piog plac:rosi plac:metf metf:rosi piog:rosi plac:acar_plac:acar:metf plac:metf_plac:acar:metf acar:sulf plac:acar Forest Plot 35 metf:sulf rosi:sulf metf:piog plac:piog plac:rosi plac:metf metf:rosi piog:rosi plac:acar_plac:acar:metf plac:metf_plac:acar:metf acar:sulf plac:acar Fixed Effect Model inconsistency multi-arm Willms1999

37 Network Meta-Analysis acar:sulf plac:acar_plac:acar:metf plac:metf plac:acar metf:rosi piog:rosi plac:metf_plac:acar:metf plac:rosi metf:piog plac:piog metf:sulf rosi:sulf Forest Plot acar:sulf plac:acar_plac:acar:metf plac:metf plac:acar metf:rosi piog:rosi plac:metf_plac:acar:metf plac:rosi metf:piog plac:piog metf:sulf rosi:sulf inconsistency 36 Random Effect Model multi-arm Willms1999

38 Arm-based Contrast-based Study Treatment1 y1 sd1 n1 Treatment2 y2 sd2 n2 Treatment3 y3 sd3 n NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA TE sete treat1 treat2 studlab contrast-based format arm-based format 37 Dias et. al. (2013) Parkinson

39 Arm-based Contrast-based pairwise(treat, event, n, time, mean, sd, TE, sete, data, studlab,...) treat 1 2 event n time mean sd TE sete data studlab 38 treat n mean sd

40 Arm-based Contrast-based Parkinson > data(parkinson) > p <- pairwise(list(treatment1, Treatment2, Treatment3), + n=list(n1, n2, n3), + mean=list(y1, y2, y3), + sd=list(sd1, sd2, sd3), + data=parkinson, studlab=study) > p TE sete studlab treat1 treat2 n1 mean1 sd1 n2 mean2 sd

41 Arm-based Contrast-based 2 data(smokingcessation) p1 <- pairwise(list(treat1, treat2, treat3), event=list(event1, event2, event3), n=list(n1, n2, n3), data=smokingcessation, sm="rr") 2 p2 <- pairwise(list(treat1, treat2, treat3), event=list(event1, event2, event3), n=list(n1, n2, n3), data=smokingcessation, sm="or") 2 treat event n 40 sm="rd"

42 Arm-based Contrast-based data(dietaryfat) p3 <- pairwise(list(treat1, treat2, treat3), list(d1, d2, d3), time=list(years1, years2, years3), studlab=id, data=dietaryfat, sm="ird") p4 <- pairwise(list(treat1, treat2, treat3), list(d1, d2, d3), time=list(years1, years2, years3), studlab=id, data=dietaryfat, sm="irr") treat time event 41

43 Network Meta-Analysis Network Meta-Analysis P Ranking R Bayesian Network Meta-Analysis 42

44 Bayesian Network Meta-Analysis Bayesian Network Meta-Analysis gemtc Network Meta-Analysis gemtc 1. JAGS JAGS exe 2. R LANGUAGE=en 3. install.packages(c("rjags","gemtc"), dep=t) gemtc > install.packages(c("rjags","gemtc"), dep=t) # > Sys.setlocale(locale="C") # locale US [1] "C" > # Sys.setlocale(locale="Japanese_Japan.932") # locale > plot.new() # History Recording ON > library(gemtc) WinBUGS 43

45 Bayesian Network Meta-Analysis > data <- read.table(textconnection(' + study treatment diff std.err + Alex1998 metf Alex1998 sulf NA NA... + Zhu2003 rosi Zhu2003 plac NA NA'), head=t) # Senn2013 > network <- mtc.network(data.re=data) > plot(network) > summary(network) $Description [1] "MTC dataset: Network" $`Studies per treatment` acar benf metf migl piog plac rosi sita sulf vild $`Number of n-arm studies` 2-arm 3-arm 25 1 $`Studies per treatment comparison` t1 t2 nr 1 acar metf 1 2 acar plac 2 3 acar sulf rosi sulf 1 44

46 Bayesian Network Meta-Analysis > # > model.fe <- mtc.model(network, linearmodel='fixed', likelihood='normal',link='identity') > result.fe <- mtc.run(model.fe, n.adapt=1000, n.iter=5000) > plot(result.fe) > gelman.diag(result.fe) # > summary(result.fe) Results on the Mean Difference scale Iterations = 1:5000 Thinning interval = 1 Number of chains = 4 Sample size per chain = Empirical mean and standard deviation for each variable, plus standard error of the mean: Mean SD Naive SE Time-series SE d.acar.metf d.acar.plac d.acar.sulf d.plac.benf d.plac.migl d.plac.piog d.plac.rosi d.plac.sita d.plac.vild

47 Bayesian Network Meta-Analysis > # > model.fe <- mtc.model(network, linearmodel='fixed', likelihood='normal',link='identity') > result.fe <- mtc.run(model.fe, n.adapt=1000, n.iter=5000) > plot(result.fe) > gelman.diag(result.fe) # > summary(result.fe) 2. Quantiles for each variable: 2.5% 25% 50% 75% 97.5% d.acar.metf d.acar.plac d.acar.sulf d.plac.benf d.plac.migl d.plac.piog d.plac.rosi d.plac.sita d.plac.vild Model fit (residual deviance): Dbar pd DIC data points, ratio 3.926, I^2 = 75% > forest(relative.effect(result.fe, t1="plac")) > ranks <- rank.probability(result.fe, preferreddirection=-1) > plot(ranks) 46

48 Bayesian Network Meta-Analysis > ranks Rank probability; preferred direction = -1 [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] acar benf metf migl piog plac rosi sita sulf vild

49 Bayesian Network Meta-Analysis Bayesian Network Meta-Analysis 48

50 Bayesian Network Meta-Analysis > # > model.re <- mtc.model(network, linearmodel='random', likelihood='normal',link='identity') > result.re <- mtc.run(model.re, n.adapt=1000, n.iter=5000) > plot(result.re) > gelman.diag(result.re) > summary(result.re) > forest(relative.effect(result.re, t1="plac")) > # DIC > result.fe$deviance$dic [1] > result.re$deviance$dic [1] > # Node-splitting Inconsistency > result.ns <- mtc.nodesplit(network, linearmodel='random', n.adapt=500, n.iter=2000) > plot(result.ns) > summary.ns <- summary(result.ns) > plot(summary.ns) > summary.ns 49

51 Bayesian Network Meta-Analysis Bayesian Network Meta-Analysis 50

52 Network Meta-Analysis Network Meta-Analysis P Ranking R Bayesian Network Meta-Analysis 51

53 Network Meta-Analysis Network Meta-Analysis R homogeneity consistency Ranking Bayesian Network Meta-Analysis 52

54 Caldwell DM (2014) An overview of conducting systematic reviews with network meta-analysis. Systematic Reviews, 3:109. Dias S, Sutton AJ, Ades AE and Welton NJ (2013) Evidence synthesis for decision making 2: A generalized linear modeling framework for pairwise and network meta-analysis of randomized controlled trials. Medical Decision Making, 33, Higgins JPT, Jackson D, Barrett JK, Lu G, et. al. (2012) Consistency and inconsistency in network meta-analysis: concepts and models for multi-arm studies. Research Synthesis Methods, 3(2), pages Hutton B, Salanti G, Caldwell DM, Chaimani A, et. al. (2015) The PRISMA extension statement for reporting of systematic reviews incorporating network meta-analyses of health care interventions: checklist and explanations. Annals of Internal Medicine, 162(11): Krahn U, Binder H and Konig J (2013) A graphical tool for locating inconsistency in network meta analyses. BMC Medical Research Methodology, 13:35. 53

55 Mills EJ, Thorlund K and Ioannidis JP (2013) Demystifying trial networks and network meta-analysis. BMJ. May14; 346:f2914. DerSimonian R and Laird N (1986) Meta-analysis in clinical trials. Controlled Clinical Trials 7: Rucker G, Schwarzer G, Krahn U and Konig J (2014) Netmeta: network meta-analysis with R. R package (version 0.8-0). Rucker G and Schwarzer G (2015) Ranking treatments in frequentist network meta-analysis works without resampling methods. BMC Medical Research Methodology, 15:58. Senn S, Gavini F, Magrez D, and Scheen A (2013) Issues in performing a network meta-analysis. Statistical Methods in Medical Research, 22 (2), Gert van Valkenhoef, and Joel Kuiper (2015) Gemtc: Network Meta-Analysis Using Bayesian Methods (Version 0.7-1). 54

56 2016/8/3

57 Backup Slides 2016/8/3

58 Q homogeneity consistency

59 Windows R CRAN R R exe 58

60 Windows R OK (N) > 59

61 Windows R 1. Message translations 2. PC 64-bit 64-bit Files 32-bit Files 32-bit Files (N) > 60

62 Windows R SDI 61

63 Windows R (N) > 62

64 Windows R R or 63

65 Windows R R R install.packages("netmeta", dep=t) Japan (Tokyo) OK R 64

untitled

untitled SAS 2016 SAS Network Meta-Analysis SAS Network Meta-Analysis homogeneity consistency Ranking Network Meta-Analysis(homogeneity) (consistency) Ranking 1 Network Meta-Analysis Network Meta-Analysis P Ranking

More information

Network Meta-Analysis

Network Meta-Analysis SAS ユーザー総会 2016 SAS での Network Meta-Analysis の実施例 ~ 頻度論に基づくアプローチ ~ 武田薬品工業株式会社日本開発センター生物統計室舟尾暢男 黒田晋吾 要旨 SAS による Network Meta-Analysis の実施方法として 均質性 (homogeneity) や一致性 (consistency) の評価 各薬剤の Ranking の算出方法を含めて紹介する

More information

Microsoft Word - StatsDirectMA Web ver. 2.0.doc

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

More information

/22 R MCMC R R MCMC? 3. Gibbs sampler : kubo/

/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. ( )

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 information

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

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

More information

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

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

More information

12/1 ( ) GLM, R MCMC, WinBUGS 12/2 ( ) WinBUGS WinBUGS 12/2 ( ) : 12/3 ( ) :? ( :51 ) 2/ 71

12/1 ( ) GLM, R MCMC, WinBUGS 12/2 ( ) WinBUGS WinBUGS 12/2 ( ) : 12/3 ( ) :? ( :51 ) 2/ 71 2010-12-02 (2010 12 02 10 :51 ) 1/ 71 GCOE 2010-12-02 WinBUGS kubo@ees.hokudai.ac.jp http://goo.gl/bukrb 12/1 ( ) GLM, R MCMC, WinBUGS 12/2 ( ) WinBUGS WinBUGS 12/2 ( ) : 12/3 ( ) :? 2010-12-02 (2010 12

More information

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

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

More information

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

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

More information

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

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

More information

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

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

More information

こんにちは由美子です

こんにちは由美子です 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

DAA09

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

More information

こんにちは由美子です

こんにちは由美子です 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

食道がん化学放射線療法後のsalvage手術

食道がん化学放射線療法後のsalvage手術 2006 2 17 52 Daly JM, et al. J Am Coll Surg 2000;190:562-573 Esophageal Cancer: Results of an American College of Surgeons Patient Care Evaluation Study Daly JM, et al. J Am Coll Surg 2000;190:562-573

More information

<30332D D8CB492982D FC95E490E690B62D8FAC8FBC96BE90E690B62E706466>

<30332D D8CB492982D FC95E490E690B62D8FAC8FBC96BE90E690B62E706466> A questionnaire survey on the use of placebo in clinical practice by nurses in Japan 1 1 2 2 1 4 1964 5, 6 1998 GCP 1 School of Nursing, Faculty of Medicine, Toho University 2Faculty of Medical Technology,

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 information

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

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)

More information

特 集 conditionally essential amino acid, グルタミンジペプチド (Glutamine dipeptide), 多臓器不全 Ⅰ はじめに sepsis RE- DOXS Ⅱ 重症患

特 集 conditionally essential amino acid, グルタミンジペプチド (Glutamine dipeptide), 多臓器不全 Ⅰ はじめに sepsis RE- DOXS Ⅱ 重症患 143 特 集 conditionally essential amino acid, グルタミンジペプチド (Glutamine dipeptide), 多臓器不全 Ⅰ はじめに 1980 2013 sepsis RE- 173 8605 2 11 1 Ⅱ 重症患者における生体内でのグルタミンの役割 1950 Eagle 1 1970 2 1/3 40 50 144 30 conditionally

More information

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

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

More information

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

2 The Bulletin of Meiji University of Integrative Medicine 3, Yamashita 10 11

2 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

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

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

More information

4703ALL01

4703ALL01 473201010 *** (1999) 16 2017 20023 2 1 2 1220 (1999a),(1999b) (1998), (2002) (2003)(1999)(2003) Conjoint Analysis Conjoint Analysis Willingness to PayStandard gamble Time trade-off Rating Scale Willingness

More information

1 Tokyo Daily Rainfall (mm) Days (mm)

1 Tokyo Daily Rainfall (mm) Days (mm) ( ) r-taka@maritime.kobe-u.ac.jp 1 Tokyo Daily Rainfall (mm) 0 100 200 300 0 10000 20000 30000 40000 50000 Days (mm) 1876 1 1 2013 12 31 Tokyo, 1876 Daily Rainfall (mm) 0 50 100 150 0 100 200 300 Tokyo,

More information

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

第11回:線形回帰モデルのOLS推定 11 OLS 2018 7 13 1 / 45 1. 2. 3. 2 / 45 n 2 ((y 1, x 1 ), (y 2, x 2 ),, (y n, x n )) linear regression model y i = β 0 + β 1 x i + u i, E(u i x i ) = 0, E(u i u j x i ) = 0 (i j), V(u i x i ) = σ 2, i

More information

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

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

More information

日本統計学会誌, 第44巻, 第2号, 251頁-270頁

日本統計学会誌, 第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

EAMS 2014

EAMS 2014 2014 2 10 RCT EAMS 2014 2015. 3. 31 Evidence Reports of Anma-Massage-Shiatsu: 2 Meta-Analisys and 10 Randomized Controlled Trials of Japan 31 Mar 2015 CONTENTS 1. (prologue).. 1 2. (steps for development

More information

untitled

untitled 19 CSPOR CRC/2009.8.8-9 2009 8 8 ( ) - CSPOR CRC SEMINAR- twatanab@oncoloplan.com http://www.oncoloplan.com - Toru Watanabe MD - 2 Phase I: Dose Finding trial Phase II : Efficacy and Safety trial Phase

More information

(lm) lm AIC 2 / 1

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

More information

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

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

More information

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

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

More information

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

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

More information

II) SIR SMR)

II) SIR SMR) II) ) 6 55 95 95Confidence Interval; 95CI) confidence intervals 95 95 1 5 95CI.66 Mortality experience of cockpit crewmembers from Japan Airlines. Kaji M, Tango T, Asukata I, Tajima N, Yamamoto K, Yamamoto

More information

2014ESJ.key

2014ESJ.key http://www001.upp.so-net.ne.jp/ito-hi/stat/2014esj/ Statistical Software for State Space Models Commandeur et al. (2011) Journal of Statistical Software 41(1) State Space Models in R Petris & Petrone (2011)

More information

kubo2015ngt6 p.2 ( ( (MLE 8 y i L(q q log L(q q 0 ˆq log L(q / q = 0 q ˆq = = = * ˆq = 0.46 ( 8 y 0.46 y y y i kubo (ht

kubo2015ngt6 p.2 ( ( (MLE 8 y i L(q q log L(q q 0 ˆq log L(q / q = 0 q ˆq = = = * ˆq = 0.46 ( 8 y 0.46 y y y i kubo (ht kubo2015ngt6 p.1 2015 (6 MCMC kubo@ees.hokudai.ac.jp, @KuboBook http://goo.gl/m8hsbm 1 ( 2 3 4 5 JAGS : 2015 05 18 16:48 kubo (http://goo.gl/m8hsbm 2015 (6 1 / 70 kubo (http://goo.gl/m8hsbm 2015 (6 2 /

More information

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

R John Fox R R R Console library(rcmdr) Rcmdr R GUI Windows R R SDI *1 R Console R 1 2 Windows XP Windows * 2 R R Console R ˆ R R John Fox 2006 8 26 2008 8 28 1 R R R Console library(rcmdr) Rcmdr R GUI Windows R R SDI *1 R Console R 1 2 Windows XP Windows * 2 R R Console R ˆ R GUI R R R Console > ˆ 2 ˆ Fox(2005) jfox@mcmaster.ca

More information

EBM Reviews Ovid Ovid EBM Reviews ACP Journal Club Cochrane Database of Systematic Revie

EBM Reviews Ovid   Ovid EBM Reviews ACP Journal Club Cochrane Database of Systematic Revie EBM Reviews 2006 4 Ovid http://www.usaco.co.jp/products/ovid/ovidplatform.html Ovid EBM Reviews ACP Journal Club Cochrane Database of Systematic Reviews CDSR Database of Abstracts of Reviews of Effectiveness

More information

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

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

More information

1999 02 03 Ver. 1.00 1999 07 09 Ver. 1.01 UID, 2003 10 01 Ver. 1.02 Storage SCU DICOM 3.0 Conformance Statement Ver.1.01 1999/07/09 1 0... 3 1... 3 1.1... 3 1.2... 3 1.3... 3 2... 4 2.1... 4 2.1.1... 4

More information

自由集会時系列part2web.key

自由集会時系列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 information

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

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

More information

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

More information

untitled

untitled 72 22 30 1985 Conductance Vessels flow Vessels< - Shear stress conductance vessel Exercise and Adenosine SPECT in the Patient with Angina SPECT,PET ACC/AHA/ASNC Guideline 2003 Dicarli JNM 2007;48:783

More information

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

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

More information

SERPWatcher SERPWatcher SERP Watcher SERP Watcher,

SERPWatcher SERPWatcher SERP Watcher SERP Watcher, SERPWatcher 112-8610 2-1-1 112-8610 2-1-1 229-8558 5-10-1 E-mail: nakabe@db.is.ocha.ac.jp, chiemi@is.ocha.ac.jp SERPWatcher SERP Watcher SERP Watcher, SERP Analysis of transition of ranking in SERP Watcher

More information

untitled

untitled 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

JMP V4 による生存時間分析

JMP 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

..,,...,..,...,,.,....,,,.,.,,.,.,,,.,.,.,.,,.,,,.,,,,.,,, Becker., Becker,,,,,, Becker,.,,,,.,,.,.,,

..,,...,..,...,,.,....,,,.,.,,.,.,,,.,.,.,.,,.,,,.,,,,.,,, Becker., Becker,,,,,, Becker,.,,,,.,,.,.,, J. of Population Problems. pp.,,,.,.,,. Becker,,.,,.,,.,,.,,,,.,,,.....,,. ..,,...,..,...,,.,....,,,.,.,,.,.,,,.,.,.,.,,.,,,.,,,,.,,, Becker., Becker,,,,,, Becker,.,,,,.,,.,.,, ,,, Becker,,., Becker,

More information

02[021-046]小山・池田(責)岩.indd

02[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 information

untitled

untitled 18 1 2,000,000 2,000,000 2007 2 2 2008 3 31 (1) 6 JCOSSAR 2007pp.57-642007.6. LCC (1) (2) 2 10mm 1020 14 12 10 8 6 4 40,50,60 2 0 1998 27.5 1995 1960 40 1) 2) 3) LCC LCC LCC 1 1) Vol.42No.5pp.29-322004.5.

More information

時間割引:双曲割引と弱加法性

時間割引:双曲割引と弱加法性 Discussion aper No. 666 June 2006 The Institute of Social and Economic Research Osaka University 6-1 Mihogaoka, Ibaraki, Osaka 567-0047, Japan Time Discounting: Declining Impatience and Interval Effect

More information

44 4 I (1) ( ) (10 15 ) ( 17 ) ( 3 1 ) (2)

44 4 I (1) ( ) (10 15 ) ( 17 ) ( 3 1 ) (2) (1) I 44 II 45 III 47 IV 52 44 4 I (1) ( ) 1945 8 9 (10 15 ) ( 17 ) ( 3 1 ) (2) 45 II 1 (3) 511 ( 451 1 ) ( ) 365 1 2 512 1 2 365 1 2 363 2 ( ) 3 ( ) ( 451 2 ( 314 1 ) ( 339 1 4 ) 337 2 3 ) 363 (4) 46

More information

i ii i iii iv 1 3 3 10 14 17 17 18 22 23 28 29 31 36 37 39 40 43 48 59 70 75 75 77 90 95 102 107 109 110 118 125 128 130 132 134 48 43 43 51 52 61 61 64 62 124 70 58 3 10 17 29 78 82 85 102 95 109 iii

More information

Fig. 1 Trends of TB incidence rates for all forms and smear-positive pulmonary TB in Kawasaki City and Japan. Incidence=newly notified cases of all fo

Fig. 1 Trends of TB incidence rates for all forms and smear-positive pulmonary TB in Kawasaki City and Japan. Incidence=newly notified cases of all fo Kekkaku Vol. 79, No. 1: 17-24, 2004 17 (Received 21 Aug. 2003/Accepted 18 Nov. 2003) Fig. 1 Trends of TB incidence rates for all forms and smear-positive pulmonary TB in Kawasaki City and Japan. Incidence=newly

More information

yamadaiR(cEFA).pdf

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

More information

2

2 Programming of Terrestrial Broadcasters and BS Digital Broadcasters by Genre ABSTRACT This study examined the programming of terrestrial and BS digital broadcasters by genre from the viewpoint of economics,

More information

Journal Club 2017/05/30 Tomohiro Matsumoto PGY-2

Journal Club 2017/05/30 Tomohiro Matsumoto PGY-2 Journal Club 2017/05/30 Tomohiro Matsumoto PGY-2 本 の論 March 21, 2017 2 Introduc)on 3 : (CIN) 3 25% 44µmol/L(0.5mg/dL) CIN 4 European Society of Urogenital Radiology. ESUR guidelines on contrast media,

More information

A5 PDF.pwd

A5 PDF.pwd Average Treatment Effect; ATE attributes Randomized Factorial Survey Experiment; RFSE cues ATE ATE Hainmueller et al. 2014 Average Marginal Component Effect ATE 67 4 2017 2 845 , ;, ATE, ;, ;, W 846 67

More information

Size 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 information

1) University Group Diabetes Program: A study of hypoglycemic agents on vascular complica- in patients with adult-onset tions diabetes. I. Design, methods and baseline results. Diabetes 19 (suppl. 2):

More information

29 Short-time prediction of time series data for binary option trade

29 Short-time prediction of time series data for binary option trade 29 Short-time prediction of time series data for binary option trade 1180365 2018 2 28 RSI(Relative Strength Index) 3 USD/JPY 1 2001 1 2 4 10 2017 12 29 17 00 1 high low i Abstract Short-time prediction

More information

自殺の経済社会的要因に関する調査研究報告書

自殺の経済社会的要因に関する調査研究報告書 17 1 2 3 4 5 11 16 30,247 17 18 21,024 +2.0 6 12 13 WHO 100 14 7 15 2 5 8 16 9 10 17 11 12 13 14 15 16 17 II I 18 Durkheim(1897) Hamermesh&Soss(1974)Dixit&Pindyck(1994) Becker&Posner(2004) Rosenthal(1993)

More information

SEJulyMs更新V7

SEJulyMs更新V7 1 2 ( ) Quantitative Characteristics of Software Process (Is There any Myth, Mystery or Anomaly? No Silver Bullet?) Zenya Koono and Hui Chen A process creates a product. This paper reviews various samples

More information

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

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

More information

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

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

More information

COE-RES Discussion Paper Series Center of Excellence Project The Normative Evaluation and Social Choice of Contemporary Economic Systems Graduate Scho

COE-RES Discussion Paper Series Center of Excellence Project The Normative Evaluation and Social Choice of Contemporary Economic Systems Graduate Scho COE-RES Discussion Paper Series Center of Excellence Project The Normative Evaluation and Social Choice of Contemporary Economic Systems Graduate School of Economics and Institute of Economic Research

More information

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

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

More information

日本製薬工業協会シンポジウム 生存時間解析の評価指標に関する最近の展開ー RMST (restricted mean survival time) を理解するー 2. RMST の定義と統計的推測 2018 年 6 月 13 日医薬品評価委員会データサイエンス部会タスクフォース 4 生存時間解析チー

日本製薬工業協会シンポジウム 生存時間解析の評価指標に関する最近の展開ー 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

JHN Journal Club 手稲渓仁会病院

JHN Journal Club 手稲渓仁会病院 JHOSPITALIST network GUIDE-IT trial; The Guiding Evidence Based Therapy Using Biomarker Intensified Treatment in Heart Failure multicentre randomised clinical trial 20171010 1 70 (EF 30%), (PCI), (HbA1c

More information

医学教育研究とは 2011_10_14 v1

医学教育研究とは 2011_10_14 v1 - - how to overview - overview - - - - - Harden RM (1986). Medical Education, 20: 522-531 PBL Description: What was done? 64% Justification: Did it work? 29% Clarification: Why or how did it work 7% Cook

More information

AJACS18_ ppt

AJACS18_ ppt 1, 1, 1, 1, 1, 1,2, 1,2, 1 1 DDBJ 2 AJACS3 2010 6 414:20-15:20 2231 DDBJ DDBJ DDBJ DDBJ NCBI (GenBank) DDBJ EBI (EMBL-Bank) GEO DDBJ Omics ARchive(DOR) ArrayExpress DTA (DDBJ Trace Archive) DRA (DDBJ

More information

Ⅰ. 緒言 Suzuki, et al., Ⅱ. 研究方法 1. 対象および方法 1 6 表 1 1, 調査票の内容 図

Ⅰ. 緒言 Suzuki, et al., Ⅱ. 研究方法 1. 対象および方法 1 6 表 1 1, 調査票の内容 図 The Journal of the Japan Academy of Nursing Administration and Policies Vol 12, No 1, pp 515, 2008 原著 病棟勤務看護師の勤務状況とエラー ニアミスのリスク要因 Relationship between Risk Factors Related Medical Errors and Working Conditions

More information

\\afs001-0m0005\project02\A32\M

\\afs001-0m0005\project02\A32\M Technical Information 2004.09 2009.04 Store Request Query Request Retrieve Request DICOM Client Application Remote SCP Remote Query/Retrieve SCP Image Stored * DICOM Server Application Remote SCU Print

More information

5 Armitage x 1,, x n y i = 10x i + 3 y i = log x i {x i } {y i } 1.2 n i i x ij i j y ij, z ij i j 2 1 y = a x + b ( cm) x ij (i j )

5 Armitage x 1,, x n y i = 10x i + 3 y i = log x i {x i } {y i } 1.2 n i i x ij i j y ij, z ij i j 2 1 y = a x + b ( cm) x ij (i j ) 5 Armitage. x,, x n y i = 0x i + 3 y i = log x i x i y i.2 n i i x ij i j y ij, z ij i j 2 y = a x + b 2 2. ( cm) x ij (i j ) (i) x, x 2 σ 2 x,, σ 2 x,2 σ x,, σ x,2 t t x * (ii) (i) m y ij = x ij /00 y

More information

10_08.dvi

10_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

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

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

More information

untitled

untitled The Japanese Association for Behavior Analysis Newsletter, No. 35 The Japanese Association for Behavior Analysis Newsletter, No. 35 The Japanese Association for Behavior Analysis Newsletter, No. 35 Table

More information

Study on Throw Accuracy for Baseball Pitching Machine with Roller (Study of Seam of Ball and Roller) Shinobu SAKAI*5, Juhachi ODA, Kengo KAWATA and Yu

Study on Throw Accuracy for Baseball Pitching Machine with Roller (Study of Seam of Ball and Roller) Shinobu SAKAI*5, Juhachi ODA, Kengo KAWATA and Yu Study on Throw Accuracy for Baseball Pitching Machine with Roller (Study of Seam of Ball and Roller) Shinobu SAKAI*5, Juhachi ODA, Kengo KAWATA and Yuichiro KITAGAWA Department of Human and Mechanical

More information

181 第 54 回土木計画学研究発表会 講演集 GPS S-502W S-502W

181 第 54 回土木計画学研究発表会 講演集 GPS S-502W S-502W 181 GPS 1 2 3 4 1 98-845 468-1 S-52W E-mail: h-ymgc@plan.civil.tohoku.ac.jp 2 98-845 468-1 S-52W E-mail: mokmr@m.tohoku.ac.jp 3 18-626 2 15-3 C 6F E-mail: h kaneda@zenrin-datacom.net 4 18-626 2 15-3 C

More information

生命倫理100_資料4-7

生命倫理100_資料4-7 26 27 Single cell clones picked from swaps12p36 (18) Single cell clones picked from swaps12p18 (24) 1.3% (p0) 1.1% (p58) 28 29 Outline 30 31 32 33 34 Columbia, The New York Stem Cell Foundation Oocyte

More information

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

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

More information

yasi10.dvi

yasi10.dvi 2002 50 2 259 278 c 2002 1 2 2002 2 14 2002 6 17 73 PML 1. 1997 1998 Swiss Re 2001 Canabarro et al. 1998 2001 1 : 651 0073 1 5 1 IHD 3 2 110 0015 3 3 3 260 50 2 2002, 2. 1 1 2 10 1 1. 261 1. 3. 3.1 2 1

More information

デジタルオルソ作成マニュアルに関する調査研究作業

デジタルオルソ作成マニュアルに関する調査研究作業 A1-No.289 16 1 ...1....1....1....4...5...11... 11...16...16...16...18...22...24...26...29...31...31....31. (TIN)...31....31....31...41....41....43 . 14 20 IT. GZ-1 OR-1 1 1 1 2 2 2 3 A A B A a a a a 3

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

橡

橡 TP T-13i AST ALT ALP Y CTP LOH CHE CRP 8.7 id 0.4 mg/d 21 AWL 131U/L 257 LU/L 12 111/L 234 11J/L 344 IU/L 0.30 mid No K a BUN ORE UA 142 meq/1 4.3 meqa 108 meqa 15.8 met! 0.71 mg/di 5.0 mg/d CEA 5_5 n8/m1

More information

IPSJ 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

IPSJ 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 information

Presentation Title Goes Here

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

More information

66-1 田中健吾・松浦紗織.pwd

66-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 information

untitled

untitled twatanab@oncoloplan.com http://www.oncoloplan.com I II - III IV Fig 3. Survival curves overall and according to response Bruzzi, P. et al. J Clin Oncol; 23:5117-5125 25 Copyright merican Society of Clinical

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

Kaplan-Meierプロットに付加情報を追加するマクロの作成

Kaplan-Meierプロットに付加情報を追加するマクロの作成 Kaplan-Meier 1, 2,3 1 2 3 A SAS macro for extended Kaplan-Meier plots Kengo Nagashima 1, Yasunori Sato 2,3 1 Department of Parmaceutical Technochemistry, Josai University 2 School of Medicine, Chiba University

More information

36:378 第 38 回日本脳卒中学会講演シンポジウム 原著 36: , 要旨 TIA 2 t-pa Key words: stroke registry, stroke subtype, onset-visi

36:378 第 38 回日本脳卒中学会講演シンポジウム 原著 36: , 要旨 TIA 2 t-pa Key words: stroke registry, stroke subtype, onset-visi 36:378 第 38 回日本脳卒中学会講演シンポジウム 原著 36: 378 384, 2014 1 2 要旨 1999 2012 10 31 29 26 80 30 TIA 2 t-pa Key words: stroke registry, stroke subtype, onset-visit time, chronological change はじめに 4 12 23 27 1 Japan

More information

講義のーと : データ解析のための統計モデリング. 第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

More information

kubostat7f p GLM! logistic regression as usual? N? GLM GLM doesn t work! GLM!! probabilit distribution binomial distribution : : β + β x i link functi

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

More information