物理学者(ポスドク)による日々の研究生活のメモ書きです ( python/emacs/html/Japascript/シェルスクリプト/TeX/Mac/C言語/Linux/git/tmux/R/ポケモンGO)
■ 2つのデータに対してKS検定をする場合
■ 1つのデータに対した正規分布のKS検定をする場合
2つの引数は、比較したい分布を指定する文字列、正規分布なら pnorm、ポアソン分布なら ppois、一様分布の場合 punif を指定できる
■ 正規分布の場合だと、平均と分散も指定可能
■ 片側検定/両側検定の切り替えはalternativaで行う
デフォルトでは両側検定になっている
■ 検定結果からp値を取り出したいときは、
結果をなにかに代入して、namesで中身を調べる→ p値のみを取り出す
ツイート
% x=rnorm(10, mean = 1)
% y=rnorm(10, mean = 1)
% x
[1] 0.91430534 1.59262261 1.52565219 1.68065990 2.35755910 1.91386531
[7] -2.26359817 0.08465784 0.06163455 -0.10943752
% y
[1] 2.2200235 2.2241460 1.0733582 0.8002340 1.6228854 0.2652880
[7] 1.4286084 -0.6869875 0.1918104 2.2083442
% ks.test(x, y)
Two-sample Kolmogorov-Smirnov test
data: x and y
D = 0.3, p-value = 0.7869
alternative hypothesis: two-sided
% y=rnorm(10, mean = 1)
% x
[1] 0.91430534 1.59262261 1.52565219 1.68065990 2.35755910 1.91386531
[7] -2.26359817 0.08465784 0.06163455 -0.10943752
% y
[1] 2.2200235 2.2241460 1.0733582 0.8002340 1.6228854 0.2652880
[7] 1.4286084 -0.6869875 0.1918104 2.2083442
% ks.test(x, y)
Two-sample Kolmogorov-Smirnov test
data: x and y
D = 0.3, p-value = 0.7869
alternative hypothesis: two-sided
■ 1つのデータに対した正規分布のKS検定をする場合
% x=rnorm(10, mean = 1)
% ks.test(x, pnorm)
One-sample Kolmogorov-Smirnov test
data: x
D = 0.43645, p-value = 0.02996
alternative hypothesis: two-sided
% ks.test(x, pnorm)
One-sample Kolmogorov-Smirnov test
data: x
D = 0.43645, p-value = 0.02996
alternative hypothesis: two-sided
2つの引数は、比較したい分布を指定する文字列、正規分布なら pnorm、ポアソン分布なら ppois、一様分布の場合 punif を指定できる
■ 正規分布の場合だと、平均と分散も指定可能
% ks.test(x=vx,y="pnorm",mean=mean(vx),sd=sd(vx))
One-sample Kolmogorov-Smirnov test
data: vx
D = 0.13793, p-value = 0.9196
alternative hypothesis: two-sided
One-sample Kolmogorov-Smirnov test
data: vx
D = 0.13793, p-value = 0.9196
alternative hypothesis: two-sided
■ 片側検定/両側検定の切り替えはalternativaで行う
デフォルトでは両側検定になっている
% ks.test(x, y, , alternative="t")
Two-sample Kolmogorov-Smirnov test
data: x and y
D = 0.3, p-value = 0.7869
alternative hypothesis: two-sided
% ks.test(x, y, , alternative="l")
Two-sample Kolmogorov-Smirnov test
data: x and y
D^- = 0.1, p-value = 0.9048
alternative hypothesis: the CDF of x lies below that of y
% ks.test(x, y, , alternative="g")
Two-sample Kolmogorov-Smirnov test
data: x and y
D^+ = 0.3, p-value = 0.4066
alternative hypothesis: the CDF of x lies above that of y
Two-sample Kolmogorov-Smirnov test
data: x and y
D = 0.3, p-value = 0.7869
alternative hypothesis: two-sided
% ks.test(x, y, , alternative="l")
Two-sample Kolmogorov-Smirnov test
data: x and y
D^- = 0.1, p-value = 0.9048
alternative hypothesis: the CDF of x lies below that of y
% ks.test(x, y, , alternative="g")
Two-sample Kolmogorov-Smirnov test
data: x and y
D^+ = 0.3, p-value = 0.4066
alternative hypothesis: the CDF of x lies above that of y
■ 検定結果からp値を取り出したいときは、
結果をなにかに代入して、namesで中身を調べる→ p値のみを取り出す
% out <- ke.test(x, y)
% names(out)
[1] "statistic" "p.value" "alternative" "method" "data.name"
% out$p.value
[1] 0.7869298
% names(out)
[1] "statistic" "p.value" "alternative" "method" "data.name"
% out$p.value
[1] 0.7869298
PR
この記事にコメントする
プロフィール
HN:coffee
職業:物理屋(自称)
趣味:映画鑑賞、登山
出身:大阪府の南の田舎
自己紹介:
import MyProfile
import coffee_pote from TWITTER
import amazonのほしい物リスト from WISH_LIST
print "先月子供が産まれました!"
# 最終更新 2022/10/25
職業:物理屋
趣味:映画鑑賞、登山
出身:大阪府の南の田舎
自己紹介:
import MyProfile
import coffee_pote from TWITTER
import amazonのほしい物リスト from WISH_LIST
print "先月子供が産まれました!"
# 最終更新 2022/10/25
カテゴリー
ブログ内検索
リンク
相互リンク募集中です
(Twitterにてお知らせください)
Demo scripts for gnuplot version 5 (gnuplotのさまざまなデモ画像と作り方がまとめられている、眺めているだけでできるようになった気分になれる)
gnuplotスクリプトの解説 (米澤進吾さんの個人ページ、gnuplotと言えばこのかた)
gnuplot のページ (Takeno Lab、うちのブログがリンクされていたのでリンク返し)
Twitterから映画の評価が分かる & 映画の鑑賞記録が残せる coco (映画の感想をまとめられるサイト、いつもお世話になっています)
Astronomy Picture of the Day Archive (天文や宇宙関連の最新の話題について画像とともにNASAが説明しているページ)
今日のほしぞら (任意の時刻の空で見える星を表示してくれる、国立天文台が管理している)
GNUPLOTとアニメーション (応用の項目の「見せてあげよう!ラピュタの雷を!!」あたりからすごすぎる)
読書メーター (読んだ本をリストできる便利なサイト)
flickr難民の写真置き場 (20XX年、flickrは有料化の炎に包まれた。あらゆるflickr無料ユーザーは絶滅したかに見えた。 しかし、tumblr移住民は死に絶えてはいなかった。)
教授でもできるMac OS X へのLaTeX, X11, gccのインストレーションと環境設定 (阪大の山中卓さんのwebページ、タイトルにセンスが溢れている、内容は超充実してる、特にTeX関連、学振DCとかPDの申請書類作成時にはお世話になっております)
英語論文執筆用の例文検索サービス (とんでもないものを見つけてしまった・・・・ arXivに収録されている 811,761報の 英語論文から,例文を検索するための検索エンジン)
Demo scripts for gnuplot version 5 (gnuplotのさまざまなデモ画像と作り方がまとめられている、眺めているだけでできるようになった気分になれる)
gnuplotスクリプトの解説 (米澤進吾さんの個人ページ、gnuplotと言えばこのかた)
gnuplot のページ (Takeno Lab、うちのブログがリンクされていたのでリンク返し)
Twitterから映画の評価が分かる & 映画の鑑賞記録が残せる coco (映画の感想をまとめられるサイト、いつもお世話になっています)
Astronomy Picture of the Day Archive (天文や宇宙関連の最新の話題について画像とともにNASAが説明しているページ)
今日のほしぞら (任意の時刻の空で見える星を表示してくれる、国立天文台が管理している)
GNUPLOTとアニメーション (応用の項目の「見せてあげよう!ラピュタの雷を!!」あたりからすごすぎる)
読書メーター (読んだ本をリストできる便利なサイト)
flickr難民の写真置き場 (20XX年、flickrは有料化の炎に包まれた。あらゆるflickr無料ユーザーは絶滅したかに見えた。 しかし、tumblr移住民は死に絶えてはいなかった。)
教授でもできるMac OS X へのLaTeX, X11, gccのインストレーションと環境設定 (阪大の山中卓さんのwebページ、タイトルにセンスが溢れている、内容は超充実してる、特にTeX関連、学振DCとかPDの申請書類作成時にはお世話になっております)
英語論文執筆用の例文検索サービス (とんでもないものを見つけてしまった・・・・ arXivに収録されている 811,761報の 英語論文から,例文を検索するための検索エンジン)
最新記事
(11/20)
(05/17)
(05/17)
(05/09)
(03/05)
(02/29)
(02/21)
(02/21)
(02/21)
(02/21)