底辺私立医大を卒業した医者って頭悪いよね? Part3 at HOSP
[2ch|▼Menu]
[前50を表示]
550:卵の名無しさん
18/05/02 16:00:52.37 HRauQnuh.net
>>548
命題「馬鹿は、死ななきゃ治らない」の対偶を教えてくれ仲間はいないみたいだな。
ド底辺シリツ医大卒の周りって裏口バカだらけなのか?
それとも馬鹿すぎて相手にされないのか?
一次方程式すらまともに立てられない馬鹿であることは実例を挙げた。
URLリンク(imagizer.imageshack.com)

551:卵の名無しさん
18/05/02 16:01:56.12 HRauQnuh.net
一次方程式で解決できる問題ができないバカはまあ、普通にいるだろう。
ところが、ド底辺シリツ医大卒は
付与の条件では計算できないと断言しちゃうんだよね。
馬鹿すぎる。
計算できないのはド底辺シリツ医大卒のオツムの問題なんだけどね。
URLリンク(imagizer.imageshack.com)
こんな算数もできないなんて馬鹿なのだろうか?
馬鹿なのにどうして大学に入れるの?
中学数学の問題だから高校入試ですら落ちても不思議ないのだが。
こんな馬鹿なのでどうしてド底辺シリツとはいえ、
大学に入れるわけ?
やっぱり、あれ?
,二     |  ̄ヽ    . /〃   -┼─   
  )    |___.     (         レ--、 
 .ノ      __)      \.      _ノ . 

552:卵の名無しさん
18/05/02 16:25:24.20 HRauQnuh.net
医師の世界では
 頭悪すぎ
に続く言葉は、やはり
 底辺私立
スレリンク(hosp板:428番)
428 卵の名無しさん sage 2018/05/02(水) 16:01:01.76 ID:V4E7oyFY
リベラルと言えばいけてると思ってんの頭悪すぎ。底辺私立か?

553:卵の名無しさん
18/05/02 16:27:43.48 HRauQnuh.net
ちはやふる は 神 の
頭悪すぎ は 底辺私立 の
が枕言葉なんだな、φ(..)メモメモ

554:卵の名無しさん
18/05/02 20:48:24.46 n+hyj4Vz.net
呪文:「ド底辺シリツ医大が悪いのではない、本人の頭が悪いんだ」を唱えると甲状腺機能が正常化するという統計処理の捏造をやってみる。
TSH : 0.34〜4.0 TSH:μIU/mlを基準値として(これが平均±2×標準偏差で計算されているとして)呪文前と呪文後のデータを正規分布で各々50個つくる。
負になった場合は検出限界以下として0にする。
URLリンク(i.imgur.com)
同じ平均値と標準偏差で乱数発生させただけなので両群には有意差はない。
横軸に呪文前のTSHの値、縦軸にTSHの変化(呪文後−呪文前)をグラフしてみると
URLリンク(i.imgur.com)
つまり、呪文前のTSHが高いほど呪文後はTSHは下がり、呪文前のTSHが低いほど呪文後のTSHは上がるという傾向がみてとれる。
これを線形回帰して確かめてみる。
回帰直線のパラメータは以下の通り。
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 2.3812 0.3323 7.165 4.10e-09 ***
before -1.0351 0.1411 -7.338 2.24e-09 ***
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
Residual standard error: 0.8365 on 48 degrees of freedom
Multiple R-squared: 0.5287, Adjusted R-squared: 0.5189
F-statistic: 53.84 on 1 and 48 DF, p-value: 2.237e-09
p = 2.24e-09 なので有意差ありとしてよい。
ゆえに、「ド底辺シリツ医大が悪いのではない、本人の頭が悪いんだ」という呪文は甲状腺機能を正常化させる。
ここで問題、この統計処理のどこが誤っているか?

555:卵の名無しさん
18/05/02 21:02:45.39 n3JqitZ1.net

て、あんたお医者さん違うやん

556:卵の名無しさん
18/05/02 21:03:44.94 n+hyj4Vz.net
乱数を変えて再検しても、同様のグラフになる。
URLリンク(i.imgur.com)

557:卵の名無しさん
18/05/02 21:09:44.36 n+hyj4Vz.net
>>555
第三法則実証乙。
>483 および >486参照
統計処理くらい理解していれば、>554に答えてみ!
ひょっとして、底辺シリツ医大卒の問題そのものがわからない馬鹿なのか?
馬鹿なのにどうして大学に入れたのだろうか?
やっぱり、あれ?

,二     |  ̄ヽ    . /〃   -┼─   
  )    |___.     (         レ--、 
 .ノ      __)      \.      _ノ .

558:卵の名無しさん
18/05/02 21:11:54.95 n+hyj4Vz.net
ド底辺シリツ医大卒には数式や英文のある問題は無理みたいだなぁ。
これなら答えられるかな?
URLリンク(imagizer.imageshack.com)

Xの先輩Yは医師である。
Yは真正の裏口入学(ガチの裏)と評価される学生Zを知っている。
Zの入学した大学をガチ裏医大と呼ぶことにする。
Zが1回生から留年を繰り返すことを知っているのでYとZは同期入学であると思われる。
問題:Xの入学した大学はどこか?

559:卵の名無しさん
18/05/02 21:24:05.97 jAwYPU+x.net
国試勉強休憩中。下には下がいるのだなとなぜかホッとする自分がいる。なんとか国試合格したい。

560:卵の名無しさん
18/05/02 21:39:35.00 jAwYPU+x.net
現時点で彼女の方が勉強が進んでいる。なんとか追いつきたい。

561:卵の名無しさん
18/05/03 00:32:28.86 0JFjpoYv.net
Mistake: Over interpreting a small P valueという例示
回帰直線の傾き=0の帰無仮説はp値は0.03で棄却されるが、
相関係数R2は 0.009407で1%未満であるとの説明。
横軸:一様分布、縦軸:正規分布はみてとれたが
nを増やしただけでは有意差が簡単にはでてこない。試行錯誤で乱数の種を虱潰しに当たらせて完成。
グラフがようやく再現できた。
URLリンク(i.imgur.com)
Call:
lm(formula = y ~ x)
Residuals:
Min 1Q Median 3Q Max
-7.2026 -1.4161 -0.0427 1.2997 9.7284
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 99.66681 0.18517 538.232 <2e-16 ***
x 0.01374 0.00632 2.175 0.0301 *
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
Residual standard error: 2.039 on 498 degrees of freedom
Multiple R-squared: 0.009407, Adjusted R-squared: 0.007417
F-statistic: 4.729 on 1 and 498 DF, p-value: 0.03013
2.5 % 97.5 %
(Intercept) 99.302994240 100.03063277
x 0.001326443 0.02616038

562:卵の名無しさん
18/05/03 00:41:57.72 0JFjpoYv.net
> cor.test(x,y)
Pearson's product-moment correlation
data: x and y
t = 2.1746, df = 498, p-value = 0.03013
alternative hypothesis: true correlation is not equal to 0
95 percent confidence interval:
0.00937691 0.18312091
sample estimates:
cor
0.09698775
URLリンク(i.imgur.com)
p値が小さくても回帰直線はほぼ水平。

563:卵の名無しさん
18/05/03 00:46:57.83 0JFjpoYv.net
裏口馬鹿よりlinear regressionの方が頼れると思うぞwww
URLリンク(www.explainxkcd.com)
URLリンク(www.explainxkcd.com)

564:卵の名無しさん
18/05/03 07:29:51.51 9nPOoDoJ.net
底辺私立医は自分がバカという自覚すらないだろ。
家族にもバカにされてるくせに。

565:卵の名無しさん
18/05/03 08:42:56.61 0JFjpoYv.net
URLリンク(i.imgur.com)
この2つの分布の相関係数とp値は
cor p
0.197477921301 0.000008644366
cor p
-0.000005514402 0.999901862465

線形回帰を使うと捏造が容易だな。
例えば、統計的に有意な方を毎日オナニーをするというド底辺シリツ医大卒10人を50か月に亘って知能指数を調べた結果のグラフとすると
相関係数が正なので p=0.000008644366でド底辺シリツ医大卒のオナニーは知能指数を上昇させるという結論が導きだせる。

格言 : 統計と女の涙は信じるな
このスレ的には 
 ド底辺シリツ医大卒と裏口バカは信じるな
というのが正しいwww

566:卵の名無しさん
18/05/03 08:45:46.33 0JFjpoYv.net
>>564
誇るべき母校に
明らかにおかしなバカ学生がいることを自慢するのが、ド底辺シリツ医大卒。
しかもガチの裏と裏口入学の存在を前提としている。
URLリンク(imagizer.imageshack.com)

567:卵の名無しさん
18/05/03 08:52:22.84 0JFjpoYv.net
>>559
>下には下がいる
学力で判定するとド底辺シリツ医大卒の下には下はないんじゃね?
ガチの裏ですって
URLリンク(imagizer.imageshack.com)

568:卵の名無しさん
18/05/03 09:10:22.59 tiEdtl/y.net
これだけまとまった勉強をするのは久しぶりだな。充実している。

569:卵の名無しさん
18/05/03 09:18:19.00 0JFjpoYv.net
>>547
> まずは基本的なところをじっくり確認
このスレのタイトルを確認すればいいだけだろ。
それすらわからんとは馬鹿なのだろうか?
馬鹿なのにどうして大学に入れたのだろうか?
やっぱり、あれ?

,二     |  ̄ヽ    . /〃   -┼─   
  )    |___.     (         レ--、 
 .ノ      __)      \.      _ノ .

570:卵の名無しさん
18/05/03 09:21:32.62 0JFjpoYv.net
当直医スレから確認させている基本法則はこれね!

Three rudimentary 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.
ド底辺特殊シリツ医大卒は裏口入学の負い目から裏口馬鹿を暴く人間を偽医者扱いしたがる。

571:卵の名無しさん
18/05/03 09:24:56.13 IqRiigyv.net
帝国男児ヨ! 括目セヨ!!
2013年03月02日20:36
医学部の学費で打線組んだwwwww
(中)金沢医科大学医学部 3950万円
(二)福岡大学医学部   3770万円
(左)川崎医科大学医学部 4565万円
(投)帝京大学医学部   4920万円
(一)東海大学医学部   4184万円
(右)北里大学医学部   3890万円
(捕)兵庫医科大学医学部 3880万円
(左)埼玉医科大学医学部 3800万円
(遊)愛知医科大学医学部 3800万円

572:卵の名無しさん
18/05/03 09:50:52.94 0JFjpoYv.net
>>565
Intelligence and semen quality are positively correlated
URLリンク(www.sciencedirect.com)
Why is intelligence correlated with semen quality?
URLリンク(www.ncbi.nlm.nih.gov)
という論文があったりするので案外、偶然以上の有意差がでるかも知れないが
たとえそうであってもド底辺シリツ医大卒の裏口バカがサンプルだと、バイアスがかかりすぎだな。
なんせ、一次方程式すら立てられない馬鹿でしたからね。

573:卵の名無しさん
18/05/03 09:51:56.15 0JFjpoYv.net
タイプミス修正
当直医スレから確認されている基本法則はこれね!

Three rudimentary 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.
ド底辺特殊シリツ医大卒は裏口入学の負い目から裏口馬鹿を暴く人間を偽医者扱いしたがる。

574:卵の名無しさん
18/05/03 09:58:55.81 uUyq0HNp.net
私立を差別するつもりは学生時代なかったが、
あまりのバカさと仕事も社会も舐めた態度、
患者殺しの実態を多々経験すると、
差別どころか、同じ空気を吸うのも悍ましい。

575:卵の名無しさん
18/05/03 10:06:33.94 g3ZJipUc.net

で?あんたは?
医学生というか
国試失敗連続23回だろ?
あのね、そーゆーことはお医者さんになってから書いた方がいいよ

576:卵の名無しさん
18/05/03 10:18:42.63 tiEdtl/y.net
勉強会はペースをつかめるし勉強が楽しくなるからいいな。

577:卵の名無しさん
18/05/03 10:21:25.22 ft8gnQDo.net
抜ける看護婦揃い
URLリンク(www.tsukubahigashi-hp.com)

578:卵の名無しさん
18/05/03 11:49:16.37 tiEdtl/y.net
さて、みんなと寿司でも食いにいくかな。

579:卵の名無しさん
18/05/03 11:54:44.32 0JFjpoYv.net
>>572
Intelligence and semen quality are positively correlated
URLリンク(www.sciencedirect.com)
のアブストラクトにはサンプルサイズの記載がないという不親切設計。しかもp値が有効数字一桁。
We found significant positive correlations between intelligence and 3 key indices of semen quality:
log sperm concentration (r = .15, p = .002),
log sperm count (r = .19, p < .001), and
sperm motility (r = .14, p = .002) in a large sample of US Army Veterans.
どうも統計処理が腑に落ちないので検討してみた。
concentrationの相関係数0.15、p値0.002からサンプルサイズを逆算すると422人
有効数字1桁のp値は0.0015〜0.0025の間とすると404人から445人ということになる。
サンプルサイズがこの範囲にあるときにmotilityの相関係数が0.14なら
そのp値は0.00308 〜 0.00482になるはず。
こちらのp値0.002が正しいとするとサンプルサイズは464〜511人になって
concentrationから逆算した区間と重なりはないので、どちらかの記述が間違っていることがわかる。
検索したら
URLリンク(www1.udel.edu)
で全文がよめてサンプルサイズは425人であった。
すると、sperm motilityのp値は0.0038が正しい。
log sperm count (r = .19)でのp値は0.000081なのでp < .001はp < .0001と書いてもよいくらいである。
Rのコードはこれ。
スレリンク(math板:194番)
統計の衣を纏っていても検算は必要ってこったね。

580:卵の名無しさん
18/05/03 11:59:53.76 0JFjpoYv.net
>私立を差別する
人種や民族など本人に選択の余地のないことを根拠に区別するのは差別であるが、
ド底辺シリツ医大進学は本人が選択したのだから、それを根拠に区別することは差別ではない。
シリツ医大の法外な学費が払える経済状況なら教育投資して学力を高めることができたはず。
金をかけても学力が身につかなかった馬鹿なのだから蔑まれても当然。
これは昔から不変。
(quote)
私は昭和の時代に大学受験したけど、昔は今よりも差別感が凄く、慶応以外の私立医は特殊民のための特殊学校というイメージで開業医のバカ息子以外は誰も受験しようとすらしなかった。
常識的に考えて、数千万という法外な金を払って、しかも同業者からも患者からもバカだの裏口だのと散々罵られるのをわかって好き好んで私立医に行く同級生は一人もいませんでした。
本人には面と向かっては言わないけれど、俺くらいの年代の人間は、おそらくは8−9割は私立卒を今でも「何偉そうなこと抜かしてるんだ、この裏口バカが」と心の底で軽蔑し、嘲笑しているよ。当の本人には面と向かっては絶対にそんなことは言わないけどね。
(unquote)
それにはこういうレスが続いた。
(quote)
それは昭和の時代だけでなく平成の今も同じです。
みんな口に出さないだけで思ってることは同じです。
そういう私も皆の前では本心とは裏腹に「国家試験に受かれば
私立も国立も変わらないよ。」と言ってます。
在日や帰化人を連想させる名前をした医師に「******(伏字は引用者)」と内心では思ってても、
ニコニコしながら話を合わせてるのと同じです。
(unquote)
>でも一生、卒業大学を隠すのは本当に辛い      
URLリンク(imagizer.imageshack.com)

581:卵の名無しさん
18/05/03 12:06:20.77 tiEdtl/y.net
w

582:卵の名無しさん
18/05/03 12:33:00.95 tiEdtl/y.net
少しハイペースで勉強してるので、4人で気分転換に海でも見にいこうとなった。

583:卵の名無しさん
18/05/03 12:45:14.57 g3ZJipUc.net
寿司に海か
たまいいなあ
算数やら2チャンネル貼ったりしてる
医師真似事務員バイトは
勉強しねえから気分転換に行く資格ないな
たまには国試の勉強しろよ
無職のくせに

584:卵の名無しさん
18/05/03 12:52:47.21 tiEdtl/y.net
ありがとうございます。みんなで国試合格を一発で決めたいです。

585:卵の名無しさん
18/05/03 13:24:59.32 0JFjpoYv.net
FisherのZ変換をつかって決定係数(アールスクエア)を求めると
> cor2ci(rs,425)^2
lower cor upper
conc 0.003099511 0.0225 0.05840578
count 0.009336431 0.0361 0.07843337
Motility 0.002068850 0.0196 0.05383846
精子の健全性は知性の数パーセントを説明することしかできない、という結果であった。
ド底辺シリツ医大卒でないこと が知性をもつ人間の必要条件であることは統計処理すら必要ない。
小学生からナマポも知っている常識であることは既述。

586:卵の名無しさん
18/05/03 13:26:29.15 g3ZJipUc.net
頑張ってね
ここにはウリュウひろゆきジイさんという
算数好きの国試浪人の成れの果てがいるからな
君たちの反面教師にはもってこいだよ

587:卵の名無しさん
18/05/03 13:31:25.65 0JFjpoYv.net
>>584
医師国家試験は落ちる方が難しい試験だが、
ド底辺シリツ医大に入学難関である。
なぜならば
こういうことだから。
卓見投稿を再引用
(quote)
私は昭和の時代に大学受験したけど、昔は今よりも差別感が凄く、慶応以外の私立医は特殊民のための特殊学校というイメージで開業医のバカ息子以外は誰も受験しようとすらしなかった。
常識的に考えて、数千万という法外な金を払って、しかも同業者からも患者からもバカだの裏口だのと散々罵られるのをわかって好き好んで私立医に行く同級生は一人もいませんでした。
(unquote)

つまり、一生、他人から蔑まれることを受け入れた人間でなければ選択できないのだからね。

当然の結果がこれね。
>でも一生、卒業大学を隠すのは本当に辛い      

URLリンク(imagizer.imageshack.com)

588:卵の名無しさん
18/05/03 13:42:17.40 tiEdtl/y.net
下には下がいると思うとどこかホッとする。

589:卵の名無しさん
18/05/03 13:53:39.55 tiEdtl/y.net
と仲間に話すと笑いながら悪趣味ねぇと言われた。

590:卵の名無しさん
18/05/03 14:17:16.90 nEmJkkk2.net
解釈の原因は解釈者自身の固定観念。解釈の自由には責任が伴う
言葉風紀世相の乱れはそう感じる人の心の乱れの自己投影。人は鏡
憤怒は一時の狂気、無知無能の自己証明。中途半端な知識主ほど激昂
「真実は一つ」は錯誤。執着する者ほど矛盾を体験(争い煩悩)
他人に不自由(制約)を与えれば己も不自由(不快)を得る
問題解決力の乏しい者ほど自己防衛の為に礼儀作法マナーを要求
情報分析力の低い者ほどデマ宗教フェイク疑似科学に感化洗脳
自己肯定感の欠けた者ほど「己の知見こそ全で真」に自己陶酔
人生経験の少ない者ほど嫌いキモイ怖いウザイ憎い想定外を体験
キリスト教は世界最大のカルト。聖書は史上最も売れているト本
全ては必然。偶然 奇跡 理不尽 不条理は思考停止 視野狭窄の産物
人生存在現象に元々意味価値理由目的義務使命はない
宗教民族領土貧困は争いの「原因」ではなく「口実動機言訳」
虐め差別犯罪テロ紛争は根絶可能。必要なのは適切十分な高度教育
体罰は指導力問題解決力の乏しい教育素人の独善甘え怠慢責任転嫁
死刑は民度の低い排他的集団リンチ殺人。「死ねば償える」は偽善
核武装論は人間不信と劣等感に苛まれた臆病な外交素人の精神安定剤
投票率低下は社会成熟の徴候。奇人変人の当選は議員数過多の証左
感情自己責任論 〜学校では教えない合理主義哲学〜 m9`・ω・)

591:卵の名無しさん
18/05/03 14:22:54.26 ffi6hOtq.net
金で医者になる私立医をボコろー
患者殺したら内部告発しょう
私立医をオリの中に入れて吊るそうぜ

592:卵の名無しさん
18/05/03 14:51:29.53 0JFjpoYv.net
>>588
一次方程式すら立てられないド底辺シリツ医大卒の裏口バカでも医者になれるらしいしねwwww

593:卵の名無しさん
18/05/03 15:02:09.20 0JFjpoYv.net
けがで死亡はO型多い、医科歯科大
URLリンク(mm.m3.com)
は無料で全文が読める。
URLリンク(ccforum.biomedcentral.com)
多重比較なので他の血液型とのpairwiseでのp値があるかと思ったらなかったのでがっかり。
与えられた数値で自分で再計算してみると
> riskratio(round(n.O*.28),round(n.other*.11),n.O,n.other)
Disease Nondisease Total
Exposed 80 204 284
Nonexposed 68 549 617
Risk ratio estimate and its significance probability
data: round(n.O * 0.28) round(n.other * 0.11) n.O n.other
p-value = 1.113e-10
95 percent confidence interval:
1.910340 3.419677
sample estimates:
[1] 2.555924
だから、多重比較によるバラつき検出ではないな。確かにp値が小さすぎる。
Bonnferonni補正しても結論は覆りそうにない。
うお座に心不全が有意に多いとう統計処理の誤謬指摘の論文を思いだした。
Pisces did not have increased heart failure: data-driven comparisons of binary proportions between levels of a categorical variable can result in incorrect statistical significance levels
URLリンク(www.jclinepi.com)(07)00188-6/abstract?code=jce-site

594:卵の名無しさん
18/05/03 15:03:33.35 0JFjpoYv.net
>>589
仲間ってド底辺シリツ医大卒の裏口バカ仲間かな?
馬鹿は死ななきゃ治らない の対偶を即答できない、馬鹿集団??

595:卵の名無しさん
18/05/03 15:26:11.82 tiEdtl/y.net
w

596:卵の名無しさん
18/05/03 15:36:42.26 WLvmaBgZ.net
対偶?
>>594
即答してやるよ
簡単だからみんなスルーしてるの気づけよ
答え:
おまえは国試浪人の
ウリュウひろゆきジイさんですね

どうだ?
大正解だろう!

597:卵の名無しさん
18/05/03 15:47:58.00 Jz8cV8mL.net
>>594
俺も、即答してやるよ
簡単だからみんなスルーしてるの気づけよ な
ほんま
答え:
おまえは大学に12年間かよった伝説の
ウリュウひろゆきジイさんですね

どうだ?
これも大正解だろう!

598:卵の名無しさん
18/05/03 15:54:21.30 kxDx/PCD.net
>>594
俺も、即答だ
簡単だからみんなスルーしてるの気づけよ!
そうだね!
答え:
おまえが3年高卒浪人してやっとこさ
まぎれこめた長崎は
おまえがはいった時代の偏差値45な
一応国立やけど、40年位前やしな

どうだ?
これも大正解だろ
あんた、これで私立の学校をあかんつーのもどうかとおもうぜ
ウリュウひろゆきさんよ
しかも、あんたまだ国家試験通ってないんやろ?
23年も
なにしてたん?
てか、今だになにしてんの?

599:卵の名無しさん
18/05/03 17:15:19.99 tiEdtl/y.net
海は広いな大きいなー

600:卵の名無しさん
18/05/03 17:49:21.68 kxDx/PCD.net
>>594
俺も即答したくなったよ
簡単だよほんま
答え:
そもそもウリュウひろゆきジイさんは
医師真似で、病院よ事務員バイトしている国試浪人である!
おんとし63歳!
どうだ?
これ、もっと大正解だろ

601:卵の名無しさん
18/05/03 17:54:04.77 Drs+1rKs.net
>>594
うんにゃ、俺が大正解を即答してやるよ
朝飯前だ
答え: 医科歯科卒業を自称しているが、競争率1倍の偏差値40台の長崎を12年かけて卒業した、いや、温情で卒業させてもらえた、
てかそもそも彼は医師ではない!
そう、それは!ウリュウひろゆきジイさんですね

どうだ? 大正解だ!

602:卵の名無しさん
18/05/04 09:10:17.72 xAGSA6fq.net
帝国男児ヨ! 括目セヨ!!
2013年03月02日20:36
医学部の学費で打線組んだwwwww
(中)金沢医科大学医学部 3950万円
(二)福岡大学医学部   3770万円
(左)川崎医科大学医学部 4565万円
(投)帝京大学医学部   4920万円
(一)東海大学医学部   4184万円
(右)北里大学医学部   3890万円
(捕)兵庫医科大学医学部 3880万円
(左)埼玉医科大学医学部 3800万円
(遊)愛知医科大学医学部 3800万円
私立医経営は美味しい美味しい商売

603:卵の名無しさん
18/05/04 10:28:39.67 6BNcb2/S.net
非線形回帰をRのnlsでやってみたが、決定係数は表示されない。
残差平方和から自分で計算しろってことか。
面倒なので関数化してみた。ついでに信頼区間の出し方を勉強して組み込んでみた。
間違いの指摘歓迎。
# To calculate R2, its CI and modified R2 from output of nls() and y
nls2R2 <- function(re.nls, y, conf.level=0.95){
sm = summary(re.nls) # re.nls :output of nonlinear least square function nls
SSalt = sum(sm$residuals^2) # Squared Sum(SS) of residuals
n = length(y) # y : response varialbles (dependent valuables)
SSnull = var(y)*(n-1) # total SS
R2 = 1-SSalt/SSnull
k = sm$df[1] # how many parameters
seR2 = sqrt((4*R2*(1-R2)^2*(n-k-1)^2)/((n^2-1)*(n+3)))
R2.lower = R2 + seR2*qnorm((1-conf.level)/2)
R2.upper = R2 - seR2*qnorm((1-conf.level)/2)
modified.R2 = 1-(n-1)/(n-k-1)*(1-R2)
data.frame(R2,R2.lower,R2.upper,modified.R2)
}
ド底辺シリツ医大卒の方のコメントはまず
>558に答えてからにしていただきたい。

604:卵の名無しさん
18/05/04 13:22:03.42 6BNcb2/S.net
同じような疑問をもつ人は他にもいた。
I googled and searched on stats.stackexchange but I cannot find the formula to calculate a 95% confidence interval for an R2 value for a linear regression. Can anyone provide it?
URLリンク(stats.stackexchange.com)
bootstrapしたら、という提言があったので
やってみた
> quantile(re$t,c(.025,.975))
2.5% 97.5%
0.9953611 0.9997229
> nls2R2(re.nls,y)
R2 R2.lower R2.upper Adjusted.R2
1 0.9971203 0.9951952 0.9990453 0.9958861
近似式と合致している。
まあ、ド底辺シリツ医大卒と裏口バカの合致率ほどではないかもしれないが。

605:卵の名無しさん
18/05/04 16:47:52.44 6BNcb2/S.net
>>604
Rのコードはこれ。
スレリンク(hosp板:477番)

606:卵の名無しさん
18/05/04 20:42:01.42 6BNcb2/S.net
重回帰分析の部分を読み直していたら
A Note on Screening Regression Equations
URLリンク(amstat.tandfonline.com)
に触れてあったので原文を読みたいなと思って、検索したら全文がここで読めた。
URLリンク(www.public.asu.edu)(Screening%20Regression%20Equations).pdf
ノイズデータから重回帰分析に使うパラメータを絞っていって有意差を出すシミュレーションが載っていたので
Rで追試してみた。
R2 (1st) = 0.6183043
p.value (1st) = 0.05389239
coef < 0.25 = 21
coef < 0.05 = 8
R2 (2nd) = 0.4742404
p.value(2nd) = 5.409517e-05
coef < 0.25 (2nd) = 20
coef < 0.05 (2nd) = 12
最初のスクリーンニンぐで回帰係数が0.25未満の説明変数だけを残して重回帰分析すると
p値0.05強のデータがp<0.0001に早変わり。
まるで、裏口バカを留年ビジネスでスクリーニングして医者にしてしまう、ド底辺シリツ医大さながらではないかwwww
Rのコードはこれスレリンク(math板:195番) 実行には全角→半角変換が必要。

607:卵の名無しさん
18/05/05 01:43:36.89 pCZ7kmcH.net
# ド底辺シリツ医大の合格者120人の入試の結果を分析したとする。
# 120人全員の成績の平均値は45、標準偏差は10であった。
# そのうち入学を辞退した20人の成績の平均値は60、標準偏差3であった。
# 100人の入学者のうち5人は特待生で成績の平均値は65、標準偏差は2であった。
# 残り95人を裏口バカと呼ぶことにする
# 入学辞退者、特待生、裏口バカの成績に有意差はあるかを分散分析で答えよ。
total.Mean <- function(Mean,N) sum(N*Mean)/sum(N)
total.SD <- function(Mean,N,SD) sqrt(sum((N-1)*SD^2 + N*(Mean-sum(Mean*N)/sum(N))^2)/(sum(N)-1))
(m.u=uniroot(function(x) total.Mean(c(60,65,x),c(20,5,95)) - 45,c(0,100))$root)
total.Mean(c(60,65,m.u),c(20,5,95))
(m.sd=uniroot(function(x) total.SD(c(60,65,m.u),c(20,5,95),c(3,2,x)) - 10 , c(0,100))$root)
total.SD(c(60,65,m.u),c(20,5,95),c(3,2,m.sd))
辞退者=c(60,20,3)
特待生=c(65, 5, 2)
裏口バカ=c(m.u,95,m.sd) # c(40.79, 95, 6.13)
# Totalc(45 120 10)
dat=rbind(辞退者,特待生,裏口バカ)
colnames(dat)=c("Mean","N","SD")
mean.G=sum(dat[,"Mean"]*dat[,"N"])/sum(dat[,"N"])
SS.bit=sum((dat[,"Mean"]-mean.G)^2*dat[,"N"])
SS.wit=sum(dat[,"SD"]^2*(dat[,"N"]-1))
df.bit=nrow(dat)-1
df.wit=sum(dat[,"N"]-1)
MS.bit=SS.bit/df.bit
MS.wit=SS.wit/df.wit
(F.ratio=MS.bit/MS.wit)
pf(F.ratio,df.bit,df.wit,lower.tail=FALSE)
(η2=(SS.bit)/(SS.bit+SS.wit)) # eta^2

608:卵の名無しさん
18/05/05 07:20:00.69 hBjvmFP3.net

おじいちゃん
国家試験通ってないんやろ
なんなんだ?それ?毎日同じのばっかり
勉強しろよ

609:卵の名無しさん
18/05/05 07:47:50.69 TePPONvo.net
偽善者と偽装社会 とパナマ文書
■涙の告発「子どもが生贄に…」
動画:URLリンク(youtu.be)
フリン補佐官がスキャンダルに巻き込まれたのは、実はフリンが小児性愛ネットワークを捜査していたから。
フリン補佐官が追いかけていた小児性愛ネットワークの一人がペンス副大統領やプリーバス首席補佐官ともつながった。
URLリンク(nqj17437.wor...press.com)

610:卵の名無しさん
18/05/05 08:57:31.90 WOhQQmHo.net
卒後20年のオナニー中毒よ、学生にも敗北したようだな。もう国試は無理だ。一生妄想の中で生きて行け。

611:卵の名無しさん
18/05/05 09:32:33.76 pCZ7kmcH.net
分散分析で有意差が確認できたのでpairwiseに検定する。
> Tk=TukeyHSD(aov(score~group,data=data)) ; Tk[[1]]
diff lwr upr p adj
辞退-特待 -5.00 -11.66109 1.661085 1.801228e-01
裏口-特待 -24.21 -30.32263 -18.097369 4.329870e-14
裏口-辞退 -19.21 -22.48753 -15.932467 1.310063e-14
> pairwise.t.test(score,group)
Pairwise comparisons using t tests with pooled SD
data: score and group
特待 辞退
辞退 0.077 -
裏口 1.1e-15 < 2e-16
P value adjustment method: holm
特待生は入学辞退すべきであったことが示唆されるwww

612:卵の名無しさん
18/05/05 09:36:27.83 pCZ7kmcH.net
ド底辺シリツ医大卒って英文や数字が出ると全くコメントできずに第三法則を実証するだけ。
馬鹿なのだろうか?

URLリンク(i.imgur.com)
統計は教養以前の問題。
それすら扱えないド底辺シリツ医ってなんと表現すればいいのだろう?
当事者の間では ガチの裏 と表現されてらしいw
URLリンク(imagizer.imageshack.com)

馬鹿なのにどうして大学に入れただろう?
やっぱり、あれ??

,二     |  ̄ヽ    . /〃   -┼─   
  )    |___.     (         レ--、 
 .ノ      __)      \.      _ノ .

613:卵の名無しさん
18/05/05 09:44:35.04 pCZ7kmcH.net
>>608
>なんなんだ?それ?
やはり、ド底辺シリツ医大卒には理解できないらしな。
新薬説明の統計を裏口バカが理解できているはずもないなぁ。

>毎日同じのばっかり
傾向スコアの話から多重回帰に話が移っているんだけどな
区別すらできないらしいな。
ド底辺シリツ医大を話題にしたクイズで興味が沸くかと思ったが、そもそも、向学心がない。
ひょっとして馬鹿なの?
馬鹿なのにどうして大学に入れたの?
ド底辺シリツ医大でも一様、大学だろ?
やっぱり、あれ??

,二     |  ̄ヽ    . /〃   -┼─   
  )    |___.     (         レ--、 
 .ノ      __)      \.      _ノ .

614:卵の名無しさん
18/05/05 09:54:11.66 pCZ7kmcH.net
□□□□□□□■□□□□□■□■□
□□□□□□□■□□□□□■□■□
□□□□□□□■□□□□□□□□□
□□□□□□□■□□□□□□□□□
□□■■■■■■■■■■■■□□□
□□□□□□□■□□□□□■□□□
□□□□□□□■□□□□□■□□□
□□□□□□■■□□□□□■□□□
□□□□□□■□□□□□□■□□□
□□□□□□■□□□□□□■□□□
□□□□□■■□□□□□□■□□□
□□□□□■□□□□□□■■□□□
□□□□■■□□□□□□■□□□□
□□□■■□□□□□□□■□□□□
□□■■□□□□□■■■■□□□□
□□□□□□□□□□□□□□□□□
□□□□□□□□□□□□□□□□□
□□□□□□□□□□□■■□□□□
□□□□□□□□■■■■□□□□□
□□□■■■■■■□□□□□□□□
□□□□□□□□■□□□□□□□□
□□□□□□□□■□□□□□□□□
□□□□□□□□■□□□□□□□□
□■■■■■■■■■■■■■■□□
□□□□□□□□■□□□□□□□□
□□□□□□□□■□□□□□□□□
□□□□□□□■■□□□□□□□□
□□□□□□□■□□□□□□□□□
□□□□□□■■□□□□□□□□□
□□□□□■■□□□□□□□□□□
□□□■■■□□□□□□□□□□□ 👀
Rock54: Caution(BBR-MD5:1341adc37120578f18dba9451e6c8c3b)


615:卵の名無しさん
18/05/05 09:54:20.78 pCZ7kmcH.net
□□□□□□□■□□□□□□□□□
■■■■■■■■■■■■■■■■□
□□□□□□□□□□□□□□□□□
□□□■■■■■■■■■■□□□□
□□□■□□□■□□□□■□□□□
□□□■■■■■■■■■■□□□□
□□□■□□□■□□□□■□□□□
□□□■■■■■■■■■■□□□□
□□□□□□□■□□□□□□□□□
□□■■■■■■■■■■■■□□□
□□□□□□□■□□□□□□□□□
■■■■■■■■■■■■■■■■□
□□□□□■□□■□□□□□■□□
■■■■■□□□□■■□□■□□□
□□□□■■■■□□■■■□□□□
□□■■■□□□□□□□■■■■□
□□□□□□□□□□□□□□□□□
□□□□□□□□□□□□□□□□□
□□■■■■■■■■■■■■□□□
□□■□□□□□□□□□□■□□□
□□■□□□□□□□□□□■□□□
□□■□□□□□□□□□□■□□□
□□■□□□□□□□□□□■□□□
□□■□□□□□□□□□□■□□□
□□■□□□□□□□□□□■□□□
□□■□□□□□□□□□□■□□□
□□■□□□□□□□□□□■□□□
□□■□□□□□□□□□□■□□□
□□■□□□□□□□□□□■□□□
□□■■■■■■■■■■■■□□□
□□□□□□□□□□□□□□□□□
□□□□□□□□□□□□□□□□□ 👀
Rock54: Caution(BBR-MD5:1341adc37120578f18dba9451e6c8c3b)


616:卵の名無しさん
18/05/05 09:55:18.98 pCZ7kmcH.net
ガチの裏口同窓生当事者の投稿です。
URLリンク(imagizer.imageshack.com)

617:卵の名無しさん
18/05/05 09:56:53.26 WOhQQmHo.net
ふん。ぶざまだな。

618:卵の名無しさん
18/05/05 10:08:07.99 rMcPOeHp.net
基本自分は私立中堅以下卒の医者にはかからないようにしている。ホームページ見たら経歴書いてたりするから事前に調べて。

619:卵の名無しさん
18/05/05 10:09:13.05 pCZ7kmcH.net
??????????????????????????????????????????????????????????????????????????????????????????????????
????■????■?????????????????????????????????■?????????????????????????????????????????????????????????
????■???■?????????■■??????????????????■■■■■■■■■■■????????????????????????????????????????????????????
????■?????????■■■■■?????????■■■■■??????■■■■■■■■■?????????????????????????????????????????????????????
■■■■■■■■■????????■?????????■??■??■■????■???■???■?????????????????????????????????????????????????????
????■???■????????■????????■???■???■■???■■■■■■■■■?????????????????????????????????????????????????????
????■???■???■■■■■■■■■■■??■????■????■???■■■■■■■■■?????????????????????????????????????????????????????
????■???■????????■???????■???■?????■??■■■■■■■■■■■????????????????????????????????????????????????????
???■????■????????■???????■???■?????■?■■■■■■■■■■■■■???????????????????????????????????????????????????
???■????■???????■■???????■■?■?????■■?????■■?■??■?????????????????????????????????????????????????????
??■????■■??????■■?????????■■■????■■???■■■■???■■??????????????????????????????????????????????????????
■■???■■■??????■■???????????????■■????????■■■■?■■■????????????????????????????????????????????????????
???????????????????????????????????????■■???????■?

620:卵の名無しさん
18/05/05 10:10:14.30 pCZ7kmcH.net
>>617
ド底辺シリツ医大卒の汚名より、ぶざまな事実は存在しない。
ちょっと、第二法則を確認してみましょう。
あんた、どこ卒?

621:卵の名無しさん
18/05/05 10:11:39.94 pCZ7kmcH.net
>>618
手先が不器用な奴は国立卒にもいるが、
ド底辺シリツ医大卒だと頭が不器用である確率が高い。
統計どころか、中学数学すら怪しい。

622:卵の名無しさん
18/05/05 10:13:29.85 WOhQQmHo.net
女医が見て嘲笑ってるぞw

623:卵の名無しさん
18/05/05 10:33:27.25 iiQZU4qb.net
ウリュウジイさん
事務員バイトで満足だろ
最近はむつかしいし、
ウリュウ先輩には国試むりだよ

624:卵の名無しさん
18/05/05 10:37:58.26 pCZ7kmcH.net
スキャッチャードプロットの手法
Wikipediaにはきちんと問題点が指摘されている。
URLリンク(ja.wikipedia.org)
呪文:「ド底辺シリツ医大が悪いのではない、本人の頭が悪いんだ」を唱えると甲状腺機能が正常化するという統計処理の捏造をやってみる。
TSH : 0.34〜4.0 TSH:μIU/mlを基準値として(これが平均±2×標準偏差で計算されているとして)呪文前と呪文後のデータを正規分布で各々50個つくる。
負になった場合は検出限界以下として0にする。
URLリンク(i.imgur.com)
同じ平均値と標準偏差で乱数発生させただけなので両群には有意差はない。
横軸に呪文前のTSHの値、縦軸にTSHの変化(呪文後−呪文前)をグラフしてみると
URLリンク(i.imgur.com)
つまり、呪文前のTSHが高いほど呪文後はTSHは下がり、呪文前のTSHが低いほど呪文後のTSHは上がるという傾向がみてとれる。これを線形回帰して確かめてみる。
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 2.3812 0.3323 7.165 4.10e-09 ***
before -1.0351 0.1411 -7.338 2.24e-09 ***
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
Residual standard error: 0.8365 on 48 degrees of freedom
Multiple R-squared: 0.5287, Adjusted R-squared: 0.5189
F-statistic: 53.84 on 1 and 48 DF, p-value: 2.237e-09
p = 2.24e-09 なので有意差ありとしてよい。
ゆえに、「ド底辺シリツ医大が悪いのではない、本人の頭が悪いんだ」という呪文は甲状腺機能を正常化させる。
ここで問題、この統計処理のどこが誤っているか?

625:卵の名無しさん
18/05/05 10:54:16.04 WOhQQmHo.net
くくく!バカ丸出し!

626:卵の名無しさん
18/05/05 11:14:23.69 r+HJWwEoO
医科歯科君(今はウリュウなの?)、以前はもう少しまともな書き込みをしてたよね。医者っぽい発言も時々あったし。
今や完全に発症しちゃったね。病識ないだろうけど。
絶対医者じゃないね。まぁ全く信用されてないだろうけど。
本当に何者?逆に興味が出てきた。

627:卵の名無しさん
18/05/05 11:38:02.77 h9pFu5Lrh
必ず誰かに肘肩まわりをやさしくさすってもらいながら
一日3回 2週間
1テーブルや安定した椅子の前に脚を軽く開いて立ち、テーブルか椅子の背もたれに両手をつく。半歩下がり上半身の体重の半分を両手に預ける
2モゾモゾといいながら腰を左右前後に動かしてほぐしてゆるめていく。ひじは伸ばしたまま、全身の力を抜いて行う
肘は伸ばしたまま

628:卵の名無しさん
18/05/05 11:56:26.89 w175uM/J.net
>>622
ド底辺コルサコフ症候群は益々、悪化しているようだな。

629:卵の名無しさん
18/05/05 12:39:36.29 M4ESiG4T.net
地方国立卒の勤務医だが当直中だ。
貧乏に生まれると人が休んでる時でさえ
俺みたいに当直しないといけない。
医科歯科爺さんは底辺私立をクソミソに
貶してるけど、俺はどっちかというと羨ましいわ。
だって親がお金持ちじゃ無かったら
本来は医者になれなかった連中なんだぜ。
それが親のお陰だとは言っても医者やってるんだからさ。
受験の時に嬉々として私立医大に進学する同級生を
羨望の眼差しで見つめながら、諦めて旧帝理や旧帝工に
進学した同級生をたくさん知ってるよ。
あの時思ったな〜。
やっぱり金持ちにはかなわねえやと。
どう考えても旧帝の理や工に進学した
同級生の方が成績が良いんだから。
やっぱり人生はお金だよ。

630:卵の名無しさん
18/05/05 12:56:31.88 DcBecp9b.net
>>629
医師真似国試浪人ジイさん
キャラかえなくてもいいよ
しかしあんた、てか、ウリュウジジイ、つまんねえ奴だな

631:卵の名無しさん
18/05/05 13:21:40.97 WOhQQmHo.net
大当たりー カッコわりーオナ中だな。

632:卵の名無しさん
18/05/05 14:00:50.74 WqhhrHYt.net
カネで医者になる私立医は全員似非医者。
バカのまま入学し、
低レベルの国試に通るだけ。
こんなインチキ国家資格は医者だけだろが

633:卵の名無しさん
18/05/05 14:12:08.00 w175uM/J.net
>>629
接客業って賤業だと気付くいたときには遅しだったな。
俺は明日が日当直、月曜日は半日内視鏡バイトしてステーキハウスに行く予定。
自宅から呼ぶ出される方が当直室から呼び出されるより苦痛。

>>632
だからこそ、ド底辺シリツ医大卒は卒業大学をひた隠し。

634:卵の名無しさん
18/05/05 16:29:39.22 U4l9hcPD.net
内視鏡バイトの後にステーキハウス行く自慢する医師初めて見たわな
本物の医者だったらそんなバイト自慢しないし、ステーキハウスも自慢しなくね?

635:卵の名無しさん
18/05/05 16:31:50.51 U4l9hcPD.net
消化器内科か、、人の糞をほじくり出す賤業だな
その汚い手でステーキ食って幸せか?w

636:卵の名無しさん
18/05/05 18:27:14.11 pCZ7kmcH.net
>>635
午前中に上部のみ。

637:卵の名無しさん
18/05/05 18:30:24.93 pCZ7kmcH.net
>>634
某ド底辺医大生は定期試験後にフェリーで高松のソープに行くらしい。

638:卵の名無しさん
18/05/05 22:43:21.99 hBjvmFP3.net
ジイさん、ソープってなに?
初耳だよ、僕ら
ソープは医師真似事務員の仲間たちで流行ってるの?
最近はそんなソープ症例が国試にでるの?
ソープって、まさか長崎のめいぶつなの?

639:卵の名無しさん
18/05/05 23:24:40.36 2KIi+pnL.net
バカ私立医なら生涯、300人は患者殺してるだろ!
全員腹を切れ

640:卵の名無しさん
18/05/06 03:01:38.30 ounCZmr/.net
私立医大卒の医者って良いよな〜。
うらやましいよ。
だってお金が無かったら
本来は医者じゃ無い連中なんだぜ。

641:卵の名無しさん
18/05/06 07:38:40.97 ZBKvHzPJ.net
>>640
あんたも医者になれたら
そんな事は気にならなくなるよ
さあ、34回目の国試の為に勉強しなきゃ

642:卵の名無しさん
18/05/06 07:53:17.14 LrYAyHYX.net
ド底辺シリツ医大卒って英文や数字が出ると全くコメントできずに第三法則を実証するだけ。
馬鹿なのだろうか?

URLリンク(i.imgur.com)
統計は教養以前の問題。
それすら扱えないド底辺シリツ医ってなんと表現すればいいのだろう?
当事者の間では ガチの裏 と表現されてらしいw
URLリンク(imagizer.imageshack.com)

馬鹿なのにどうして大学に入れただろう?
やっぱり、あれ??

,二     |  ̄ヽ    . /〃   -┼─   
  )    |___.     (         レ--、 
 .ノ      __)      \.      _ノ .

643:卵の名無しさん
18/05/06 13:23:55.05 itQRaw6a.net
>>631
ド底辺シリツ医大卒よりカッコー悪いことは医療の世界には存在しない。
ド底辺シリツ医大卒の汚名より、ぶざまな事実は存在しない。
ちょっと、第二法則を確認してみましょう。
あんた、どこ卒?

644:卵の名無しさん
18/05/06 15:13:15.10 4gENZ/hC.net
順位和での傾向検定
jonckheere <- function(L){ # L : list of vectors A1,A2,...,Ak, with assumed tendency
How.Many.Greater.Pairs <- function(A,B){ # How many pairs of A[i] > B[j], count as 0.5 when equal,
n.a = length(A)
n.b = length(B)
how.many.greater.pairs = 0
for(i in 1:n.a){
for(j in 1:n.b){
how.many.greater.pairs = how.many.greater.pairs+ifelse(A[i]==B[j],0.5,A[i]>B[j])
}
}
return(how.many.greater.pairs)
}
Sum.of.Greater.Pairs <- function(L){ #L=list(A1,,,,Ak),A1 < A2 < A3,..,< Ak : vector
k = length(L)
comb = combn(1:k,2) # possible combinaition of pairs to compare
n.comb = ncol(comb) # how many combinations
J = 0 # sum of greater pairs
for(i in 1:n.comb){
J = J + How.Many.Greater.Pairs(L[[comb[1,i]]],L[[comb[2,i]]])
}
return(J)
}
J = Sum.of.Greater.Pairs(L)
n = sapply(L,length)
N = sum(n)
EJ = (N^2-sum(n^2))/4
VJ = (N^2*(2*N+3)-sum(n^2*(2*n+3)))/72
Z = (J-EJ)/sqrt(VJ)
p.value = pnorm(-abs(Z))*2
return(p.value)
}

645:卵の名無しさん
18/05/06 15:51:37.60 4gENZ/hC.net
パッケージの結果と一致して気持ちが(・∀・)イイ!!
> A1<-c(205,206,164,190,194,203)
> A2<-c(201,221,197,185)
> A3<-c(248,265,197,220,212,281)
> A4<-c(202,276,237,254,230)
> L<-list(A1,A2,A3,A4)
> g=rep(1:4,sapply(L,length))
> jonckheere(L)
[1] 0.004279366
> clinfun::jonckheere.test(c(A1,A2,A3,A4),g, alternative ='two.sided')
Jonckheere-Terpstra test
data:
JT = 127.5, p-value = 0.004279
alternative hypothesis: two.sided

646:卵の名無しさん
18/05/06 17:40:41.16 ToTSIDGj.net
>>640
だよね〜
本当に私立医大卒の医者って良いよな〜。
うらやましいよ。
だってお金が無かったら
本来は医者じゃ無い連中なんだよな〜。

647:卵の名無しさん
18/05/06 17:48:37.68 4gENZ/hC.net
Jonckheere's trend test
URLリンク(en.wikipedia.org)
での統計量が正規分布で近似できることをシミュレーションで確認してみる。
統計量○○は××分布に従う、とか○○の従う分布を××分布という、などの記載を読むとそんなものかと「納得」して次に進むことになるのだが、これがなんとも気持ち悪い。
まあ、ド底辺シリツ卒を裏口バカと呼ぶのはπ≒3.14程度のことで万人が受け入れることができるんだが(笑)
# simulation
How.Many.Greater.Pairs <- function(A,B){ # How many pairs of A[i] > B[j], count as 0.5 when equal.
n.a = length(A)
n.b = length(B)
how.many.greater.pairs = 0
for(i in 1:n.a){
for(j in 1:n.b){
how.many.greater.pairs = how.many.greater.pairs+ifelse(A[i]==B[j],0.5,A[i]>B[j])
}
}
return(how.many.greater.pairs)
}
Sum.of.Greater.Pairs <- function(L){ #L=list(A1,,,,Ak),A1 < A2 < A3,..,< Ak : vector
k = length(L)
comb = combn(1:k,2) # possible combinaition of pairs to compare
n.comb = ncol(comb) # how many combinations
J = 0 # sum of greater pairs
for(i in 1:n.comb){
J = J + How.Many.Greater.Pairs(L[[comb[1,i]]],L[[comb[2,i]]])
}
return(J)
}

648:卵の名無しさん
18/05/06 17:49:08.82 4gENZ/hC.net
repeat.jonckheere.simulation <- function(n.sim=10^3){
Max.score=100
Max.sample.size=10
N.group=sample(5:10,1)
N.sample=sample(3:Max.sample.size, N.group, replace=TRUE)
jonck.sim <- function(){
L=list()
for(i in 1:N.group) L=append(L,list(sample(1:Max.score,N.sample[i],replace=TRUE)))
Sum.of.Greater.Pairs(L)
}
rep.jon = replicate(n.sim,jonck.sim())
hist(rep.jon,freq=FALSE,col=sample(colours(),1),main='Gaussian on Jonckheere',xlab='statistic J',breaks=20)
n = N.sample
N = sum(n)
EJ = (N^2-sum(n^2))/4
VJ = (N^2*(2*N+3)-sum(n^2*(2*n+3)))/72
curve(dnorm(x,EJ,sqrt(VJ)), add=TRUE,lwd=2)
invisible(list(n=n,EJ=EJ,VJ=VJ,J=rep.jon))
}
re=repeat.jonckheere.simulation() ; names(re)
qqnorm(re$J) ; qqline(re$J)
shapiro.test(re$J)
ks.test(re$J, "pnorm", mean=re$EJ, sd=sqrt(re$VJ))

649:卵の名無しさん
18/05/06 17:53:49.64 4gENZ/hC.net
正規性を検定して、グラフ化してみる。

> shapiro.test(re$J)
Shapiro-Wilk normality test
data: re$J
W = 0.99875, p-value = 0.723
> ks.test(re$J, "pnorm", mean=re$EJ, sd=sqrt(re$VJ))
One-sample Kolmogorov-Smirnov test
data: re$J
D = 0.040292, p-value = 0.07779
alternative hypothesis: two-sided
URLリンク(i.imgur.com)


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

2182日前に更新/392 KB
担当:undef