[過去ﾛｸﾞ] 当直医のスレ Part24 [無断転載禁止]©2ch.net (974ﾚｽ)
上下前次1-新

---

このｽﾚｯﾄﾞは過去ﾛｸﾞ倉庫に格納されています｡
次ｽﾚ検索 歴削→次ｽﾚ 栞削→次ｽﾚ 過去ﾛｸﾞﾒﾆｭｰ

---

94: 2017/09/24(日)18:32 ID:vQO8g+Lx(17/20) AA×

外部ﾘﾝｸ[html]:oku.edu.mie-u.ac.jp
画像ﾘﾝｸ[png]:i.imgur.com

[240|320|480|600|100%|JPG|べ|ﾚｽ栞|ﾚｽ消]
94: [sage] 2017/09/24(日) 18:32:14.64 ID:vQO8g+Lx https://oku.edu.mie-u.ac.jp/~okumura/stat/brunner-munzel.htmlのBrunner-Munzel検定の解説に 漸近的にn1やn2 が大きくなると），統計量WBFNは、標準正規分布 N(0,1) に近づきます、と書いてあるのでスクリプトを書いて検証してみる。 bm <- function(X1,X2){ # two.sided X1=na.omit(X1) X2=na.omit(X2) n1=length(X1) n2=length(X2) n12=n1+n2 R=rank(c(X1,X2)) R1=R[1:n1] R2=R[(n1+1):n12] R1_=mean(R1) R2_=mean(R2) r1=rank(X1) r2=rank(X2) v1=sum((R1-r1-R1_+(n1+1)/2)^2)/(n1-1) v2=sum((R2-r2-R2_+(n2+1)/2)^2)/(n2-1) sig1=v1/((n12-n1)^2) sig2=v2/((n12-n2)^2) sigN=n12*(sig1/n1+sig2/n2) (R2_-R1_)/(sqrt(sigN*n12)) } N=10^4 ; x1=rnorm(N) ; x2=rnorm(N,0.05,5) n1=100 ; n2=200 foo <- function() bm(sample(x1,n1),sample(x2,n2)) BM=replicate(10^4,foo()[1]) hist(BM,freq=FALSE,col='lightblue',main='statistics of Brunner-Munzel') #統計量WBFNのヒストグラムに curve(dnorm(x),add=TRUE,lwd=2) # N(0,1)の正規分布を重ねる http://i.imgur.com/GqJlZET.png　記述どおりにグラフが書けた。 http://egg.5ch.net/test/read.cgi/hosp/1506113587/94
の検定の解説に 漸近的にや が大きくなると統計量は標準正規分布 に近づきますと書いてあるのでスクリプトを書いて検証してみる 統計量のヒストグラムに の正規分布を重ねる 記述どおりにグラフが書けた

---

上下前次1-新書関写板覧索設栞歴

---

あと 880 ﾚｽあります
ｽﾚ情報 赤ﾚｽ抽出 画像ﾚｽ抽出 歴の未読ｽﾚ AAｻﾑﾈｲﾙ

---

ぬこの手 ぬこTOP 0.040s