底辺私立医大を卒業し ..
[2ch|▼Menu]
2:卵の名無しさん
18/05/19 06:10:57.97 rrC4yXIM.net
これが現実
1次方程式もできないド底辺特殊シリツ医大卒の記録
URLリンク(imagizer.imageshack.com)

何度読んでも馬鹿すぎる。
男女別の割合と全体での割合から男女比が計算できるとも思わないとは。
なんでこんなのが大学に入れるわけよ?
裏口入学以外に説明がつく?
中学生でも解ける一次方程式の問題だろ。
それすらできない馬鹿が自信を持って発言。
>患者の男女比が必要なのもわからないのか?
だとさ。
URLリンク(imagizer.imageshack.com)

0.2575549
と答を書いてやったら
>単位も書かずに答えだとか…
ド底辺シリツ医大では確率に単位があるらしいぞwww
何でこんな馬鹿が大学に入れるわけ?
裏口入学以外に説明がつく?
URLリンク(imagizer.imageshack.com)

3:卵の名無しさん
18/05/19 06:12:14.77 rrC4yXIM.net
このような国民世論が投稿されております。
医者が信用ならない理由 [転載禁止](c)2ch.net
スレリンク(hosp板:47番)
47 名前:卵の名無しさん[] 投稿日:2016/06/21(火) 05:34:59.81 ID:Eoye1u41
出身大学の明示を義務化して欲しい
帝京、金沢医、埼玉医、川崎医出身の医者はお断りしたい
できれば、受験歴も知りたい

底辺シリツ医大卒は高卒の学力すら欠いている。
馬鹿が答えて自爆した資料(その一)
URLリンク(i.imgur.com)
(その二)全く進歩の跡がみられません!
URLリンク(i.imgur.com)

4:卵の名無しさん
18/05/20 14:33:32.61 gHqbdBVY.net
「帝京大医学部の女の子がテラスハウス出てたけど、
その肩書きがテロップに出てる限り、「私はめちゃめちゃバカだけど、
親に大金払ってもらって医者になるよ!」って
吹き出しに書いてあるようなもんじゃん」
東大卒・高橋まつりさまの至言です。

5:卵の名無しさん
18/05/20 16:56:36.42 n2fbjQMc.net
p<0.05で有意差ありと判断した試験の再現性のなさで
有意水準を0.005あるいは0.001にすべきとか、p値を判断に用いるべきでないという議論が盛んなようである。
p値が0.05をわずかに切った程度ではFalse Positive Risk(棄却された帰無仮説が正しい割合が高い。
FPRは帰無仮説が正しいときの棄却率(危険率)ではない!
URLリンク(arxiv.org)
にはt検定での平均値の比較が論じられている。
同じロジックで比率の比較をするスクリプトを組んでみた。
p=0.0434のEMPA-REG OUTCOMEのデータ
URLリンク(square.umin.ac.jp)を用いてシミュレーションしてみた。
帰無仮説が正しい比較対照としてプラセボ通しの比較でやってみた。
(理論上はp-値が一様分布するはずなのでこれはシミュレーションしなくてもいいが、
スクリプトの動作確認も含めて試行。)
URLリンク(i.imgur.com)
p<0.05でのFPR
> s0/(s0+s1)
[1] 0.07757254
0.040< p < 0.045でのFPR
> S0/(S0+S1)
[1] 0.1909091
確率密度関数で線分比(FPR.equal)および面積比(FPR.less)とほぼ一致した。
> calc.FPR.p2(rp,rb,np,nb,Print=FALSE)
power p.value FPR.equal LH10 FPR.less
0.53163451 0.04342913 0.20648690 3.84292218 0.07877384
スクリプトに大きなバグはないと確信できて、気持ちが(・∀・)イイ!!
理論値とシミュレーションの一致はまるで
ド底辺シリツ医大卒と裏口バカの一致をみるようだ。

6:卵の名無しさん
18/05/20 17:17:59.33 n2fbjQMc.net
心血管イベントの発生がポワソン分布にしたがうモデルでシミュレーションしてみた。

> s0/(s0+s1)
[1] 0.07810437
> S0/(S0+S1)
[1] 0.1952311

> calc.FPR.p2(rp,rb,np,nb,Print=FALSE)
power p.value FPR.equal LH10 FPR.less
0.53163451 0.04342913 0.20648690 3.84292218 0.07877384
ほぼ、同様の値になった。
英文で学力を試そうが、中学数学で学力を試そうが
ド底辺シリツ医大卒は低学力の裏口バカであることが示されるのに似ている。

7:卵の名無しさん
18/05/20 18:19:26.39 n2fbjQMc.net
これくらい差がつくと
Apixaban for Extended Treatment of Venous Thromboembolism
URLリンク(www.nejm.org)
FPR(False Positive Risk)がp値の20倍大きいけど全く問題ないな。
一応、計算してみると
r1=14
n1=840
r2=73 # placebo
n2=829 # placebo
> calc.FPR.p2(r1,r2,n1,n2)
power p.value FPR.equal LH10 FPR.less
9.999998e-01 1.116811e-10 2.171983e-09 4.604087e+08 1.561171e-10

8:卵の名無しさん
18/05/20 21:18:47.27 ZIRLy3z0.net
中国の「シェア看護師」、高齢化社会で不可欠 不安要素も
URLリンク(headlines.yahoo.co.jp)

9:卵の名無しさん
18/05/21 00:12:50.70 r/Pmsrg8.net
p値にreproducibilityがないが、あるp値が算出された試験をサンプル数を維持して再現したときに得られるp値の95%信頼区間を算出するという触れ込みで
紹介されている論文: 
P-values in genomics: Apparent precision masks high uncertainty
URLリンク(www.ncbi.nlm.nih.gov)

数理展開が納得できない(理解できない)のでシミュレーションしてみた。
URLリンク(i.imgur.com)
95%の範囲には収まらないでもっと拡散する。
high uncertaintyどころかcompletely unpredictableだと思う。

10:卵の名無しさん
18/05/21 00:13:47.86 r/Pmsrg8.net
>>9
シミュレーションのコードはここに置いた。
スレリンク(hosp板:503番)

11:卵の名無しさん
18/05/21 00:34:25.67 81lS21Ih.net
底辺私立医大を卒業した医者って頭悪いよね?
って、、、、、
何を今更、、、、、

12:卵の名無しさん
18/05/21 06:59:08.06 X6Stfedi.net
医師、患者の共通認識を確認するスレね。
名投稿を再度、引用
(quote)
私は昭和の時代に大学受験したけど、昔は今よりも差別感が凄く、慶応以外の私立医は特殊民のための特殊学校というイメージで開業医のバカ息子以外は誰も受験しようとすらしなかった。
常識的に考えて、数千万という法外な金を払って、しかも同業者からも患者からもバカだの裏口だのと散々罵られるのをわかって好き好んで私立医に行く同級生は一人もいませんでした。
本人には面と向かっては言わないけれど、俺くらいの年代の人間は、おそらくは8−9割は私立卒を今でも「何偉そうなこと抜かしてるんだ、この裏口バカが」と心の底で軽蔑し、嘲笑しているよ。当の本人には面と向かっては絶対にそんなことは言わないけどね。
(unquote)
それにはこういうレスが続いた。
(quote)
それは昭和の時代だけでなく平成の今も同じです。
みんな口に出さないだけで思ってることは同じです。
そういう私も皆の前では本心とは裏腹に「国家試験に受かれば
私立も国立も変わらないよ。」と言ってます。
在日や帰化人を連想させる名前をした医師に「******(伏字は引用者)」と内心では思ってても、
ニコニコしながら話を合わせてるのと同じです。
(unquote)
>でも一


13:カ、卒業大学を隠すのは本当に辛い       http://imagizer.imageshack.com/img922/6510/c71CG3.jpg は 昭和の時代でなく平成の今も、そして次の御世でも同じであることを示していると言えましょう。



14:卵の名無しさん
18/05/21 07:01:13.62 X6Stfedi.net
頭が悪いのを自覚していない医者っているよなぁ。
ICU Bookの最終章の冒頭で著者がこう書いている。
In clinical matters, ignorance can be dangerous,
but ignorance of ignorance can be fatal.
経験上、これが当てはまるのがシリツ医大卒。

15:卵の名無しさん
18/05/21 07:33:00.62 X6Stfedi.net
Record of a DoTeihen lowerest tiered school graduate incompetent to solve linear equation
URLリンク(imagizer.imageshack.com)
We are taken aback by the daunting stupidity of the DoTeihen that he is not able to calculate the male/female proportion when each ratio and overall ratio are given
How come could such a retard go to the DoTeihen medical school?
Is there any better accounting for it than his bribing his way into the DoTeihen medical school, in other word Uraguchi?
Mr. Kazuki Enari, who insists DoTeihen medical school accommodates such feebs who cannot solve quadratic equations, would be flabbergasted .
It is as easy as anything even for junior high school students to solve the problem with a linear equation.
That Doteihen dumb contends unwisely that
male/female proportion is necessary to figure out.
URLリンク(imagizer.imageshack.com)
After I posted the answer of its probability 0.2575549,
he responded that the answer needed unit.
There seems to be some unit to denote probability in DoTeihen medical school , LOL
How come could such a retard go to the DoTeihen medical school?
Is there any better accounting for it than his bribing his way into the DoTeihen medical school, in other word Uraguchi?
URLリンク(imagizer.imageshack.com)

16:卵の名無しさん
18/05/21 09:26:33.03 J8X+8OTV.net
帝京大医学部の女の子がテラスハウス出てたけど、
その肩書きがテロップに出てる限り、「私はめちゃめちゃバカだけど、
親に大金払ってもらって医者になるよ!」って
吹き出しに書いてあるようなもんじゃん」
東大卒・高橋まつりさまの至言です。

17:卵の名無しさん
18/05/21 10:28:18.46 DAkaDn2l.net
よしよしもっとやれもっとやれ。他人がどんどん落ちていくのを見るのは楽しいもんだな。世の中厳しいからなw

18:卵の名無しさん
18/05/21 11:39:04.86 BarEo4sL.net
>>16
ド底辺シリツ医大卒より恥ずかしい資格に墜ちることはない。
Last but not least, three laws of Do-Teihen(lowest-tier) Medical School, currently called Gachi'Ura by its graduates.
最後にド底辺医大の三法則を掲げましょう。
1: It is not the lowest-tier medical school but its enrollee that is despicable, which deserves to be called a bona fide moron beyond redemption.
ド底辺シリツ医大が悪いのではない、本人の頭が悪いんだ。
2: The graduates of Do-Teihen(lowest-tier medicalschool) are so ashamed that none of them dare to mention their own alma mater they had gone through.
ド底辺シリツ医大卒は恥ずかしくて、学校名を皆さま言いません。
3: The Do-Teihen(lowest-tier medicalschool) graduates are so ashamed of having bought their way into the exclusively bottom-leveled medical school that they tend to call a genuine doctor a charlatan who elucidates their imbecility.
ド底辺特殊シリツ医大卒は裏口入学の負い目から裏口馬鹿を暴く人間を偽医者扱いしたがる。

19:卵の名無しさん
18/05/21 11:39:46.68 BarEo4sL.net
>>16
さてー、第二法則を確認しましょう。
あんたどこ卒?

20:卵の名無しさん
18/05/21 12:43:48.51 U4ocjQOb.net
>>16
確かに、
ド底辺シリツ医大卒という烙印より墜ちることはないから
「楽しい」と言うのは一理あるなwww
患者に「失礼ですがどこの大学卒業ですか?」と聞かれて医科歯科卒と答えたら、
「よかった、自分より馬鹿に診てもらいたくなかったので」と言われたな。
懇意なw看護師から「私達も今の先生はどこ卒ですかとよく聞かれるんですよと」言われた。
シリツ医大卒の時は「私達はそういう個人情報は教わってないんです」といって逃げるそうです。
国立卒の医師の時は正しく教えて患者が安心する顔をみてほっとする、と言ってたな。

21:卵の名無しさん
18/05/21 12:45:10.35 U4ocjQOb.net
一般的には落ちた人間と評価されるであろうナマポからも蔑まれるのがド底辺シリツ医大卒。

517 卵の名無しさん 2018/02/25(日) 11:36:00.56 ID:gq76tAvs
福岡のあの歯科大かな?
歯科口腔外科で抜歯依頼したら爺ちゃんが〇〇歯科大卒の先生は避けて下さいねがあった
けど。
この爺ちゃん、聖マリ卒の先生もよけて下さいと初診時に言ってた札付き爺さん。
生保受給者のくせにね。

22:卵の名無しさん
18/05/21 12:46:15.48 U4ocjQOb.net
おーい、ド底辺!
前々スレの宿題が残っているぞ!!

ド底辺シリツ医はpulseless VTにcardioversionしないで何するんだい?
私立医三年生の正答率ほぼ100%なんだろ?

ネットにド底辺シリツ医大卒の馬鹿の証を保存しておきました。
URLリンク(imagizer.imageshack.com)

pulseless VTにド底辺シリツ医大はcardioversionせずに何をするんだ?
PEAをpulseless VTと混同するような馬鹿は医師ではないだろう。
いや、ド底辺の裏口バカならありうるかもしれん。
だって、一次方程式すらまともに立てれない馬鹿ですもん。

23:卵の名無しさん
18/05/21 12:47:15.58 U4ocjQOb.net
これには、ド底辺シリツ医大卒の裏口バカはダンマリだなぁ
同窓会名簿は金払っているから捨ててないよ。

URLリンク(imagizer.imageshack.com)
このスレを表示したタブレット画面だから俺の手元にあるのが納得できた?
しかし医科歯科卒程度の学歴が羨ましいかねぇ?
医科歯科蹴る合格者もいるし。

24:卵の名無しさん
18/05/21 14:14:24.16 HybFwdZX.net
開業医スレにこんな投稿が
>>
【ウハも】 開業医達の集い 15診 【粒も】
スレリンク(hosp板:941番)
941 名前:卵の名無しさん[] 投稿日:2018/05/21(月) 13:55:59.02 ID:CoezPpXl
接客業だ。いかに患者様に安心していただくか。患者様の笑顔のために栗は心を尽くします
<<
ナマポからも忌避される、ド底辺シリツ医大卒ってどうすりゃいんだろね?
>>
517 卵の名無しさん 2018/02/25(日) 11:36:00.56 ID:gq76tAvs
福岡のあの歯科大かな?
歯科口腔外科で抜歯依頼したら爺ちゃんが〇〇歯科大卒の先生は避けて下さいねがあった
けど。
この爺ちゃん、聖マリ卒の先生もよけて下さいと初診時に言ってた札付き爺さん。
生保受給者のくせにね。
<<
ナマポにも安心感を与えられない「社会的な」資格が、ド底辺シリツ医大卒ってことだな。
なんでだろう?
ド底辺シリツ医大卒=裏口バカという公式が上級国民から生活保護受給者までの共通認識だからではないだろうか?

25:卵の名無しさん
18/05/22 01:24:04.14 pgTq+n72.net
Intuitive Biostatisticsの最新版で加筆された章
Statistics and Reproducibilityは実に面白い。
この本は数式は殆ど皆無で結論が数値や図表として書かれてとその根拠論文が掲載してあるだけ
算出方法を知りたければ原著にあた


26:って自分で勉強することになる。 P values are not very reproducible.を例示した表が整理引用されていた。 全文がfree accessであったので Explorations in statistics: statistical facets of reproducibility https://www.physiology.org/doi/pdf/10.1152/advan.00042.2016 に目を通した。算出根拠となる理論さえわかればあとはプログラムを組めばいい。 みつかった。 If we assume the magnitude of the effect observed in our initial experiment equals the magnitude of the true effect, then the probability that the P value from a second experiment, P2, will be less than alpha(=0.05) is Pr{z > z* - z1} = Pr{P2<alpha}. 但し、理論説明のためか母集団の分散が既知であるのを前提にしているz検定。 母集団の平均値に差があるかを知りたいのに母集団の平均値がわからずに分散がわかっているというケースはあまりないと思う。このスレ的に言えばこういう例示がわかりやすい。 ド底辺シリツ医大の入学者の偏差値は60未満であると容易に推測できるが、分散の推測は困難である。ガチの裏から特待生までいるが裏口バカも多数派だろうから。 z検定とともにt検定で算出できるようにした。



27:卵の名無しさん
18/05/22 01:24:35.65 pgTq+n72.net
Even if the P value from the first experiment is 0.001, you can be far from sure that the repeat experiment will have a P value less than 0.05.
It turns out that that probability is only 78% or 91%(the two values differ because two methods were used to compute.)
前者のGoodmanの論文はabstractしか入手できず、後者の上記論文でプログラムした。
このスレの趣旨の問題にするとこんな感じ。
帰無仮説:ド底辺シリツ医大の合格者と入学者の偏差値の平均値は等しいが、
無作為抽出比較でp=0.01で棄却された。
この実験での平均値の差が母集団の平均値の差を表していると仮定して同じ実験を繰り返したときに得られるp値が0.05未満になる確率はいくらか?
答. t検定だと77%、(z検定だと73%)
URLリンク(i.imgur.com)
URLリンク(i.imgur.com)
# calculate probability of repeat experiment is significant
p2p2 <- function(p.value,n1=10,n2=10,sd1=1,sd2=1,alpha=0.05){
p=p.value/2 # two.sided comparison
z=qnorm(p)
d.z=z*sqrt(sd1^2/n1+sd2^2/n2) # difference of means by z.test
p2.z=pnorm(-qnorm(alpha/2)+z,lower=FALSE)
df=n1-1+n2-1
t=qt(p,df)
d.t=t*sqrt((1/n1+1/n2)*((n1-1)*sd1^2+(n2-1)*sd2^2)/((n1-1)+(n2-1)))
p2.t=pt(-qt(alpha/2,df)+t,df,lower=FALSE)
data.frame(p2.t,p2.z,d.t,d.z)
}

28:卵の名無しさん
18/05/22 01:27:24.16 pgTq+n72.net
再検でも有意水準0.05で有意になる確率が95%以上あるためには最初のt検定で必要なp値は
> uniroot(function(x)p2p2(x)$p2.t - 0.95,c(0,1))$root
[1] 0.001213122
と算出された。
z検定だと
> uniroot(function(x)p2p2(x)$p2.z - 0.95,c(0,1))$root
[1] 0.0003217581
これはテキストの数値と一致して気持ちが( ・∀・)イイ!!

29:卵の名無しさん
18/05/22 03:08:44.00 p6I2zHsz.net
>>3
最底辺聖マがぬけてる

30:卵の名無しさん
18/05/22 06:37:09.14 pgTq+n72.net
>>27
そこはナマポからも蔑まれてているよな。
517 卵の名無しさん 2018/02/25(日) 11:36:00.56 ID:gq76tAvs
福岡のあの歯科大かな?
歯科口腔外科で抜歯依頼したら爺ちゃんが〇〇歯科大卒の先生は避けて下さいねがあった
けど。
この爺ちゃん、聖マリ卒の先生もよけて下さいと初診時に言ってた札付き爺さん。
生保受給者のくせにね。

31:卵の名無しさん
18/05/22 06:51:45.14 pgTq+n72.net
合格者の偏差値と入学者の偏差値はド底辺シリツ医大では乖離があるとされているが、
一応、この偏差値が学力を反映する。
データは各シリツ医大の合格者から無作為抽出した値であるとして
都内シリツと都外シリツで有意差があるか検定すると
> t.test(intelligence~tokyo,dat)
Welch Two Sample t-test
data: intelligence by tokyo
t = -3.4885, df = 37.567, p-value = 0.001256
以上から、このスレに投稿している裏口バカは
都外のド底辺シリツ医大卒と推測される。

32:卵の名無しさん
18/05/22 08:54:31.17 C2xlYb8n.net
帝京大医学部の女の子がテラスハウス出てたけど、
その肩書きがテロップに出てる限り、「私はめちゃめちゃバカだけど、
親に大金払ってもらって医者になるよ!」って
吹き出しに書いてあるようなもんじゃん」
東大卒・高橋まつりさまの至言です。

33:卵の名無しさん
18/05/22 11:02:40.90 pgTq+n72.net
このスレに投稿している裏口バカが都内シリツかを検討するのに慶應を外すのは最初から慶應卒でないことが前提になっているからバイアスといえる。
全データを使って解析をやり直す。
URLリンク(2chreport.net)
慶應義塾学医学部 2051 70
順天堂大学医学部 2090 66
昭和大学医学部 2200 65
東京慈恵会医科大学医学部 2250 69
自治医科大学医学部 2260 66
産業医科大学医学部 2749 65
日本医科大学医学部 2813 66
東京医科大学医学部 2995 64
関西医科大学医学部 3014 65
大阪医科大学医学部 3141 66
東邦大学医学部 3180 62
久留米大学医学部 3237 61
東京女子医科大学医学部 3284 60
日本大学医学部 3310 61
岩手医科大学医学部 3400 60
聖マリアンナ医科大学医学部 3440 58
近畿大学医学部 3580 62
藤田保健衛生大学医学部 3650 60
獨協医科大学医学部 3660 59
帝京大学医学部 3750 59
杏林大学医学部 3755 61
東海大学医学部 3760 61
福岡大学医学部 3770 60
愛知医科大学医学部 3800 59
埼玉医科大学医学部 3800 57
兵庫医科大学医学部 3880 60
北里大学医学部 3890 60
金沢医科大学医学部 3950 58
川崎医科大学医学部 4565 57

34:卵の名無しさん
18/05/22 11:06:05.43 pgTq+n72.net
p値は0.0268で頻度主義統計では有意差(両側検定)ありと判断された。
> t.test(intelligence~tokyo,dat)
Welch Two Sample t-test
data: intelligence by tokyo
t = 2.4141, df = 17.779, p-value = 0.02679、
alternative hypothesis: true difference in means is not equal to 0
95 percent confidence interval:
0.4037647 5.8588616
sample estimates:
mean in group 1 mean in group 2
63.90909 60.77778

35:卵の名無しさん
18/05/22 11:10:30.94 pgTq+n72.net
上記の偏差値がある試験でのサンプリングによるとして、それは母集団の平均や分散を反映しているとの仮定をおくと
別の試験で調査したときのp<0.05となる確率はいくらになるか?
等分散を仮定しないWelch法でのt検定からでも算出できるように
スクリプトを変更。
p2p2w <- function(p.value,n1=10,n2=10,sd1=1,sd2=1,alpha=0.05){
p=p.value/2 # two.sided comparison
# Z-test
z=qnorm(p)
d.z=z*sqrt(sd1^2/n1+sd2^2/n2) # estimated difference of means by z.test
p2.z=pnorm(-qnorm(alpha/2)+z,lower=FALSE)
# Student
df=n1-1+n2-1
t=qt(p,df)
d.t=t*sqrt((1/n1+1/n2)*((n1-1)*sd1^2+(n2-1)*sd2^2)/((n1-1)+(n2-1)))
p2.t=pt(-qt(alpha/2,df)+t,df,lower=FALSE)
# Welch
var1=sd1^2 ; var2=sd2^2
dfw=(var1/n1+var2/n2)^2 / (var1^2/n1^2/(n1-1)+var2^2/n2^2/(n2-1))
t.w=qt(p,dfw)
d.w=t.w*sqrt(var1/n1+var2/n2)
p2.w=pt(-qt(alpha/2,dfw)+t.w,dfw,lower=FALSE)
data.frame(p2.w, p2.t, p2.z)
}

36:卵の名無しさん
18/05/22 11:52:57.04 pgTq+n72.net
>>31
こんな計算も面白い。
このデータから偏差値1点をいくらで買うことになるか、95%信頼区間とともに述べよ。
線形回帰の係数で求めればいいので
1点163万円(95%信頼区間は130〜196万円)になる。
これには「任意の」寄付や学債は含まない。

37:卵の名無しさん
18/05/22 15:25:20.89 pgTq+n72.net
>>32
この試験の再現性を計算すると
> (n1=sum(dat$group==1)) ; (n2=sum(dat$group==2))
[1] 11
[1] 18
> (p.siri=t.test(score~group,dat)$p.value)
[1] 0.0267907
> (var=tapply(dat$score,dat$group,var))
1 2
13.290909 8.535948
> (sd1=sqrt(var[1])) ; (sd2=sqrt(var[2]))
1
3.64567
2
2.921634
> p2p2w(p.siri,n1,n2,sd1,sd2)
p2.w p2.t p2.z d.w d.t d.z
1 0.6203964 0.6131503 0.6004801 -3.131313 -2.876468 -2.872511
再現(次回もp<0.05になる確率)確率は0.62と算出された。
どうやって検証するか?
ブートストラップ法でやってみる。

38:卵の名無しさん
18/05/22 15:36:53.93 pgTq+n72.net
library(boot)
data2p2.w <- function(data,i){
dat=data[i,]
(n1=sum(dat$group==1)) ; (n2=sum(dat$group==2))
(p.siri=t.test(score~group,dat)$p.value)
(var=tapply(dat$score,dat$group,var))
(sd1=sqrt(var[1])) ; (sd2=sqrt(var[2]))
c(p=p.siri,p2=p2p2w(p.siri,n1,n2,sd1,sd2)$p2.w)
}
data2p2.w(dat,1:29)
re.boot=boot(dat,function(x,i)data2p2.w(x,i),R=10000)
re.boot$t0
> mean(re.boot$t[,1]<0.05)
[1] 0.6284
5000回のブートストラップでt検定するとp値が0.05未満になるの割合は0.628で上記の計算とほぼ一致した。
やはり、都内私立医大と都外私立医大の合格者の偏差値の平均値は等しくないことが確認できた。
このスレに投稿しているド底辺シリツ医大は都外のド底辺シリツ医大卒と推測される。
聖魔か不治多ではないかと思う。
皮裂は学生の不祥事は多いが医師の不祥事は聞かないように思う。

39:卵の名無しさん
18/05/22 15:47:51.89 pgTq+n72.net
ブートストラップでのp値をグラフにしてみた。
URLリンク(i.imgur.com)
単発データから重複サンプリングして分布を推測するブートストラップって便利だな。

インフルエンザ新薬ゾフルーザの解析でも使われていたなぁ。
対称な分布にはならないだろうから信頼区間はパーセンタイル法でなくて
Highest Density Intervalにすればいいのにと思ったな。

40:卵の名無しさん
18/05/22 19:08:19.80 pgTq+n72.net
 気仙沼市は、市立病院で受診し、同病院に入院中に死亡した30代男性の遺族に損害賠償金5000万円を支払って和解することを決めた。
 同病院によると、男性は2014年3月、喉の強い痛みなどを訴えて受診した。当直医の診察を受け入院したが、翌日に容体が急変し死亡した。
遺族は15年10月、医師と看護師の連携など救急医療体制の不備で病院側に過失があったとして、約9330万円の賠償を求め、仙台地裁に提訴。今年2月に地裁から和解勧告を受けた。
URLリンク(mainichi.jp)
急性喉頭蓋炎かな?

41:卵の名無しさん
18/05/22 21:22:34.79 pgTq+n72.net
>>25
Goodman(1992) URLリンク(www.ncbi.nlm.nih.gov) の78%と
Curran-Everett(2016) URLリンク(www.ncbi.nlm.nih.gov) の91%では
計算法の違いで説明するには差が大きすぎるように感じる。
後者のみFree Accessだったので>33のスクリプトにできた。
Even if the P value from the first experiment is 0.001,のデータを作ってブートストラップでシミュレーションしてみることにした。
母集団の平均=0の帰無仮説がp≒0.001で棄却される標本を乱数発生させて作成
> dat
[1] -1.03721684 0.08464184


42:-0.76377658 -1.87636634 0.10920118 0.25192097 [7] -0.96821719 -1.21244950 -1.49732560 -0.14440373 -1.62905565 -0.24083969 [13] -0.98561242 0.47706931 -1.69032737 -1.41328417 > t.test(dat) One Sample t-test data: dat t = -4.0741, df = 15, p-value = 0.0009973 alternative hypothesis: true mean is not equal to 0 95 percent confidence interval: -1.1934074 -0.3735978 sample estimates: mean of x -0.7835026 0.0009973なので、これを使う。



43:卵の名無しさん
18/05/22 21:32:26.29 pgTq+n72.net
青〇がブートストラップでのシミュレーション、黒線がt分布の累積密度関数からの作成した曲線。
URLリンク(i.imgur.com)
0.78より、0.91の方が近


44:似している 但し、シミュレーションでは > mean(p<0.05) [1] 0.9796 98%という結果になった。



45:卵の名無しさん
18/05/23 11:11:11.55 dY2QXaHx.net
These are the very benificial laws for not only doctors but also everyone who has common sense,
Three laws of Do-Teihen lowest tier Medical School, currently called Gachi'Ura by its graduates.
1: It is not bottom medical school but its enrollee that is despicable, which deserves to be called a bona fide moron beyond redemption.
2: The graduates of Do-Teihen are so ashamed that none of them dare to mention their own alma mater they have gone through.
3: The Do-Teihen graduates are so ashamed of having bought their way into the privately-founded lowest tier medical school that they tend to call a genuine doctor a charlatan who elucidates their imbecility.
Glossary for better understanding
bona fide : neither specious nor counterfeit : genuine
despicable : very bad or unpleasant : deserving to be despised
matriculation : the formal process of entering a university, or of becoming eligible to enter by fulfilling certain academic requirements such as a matriculation examination.
alma mater : a school, college, or university which one has attended or from which one has graduated
charlatan : someone who professes knowledge or expertise, esp in medicine, that he or she does not have; quack.

46:卵の名無しさん
18/05/23 11:11:39.60 dY2QXaHx.net
>>40
シミュレーションと理論グラフが乖離している。
しかも理論グラフの上下にばらつくのでなくて一定の傾向(シミュレーション>理論)をもっているので乖離には理由があるはず。
理論曲線は正規分布で計算しているのに
シミュレーションに使った値はt検定でp=0.001になった数値だからだろう。
母集団の分散既知のz.testでp=0.001となる数値でシミュレーションすればいいんだろうな。

47:卵の名無しさん
18/05/23 12:24:25.43 k6Vtmkg0.net
帝京大医学部の女の子がテラスハウス出てたけど、
その肩書きがテロップに出てる限り、「私はめちゃめちゃバカだけど、
親に大金払ってもらって医者になるよ!」って
吹き出しに書いてあるようなもんじゃん」
東大卒・高橋まつりさまの至言です。

48:卵の名無しさん
18/05/23 12:42:04.64 fmIpyRoK.net
分散1の母集団から16個サンプリングすることにして
母集団の平均値=0を帰無仮説とした検定でp=0.001となるようなサンプルを見出す。
n=16
sigma=1
z.test <- function(x){
z=sqrt(n)*mean(x)/sigma
2*pnorm(abs(z),lower=FALSE)
}
> (dat=rnorm(n))
[1] 0.1152750 0.1132882 0.7574643 1.1997442 -0.4105303 0.5980804
[7] 1.7610516 -0.7704866 0.3171577 1.0221521 1.3038936 1.9375013
[13] 0.5715436 1.2104586 1.2727820 2.1627381
> z.test(dat)
[1] 0.0009999938
ほぼp値0.001のサンプルがみいだせた。
これを使ってブートストラップ法でシミュレーションする。
URLリンク(i.imgur.com)
シミュレーションと理論曲線が一致して気持ちが( ・∀・)イイ!!
まるで、ド底辺シリツ医大卒と裏口バカの一致のようではないか(笑)

49:卵の名無しさん
18/05/23 18:08:51.15 fmIpyRoK.net
>>44
16個のサンプルを64個にしてt検定でp=0.001となるデータを作成して>40のグラフを再作成してみた。
URLリンク(i.imgur.com)
近似してきた。
自由度無限大のt分布が正規分布だから当然といえば当然。

50:卵の名無しさん
18/05/23 18:33:55.04 fmIpyRoK.net
Even if the P value from the first experiment is 0.001, you can be far from sure that the repeat experiment will have a P value less than 0.05.
It turns out that that probability is only 78% or 91%(the two values differ because two methods were used to compute.)
の結論はシミュレーションでは後者が正しそうであるのでが
Explorations in statistics: statistical facets of reproducibility
URLリンク(www.physiology.org)
の例示では平均差で説明してあるが、比率の差でも成立するだろうか?
サンプル数が増えれば正規分布で近似でp値を判断できるのでたぶん成立するんだろうな。たぶんだと気持ちが悪いのでシミュレーションしたい。
まず、p値=0.001になるデータ作成から始める。
EMPA-REG試験の数字を種に総当たりで作成。
> M.seed=matrix(m[954,],ncol=2) ; M.seed
[,1] [,2]
[1,] 295 2038
[2,] 471 4231
> chisq.test(M.seed)
Pearson's Chi-squared test with Yates' continuity correction
data: M.seed
X-squared = 10.827, df = 1, p-value = 0.001
周辺度


51:数を固定して二項分布でイベント発生するとしてこれからMCMCでシミュレーションデータ群をつくってそのデータでのχ2乗検定で何%がp<0.05になるか計算すればいい。



52:卵の名無しさん
18/05/23 19:05:28.77 fmIpyRoK.net
>>46
シミュレーションでのp値の分布はこんなかんじ。
URLリンク(i.imgur.com)
p<0.05で有意差が再現される確率は90.1%となり、理論値91%に近似した結果となった。

53:卵の名無しさん
18/05/23 22:10:40.90 fmIpyRoK.net
サパイラに新版がでてるな。
URLリンク(www.amazon.co.jp)
kappa値の計算方法が丁寧に解説してあったのを思い出した。
10-finger testは感度が高いという記述も面白いと思ったが、その後に読んだ
Modern EpidemiologyにRothmanが
頭があるということは髄膜炎の診断には感度100%である
ということを書いていてこっちの方が極端で感銘したなぁ。
ド底辺シリツ医大卒にもわかる説明だと思うが、やつらは 「馬鹿は、死ななきゃ治らない」の対偶すら即答できない底抜け馬鹿なのでどうだろうね。

54:卵の名無しさん
18/05/24 10:01:59.03 o31IrGHvu
すごい。完全に病気だ。おぞましいコンプレックス。

55:卵の名無しさん
18/05/24 13:11:33.62 a1LNB/B6.net
バカ私立は使えない上に性格最悪
底辺駅弁にも似たようなのがいる
ひねくれ根性がハンパない
クビにしないと他の医者が迷惑この上ない

56:卵の名無しさん
18/05/24 13:50:56.92 R1FtDNpg.net
性格の悪いのはどんな組織にもいるだろうが、
頭が悪い奴が進学するのがド底辺シリツ医大であることは申すまでもない。

57:卵の名無しさん
18/05/24 14:38:15.29 R1FtDNpg.net
一次方程式で解ける比率計算すらできない馬鹿はどこにでもいるが、
ド底辺シリツ医大卒の裏口は自分が馬鹿だという自覚がないから、
与えられた数値から計算できない、問題が悪いと言い張るんだよなぁ。
悪いのはそいつの頭なのに。
これが現実
1次方程式もできないド底辺特殊シリツ医大卒の記録
URLリンク(imagizer.imageshack.com)

58:卵の名無しさん
18/05/24 18:21:18.09 QkjtOtjO.net
確かに馬鹿の自覚がない。
寧ろ国立よりどーでもいい枝葉末節説ひけらかし、
威張って堂々としてる馬鹿ボッチャンが多い。
はっきり言うが、命を預かるこの仕事に
底辺私立を就かせてはならない。
これまで、そしてこれから、
どれだけの患者を奴らに殺させるのか。

59:卵の名無しさん
18/05/24 21:01:35.04 R1FtDNpg.net
COMPASS試験のデータで遊んでみる。
URLリンク(www.nejm.org)
Table.2から、脳出血の発生率は
リバロキサバン+アスピリン(r1/n1)、リバロキサバン単独(r2/n2)、アスピリン単独(r3/n3)
# R+A ; R alone ; A alone
r1=15 ; r2=27 ; r3=10
n1=9152 ; n2=9117 ; n3=9126
これは多重比較なので補正が必要なような気がしないでもないのでbonrerroniの補正でやってみる
Pairwise comparisons using Pairwise comparison of proportions
data: c(`Rivaloxban+Aspirin` = r1, `Rivaroxavan alone` = r2, `Aspirin alone` = r3) out of c(n1, n2, n3)
Rivaroxban+Aspirin Rivaroxavan alone
Rivaroxavan alone 0.261 -
Aspirin alone 1.000 0.025
P value adjustment method: bonferroni
結論は変わらない。Rivaloxban単独はアスピリン単独より脳出血が多いという結論は変わらず。
補正なしだと
Rivaloxban+Aspirin Rivaroxavan alone
Rivaroxavan alone 0.0869 -
Aspirin alone 0.4275 0.0084
P value adjustment method: no


60:ne 脳出血の発生がポワソン分布に従うとしてこの試験を繰り返したときに アスピリン単独とRivaroxbanで有意差がでる確率をシミュレーションしてみた。 Pr(p.value < 0.05 ) = 0.7468 Original p.value = 0.008387501 になった。多重比較補正をしてもPr(p.value < 0.01666667 ) = 0.5754 の確率で で脳出血が多いという結果が再現される。



61:卵の名無しさん
18/05/24 21:24:22.16 R1FtDNpg.net
stanを含むRのコードは スレリンク(hosp板:505番)-506
論文では計算されていないが、
リバロキサバン+アスピリン併用群とリバロキサバン単独群の比較をやってみる(原著と同じく多重比較補正なしとする)
Original p.value = 0.08694743で有意差なし。
有意差があれば、アスピリン併用はリバロキサバンによる脳出血を抑制するという解釈でいいのだろうかwww
面白いのでこの結論がp<0.05で再現される確率を、脳出血の発生数はポアソン分布に従うとしてシミュレーションすると
Pr(p.value < 0.05 ) = 0.3868
という結果になった。
p値が0.05を少し下回るくらいで、統計的有意差があったなんて結論されているときには
シミュレーションで再現確率を出すと( ・∀・)イイ!!

62:卵の名無しさん
18/05/24 21:37:10.91 R1FtDNpg.net
>>53
性別や年齢で入院日数や手術成績に相関があるかという論文がでてたよなぁ。
ド底辺シリツ医大卒か否かでわけて相関がだせたら面白いのになぁと思う。
大学の偏差値を傾向スコアにして解析できそうだが、データ集積の方が困難かもな。

63:卵の名無しさん
18/05/25 06:54:26.07 8Ot6XASs.net
叱られないと勉強しない子供が勉強しているのをみたら
叱られたから勉強しているんだな、と考えるのが一般常識人の論理だよね。
ところが、底辺私立医大出身者は
勉強すると叱られる
と主張するんだよ。
しかも対偶をとると意味が逆になる例、とか底抜け馬鹿だぞ。
スレリンク(hosp板:499番)
URLリンク(i.imgur.com)

底辺シリツ医大出身者には通読できないICU Bookの最終章の冒頭で著者がこう書いている。
In clinical matters, ignorance can be dangerous,
but ignorance of ignorance can be fatal.
「叱られないと勉強しない」の対偶を「勉強すると叱られる」
と答えるのはignorance can be dangerousの範疇だが、
ドヤ顔で
>対偶をとれば意味が逆になる例文。
というのは、まさに
ignorance of ignorance can be fatal.

64:卵の名無しさん
18/05/25 07:44:03.59 8Ot6XASs.net
問題:
命題とその対偶は真偽が一致する。
命題「馬鹿でないと裏口医大には入学しない」
の対偶を述べよ。
国立大学卒の答:
「裏口医大に入学するのは馬鹿である」
裏口医大出身の答:
「裏口医大に入学すると馬鹿になる」
URLリンク(i.imgur.com)
last but not least
底辺私立医大が悪いんじゃない
本人の頭が悪いんだ。

65:卵の名無しさん
18/05/25 10:32:09.02 8Ot6XASs.net
私立医大と裏口の相関のデータにグーグル検索を用いてみる。
検索単語組み合わせと ヒットした件数は以下の通りである。
医学部 私立 裏口 29,700件
医学部 国立 裏口 25,500件
医学部 私立  2,110,000件
医学部 国立  5,580,000件
ド底辺私立医 裏口33,200 件
ド底辺私立医  208,000 件
ド底辺私立医大と裏口の相関関係は他と比べて有意に高いといえるか?
多重検定になるのでまず、familywise のp値を求めてみる。
χ2乗検定で両側検定すると、
3-sample test for equality of proportions without continuity correction
data: c(r1, r2, r3) out of c(n1, n2, n3)
X-squared = 437710, df = 2, p-value < 2.2e-16
小さすぎてRでは表示してくれないのでMathematicaにお願いすると 1-CDF[ChiSquareDistribution[2],437710]
p値は 3.0280501516×10^-95048
で偶然にこのような分布になる確率は0といってよい。 
では、どの二群間で有意差があるかやってみよう。

66:卵の名無しさん
18/05/25 10:36:49.56 8Ot6XASs.net
Pairwise comparisons using Pairwise comparison of proportions
data: c(私立 = r1, 国立 = r2, 底辺 = r3) out of c(n1, n2, n3)
私立 国立
国立 <2e-16 -
底辺 <2e-16 <2e-16
P value adjustment method: bonferroni
シリツと裏口の相関がありすぎて、bonferroni補正してもp値が小さすぎ!!

67:卵の名無しさん
18/05/25 10:59:34.16 8Ot6XASs.net
上記は、サンプルサイズを大きくすればp値はいくらでも小さくできるという実例。
こういう手法は、10mg投与と25mg投与を合算してサンプルサイズを増やして有意差を捏造したのがEMPA-REG OUTCOME試験の解析でも使われている。
n増加でのp-hackingは面白くないので裏口との関係リスク比を底辺私立医大と私立医大で比べて、信頼区間もだしてみる。
Incidence rate ratio estimate and its significance probability
data: r3 r1 n3 n1
p-value < 2.2e-16
99 percent confidence interval:
11.10877 11.57538
sample estimates:
[1] 11.33968
つまり、底辺私立医大は私立医大一般よりも11倍裏口と関係が深いことが確認された。
その99%信頼区間は11.1-11.6と信頼性の高い数値であることも判明した。

68:卵の名無しさん
18/05/25 11:04:33.11 8Ot6XASs.net
底辺私立医と国立で比べると
裏口関係がリスク比は
p-value < 2.2e-16
99 percent confidence interval:
214.1801 223.5677
sample estimates:
[1] 218.8235
実に200倍以上の数値を得る。

69:卵の名無しさん
18/05/25 11:05:22.23 8Ot6XASs.net
底辺私立医と国立で比べると
裏口との関係がある、リスク比は
p-value < 2.2e-16
99 percent confidence interval:
214.1801 223.5677
sample estimates:
[1] 218.8235
実に200倍以上の数値を得る

70:卵の名無しさん
18/05/25 13:04:02.31 XgMJ+HGa.net
実は国立でも昔から裏口がある。
阪奈香だけではない。
鹿饂飩あたりは常態的にやってただろう、
そういう土地柄だ。
今は推薦、AO、地域枠で
どこでも堂々とやってるな。

71:卵の名無しさん
18/05/25 13:26:53.01 8Ot6XASs.net
ド底辺シリツ医大にも地域枠があってびっくりした。
額に地底と刺青してほしいねw

72:卵の名無しさん
18/05/25 13:28:31.50 8Ot6XASs.net
ド底辺特殊シリツ医大は受験しただけでも人生の汚点。
進学した方の悲劇をご報告いたします。
親の期待に反して底辺シリツにしか進学できない子供をもつとこうなるそうです。
赤裸々な告白より引用
(quote)
恥ずかしながら実はうちの親がそうなんです
私立洗顔(ママ)でした
本人は隠したがっていたけど祖父が公言してたからね
祖父は息子(うちの父親)の学力など知らずに進路指導で当たり前に自分の出身医大(国立)を希望させ
父親の担任に無理だと断言されこれが一生で最大の屈辱だったそうです
祖父が何十年もたってるのに死ぬまで祖母をその事で叱りつけてました
お前がしっかり勉強させないからこうなったっとね
自分は教育に失敗したが祖父の口癖
これが私大医学部OBの悲しい現実

73:卵の名無しさん
18/05/25 13:54:47.65 8Ot6XASs.net
URLリンク(www.kawasaki-m.ac.jp)
特別推薦入試専願
特別
推薦
入試
専願
ですって
ド底辺
特殊
シリツ医大
裏口枠
限定
と条件反射で脳内変換してしまうな。

74:卵の名無しさん
18/05/25 18:35:37.49 8Ot6XASs.net
こんなスレが立ってるな。
【日大】学生に「大学のイメージ下げる行動は慎んで」と通達 現役学生「イメージ下げているのは大学。むしろ学生に謝罪して」と憤慨★3

ド底辺シリツ医大卒の学歴ロンダは歯止めをかけてほしいなぁ。

75:卵の名無しさん
18/05/25 18:59:41.78 8Ot6XASs.net
# 医学部 私立 裏口 24,800 件
# 医学部 国立 裏口 23,800 件
# 医学部 私立  2,430,000 件
# 医学部 国立  5,370,000 件
#
# ド底辺私立医 裏口3,480 件
# ド底辺私立医 12,200 件
r1=24800
r2=23800
r3=3480
n1=2430000 # 私立
n2=5370000 # 国立
n3=12200 # 底辺
prop.test(c(r1,r2,r3),c(n1,n2,n3))
prop.test(c(r1,r2),c(n1,n2)) # 私立vs国立
prop.test(c(r1,r3),c(n1,n3)) # 私立vs底辺
prop.test(c(r3,r2),c(n3,n2)) # 底辺vs国立
pairwise.prop.test(c(私立=r1,国立=r2,底辺=r3),c(n1,n2,n3),p.adjust='bonferroni')

76:卵の名無しさん
18/05/25 19:33:16.06 VgWUy42F.net
>>68
底辺バカ私立の東大ロンダ多いな
糞バカで恥知らずだが、
野心だけは凄いんだよな
国立の教授なったり
政治家になってるよ

77:卵の名無しさん
18/05/26 15:02:31.72 S2P1Zub3.net
統計の知識は臨床医に必須の教養だと思う。
寄付の有無が仮想大学の合否に関係するか2つの調査があったとする。
調査A 合格 不合格
寄付有  30  170
寄付無  10   90
合格率は15%と10%で寄付有が高い。p-value = 0.28と有意差なし(fisher.test)。
調査B 合格 不合格
寄付有  90   10
寄付無 170   30
合格率は90%と85%で寄付有が高い。p-value = 0.28と有意差なし。
単純に足し算すると
    合格 不合格
寄付有 120  180
寄付無 180  120
合格率は40%と60%で寄付無が高い。p-value = 1.349e-06と有無を言わさぬ有意差あり。
これと同じ論理が展開されたのがEMPA-REG OUTCOME試験
#             Subjects   CV-events   %   HR   p value 95%CI
# Placebo            2333  282      12.1  
   
# Empagliflozin 10mg     2345  243      10.4   0.85   0.07    0.72-1.01
# Empagliflozin 25mg     2342  247      10.5   0.86   0.09    0.73-1.02
# Empagliflozin 10mg+25mg 4687  490      10.5   0.86   0.04    0.74-0.99
ド底辺シリツ医大卒の医者なら簡単に騙されるだろうな。有意差がでなかったから別剤型の症例を加えてp<0.05にしただけ。
どうせなら、メタアナリスを使ってこういう風に解析すればいいのに URLリンク(i.imgur.com)

78:卵の名無しさん
18/05/27 14:06:38.44 YM9KEdJB.net
>>5
>同じロジックで比率の比較をするスクリプト
スレリンク(hosp板:498番)
を少し、refineして著者のColquhun氏に送付したら、
比率については考えていなかったので、ありがとさんという趣旨の返事がきた。
underlying mathsを書いてくれというのだが、専門家に非心χ二乗分布の講義をするわけにもいかず、
どうしようかなと思っている。自分自身日本語でかかれた統計数理を学習しながら書いたスクリプトなんだよなぁ。
帰無仮説が正しいときの非心χ二乗分布の非心パラメータ設定を書いた適当な英語文献を検索中。

教育者ってきちんと対応してくれて好感がもてる、Modern EpidemiologyやIntuitive Biostatisticsの著者からも返事をもらえたし。

79:卵の名無しさん
18/05/27 14:59:25.69 F/qyBnY/.net
でも私立医大卒の医者って羨ましいよな。
だって親がお金持ちじゃ無かったら
本来は医者じゃない連中なんだぜ。

80:卵の名無しさん
18/05/27 15:07:54.55 JRoGGzmX.net
逆正弦法近似でないが
非心Χ二乗分布の解説がみつかったので
それに合わせてスクリプトを書き直して著者に送付してみた。
URLリンク(digitalcommons.calpoly.edu)
比率のFPR計算のweb calculatorのプロトタイプになればいいなぁ。
GUIのプログラムは俺には作れないから。

81:卵の名無しさん
18/05/27 15:14:43.96 6yKK2V+Q.net
いえいえ全ての医者は頭悪いです。
俺医師の言うこと聞かない患者として有名wwwww
何で俺よりアホな医者の言うこと聞かないとあかんの?
今の医者特に総合病院の勤務医はサラリーマンよりアホwwwwww
そらな、検査漬けと薬漬けしかしない医師に賢い医者がいないのは当然。
向上


82:心のある医者なら過去の知識を忘れず常に新しい知見を吸収してるが 今の医者は基礎的な医学知識もないもんなwwwwww そらバカにされて当然だろwwwww



83:卵の名無しさん
18/05/27 15:17:32.28 YM9KEdJB.net
>>75
本当に頭のいい奴は医学部にいかずに理学部か工学部に行く。
本当に頭の悪い奴が行くのがド底辺シリツ医大であることは申し上げるまでもない。

84:卵の名無しさん
18/05/27 15:20:41.66 YM9KEdJB.net
>71のような記述があると
調査Aと調査Bを統合した結果を出してみたくなるよね?
寄付金無で不合格になるリスクで考えて、この2つの調査を統合してみる。
調査A 不合格 合格
寄付無  90  10
寄付有 170  30
寄付無での不合格リスク比0.90/0.85= 1.06
調査B 不合格 合格
寄付無 30  170 
寄付有 10  90 
寄付無での不合格リスク比0.15/0.10= 1.50
ド底辺シリツ裏口バカの計算(単純な足し算、流石に足し算くらいできるだろ?)
   不合格 合格
寄付無 120  180
寄付有 180  120
調査A、Bの統合リスク比(ド底辺法w)0.40/0.60=0.67
問題:メタアナリシスの手法を用いて調査Aと調査Bを統合すると
寄付有に対して寄付無で不合格になるリスク比を求めよ。

85:卵の名無しさん
18/05/27 15:23:16.72 JRoGGzmX.net
>>75
勤務医って基本中サラリーマンですが何か?

86:卵の名無しさん
18/05/27 15:27:54.48 YM9KEdJB.net
>>75
自分より馬鹿な医者に診てもらいたくない、という患者はいるよ。
その点、ド底辺シリツ医大卒の医者は最優良患者といえる(笑)

87:卵の名無しさん
18/05/27 19:07:22.46 YM9KEdJB.net
>>77
原理を忘れないようにパッケージを使わずに関数を書いてみた。
a=c(90,30)
b=c(10,170)
c=c(170,10)
d=c(30,90)
RR.FE <- function(a,b,c,d){ # contigency table
lrr=log((a/(a+b))/(c/(c+d))) # effect=log(rr)
v=1/a-1/(a+b)+1/c-1/(c+d)
se=sqrt(v)
rr=exp(lrr)
z=qnorm(.975)
lwr=exp(lrr-z*se)
upr=exp(lrr+z*se)
study=cbind(rr,lwr,upr)
LRR=weighted.mean(lrr,1/v) # weight=1/v
SE=sqrt(1/sum(1/v))
list(study=study,
FE=c(estimate=exp(LRR),
LWR=exp(LRR-z*SE),
UPR=exp(LRR+z*SE)))
}
> print(RR.FE(a,b,c,d),digits=3)
$study
rr lwr upr
[1,] 1.06 0.970 1.16
[2,] 1.50 0.764 2.94
$FE
estimate LWR UPR
1.065 0.976 1.161

88:卵の名無しさん
18/05/27 21:17:26.34 YM9KEdJB.net
     不合格 合格
寄付無 a    b     
寄付有 c    d
このリスク比a/(a+b) ÷ c/(c+d)をRRとするとlog(RR)の分散vが 1/a -1/(a+b) +1/c -1/(c+d) になるとして信頼区間計算をしている。
調査の数だけその分散がでるわけでこれをViとすると
その複数の分散の逆数1/Viでの重み付き平均を固定モデルでの統合対数リスク比として算出する。
まあ、そういうものとして計算したというのでいいんだが、
その統合対数リスク比の分散は1/(Σ1/Vi)である、との記述がどうもよくわからん。
統計数理が追えないときは、シミュレーションで納得することにしよう。 困ったときのシミュレーション。
>77のデータをつかってのシミュレーションが>80の値になるかやってみる。
> RR.FE(a,b,c,d)
$study
rr lwr upr
[1,] 1.058824 0.9701055 1.155655
[2,] 1.500000 0.7643306 2.943752
$FE
estimate LWR UPR
1.0649511 0.9764301 1.1614972
シミュレーション結果はこれ、Bの近似に難があるにもかかわらず、メタアナリシスの信頼区間はほぼ一致している。
URLリンク(i.imgur.com)

89:卵の名無しさん
18/05/27 23:06:10.96 YM9KEdJB.net
>>81
Bの分布をみると正規分布で近似するには何があると思う、
むしろシミュレーション結果の方が現実に近いのかもしれん。

90:卵の名無しさん
18/05/28 07:20:10.79 c2RA2ZIc.net
ベイジアンネットワークメタアナリシスという手法を用いたという論文の


91:。評がメール配信されてきた。 ネットワークメタアナリシスって こういうことらしい。 アホ>ドアホ バカ=アホ ドアホ>>シリツ から バカ>>シリツという結論を導く手法。



92:卵の名無しさん
18/05/28 11:44:46.02 c2RA2ZIc.net
やはり、医師の世界は学歴で評価が決まるようだな。
スレリンク(hosp板:86-88番)


次ページ
最新レス表示
スレッドの検索
類似スレ一覧
話題のニュース
おまかせリスト
▼オプションを表示
暇つぶし2ch

620日前に更新/511 KB
担当:undef