Upgrade to Pro
— share decks privately, control downloads, hide ads and more …
Speaker Deck
Features
Speaker Deck
PRO
Sign in
Sign up for free
Search
Search
R を用いた検定(補講) (1) — Welch 検定 / Other Tests Using...
Search
Kenji Saito
PRO
January 18, 2024
Business
0
100
R を用いた検定(補講) (1) — Welch 検定 / Other Tests Using R (1) - Welch Test
早稲田大学大学院経営管理研究科「企業データ分析」2023 冬のオンデマンド教材 第9回で使用したスライドです。
Kenji Saito
PRO
January 18, 2024
Tweet
Share
More Decks by Kenji Saito
See All by Kenji Saito
続・インクルーシブな社会へ / Continuing Towards an Inclusive Society
ks91
PRO
0
17
AGI (人工一般知能) と創る新しく奇妙な社会 / New and Stranger Society built with AGI
ks91
PRO
0
64
回帰分析/大規模言語モデルと統計 / Regression Analysis, Large Language Models and Statistics
ks91
PRO
0
69
多重比較/相関分析 / Multiple Comparison and Correlation Analysis
ks91
PRO
0
65
アカデミーキャンプ 2025冬「考えるのは奴らだ」 / Academy Camp 2025 Winter - Live and Let Think DAY 3
ks91
PRO
0
61
アカデミーキャンプ 2025冬「考えるのは奴らだ」 / Academy Camp 2025 Winter - Live and Let Think DAY 2
ks91
PRO
0
48
アカデミーキャンプ 2025冬「考えるのは奴らだ」 / Academy Camp 2025 Winter - Live and Let Think DAY 1
ks91
PRO
1
74
インクルーシブな社会へ / Toward an Inclusive Society
ks91
PRO
0
22
P 値と有意差/分散分析 / P-value, Significant Difference and Analysis of Variance
ks91
PRO
0
71
Other Decks in Business
See All in Business
株式会社リブセンス 会社説明資料(報道関係者様向け)
livesense
PRO
0
1k
Spice Factory Co., Ltd. Culture Deck
spicefactory
0
470
#TRG24 / Daniel Sanchez-Crespo / Cómo cargarte una industria...
tarugoconf
0
1.1k
RAKSUL会社紹介資料 / RAKSUL Introduction
raksulrecruiting
9
83k
株式会社shizai - Recruit Deck
shizai
3
56k
株式会社AbemaTV 会社説明資料
abematv
2
1.6k
因果推論が浸透した組織の現状と未来 / The Present and Future of Organizations Embracing Causal Inference
yusukekayahara
0
680
エンジニア→PM進化論
natty_natty254
2
200
株式会社B4A 会社紹介
b4a
0
6.7k
メドピアグループ紹介資料
medpeer_recruit
10
120k
SHIFT ASIA 会社説明資料 V2.1
shiftasiarec
0
250
【業界・業種別】副業・兼業トラブルに関する実態調査
fkske
0
210
Featured
See All Featured
XXLCSS - How to scale CSS and keep your sanity
sugarenia
248
1.3M
VelocityConf: Rendering Performance Case Studies
addyosmani
328
24k
Bootstrapping a Software Product
garrettdimon
PRO
306
110k
Building Applications with DynamoDB
mza
93
6.2k
Testing 201, or: Great Expectations
jmmastey
42
7.2k
Templates, Plugins, & Blocks: Oh My! Creating the theme that thinks of everything
marktimemedia
30
2.2k
Product Roadmaps are Hard
iamctodd
PRO
50
11k
The Web Performance Landscape in 2024 [PerfNow 2024]
tammyeverts
4
430
Why You Should Never Use an ORM
jnunemaker
PRO
55
9.2k
Building Adaptive Systems
keathley
40
2.4k
Building a Modern Day E-commerce SEO Strategy
aleyda
38
7.1k
RailsConf & Balkan Ruby 2019: The Past, Present, and Future of Rails at GitHub
eileencodes
133
33k
Transcript
generated by Stable Diffusion XL v1.0 2023 9 R (
) (1) — Welch (WBS) 2023 9 R ( ) (1) — Welch — 2024-01 – p.1/10
https://speakerdeck.com/ks91/collections/corporate-data-analysis-2023-winter 2023 9 R ( ) (1) — Welch —
2024-01 – p.2/10
( 20 ) 1 • 2 R • 3 •
4 • 5 • 6 ( ) • 7 (1) • 8 (2) • 9 R ( ) (1) — Welch • 10 R ( ) (2) — χ2 11 R ( ) (1) — 12 R ( ) (2) — 13 GPT-4 14 GPT-4 15 ( ) LaTeX Overleaf 8 (12/21 ) / (2 ) OK / 2023 9 R ( ) (1) — Welch — 2024-01 – p.3/10
t Welch R t.test() Welch Welch 2023 9 R (
) (1) — Welch — 2024-01 – p.4/10
2 t ( ) (1/2) 2 ( ) xA −
xB (1) : (2) : σ ( ) σ sp sp = s2 A (nA − 1) + s2 B (nB − 1) nA + nB − 2 (R var() ) nA + nB − 2 t Welch A B (µA = µB ) A B (µA = µB ) 2023 9 R ( ) (1) — Welch — 2024-01 – p.5/10
2 t ( ) (2/2) xA − xB Student µA
= µB t = (xA − xB ) − (µA − µB ) sp 1 nA + 1 nB = xA − xB sp 1 nA + 1 nB (t ) t dfp = nA + nB − 2 t ( ) t0.05 (dfp ) t0.05 (dfp ) < |t| (P < 0.05) 2023 9 R ( ) (1) — Welch — 2024-01 – p.6/10
Welch t t = xA − xB s2 A nA
+ s2 B nB ( ) v . . . v ≈ ( s2 A nA + s2 B nB )2 s4 A n2 A (nA−1) + s4 B n2 B (nB−1) R 2023 9 R ( ) (1) — Welch — 2024-01 – p.7/10
( ) - (Shapiro-Wilk test) - (Anderson-Darling test for normality)
- (Kolmogorov-Smirnov test for normality) ( ) ( ) (Bartlett’s test for homogeneity of variances) 2023 9 R ( ) (1) — Welch — 2024-01 – p.8/10
.txt A /B g <- read.table(" .txt", header=T) colnames(g) <-
c(" ", " ") sampleA <- g$ sampleB <- g$ # ( ) shapiro.test(x=sampleA) shapiro.test(x=sampleB) # ( ) samples <- c(sampleA, sampleB) group_factor <- factor(rep(c("A", "B"), c(length(sampleA), length(sampleB)))) bartlett.test(formula=samples~group_factor) # Welch (t.test() ) ( Welch ) t.test(sampleA, sampleB) 2023 9 R ( ) (1) — Welch — 2024-01 – p.9/10
2023 9 R ( ) (1) — Welch — 2024-01
– p.10/10