[過去ﾛｸﾞ] 当直医のスレ Part24 [無断転載禁止]©2ch.net (1002ﾚｽ)
上下前次1-新
抽出解除 ﾚｽ栞
このｽﾚｯﾄﾞは過去ﾛｸﾞ倉庫に格納されています｡
次ｽﾚ検索 歴削→次ｽﾚ 栞削→次ｽﾚ 過去ﾛｸﾞﾒﾆｭｰ
1: 2017/09/23(土)05:53 ID:LzeuIIdY(1/7) AAS
extend:checked:vvvvv:1000:512
こんばんは、当直医さん
ご苦労さまです。
ゆっくり寝てくださいね
問題は過去を克服することではありません。
さようなことができるわけはありません。
後になって過去を変えたり、
起こらなかったことにするわけにはまいりません。
しかし過去に目を閉ざす者は結局のところ現在にも盲目となります。(ヴァイツゼッカー)
画像ﾘﾝｸ[jpg]:imgur.com
省19
3(7): 2017/09/23(土)05:53 ID:LzeuIIdY(3/7) AAS
当直の時間帯は飲酒もしないので勉強時間がとれていい。
ド底辺シリツ医大卒様も当直でのイベント報告すればいいのに。
繰り返しになりますが、当直でのイベント報告です。
当直バイト医のド底辺シリツ医の顛末ご報告。敬称ｗ略。
起： 当直室のパソコンにmalwareの感染が発見される。画像ﾘﾝｸ[jpg]:i.imgur.com
承： 感染した日が表示されているので誰が当直だったか調べてみた。ド底辺シリツ医のバイトの日であった。
省3
8: 2017/09/23(土)18:52 ID:LzeuIIdY(5/7) AAS
法則を忘れている向きがあるので復習しましょう。
ド底辺シリツ医大卒の第一法則：
ド底辺シリツ医大が悪いのではない、本人の頭が悪いんだ。
ド底辺シリツ医大卒の第二法則：
ド底辺シリツ医大卒は恥ずかしくて、学校名を皆さま言いません。
ド底辺シリツ医大卒の第三法則
ド底辺シリツ医大に入学すると馬鹿になる。
第三法則の参考資料
画像ﾘﾝｸ[jpg]:i.imgur.com
省12
59(1): 2017/09/24(日)11:32 ID:vQO8g+Lx(2/20) AAS
熱傷診療ガイドライン(www.jsbi-burn.org/members/guideline/pdf/guideline2.pdf)に引用された論文としてこんなのが商用パンフに掲載されていた。
画像ﾘﾝｸ[jpg]:imagizer.imageshack.com
早期群: n1=17 m1=36.2 sd1=21.5
後期群: n2=12 m2=56.1 sd2=30.1
で早期群の方が上皮化期間が有意(t検定でp<0.05)に短いと主張している。
生データがないので標本数・平均・標準偏差から検定してみる。
ttest2 = function(n1,m1,s1, n2,m2,s2){ # ni:標本数 mi:平均 si:標準偏差
n=n1+n2-2
u=((n1-1)*s1^2+(n2-1)*s2^2)/n
t=(m1-m2)/sqrt(u/n1+u/n2)
省17
91: 2017/09/24(日)18:30 ID:vQO8g+Lx(15/20) AAS
外部ﾘﾝｸ[html]:oku.edu.mie-u.ac.jpの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]
省19
92: 2017/09/24(日)18:31 ID:rkMoOsfQ(1) AAS
どっちのほうが美人だと思う？
画像ﾘﾝｸ[jpg]:scontent.cdninstagram.com
マジレスお願い
93: 2017/09/24(日)18:31 ID:vQO8g+Lx(16/20) AAS
画像ﾘﾝｸ[png]:i.imgur.com
94: 2017/09/24(日)18:32 ID:vQO8g+Lx(17/20) AAS
外部ﾘﾝｸ[html]:oku.edu.mie-u.ac.jpの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]
省19
95: 2017/09/24(日)19:20 ID:vQO8g+Lx(18/20) AAS
小標本の場合には，より良い近似として，(n1*v1+n2*v2)^2/((n1*v1)^2/(n1-1)+(n2*v2)^2/(n2-1))t 分布で検定することをBrunnerとMunzelは提案しています。n1=10, n2=20の小標本としてこれを検証してみる
bm2 <- 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]
省21
96: 2017/09/24(日)19:21 ID:vQO8g+Lx(19/20) AAS
小標本の場合には，より良い近似として，自由度(n1*v1+n2*v2)^2/((n1*v1)^2/(n1-1)+(n2*v2)^2/(n2-1))のt 分布で検定することをBrunnerとMunzelは提案しています。n1=10, n2=20の小標本としてこれを検証してみる
bm2 <- 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]
省21
97: 2017/09/24(日)22:51 ID:vQO8g+Lx(20/20) AAS
ド底辺特殊シリツ医大は1学年100人として全学生600人とする。
このうち無作為に20人を抽出調査したところ全員裏口入学であった。
全学年での正規入学者数の期待値とその95%信頼区間を求めよ。
理論期待値＝26.36人
信頼区間は算出方式により下記の通り
method lower upper
agresti-coull 0 113.89
asymptotic 0 0.00
bayes 0 54.29
cloglog 0 101.06
省16
102: 2017/09/25(月)06:50 ID:1DSpMtz1(2/7) AAS
Ｒでスクリプトを組んでシミュレーションしてみた。
> atari(600,20,1,col='red')
Min. 1st Qu. Median Mean 3rd Qu. Max.
1.00 27.00 46.00 53.63 73.00 241.00
2.5% 5% 50% 95% 97.5%
7 10 46 122 140
53.63で理論値とほぼ一致している。
画像ﾘﾝｸ[png]:i.imgur.com
省4
105: 2017/09/25(月)07:37 ID:1DSpMtz1(5/7) AAS
法則を忘れている向きがあるので復習しましょう。
ド底辺シリツ医大卒の第一法則：
ド底辺シリツ医大が悪いのではない、本人の頭が悪いんだ。
ド底辺シリツ医大卒の第二法則：
ド底辺シリツ医大卒は恥ずかしくて、学校名を皆さま言いません。
ド底辺シリツ医大卒の第三法則
ド底辺シリツ医大に入学すると馬鹿になる。
第三法則の参考資料
画像ﾘﾝｸ[jpg]:i.imgur.com
省7
126(1): 2017/09/25(月)18:01 ID:v6ZqxKdw(2/2) AAS
お尻のことを「むきたてのゆで玉子のような」て形容してたのは
医師真似事務員おじいちゃんなんだから！
画像ﾘﾝｸ[jpg]:blog-imgs-87....
132: 2017/09/25(月)21:45 ID:uD79ucO8(11/15) AAS
>103の尤度対数関数をグラフで書いてやったぞ。
６００×（１/２０） = ３０ で最大値になっているのがわかるだろ。
ド底辺特殊シリツ医大卒の裏口バカには尤度対数関数の意味もわからんだろう。
画像ﾘﾝｸ[png]:i.imgur.com
155: 2017/09/26(火)10:52 ID:dTnYYVrv(5/9) AAS
シミュレーション結果をグラフにしてみた。
画像ﾘﾝｸ[jpg]:i.imgur.com
> ottery(100,50,10^4)
Min. 1st Qu. Median Mean 3rd Qu. Max.
1.000 1.000 2.000 2.838 4.000 9.000
2.5% 5% 50% 95% 97.5%
1.000 1.000 2.000 7.000 7.975
理論値が2.923077だから、概ね合致している。
160: 2017/09/26(火)13:07 ID:dTnYYVrv(7/9) AAS
>>156
何人めに発見されたかと、100人中の該当者数の期待値の関係をグラフにしてみた。
画像ﾘﾝｸ[jpg]:i.imgur.com
昼休みの暇つぶしになった。
162(1): 2017/09/26(火)14:57 ID:dTnYYVrv(8/9) AAS
さて、午後のCF3件終了。
３件目の開始前にRを走らせておいたら計算終了していた。
何人めに発見されたかと、100人中の該当者数の期待値の関係のグラフに９５％信頼区間を重ねてみた。
画像ﾘﾝｸ[jpg]:i.imgur.com
早いうちに該当者が発見されると期待値の信頼区間が広いのがみてとれる。
省2
188(1): 2017/09/27(水)11:39 ID:By0OihfH(5/28) AAS
>>179
両女子大の確率密度関数を書いてみた。
画像ﾘﾝｸ[png]:i.imgur.com
203: 2017/09/27(水)15:44 ID:By0OihfH(13/28) AAS
ド底辺特殊シリツ医大卒のために、確率密度関数のグラフを作ってやったぞ。
画像ﾘﾝｸ[png]:i.imgur.com
上下前次1-新書関写板覧索設栞歴
ｽﾚ情報 赤ﾚｽ抽出 画像ﾚｽ抽出 歴の未読ｽﾚ AAｻﾑﾈｲﾙ
ぬこの手 ぬこTOP 1.726s*