18/12/28 13:38:51.13 4e3e1x1A.net
直線解ならプログラム組めそうだけど答が円弧を含むなら
俺の能力では無理だな。
こういう問題を解ける頭脳には平伏するだけだな。
5等分の場合はこんな形が解になるそうです
URLﾘﾝｸ(imgur.com)
正三角形5等分問題だと解が完全に非対称的な形になって面白い
URLﾘﾝｸ(imgur.com)
625:卵の名無しさん
18/12/28 15:09:23.16 4e3e1x1A.net
自分がどれだけ馬鹿なのか自覚できていいな。
三角形なら直線で左右対象だろうとアプリオリに決めつけてプログラムしちゃうよなぁ。
626:卵の名無しさん
18/12/28 15:10:49.35 4e3e1x1A.net
理系で一番馬鹿なのが医学部を実感するなぁ。
627:卵の名無しさん
18/12/29 12:03:59.78 8oggYPKD.net
分岐角度が120°、対象図形と分割線のなす角が90°という
いわゆるplatoの法則で立式しないとコンピュータ解も難しそう。
628:卵の名無しさん
18/12/29 12:25:15.44 8oggYPKD.net
>>721
４分割ですらの正答出せなかった。
URLﾘﾝｸ(i.imgur.com)
URLﾘﾝｸ(i.imgur.com)
URLﾘﾝｸ(i.imgur.com)
の順に小さくなるのは計算できたけど。
629:卵の名無しさん
2018/12/30
630:(日) 15:37:28.61 ID:NpJOUauM.net
631:卵の名無しさん
18/12/31 21:09:11.51 2cJ2vbGs.net
開業するのかw大晦日にw
632:卵の名無しさん
18/12/31 23:39:12.39 anCmupmN.net
数学板の統計スレにこんな課題があったのでシミュレーション解を投稿しておいた。
例えば、母集団から50枚の答案用紙を選んで、平均が
60点だとする。母集団は2500枚の答案用紙から成り立っているとして、平均は70点とする。このとき50枚の答案用紙をランダムに選んだかの検定はどの
ようにすればいいでしょうか。（つまり60点の平均点が低いので、ランダムに選んだかどうかを疑っているわけです。）
633:卵の名無しさん
19/01/01 10:02:25.14 fQ4brTV9.net
甲乙二人がおのおの32ピストル（当時のお金の単位）の金を賭けて勝負したとする。
そしてどちらかが先に３点を得たものを勝ちとし、勝った方がかけ金の総額64ピストルをもら
えるとする。ところが甲が２点、乙が１点を得たとき、勝負が中止になってしまった。
このとき、二人のかけ金の総額64ピストルを甲と乙にどのように分配すればよいだろうか。
ただし二人の力は互角で、勝つ確率はそれぞれ1/2ずつだとする。
プログラムで計算させてみる
n=0:1
x=expand.grid(n,n,n,n,n)
y=x[apply(x[,1:3],1,sum)==2,]
64*sum(apply(y,1,sum)>2)/nrow(y)
634:卵の名無しさん
19/01/01 21:19:05.09 nLWjX1/M.net
インフルエンザ脳症成人例は死亡率が高い
URLﾘﾝｸ(www.niid.go.jp)
p0_4=c(14,202)
p5_19=c(20,408)
p20_59=c(7,72)
p60_=c(10,66)
p.mat=as.matrix(rbind(p0_4,p5_19,p20_59,p60_))
chisq.test(p.mat)
pairwise.prop.test(p.mat,p.adjust='holm')
fmsb::pairwise.fisher.test(p.mat)
p.mat
prop.test(c(34,17),c(610,138))
Fisher.test(c(34,17),c(610,138)
635:卵の名無しさん
19/01/03 18:09:59.16 6cY0SNEs.net
## m sd n
A=c(159.0625,sqrt(3924.729167),16)
B=c(240,sqrt(22027.5),17)
C=c(366.35,sqrt(5329.292105),20)
lh=rbind(A,B,C)
colnames(lh)=c("m","sd","n") ; lh
mean.G=sum(lh[,"m"]*lh[,"n"])/sum(lh[,"n"])
SS.bit=sum((lh[,"m"]-mean.G)^2*lh[,"n"])
SS.wit=sum(lh[,"sd"]^2*(lh[,"n"]-1))
df.bit=nrow(lh)-1
df.wit=sum(lh[,"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))
A=c(159.0625,sqrt(3924.729167),16)
B=c(240,sqrt(22027.5),17)
C=c(366.35,sqrt(5329.292105),20)
a=A[1]+scale(rnorm(A[3]))*A[2]
b=B[1]+scale(rnorm(B[3]))*B[2]
c=C[1]+scale(rnorm(C[3]))*C[2]
x=c(a,b,c)
g=c(rep('A',A[3]),rep('B',B[3]),rep('C',C[3]))
g=factor(g,labels=LETTERS[1:3])
pairwise.t.test(x,g,p.adjust='holm')
636:卵の名無しさん
19/01/03 18:10:26.09 6cY0SNEs.net
>>598
URLﾘﾝｸ(i.imgur.com)
637:卵の名無しさん
19/01/03 22:28:02.75 pv3ynAMm.net
It is common knowledge among doctors and patients that Do-Teihen(exclusively bottom-leveled medical school) graduates mean morons who bought their way to Gachi'Ura(currently called by themselves)
According to the experience of entrance exam to medical school in the era of Showa, when the sense of discrimination against
privately-founded medical schools were more intense than it is now,
all such schools but for Keio had been so compared to some specialized institution for educable mentally retarded kids that nobody but imbecile successors of physicians in private practice had applied for admission.
There had been NOT a single classmate who chose willingly against his/her common sense to go to the Do-Teihen(exclusively bottom-leveled medical school, currently also known as Gachi'Ura),
which would have cost outrageous money and its graduates are destined to be called Uraguchi morons who bought thier way into the Do-Teihen, by thier colleagues and even by thier own clients.
Although people won't call them names to their face,
certain 80-90% people of about my age have been yet scorning and sneering at Uraguchi graduates, speaking in the back of our mind,
" Uraguchi morons shall not behave like somebody."
We never speak out face to face in real life.
638:卵の名無しさん
19/01/04 19:52:06.17 G2Rj2t04.net
It is common knowledge among doctors and patients that Do-Teihen(exclusively bottom-leveled medical school) graduates mean morons who bought their way to Gachi'Ura(currently called by themselves)
According to the experience of entrance exam to medical school in the era of Showa, when the sense of discrimination against
privately-founded medical schools were more intense than it is now,
all such schools but for Keio had been so compared to some specialized institution for educable mentally retarded kids that nobody but imbecile successors of physicians in private practice had applied for admission.
There had been NOT a single classmate who chose willingly against his/her common sense to go to the Do-Teihen(exclusively bottom-leveled medical school, currently also known as Gachi'Ura),
which would have cost outrageous money and its graduates are destined to be called Uraguchi morons who bought thier way into the Do-Teihen, by thier colleagues and even by thier own clients.
Although people won't call them names to their face,
certain 80-90% people of about my age have been yet scorning and sneering at Uraguchi graduates, speaking in the back of our mind,
" Uraguchi morons shall not behave like somebody."
We never speak out face to face in real life.
639:卵の名無しさん
19/01/05 00:22:46.55 X56IPw/0.net
It is common knowledge among doctors and patients that Do-Teihen(exclusively bottom-leveled medical school) graduates mean morons who bought their way to Gachi'Ura(currently called by themselves)
According to the experience of entrance exam to medical school in the era of Showa, when the sense of discrimination against
privately-founded medical schools were more intense than it is now,
all such schools but for Keio had been so compared to some specialized institution for educable mentally retarded kids that nobody but imbecile successors of physicians in private practice had applied for admission.
There had been NOT a single classmate who chose willingly against his/her common sense to go to the Do-Teihen(exclusively bottom-leveled medical school, currently also known as Gachi'Ura),
which would have cost outrageous money and its graduates are destined to be called Uraguchi morons who bought thier way into the Do-Teihen, by thier colleagues and even by thier own clients.
Although people won't call them names to their face,
certain 80-90% people of about my age have been yet scorning and sneering at Uraguchi graduates, speaking in the back of our mind,
" Uraguchi morons shall not behave like somebody."
We never speak out face to face in real life.
640:卵の名無しさん
19/01/05 13:06:20.59 1KoYjjkI
641:.net
642:卵の名無しさん
19/01/05 17:38:24.01 +iYQLowq.net
ド底辺シリツ医にもわかるように原文をつけておけよ。
>>
私は昭和の時代に大学受験したけど、昔は今よりも差別感が凄く、特殊民のための特殊学校というイメージで開業医のバカ息子以外は誰も受験しようとすらしなかった。
常識的に考えて、数千万という法外な金を払って、しかも同業者からも患者からもバカだの裏口だのと散々罵られるのをわかって好き好んでド底辺医に行く同級生は一人もいませんでした。
本人には面と向かっては言わないけれど、俺くらいの年代の人間は、おそらくは８－９割はド底辺医卒を今でも「何偉そうなこと抜かしてるんだ、この裏口バカが」と心の底で軽蔑し、嘲笑しているよ。
当の本人には面と向かっては絶対にそんなことは言わないけどね。
<<>>603
643:卵の名無しさん
19/01/05 19:11:13.14 71bgktVe.net
It is common knowledge among doctors and patients that Do-Teihen(exclusively bottom-leveled medical school) graduates mean morons who bought their way to Gachi'Ura(currently called by themselves)
According to the experience of entrance exam to medical school in the era of Showa, when the sense of discrimination against
privately-founded medical schools were more intense than it is now,
all such schools but for Keio had been so compared to some specialized institution for educable mentally retarded kids that nobody but imbecile successors of physicians in private practice had applied for admission.
There had been NOT a single classmate who chose willingly against his/her common sense to go to the Do-Teihen(exclusively bottom-leveled medical school, currently also known as Gachi'Ura),
which would have cost outrageous money and its graduates are destined to be called Uraguchi morons who bought thier way into the Do-Teihen, by thier colleagues and even by thier own clients.
Although people won't call them names to their face,
certain 80-90% people of about my age have been yet scorning and sneering at Uraguchi graduates, speaking in the back of our mind,
" Uraguchi morons shall not behave like somebody."
We never speak out face to face in real life.
644:卵の名無しさん
19/01/06 08:58:44.31 g7qyfOQm.net
It is common knowledge among doctors and patients that Do-Teihen(exclusively bottom-leveled medical school) graduates mean morons who bought their way to Gachi'Ura(currently called by themselves)
According to the experience of entrance exam to medical school in the era of Showa, when the sense of discrimination against
privately-founded medical schools were more intense than it is now,
all such schools but for Keio had been so compared to some specialized institution for educable mentally retarded kids that nobody but imbecile successors of physicians in private practice had applied for admission.
There had been NOT a single classmate who chose willingly against his/her common sense to go to the Do-Teihen(exclusively bottom-leveled medical school, currently also known as Gachi'Ura),
which would have cost outrageous money and its graduates are destined to be called Uraguchi morons who bought thier way into the Do-Teihen, by thier colleagues and even by thier own clients.
Although people won't call them names to their face,
certain 80-90% people of about my age have been yet scorning and sneering at Uraguchi graduates, speaking in the back of our mind,
" Uraguchi morons shall not behave like somebody."
We never speak out face to face in real life.
645:卵の名無しさん
19/01/06 14:48:30.19 SEvsBMAu.net
Zが標準正規分布に従う時、次の値を求めよ。
１：P｛Z＞u1｝＝0.05を満たすu1の値 → u1=1.645
２：P｛Z＜u2｝＝0.005を満たすu2の値 → u2=-2.576
３：P｛-u3＜Z＜u3｝＝0.99を満たすu3の値 → u3=2.576
(u1=qnorm(0.05,lower.tail = FALSE))
(u2=qnorm(0.005,lower.tail= TRUE))
(u3=qnorm((1-0.99)/2,lower.tail=FALSE))
par(mfrow=c(3,1))
curve(dnorm(x),-4,4,bty='l')
curve(dnorm(x),u1,4,type='h',add=T,col=2)
curve(dnorm(x),-4,4,bty='l')
curve(dnorm(x),-4,u2,type='h',add=T,col=3)
curve(dnorm(x),-4,4,bty='l')
curve(dnorm(x),-u3,u3,type='h',add=T,col=4)
Xが正規分布N（10, 5^2）に従う時、次の確率を求めよ。
１：P｛X＞20｝ → 0.0228
２：P｛X＜5｝ → 0.1587
３：P｛0＜X＜20｝ → 0.9544
1-pnorm(20,10,5)
pnorm(5,10,5)
pnorm(20,10,5)-pnorm(0,10,5)
646:卵の名無しさん
19/01/06 14:48:45.04 SEvsBMAu.net
>>606
ド底辺シリツ医にもわかるように原文をつけておけよ。
>>
私は昭和の時代に大学受験したけど、昔は今よりも差別感が凄く、特殊民のための特殊学校というイメージで開業医のバカ息子以外は誰も受験しようとすらしなかった。
常識的に考えて、数千万という法外な金を払って、しかも同業者からも患者からもバカだの裏口だのと散々罵られるのをわかって好き好んでド底辺医に行く同級生は一人もいませんでした。
本人には面と向かっては言わないけれど、俺くらいの年代の人間は、おそらくは８－９割はド底辺医卒を今でも「何偉そうなこと抜かしてるんだ、この裏口バカが」と心の底で軽蔑し、嘲笑しているよ。
当の本人には面と向かっては絶対にそんなことは言わないけどね。
647:卵の名無しさん
19/01/06 23:10:04.43 PSVRj6pz.net
It is common knowledge among doctors and patients that Do-Teihen(exclusively bottom-leveled medical school) graduates mean morons who bought their way to Gachi'Ura(currently called by themselves)
According to the experience of entrance exam to medical school in the era of Showa, when the sense of discrimination against
privately-founded medical schools were more intense than it is now,
all such schools but for Keio had been so compared to some specialized institution for educable mentally retarded kids that nobody but imbecile successors of physicians in private practice had applied for admission.
There had been NOT a single classmate who chose willingly against his/her common sense to go to the Do-Teihen(exclusively bottom-leveled medical school, currently also known as Gachi'Ura),
which would have cost outrageous money and its graduates are destined to be called Uraguchi morons who bought thier way into the Do-Teihen, by thier colleagues and even by thier own clients.
Although people won't call them names to their face,
certain 80-90% people of about my age have been yet scorning and sneering at Uraguchi graduates, speaking in the back of our mind,
" Uraguchi morons shall not behave like somebody."
We never speak out face to face in real life.
648:卵の名無しさん
19/01/07 15:11:59.69 9K2vk7ik.net
なんだか疲れた
働きたくない
649:卵の名無しさん
19/01/07 15:25:03.66 31hT7o17.net
URLﾘﾝｸ(blogs.wsj.com)
URLﾘﾝｸ(www.stats.org)<)
URLﾘﾝｸ(politicalarithmetik.blogspot.com)<)
URLﾘﾝｸ(serc.carleton.edu)
URLﾘﾝｸ(www.statlit.org)
650:卵の名無しさん
19/01/07 18:20:24.95 ubVk9gDC.net
戸田Ａ哲君もアメリカで経済学助狂をやっているらしい。
バブル崩壊で安定を求めて理３に入ってしまったらしい。
無駄に齢とっちまったな。
651:卵の名無しさん
19/01/07 20:32:18.68 31hT7o17.net
>>610
正月休みで１週間ぶりに内視鏡したが違和感があるなぁ。
平日は週３やってる。
652:卵の名無しさん
19/01/07 21:02:56.59 31hT7o17.net
URLﾘﾝｸ(www.ncbi.nlm.nih.gov)
653:卵の名無しさん
19/01/07 21:09:46.14 9K2vk7ik.net
なんだかだるい
あたりが暗く見える
笑い声を聞くとムカつく
疲れた
働きたくない
654:卵の名無しさん
19/01/07 21:22:19.83 /APKkASA.net
It is common knowledge among doctors and patients that Do-Teihen(exclusively bottom-leveled medical school) graduates mean morons who bought their way to Gachi'Ura(currently called by themselves)
According to the experience of entrance exam to medical school in the era of Showa, when the sense of discrimination against
privately-founded medical schools were more intense than it is now,
all such schools but for Keio had been so compared to some specialized institution for educable mentally retarded kids that nobody but imbecile successors of physicians in private practice had applied for admission.
There had been NOT a single classmate who chose willingly against his/her common sense to go to the Do-Teihen(exclusively bottom-leveled medical school, currently also known as Gachi'Ura),
which would have cost outrageous money and its graduates are destined to be called Uraguchi morons who bought thier way into the Do-Teihen, by thier colleagues and even by thier own clients.
Although people won't call them names to their face,
certain 80-90% people of about my age have been yet scorning and sneering at Uraguchi graduates, speaking in the back of our mind,
" Uraguchi morons shall not behave like somebody."
We never speak out face to face in real life.
655:卵の名無しさん
19/01/07 21:26:07.02 31hT7o17.net
移民法でお先真っ暗だよなぁ。
656:卵の名無しさん
19/01/07 23:44:52.48 fKvhbwc+.net
窓の外で誰かが見ている
657:卵の名無しさん
19/01/08 11:14:42.05 Fx7hufSt.net
お天道様がみている、というのが日本人の道徳。
移民法でお先真っ暗。
658:卵の名無しさん
19/01/08 12:13:16.36 C0NmoK5j.net
電波が、電波の声が聞こえる
659:卵の名無しさん
19/01/08 12:44:58.49 o3Ur7/E6.net
>>620
便利でいいな。
オービスも聞こえればスピード違反で捕まることもないよな。
660:卵の名無しさん
19/01/08 13:32:28.48 o3Ur7/E6.net
rm(list=ls())
x="1-a,2-b"
y=strsplit(x,",")
z=unlist(y)
w=NULL
for(i in 1:length(z)){
w=rbind(w,unlist(strsplit(z[i],"-")))
}
data.frame(NUM=as.numeric(w[,1]),TEXT=w[,2])
661:卵の名無しさん
662:
ダルい 疲れた 働いたら負け
663:卵の名無しさん
19/01/08 14:31:07.32 o3Ur7/E6.net
>>623
正月休みで１週間内視鏡から離れたら却って体調が悪いなぁ。
金のために不愉快な仕事をしないというのが最大幸福だな。
おかげで自宅で安眠できる。
664:卵の名無しさん
19/01/08 17:02:04.26 C0NmoK5j.net
>>624
これ即答してみ！
できなきゃ、ド底辺学力と認定してあげよう。
フタナリとガチユリの相関性について
原稿用紙30枚以上で考察せよ
665:卵の名無しさん
19/01/08 17:30:33.63 5YyNsaaz.net
It is common knowledge among doctors and patients that Do-Teihen(exclusively bottom-leveled medical school) graduates mean morons who bought their way to Gachi'Ura(currently called by themselves)
According to the experience of entrance exam to medical school in the era of Showa, when the sense of discrimination against
privately-founded medical schools were more intense than it is now,
all such schools but for Keio had been so compared to some specialized institution for educable mentally retarded kids that nobody but imbecile successors of physicians in private practice had applied for admission.
There had been NOT a single classmate who chose willingly against his/her common sense to go to the Do-Teihen(exclusively bottom-leveled medical school, currently also known as Gachi'Ura),
which would have cost outrageous money and its graduates are destined to be called Uraguchi morons who bought thier way into the Do-Teihen, by thier colleagues and even by thier own clients.
Although people won't call them names to their face,
certain 80-90% people of about my age have been yet scorning and sneering at Uraguchi graduates, speaking in the back of our mind,
" Uraguchi morons shall not behave like somebody."
We never speak out face to face in real life.
666:卵の名無しさん
19/01/08 20:39:41.24 j/PKfHLe.net
>>625
原稿用紙30枚以上と書くあたり
部屋割カウントできない馬鹿とすぐ分かるぞ
667:卵の名無しさん
19/01/08 20:48:31.02 j/PKfHLe.net
こういう問題の方が楽しいね！
12345→51234のように、十進表記の末尾一桁を先頭にもってくる操作を考える。
ただし末尾は0ではないとし、全桁が同じ数字からなる数も考えないものとする。
ある数Mをこの操作でLへと置き換えた。
MがLで割り切れるMのうち、桁数が最も小さいものを全て求めよ。
668:卵の名無しさん
19/01/08 21:50:24.74 xj6Kc+5j.net
>>628
ナゾナゾは出せど
他人からの問題には全く回答不能なんだね
中卒ニートのコンビニバイトくん
669:卵の名無しさん
19/01/08 21:57:42.66 j/PKfHLe.net
>>629
数学板ででた問題、面白いだろ。
実用性は？だけど。
こっちは実用的。
Rのスレに質問あったのでレスしたけど
いつも達人がエレガントなコードを出してくる。
他人のコードを読むのは勉強になる。
URLﾘﾝｸ(rio2016.2ch.net)
670:卵の名無しさん
19/01/08 22:00:36.94 j/PKfHLe.net
>>629
では模範解答をお願いしますね。
原稿用紙３０枚以上必要らしいけど。
はよ、はよ！
671:卵の名無しさん
19/01/08 22:00:48.42 xj6Kc+5j.net
>>630
これ即答してみ！
できなきゃ、ド底辺学力と認定してあげよう。
医師板で算数クイズ出して荒らしまわってて、
いい年してそんなものが何故必要かと聞かれ、
ガソリンスタンドで客からの注文に応えるため必要だ、と答えたド底辺学力のニートくんは、
なぜコンビニバイトに転職したのでしょう？
672:卵の名無しさん
19/01/08 22:14:44.26 j/PKfHLe.net
>>632
俺にはコンビニバイトは無理だな。
内視鏡や麻酔の方がずっと楽。
673:卵の名無しさん
19/01/08 22:16:41.70 j/PKfHLe.net
>>632
算数クイズができると論文の検証ができるわけね。
内視鏡の洗浄待機時間に 胃と腸 の所見用語集２０１７を読んでいたらこんな記載があった。(602頁Fig4の説明）
white zoneの視認性と腺管密度の関係
white zoneが鮮明な癌 窩間部の長さ136±45 (p<0.05)
white zoneが不鮮明な癌 窩間部の長さ77±22 (p<0.05)
サンプルサイズが記載なしでp<0.05としか書いてない。
p<0.01ならそう記載するだろうから、0.01<p<0.05なのだろう。
±のあとの数字は標準偏差なのか標準誤差なのか判然としない。
検定はｔ検定（等分散を前提とはできないのでWelch法を使用したとする）で
サンプルサイズは等しいとして±の後の数字は標準偏差あるとすると
> n2p(5:9)
p.value p.value p.value p.value p.value
0.055446074 0.032947710 0.019850610 0.012077956 0.007403284
サンプルサイズは僅か６～８と推定される。
674:卵の名無しさん
19/01/08 22:20:06.62 j/PKfHLe.net
これ答えられるかい？
コンビニバイト
675:以下の教養ならまあ無理だな。 >643の±の後の数字は標準誤差とすると サンプルサイズはいくつと推定できるか？
676:卵の名無しさん
19/01/08 22:26:48.77 j1+ZKzaK.net
コンビニバイトくんのレスはこれだな
ｽﾚﾘﾝｸ(hosp板:63番)
677:卵の名無しさん
19/01/08 22:33:11.95 xj6Kc+5j.net
>>636
それそれ
医師はそんなレスしないだろ
678:卵の名無しさん
19/01/08 22:34:22.63 j/PKfHLe.net
>>636
そういう計算もできない馬鹿がこういう事故を起こすわけなんだな。
URLﾘﾝｸ(i.imgur.com)
統計スレなんだから>635にでも答えてみ！
できなきゃガソリンスタンドのバイト以下の頭脳な。
679:卵の名無しさん
19/01/08 22:36:19.30 j/PKfHLe.net
>>637
ガソリンスタンドのレスは理解できたらしいが
>634のレスは理解できんだろ？ド底辺頭脳クン。
680:卵の名無しさん
19/01/08 23:00:15.32 C0NmoK5j.net
>>638
>統計スレなんだから
開業医でもないのに開業医スレに書き込むコンビニバイトくんに
特大ブーメランW
681:卵の名無しさん
19/01/09 07:21:21.07 7x5Uv3Et.net
>>640
標準誤差の場合のサンプルサイズ推定の
Rのコードはこのスレに既出なんだけどね。
開業医スレに投稿したら礼を言われましたぜ！
URLﾘﾝｸ(egg.2ch.net)
682:卵の名無しさん
19/01/09 07:23:46.55 7x5Uv3Et.net
ちゃんと主張の正しさにも同意を得たぞ。
URLﾘﾝｸ(egg.2ch.net)
683:卵の名無しさん
19/01/09 08:32:19.46 e7U7282c.net
>>641
コンビニバイトくんのレスは
何の釈明にもなっていない
>632への回答 はよはよ
684:卵の名無しさん
19/01/09 09:31:40.92 q6OvMTx7.net
>>643
礼を言われるような投稿がだよ～ん。
あんたのは馬鹿を晒しているだけ。
ガソリンスタンドのレスは理解できたらしいが
>634のレスは理解できんだろ？ド底辺頭脳クン。
685:卵の名無しさん
19/01/09 13:44:05.50 e7U7282c.net
>>644
ところでコンビニバイトくん
非医師のアラシという自覚ある？
ガソリンスタンドのレスで理解できるのは
オマエがガソリンスタンドで給油の仕事をしているということのみだな
686:卵の名無しさん
19/01/09 14:24:39.06 7x5Uv3Et.net
>>645
悪いが非医師じゃないんだね。
非医師にこういう議論ができると思う？
URLﾘﾝｸ(egg.2ch.net)
こう言う方が非医師のアラシっぽいぞ。
URLﾘﾝｸ(egg.2ch.net)
687:卵の名無しさん
19/01/09 14:25:57.19 7x5Uv3Et.net
>>645
標準誤差のときのサンプルサイズの計算まだぁ？
ひょっとして、標準誤差が何かも知らない馬鹿ぁ？
688:卵の名無しさん
19/01/09 14:28:10.74 7x5Uv3Et.net
>>645
んで、あんたド底辺シリツ医大卒？
頭脳レベルはド底辺みたいだけど。
689:卵の名無しさん
19/01/09 14:52:34.20 7x5Uv3Et.net
URLﾘﾝｸ(fpr-calc.ucl.ac.uk)
の比率版を非心Χ二乗分布のRスクリプトをメールしたら返事が来たけどunderlining mathの議論が進まなかったな。
Intuitivel Biostatisticsの著者は即効で対応してErrataの追加と次版の修正の返事が来た。
690:卵の名無しさん
19/01/09 16:26:11.62 q6OvMTx7.net
>>645
どアホにもわかる計算の例示として税込みガソリン価格を上げただけということすらわからんらしいな。
こういう問題に算数が必要と書いてもオマエの頭じゃ正解だせないだろ？
診療所から病院に患者を救急搬送する。
病院から医師搭乗の救急車が診療所に向かっており１０時到着予定と連絡が入った。
患者の病態が悪化したら、診療所の普通車で病院に向かい救急車と出会ったら
救急車に患者を移して搬送し病院到着を急ぐという計画を立てた。
普通車から救急車への患者の乗り換えで１０分余分に時間がかかる。
道路事情から救急車は病院から診療所への道は
平均時速６０ｋｍで、逆方向は平均時速４５ｋｍで定速走行する。
691:診療所の普通車は信号待ちもあり平均時速３０kmで定速走行する。 何時以降の病態悪化は診療所の車を使わずに救急車の到着を待つ方が病院に早く着くか？
692:卵の名無しさん
19/01/09 16:42:16.13 PAq6vTKi.net
>>650
非正規雇用のオッさんの言い訳キモス
反面教師にしかならんわね
もうすぐ冬休みも終わって
再来年は高校受験だわ～
693:卵の名無しさん
19/01/09 18:36:24.74 7x5Uv3Et.net
>>651
答まだぁ？
阿呆じゃね？？
694:卵の名無しさん
19/01/09 19:44:34.23 e7U7282c.net
中学生にすらバカにされるコンビニバイトくんW
695:卵の名無しさん
19/01/10 14:24:50.58 9bNrFpHw.net
>>653
その頭脳すらなくて答られないのがオマエだぞ。
696:卵の名無しさん
19/01/10 14:25:28.63 9bNrFpHw.net
>>653
>>645
どアホにもわかる計算の例示として税込みガソリン価格を上げただけということすらわからんらしいな。
こういう問題に算数が必要と書いてもオマエの頭じゃ正解だせないだろ？
診療所から病院に患者を救急搬送する。
病院から医師搭乗の救急車が診療所に向かっており１０時到着予定と連絡が入った。
患者の病態が悪化したら、診療所の普通車で病院に向かい救急車と出会ったら
救急車に患者を移して搬送し病院到着を急ぐという計画を立てた。
普通車から救急車への患者の乗り換えで１０分余分に時間がかかる。
道路事情から救急車は病院から診療所への道は
平均時速６０ｋｍで、逆方向は平均時速４５ｋｍで定速走行する。診療所の普通車は信号待ちもあり平均時速３０kmで定速走行する。
何時以降の病態悪化は診療所の車を使わずに救急車の到着を待つ方が病院に早く着くか？
697:卵の名無しさん
19/01/10 14:26:28.27 9bNrFpHw.net
>>653
>>645
悪いが非医師じゃないんだね。
非医師にこういう議論ができると思う？
URLﾘﾝｸ(egg.2ch.net)
こう言う方が非医師のアラシっぽいぞ。
URLﾘﾝｸ(egg.2ch.net)
698:卵の名無しさん
19/01/10 14:44:10.09 9bNrFpHw.net
バレット食道がんって右にあるよなぁ。
URLﾘﾝｸ(www.hindawi.com)
699:卵の名無しさん
19/01/10 19:43:55.09 xCiYE3NF.net
>>654～657
これ即答してみ！
できなきゃ、ド底辺学力と認定してあげよう。
医師板で算数クイズ出して荒らしまわってて、
いい年してそんなものが何故必要かと聞かれ、
ガソリンスタンドで客からの注文に応えるため必要だ、と答えたド底辺学力のニートくんは、
なぜコンビニバイトに転職したのでしょう？
>>654
送り仮名がオカシイところをみると
日本の方ですらないのですかね
700:卵の名無しさん
19/01/10 20:04:49.50 9bNrFpHw.net
>>653
>>645
悪いが非医師じゃないんだね。
非医師にこういう議論ができると思う？
URLﾘﾝｸ(egg.2ch.net)
こう言う方が非医師のアラシっぽいぞ。
URLﾘﾝｸ(egg.2ch.net)
701:卵の名無しさん
19/01/11 07:03:12.10 eIjYcl4Y.net
>>659
非正規雇用のオッさんの言い訳キモス
反面教師にしかならんわね
もうすぐ冬休みも終わって
再来年は高校受験だわ～
702:卵の名無しさん
19/01/11 07:05:30.85 SA5KNNkQ.net
>>645
どアホにもわかる計算の例示として税込みガソリン価格を上げただけということすらわからんらしいな。
こういう問題に算数が必要と書いてもオマエの頭じゃ正解だせないだろ？
診療所から病院に患者を救急搬送する。
病院から医師搭乗の救急車が診療所に向かっており１０時到着予定と連絡が入った。
患者の病態が悪化したら、診療所の普通車で病院に向かい救急車と出会ったら
救急車に患者を移して搬送し病院到着を急ぐという計画を立てた。
普通車から救急車への患者の乗り換えで１０分余分に時間がかかる。
道路事情から救急車は病院から診療所への道は
平均時速６０ｋｍで、逆方向は平均時速４５ｋｍで定速走行する。診療所の普通車は信号待ちもあり平均時速３０kmで定速走行する。
何時以降の病態悪化は診療所の車を使わずに救急車の到着を待つ方が病院に早く着くか？
703:卵の名無しさん
19/01/11 08:08:05.29 wEewwlyz.net
>>661
>税込みガソリン価格を上げただけ ？W
違うよねぇ
ｽﾚﾘﾝｸ(hosp板:62番)-63
>税込みで5000円分給油してと言われて対応できないと仕事にならんだろ？
巧妙な話のすり替えワロス
医師は給油の仕事はしないよね
コンビニバイトくん
704:卵の名無しさん
19/01/11 08:59:52.18 SA5KNNkQ.net
>>645
どアホにもわかる計算の例示として税込みガソリン価格を上げただけということすらわからんらしいな。
こういう問題に算数が必要と書いてもオマエの頭じゃ正解だせないだろ？
診療所から病院に患者を救急搬送する。
病院から医師搭乗の救急車が診療所に向かっており１０時到着予定と連絡が入った。
患者の病態が悪化したら、診療所の普通車で病院に向かい救急車と出会ったら
救急車に患者を移して搬送し病院到着を急ぐという計画を立てた。
普通車から救急車への患者の乗り換えで１０分余分に時間がかかる。
道路事情から救急車は病院から診療所への道は
平均時速６０ｋｍで、逆方向は平均時速４５ｋｍで定速走行する。診療所の普通車は信号待ちもあり平均時速３０kmで定速走行する。
何時以降の病態悪化は診療所の車を使わずに救急車の到着を待つ方が病院に早く着くか？
705:卵の名無しさん
19/01/11 09:01:14.48 SA5KNNkQ.net
数学板の統計スレで助言したら礼を言われた。
URLﾘﾝｸ(rio2016.2ch.net)
706:卵の名無しさん
19/01/11 09:22:06.72 SA5KNNkQ.net
>>662
そのスレの続きみたぁ？
ガソリンスタンドバイトでないことは普通の頭脳があれば理解できてるぞ。
ｽﾚﾘﾝｸ(hosp板:64番)-65
ｽﾚﾘﾝｸ(hosp板:96番)
707:卵の名無しさん
19/01/11 12:10:41.16 2Jz2b/ri.net
いい年してバイトのオッさん
おおキモイキモイ
708:卵の名無しさん
19/01/11 13:28:05.64 wEewwlyz.net
>>663-665
>税込みガソリン価格を上げただけ ？W
違うよねぇ
ｽﾚﾘﾝｸ(hosp板:62番)-63
>税込みで5000円分給油してと言われて対応できないと仕事にならんだろ？
稚拙な話のすり替えワロス
医師は給油の仕事はしないよね
コンビニバイトくん
スレを辿ると
受験は遥か昔の医師の人生にとって一次方程式が何の役に立つのか？
というレスにオマエが
給油の仕事で対応できない
と答えたんだよね
普通に見れば、ナリスマシ非医師のボロが出た、としか考えられんよねW
709:卵の名無しさん
19/01/11 13:49:59.59 SA5KNNkQ.net
>>667
一次方程式は役に立つよ。
これも二次方程式は必要ないぞ。
答えてみ！
>>645
どアホにもわかる計算の例示として税込みガソリン価格を上げただけということすらわからんらしいな。
こういう問題に算数が必要と書いてもオマエの頭じゃ正解だせないだろ？
診療所から病院に患者を救急搬送する。
病院から医師搭乗の救急車が診療所に向かっており１０時到着予定と連絡が入った。
患者の病態が悪化したら、診療所の普通車で病院に向かい救急車と出会ったら
救急車に患者を移して搬送し病院到着を急ぐという計画を立てた。
普通車から救急車への患者の乗り換えで１０分余分に時間がかかる。
道路事情から救急車は病院から診療所への道は
平均時速６０ｋｍで、逆方向は平均時速４５ｋｍで定速走行する。診療所の普通車は信号待ちもあり平均時速３０kmで定速走行する。
何時以降の病態悪化は診療所の車を使わずに救急車の到着を待つ方が病院に早く着くか？
710:卵の名無しさん
19/01/11 13:51:40.72 SA5KNNkQ.net
>>667
そのスレの続きみたぁ？
シンプソンのパラドックスに感心したらしい開業医のレス読めないの？
ｽﾚﾘﾝｸ(hosp板:64番)-65
ｽﾚﾘﾝｸ(hosp板:96番)
711:卵の名無しさん
19/01/11 14:41:37.54 wEewwlyz.net
>668
>税込みガソリン価格を上げただけ ？W
違うよねぇ
ｽﾚﾘﾝｸ(hosp板:62番)-63
>税込みで5000円分給油してと言われて対応できないと仕事にならんだろ？
稚拙な話のすり替えワロス
医師は給油の仕事はしないよね
コンビニバイトくん
スレを辿ると
受験は遥か昔の医師の今後の人生にとって算数が何の役に立つのか？
というレスにオマエが
給油の仕事で対応できない
と答えたんだよね
普通に見れば、ナリスマシ非医師のボロが出た、としか考えられんよねW
712:卵の名無しさん
19/01/11 15:58:22.65 JAXb750V.net
# calculate point FPR by David Colquhoun
calc.FPR = function(pval=0.0475,nsamp=16,prior=0.5,sigma=1,delta1=1){
sdiff=sqrt(sigma^2/nsamp + sigma^2/nsamp)
df=2*(nsamp-1)
# Note FPR doesn't need calculation of power for p-equals case #
# under H0, use central t distribution
tcrit=-qt(pval/2,df,ncp=0)
x0=tcrit
y0=dt(x0,df,0)
#
# under H1 use non-central t distribution
ncp1=delta1/sdiff #non-centrality paramater
curve(dt(x,df),-3-ncp1,3+ncp1, ann=FALSE,bty='n')
curve(dt(x,df,ncp1),add=TRUE,lty=2)
abline(v=x0,lty=3)
x1=x0 #tcrit
y1=dt(x1,df,ncp=ncp1)
#
# Calc false positive rate
p0=2*y0 # p.value : two-sided verified with simulation
p1=y1
LRH1=p1/p0
FPR=((1-prior)*p0)/(((1-prior)*p0) + prior*p1)
FPR
output=c(FPR,LRH1)
return(FPR)
}
# end of function calc.FPR
calc.FPR()
713:卵の名無しさん
19/01/11 15:58:53.36 JAXb750V.net
>>670
そのスレの続きみたぁ？
シンプソンのパラドックスに感心したらしい開業医のレス読めないの？
ｽﾚﾘﾝｸ(hosp板:64番)-65
ｽﾚﾘﾝｸ(hosp板:96番)
714:卵の名無しさん
19/01/11 16:08:40.06 JAXb750V.net
内視鏡の洗浄待機時間に 胃と腸 の所見用語集２０１７を読んでいたらこんな記載があった。(602頁Fig4の説明）
white zoneの視認性と腺管密度の関係
white zoneが鮮明な癌 窩間部の長さ136±45 (p<0.05)
white zoneが不鮮明な癌 窩間部の長さ77±22 (p<0.05)
サンプルサイズが記載なしでp<0.05としか書いてない。p<0.01ならそう記載するだろうから、0.01<p<0.05なのだろう。
±のあとの数字は標準偏差なのか標準誤差なのか判然としない。
検定はｔ検定（等分散を前提とはできないのでWelch法を使用したとする）で
サンプルサイズは等しいとして±の後の数字は標準偏差あるとすると
> n2p(5:9)
p.value p.value p.value p.value p.value
0.055446074 0.032947710 0.019850610 0.012077956 0.007403284
サンプルサイズは僅か６～８と推定される。
n=7
p.value=0.019850610
としてFalse Positive Riskを計算してみる（非心ｔ分布でのコードは既出）
# calculate point FPR by David Colquhoun
を実行してみる。
> calc.FPR(0.01985,nsamp=7)
[1] 0.1245345
Web caluculator
URLﾘﾝｸ(fpr-calc.ucl.ac.uk)
の結果と一致して、気分が( ・∀・)ｲｲ!!
715:卵の名無しさん
19/01/11 16:23:30.11 JAXb750V.net
Colquhounの方法って
尤度比を確率密度比で算出してんだよなあ
reverse Bayesとか呼んでるけど。
式で書けばこうなるはず、
# FPR=1/(1+L10*PH1/(1-PH1)) where L10 likelihood ratio = P(Data|H1)/P(Data|H0) ,PH1: prior of H1
4.
The point null
Some people don’t like the assumption of a point null that’s made in this
proposed
approach to calculating the false positive risk
, but it seems natural to users who wish
to know how likely their observations would be if there were really no effect (i.e. if the
point null were true). It’s often claimed that the null is never exactly true. This isn’t
necessarily so (just giv
e identical treatments to both groups). But more importantly, it
doesn’t matter. We aren’t saying that the effect size is exactly zero: that would
obviously be impossible. We are looking at the likelihood ratio, i.e. at the probability
of the data if the t
rue effect size were not zero relative to the probability of the data if
the effect size were zero. If the latter probability is bigger than the former then
clearly we can’t be sure that there is anything but chance at work. This does
not
mean that we ar
e saying that the effect size is exactly zero. From the point of view
of the experimenter, testing the point null makes total sense. In any case, it makes
little difference if the null hypothesis is a narrow band centred on zero
(Berger and Delampady 1987)
.
716:卵の名無しさん
19/01/11 16:26:49.23 JAXb750V.net
>>674
ここにレクチャーが
URLﾘﾝｸ(www.youtube.com)
あるんだけど、
非心ｔ分布にはなんにも言及がないんだよなぁ。
非心分布の解説本ってこれを凌駕するのはないな。
サンプルサイズの決め方 (統計ライブラリー)
GPowerとかのソフト使ってクリック猿になるのは楽しくない。 <
717:卵の名無しさん
19/01/11 16:35:41.41 JAXb750V.net
比率のFalse Positive Risk算出のスクリプトを書いてみた。
calc.FPR.prop <- function(r1,r2,n1,n2,alpha=0.05,Print=TRUE,Show.chi=TRUE){
p.val=prop.test(c(r1,r2),c(n1,n2))$p.value
k=length(c(r1,r2)) # 2
df=k-1 # 1
p1=r1/n1 ; p2=r2/n2
PO=c(p1,p2) # proporition observed
# approximation with arcsine -- deprecated
# theta.i=asin(sqrt(PO)) # arcsine conversion
# delta=4*var(theta.i)*(k-1) # sum of squares
PE=(r1+r2)/(n1+n2) # PE: proportion expected
delta=sum((PO-PE)^2/PE)
N=2/(1/n1+1/n2) # harmonic mean to fit when n1 is not equal to n2
ncp1=N*delta # non-central parameter
power.chi=pchisq(qchisq(1-alpha,df),df,ncp1,lower=FALSE)
q.alpha=qchisq(1-alpha,df) # 3.841 when alpha=0.05
qcrit=qchisq(max(p.val,1-p.val),df,ncp=0) # qcrit = prop.test(c(r1,r2),c(n1,n2))$statistic
y0=dchisq(qcrit,df,0)
y1=dchisq(qcrit,df,ncp=ncp1)
FPR.equal=y0/(y0+y1) # length ratio = likelyhood ratio
if(Print){
curve(dchisq(x,df),0.5,3*qcrit,ylim=c(0,0.3),xlab=expression(paste(chi,'2')),ylab='density',bty='n',lwd=2,type='n') # H0
curve(dchisq(x,df,ncp1),qcrit,3*qcrit,col='lightgreen',type='h',add=TRUE,lwd=2)
curve(dchisq(x,df),qcrit,3*qcrit,col='orange',type='h',add=TRUE)
curve(dchisq(x,df),lwd=2,add=TRUE) # H0
curve(dchisq(x,df,ncp1),add=TRUE,lty=2,lwd=2) # H1
abline(v=qcrit)
abline(v=q.alpha,col='gray',lty=3)
718:卵の名無しさん
19/01/11 16:36:13.06 JAXb750V.net
legend('topright',bty='n',
legend=c('H0','H1',expression(paste(chi,'2:p.value')),expression(paste(chi,'2:alpha'))),
col=c('orange','lightgreen',1,'gray'),lty=c(1,2,1,3),lwd=c(2,2,1,1))
text(qcrit,y0,'y0',cex=0.75)
text(qcrit,y1,'y1',cex=0.75)
if(Show.chi){
text(qcrit,0,round(qcrit,2),cex=0.75)
text(q.alpha,0.30,round(q.alpha,2),cex=0.75)
}
}
ES=pwr::ES.h(r1/n1,r2/n2)
my.power=pwr::pwr.2p2n.test(ES,n1,n2,sig.level=p.val)$power
FPR.less=p.val/(p.val+my.power) # area ratio for 'less than' FDR
FPR.alpha=alpha/(alpha+power.chi) # FPR at alpha
output=c(power=power.chi,p.value=p.val,FPR.equal=FPR.equal,LH10=y1/y0,
FPR.less=FPR.less,FPR.alpha=FPR.alpha)
return(output)
}
##
calc.FPR0.p2 <- function(r1,r2,n1,n2,prior=0.5,alpha=0.05){
pval=prop.test(c(r1,r2),c(n1,n2))$p.value
ES=pwr::ES.h(r1/n1,r2/n2)
my.power=pwr::pwr.2p2n.test(ES,n1,n2,sig.level=pval)$power
power=pwr::pwr.2p2n.test(ES,n1,n2,sig.level=alpha)$power
LH10=my.power/pval # NOTE: sig.level=pval, not 0.05
FPR=pval*(1-prior)/(pval*(1-prior)+my.power*prior)
output=c(my.power=my.power,power=power,p.value=pval,FPR.less=FPR,LH10=LH10)
return(output)
}
719:卵の名無しさん
19/01/11 16:39:30.67 JAXb750V.net
> # EMPA-REG OUTCOME test
> # URLﾘﾝｸ(square.umin.ac.jp)
> np=2333 ; rp=282
> n10=2345 ; r10=243
> n25=2342 ; r25=247
> calc.FPR0.p2(rp,r10+r25,np,n10+n25)
my.power power p.value FPR.less LH10
0.50788496 0.53163451 0.04342913 0.07877384 11.69456891
ここに上げたけど、外部パッケージが必要なのでグラフしか出力されないな。
URLﾘﾝｸ(tpcg.io)
720:卵の名無しさん
19/01/11 16:55:13.20 JAXb750V.net
p値って（帰無仮説棄却の）結論が誤っている確率だと思っている椰子が多いんだよなぁ。
ｽﾚﾘﾝｸ(math板:701番)
721:卵の名無しさん
19/01/11 19:59:26.63 JAXb750V.net
calc.FPR = function(nsamp,pval,sigma,prior,delta1){
# sdiff=sqrt(sigma^2/nsamp + sigma^2/nsamp)
ns1=nsamp
ns2=nsamp
# CALC sdiff = sd of difference between means
sdiff=sqrt(sigma^2/ns1 + sigma^2/ns2)
df=ns1+ns2-2
# Note FPR doesn't need calculation of power for p-equals case
#
#under H0, use central t distribution
tcrit=qt((1-pval/2),df,ncp=0)
x0=tcrit
y0=dt(x0,df,0)
#
# under H1 use non-central t distribution
ncp1=delta1/sdiff #non-centrality paramater
x1=x0 #tcrit
y1=dt(x1,df,ncp=ncp1)
# check solution
# pcheck=pt(y1,df,ncp=ncp1)
# pcheck
# Calc false positive risk
p0=2*y0
p1=y1
FPR=((1-prior)*p0)/(((1-prior)*p0) + prior*p1)
cat('FPR = ',FPR,'\n')
output=c(FPR,x0,y0,x1,y1)
invisible(output)
}
# end of function calc.FPR
#
722:卵の名無しさん
19/01/11 19:59:55.21 JAXb750V.net
# calc.FPR0 gives FPR for given nsamp, for p-less-than case
calc.FPR0 = function(nsamp,pval,sigma,prior,delta1){
myp=power.t.test(n=nsamp,sd=sigma,delta=delta1,sig.level=pval,type="two.sample",alternative="two.sided",power=NULL)
power = myp$power
PH1=prior
PH0=1-PH1
FPR0=(pval*PH0/(pval*PH0 + PH1*power))
output=c(FPR_less=FPR0,power=power)
return(output)
}
#
calc.FPR(nsamp=7,pval=0.01985,sigma=1,prior=0.5,delta1=1)
calc.FPR0(nsamp=7,pval=0.01985,sigma=1,prior=0.5,delta1=1)
> calc.FPR(nsamp=7,pval=0.01985,sigma=1,prior=0.5,delta1=1)
FPR = 0.1245345
> calc.FPR0(nsamp=7,pval=0.01985,sigma=1,prior=0.5,delta1=1)
FPR_less power
0.07274575 0.25301819
>
723:卵の名無しさん
19/01/11 20:22:21.59 JAXb750V.net
>>680
そのシミュレーション
# t.test simulation when 0.045< p.value < 0.05
N=1000
n=10
x0=rnorm(N)
x1=rexp(N)
f<-function(prior=0.5){
choice=sample(1:0,1,prob=c(prior,1-prior))
x=sample(x0,n)
if(choice){
y=sample(x1,n) # choose different mean population
}else{
y=sample(x0,n) # choose the same population
}
c(t.test(x,y)$p.value,choice)
}
re=replicate(1e4,f(0.5))
res=re[2,re[1,]<0.045 & re[1,]<=0.050]
mean(res==0) # proportion of the same population = FPR
724:卵の名無しさん
19/01/14 21:11:37.39 W3rmWG0w.net
疲れた
ダルい
働いたら負け
725:卵の名無しさん
19/01/14 21:59:39.64 lm1UHN1Y.net
働きたくない
726:卵の名無しさん
19/01/15 11:09:44.07 dvPyxe58.net
待機で賃金が発生する職場が( ・∀・)ｲｲ!!
727:卵の名無しさん
19/01/15 13:02:35.43 9BKWXGSR.net
>>685
職場に出勤しなくちゃいけない時点で負け組
728:卵の名無しさん
19/01/15 20:42:00.10 LzPFbWpe.net
>>686
自宅だと飲酒するから身体に悪い
729:卵の名無しさん
19/01/16 08:14:37.24 gq6jTqnA.net
>>687
焼肉屋のオーナーなんぞ
従業員を働かせて
毎日ゴルフ三昧だぞ
730:卵の名無しさん
19/01/16 17:55:44.22 ZhNONJwp.net
>>688
別にゴルフしたくないんだが。
こういう問題が解けた方が楽しい。
URLﾘﾝｸ(rio2016.2ch.net)
まとまった時間がとれないとできないな。
731:卵の名無しさん
19/01/17 13:00:45.59 jUkqCyOy.net
>>689
ゴルフは単なる症例報告だ
俺もゴルフは全くやらん
調剤薬局もそうだ
医師は自分で働かないと無収入になるのに対して
薬局オーナーは出勤せずに毎日遊び歩いていても
病院オーナーよりはるかに高い収入を得ている
しかも薬局オーナーは薬剤師である必要すらない
732:卵の名無しさん
19/01/17 14:01:12.40 qxQFIIFv.net
>>690
オーナー理事長で役員報酬の不労所得の医師もいるぞ。
733:卵の名無しさん
19/01/17 15:18:23.26 qxQFIIFv.net
治療薬ハンドブックを読んで該当箇所のポケット医薬品集を読むと
同効薬の違いが分かって(・∀・)ｲｲ!!
防腐剤を含まない点眼抗菌薬はどれかとか実に役立つ。
734:卵の名無しさん
19/01/17 19:35:01.75 jUkqCyOy.net
>>691
そいつは働いたことが無いのかい？
でなきゃ妬ましくないな
俺の同級生の土建屋やパチンコ屋の子弟はひどいぞ
最終学歴が高卒だけど
就業経験ゼロで役員になって一生豪遊してるからな
735:卵の名無しさん
19/01/19 10:34:13.59 8uFRJt6f.net
>>693
パチンコ屋や土建屋の子弟ってシリツ医大に多いよな。
736:卵の名無しさん
19/01/19 14:00:04.35 +22CUuCd.net
>>694
知り合いの話かい？
他大学のことは知らんが
俺の出身の国立医学部にも土建屋の娘は居た
それはそうとニュースでミッキーハウス事件を見たことがあるかい？
あのニュースを見て理解できるのは
土建屋は病院なんぞ比べ物にならんくらい稼いでいるという事実だ
まあいずれにしても言えることは
働いたら負け
ということだ
737:卵の名無しさん
19/01/19 14:27:08.71 +Z1Yrlmo.net
>>695
この国には勤労・納税・教育の義務があんのよ。
738:卵の名無しさん
19/01/19 14:28:40.93 +Z1Yrlmo.net
>>696
納税したら負け、という価値観の椰子は多いだろうなぁ。
739:卵の名無しさん
19/01/19 15:20:32.33 tosZV9uO.net
働いたら負けだお
お家でカウチポテトしてコーラ飲んで
動画見てるのが最高だお
740:卵の名無しさん
19/01/19 16:57:43.58 8uFRJt6f.net
>>698
いや、学会出張に愛人を連れて行く方が楽しいんだが。
741:卵の名無しさん
19/01/19 17:22:09.90 HPdDoYmX.net
>>699
学会のついでになど下らん
平日の昼間に遊び歩いてこそ勝ち組だ
742:卵の名無しさん
19/01/19 17:46:18.37 8uFRJt6f.net
>>700
学会のついでの方が燃えるぞ！
743:卵の名無しさん
19/01/19 20:05:18.74 nd1Bx6ok.net
URLﾘﾝｸ(virtualsan-looking.jp)
三権分立が存在しない権力のゴキブリ自民レイプ山口敬之安倍官製ドワンゴとゴキブリ共は庵のは原爆で細胞ごと死滅しろw
744:卵の名無しさん
19/01/19 21:15:57.55 m2smxm4o.net
疲れた
ダルい
働きたくない
745:卵の名無しさん
19/01/19 22:13:35.46 +Z1Yrlmo.net
curve(dnorm(x,50,10),0,100,bty='l')
curve(dnorm(x,50,10),q60,100,type='h',col='yellow',add=T)
curve(dnorm(x,50,10),q30,100,type='h',col='green',add=T,lty=3)
segments(q60,0,q60,dnorm(q60,50,10),lty=3)
pnorm(c(65,68,69),50,10,lower=F)
(q30=qnorm(1-0.3,50,10))
(q60=qnorm(1-0.6,50,10))
integrate(function(x) x*dnorm(x,50,10)/0.3,q30,Inf)$value - integrate(function(x) x*dnorm(x,50,10)/0.6,q60,Inf)$value
library(BEST)
sim <- function(){
N=1e4
x=scale(rnorm(N/0.3))*10+50
x30=x[x>=q30]
x=scale(rnorm(N/0.6))*10+50
x60=x[x>=q60]
d=x30-x60
mean(d)
}
re=replicate(1e3,sim())
summary(re)
746:卵の名無しさん
19/01/19 22:32:24.95 m2smxm4o.net
働いたら負け
747:卵の名無しさん
19/01/19 22:49:34.18 NSnAx4LQ.net
>>700
盆や年末に俗世を離れてコミケに行く
アタイは勝ち組
748:卵の名無しさん
19/01/20 08:22:51.49 cvbKrpB0.net
f <- function(p,sd=10){
m=sample(100,1)
ma=integrate(function(x)x*dnorm(x,m,sd)/(1-p),-Inf,qnorm(1-p,m,sd))$value
m-ma
}
f=Vectorize(f)
pp=seq(0.001,0.999,0.001)
plot(pp,f(pp),bty='l',type='l',xlab='辞退率',ylab='虚飾値',lwd=2)
749:卵の名無しさん
19/01/20 08:42:33.97 cvbKrpB0.net
'(1)予備校の発表する偏差値は合格者の偏差値の平均値
(2)合格者の成績分布は正規分布で標準偏差は母集団と同じ
(3)成績のよい方から入学を辞退する
虚飾値＝合格者の偏差値 - 入学者の偏差値
'
f <- function(p,sd=10){
m=sample(100,1)
ma=integrate(function(x)x*dnorm(x,m,sd)/(1-p),-Inf,qnorm(1-p,m,sd))$value
m-ma
}
f=Vectorize(f)
pp=seq(0.001,0.999,0.001)
plot(pp,f(pp),bty='l',type='l',xlab='辞退率',ylab='虚飾値',lwd=2)
p=seq(0.1,0.9,0.1)
data.frame(辞退率=p,虚飾値=round(f(p),1))
750:卵の名無しさん
19/01/20 08:46:51.24 cvbKrpB0.net
>>708
虚飾値がｍの値によらないのは興味深いな。
751:卵の名無しさん
19/01/20 18:02:46.35 4pPRgzzq.net
明日は月曜
暗い日曜日
752:卵の名無しさん
19/01/20 18:31:46.86 vchoPnFo.net
>>710
明日の午後から火曜日まで代休。
753:卵の名無しさん
19/01/20 18:32:59.56 vchoPnFo.net
>>709
以下の(1)(2)(3)を前提にして計算
(1)予備校の発表する偏差値は合格者の偏差値の平均値
(2)合格者の成績分布は正規分布で標準偏差は母集団と同じ
(3)成績のよい方から入学を辞退する
虚飾値＝合格者の偏差値 - 入学者の偏差値
と定義して、辞退率と虚飾値の関係をグラフにすると
URLﾘﾝｸ(i.imgur.com)
辞退率 虚飾値
0.1 1.9
0.2 3.5
0.3 5.0
0.4 6.4
0.5 8.0
0.6 9.7
0.7 11.6
0.8 14.0
0.9 17.5
754:卵の名無しさん
19/01/20 19:47:26.59 4pPRgzzq.net
負～け負け負け負け～負け～
働いたら～負け～負け～
働い～たら～負～け～な～の～ら～♪
755:卵の名無しさん
19/01/21 06:45:51.44 20jAC8Su.net
教養がないのが負け組
756:卵の名無しさん
19/01/21 07:02:55.12 zCz5KSqt.net
ダルい
疲れた
働いたら負け
757:卵の名無しさん
19/01/21 11:08:16.93 Ed3wv7za.net
国試諦めてくれてありがとう。
758:卵の名無しさん
19/01/21 13:41:29.97 ulNFNnBT.net
先週は祝日だったので今日は内視鏡予約が多くて
この時間にステーキハウス。
眠剤のんでる老人はsedationの効きが悪いぜい。
759:卵の名無しさん
19/01/21 14:28:43.64 ulNFNnBT.net
>>709
この証明は難しいな。
数値積分での実感は簡単（ド底辺シリツ医を除く）だけどね。
760:卵の名無しさん
19/01/21 14:33:18.08 ulNFNnBT.net
オセロやルービックをプログラムで解く秀才は尊敬するな。
医学部にはいない。
761:卵の名無しさん
19/01/21 14:57:25.09 Ed3wv7za.net
良かったなー本当にありがとうございます。医者にならないで。良かったなー。
762:卵の名無しさん
19/01/21 15:04:18.32 Ed3wv7za.net
ありがとうございます。患者さんや医療従事者にとっては医者を諦めていただいて良かったです。
763:卵の名無しさん
19/01/21 15:04:35.86 ulNFNnBT.net
馬鹿でもできるのが臨床医。
算数の掛け算すら怪しいのがシリツ医がこういう事故を起こすわけなんだな。
URLﾘﾝｸ(i.imgur.com)
764:卵の名無しさん
19/01/21 15:12:33.23 Ed3wv7za.net
問題のある人物を医療の世界に入れることぐらい怖いものはないです。本当にありがとうございます。
765:卵の名無しさん
19/01/21 15:40:19.56 Ed3wv7za.net
医師の道を諦められたことに 同意いたします。その判断は正しいです。
766:卵の名無しさん
19/01/21 15:49:54.66 Ed3wv7za.net
ホッと一息 大学 厚生労働省。良かったね。
767:卵の名無しさん
19/01/21 15:58:20.70 xX0sMtrn.net
カネが欲しいから医学部にいくんや
文句あっか
768:卵の名無しさん
19/01/21 16:06:43.94 Ed3wv7za.net
危な過ぎて 許容できない。by 大学 厚生労働省
769:卵の名無しさん
19/01/21 16:34:11.82 5GFeLbpQ.net
東京医大が本来合格者に追加合格させても
本来不合格者を除籍処分しないのだから
シリツ＝裏口容疑者が確定だよ。
シリツ出身者こそ、裏口入学に厳しい処分せよを訴えるべき。
裏口入学医師の免許剥奪を！の国民運動の先頭に立てばよいぞ。
僕も裏口入学とか、言ってたら信頼の回復はない。
770:卵の名無しさん
19/01/21 16:37:46.83 5GFeLbpQ.net
>>709
証明は複雑な定積分が必要だな。
俺には無理なので数値積分で確認。
中心極限定理もそうだな。
771:卵の名無しさん
19/01/21 16:38:36.40 5GFeLbpQ.net
>>726
接客業は賤業。
772:卵の名無しさん
19/01/21 18:54:15.63 guiK2BJ8.net
>>730
仕事は何をしてるのかね？
いまどき正規職員でないものは完全に負け組だぞ
ゆっくり霊夢はFランク大学の就職課に就職したようです【第16話】
URLﾘﾝｸ(www.youtube.com)
773:卵の名無しさん
19/01/21 19:06:50.31 5GFeLbpQ.net
>>729
これを一般式でするのは大変だろうな。
f(x) = 1/{σ√2π)} exp［-(x-m)^2/ 2σ^2],
下位30% ････ 偏差値(m-0.5244σ)以下
f(x)
ma = ∫[-∞, m-0.5244σ] x･f(x) dx } / ∫[-∞, m-0.5244σ] f(x) dx{
= m + ∫[-∞, m-0.5244σ] (x-m) f(x) dx / 0.3
= m + (σ/√2π)∫[-∞, -0.5244] t･exp(-tt/2) dt / 0.3
= m + (σ/√2π) [ -exp(-tt/2) ](t=-∞, -0.5244) / 0.3
= m + (σ/√2π) [ -0.8715 / 0.3 ]
= m - 1.159σ
= m - 11.59
774:卵の名無しさん
19/01/21 19:07:33.35 5GFeLbpQ.net
累積密度関数の逆関数が必要になる:(
775:卵の名無しさん
19/01/21 23:54:21.71 zCz5KSqt.net
タバコって燃えるゴミかな？
有害物質を発生するからプラゴミと同じかな？
灰皿は不燃ゴミ
電子タバコは家電ゴミでいいのかな？
776:卵の名無しさん
19/01/22 00:16:37.39 7DkisZz5.net
>>734
吸ってる椰子がゴミ
777:卵の名無しさん
19/01/22 17:20:22.09 XyJ4p7WH.net
f <-function(a){
qa=qnorm(a,50,10,lower=F)
ma=integrate(function(x) x*dnorm(x,50,10)/a,qa,Inf)$value
va=integrate(function(x) (x-ma)^2*dnorm(x,50,10)/a,qa,Inf)$value
c(mean=ma,sd=sqrt(va))
}
m1=f(0.27)[1]
s1=f(0.27)[2]
m2=f(0.50)[1]
s2=f(0.50)[2]
g <- function(x){
y= (x-50)/s1*10+m1
50+(y-m2)/s2*10
}
x=40:70
g(x)
778:卵の名無しさん
19/01/22 17:40:34.59 XyJ4p7WH.net
f <-function(a){
qa=qnorm(a,50,10,lower=F)
ma=integrate(function(x) x*dnorm(x,50,10)/a,qa,Inf)$value
va=integrate(function(x) (x-ma)^2*dnorm(x,50,10)/a,qa,Inf)$value
c(mean=ma,sd=sqrt(va))
}
m1=f(0.27)[1]
s1=f(0.27)[2]
m2=f(0.50)[1]
s2=f(0.50)[2]
g <- function(x){
y= (x-50)/s1*10+m1
50+(y-m2)/s2*10
}
g(c(50,55,57.5))
母集団（例えば18歳人口)の成績が平均値50点標準偏差10点の正規分布に従うとする。
上位27%が進学すると
平均値62.2点標準偏差5.01点の集団Aとなる。
(成績が62.2点なら偏差値50に67.21点なら偏差値60になる)
上位50 %が進学すると
平均値58.0標準偏差6.03の集団Bとなる。
集団Aでの偏差値50, 55, 57.5は各々
集団Bでは偏差値57.0, 73.0, 81.9になる。
779:卵の名無しさん
19/01/22 18:22:22.34 XyJ4p7WH.net
>>737
バグ修正
f <-function(a){
qa=qnorm(a,50,10,lower=F)
ma=integrate(function(x) x*dnorm(x,50,10)/a,qa,Inf)$value
va=integrate(function(x) (x-ma)^2*dnorm(x,50,10)/a,qa,Inf)$value
c(mean=ma,sd=sqrt(va))
}
m1=f(0.27)[1]
s1=f(0.27)[2]
m2=f(0.50)[1]
s2=f(0.50)[2]
g <- function(x){
y= (x-50)/10 *s1+m1
50+(y-m2)/s2*10
}
g(c(50,55,57.5))
母集団（例えば18歳人口)の成績が平均値50点標準偏差10点の正規分布に従うとする。
上位27%が進学すると
平均値62.2点標準偏差5.01点の集団Aとなる。
(成績が62.2点なら偏差値50に67.21点なら偏差値60になる)
上位50 %が進学すると
平均値58.0標準偏差6.03の集団Bとなる。
集団Aでの偏差値50, 55, 57.5は各々
集団Bでは偏差値57.0, 61.2 63.3になる。
780:卵の名無しさん
19/01/22 18:32:41.55 XyJ4p7WH.net
f <-function(a){
qa=qnorm(a,50,10,lower=F)
ma=integrate(function(x) x*dnorm(x,50,10)/a,qa,Inf)$value
va=integrate(function(x) (x-ma)^2*dnorm(x,50,10)/a,qa,Inf)$value
c(mean=ma,sd=sqrt(va))
}
p1=0.27
p2=0.50
m1=f(p1)[1]
s1=f(p1)[2]
m2=f(p2)[1]
s2=f(p2)[2]
g <- function(x){
y= (x-50)/10 *s1+m1
50+(y-m2)/s2*10
}
g(c(50,55,57.5))
xx=seq(40,80,0.1)
plot(xx,g(xx))
781:卵の名無しさん
19/01/22 18:47:08.16 XyJ4p7WH.net
f <-function(a){
qa=qnorm(a,50,10,lower=F)
ma=integrate(function(x) x*dnorm(x,50,10)/a,qa,Inf)$value
va=integrate(function(x) (x-ma)^2*dnorm(x,50,10)/a,qa,Inf)$value
c(mean=ma,sd=sqrt(va))
}
sim <- function(x=50,p1=0.3,p2=0.5){
m1=f(p1)[1]
s1=f(p1)[2]
m2=f(p2)[1]
s2=f(p2)[2]
y= (x-50)/10 *s1+m1
50+(y-m2)/s2*10
}
sim()
782:卵の名無しさん
19/01/22 23:57:18.92 poZn4cZ2.net
>>740
仕事は何をしてるのかね？
いまどき正規職員でないものは完全に負け組だぞ
783:卵の名無しさん
19/01/23 00:04:12.53 p2qfG+zh.net
>>741
何のプログラムかわかる？
784:卵の名無しさん
19/01/24 01:16:27.12 sxN7Of5+.net
テロ工作機関
福山友愛病院
なんと朝木明代市議のような他殺だと日本人に騒がれるから
日本国警察が、
朝鮮殺戮殺人教団に乗っとられたことにより警察から与えられていたのは、
犯罪ライセンス！！！！！！
病院で
785:医師から処方された薬は、何ら疑いを持つことなく飲む人も少なくないだろう。 医学や薬学の知識がある人や、持病で何年も投薬治療を行っている人などは、 また違った見方をしているかもしれないが、一般的に「医師からの処方薬」に対する信頼は篤い。 しかし、そんな「信頼」を揺るがすとんでもない事件が話題となっている。 ■「在庫処理」で不要な治療薬を投与広島県の福山友愛病院は、昨年11～12月の間に、 統合失調症などの患者6名に、本来必要のないパーキンソン病の治療薬を投与していたことが判明。 これは、病院を運営する医療法人会長の指示によるもので、病院側は「使用期限が迫った薬の在庫処理がきっかけのひとつ。」と説明しているとのこと。 会長は、同病院で精神科としても勤務しており、パーキンソン病の治療薬である 「レキップ」の錠剤を統合失調症患者など6人に投与するよう看護師らに指示。 投薬は複数回行われ、末丸会長は通常の 8倍の量を 指示していたという。
786:卵の名無しさん
19/01/24 07:39:17.09 YHL/Shf0.net
ところで、タバコって燃えるゴミですか？
有害物質を発生するのでプラゴミと同じ扱いでしょうか？
電子タバコは家電ゴミですよね
今度、小学校に教育実習に行くので
お家のタバコ用具を全てゴミに出すように指導してきます
787:卵の名無しさん
19/01/24 11:08:18.87 EfYrEwvB.net
URLﾘﾝｸ(www.niid.go.jp)
H3N2で 2/21vs0/21で有意差あるか？
n<5だとΧ二乗検定はあてにならない。
788:卵の名無しさん
19/01/24 15:16:53.44 7fT+1bS6.net
It is common knowledge among doctors and patients that Do-Teihen(exclusively bottom-leveled medical school) graduates mean morons who bought their way to Gachi'Ura(currently called by themselves)
According to the experience of entrance exam to medical school in the era of Showa, when the sense of discrimination against
privately-founded medical schools were more intense than it is now,
all such schools but for Keio had been so compared to some specialized institution for educable mentally retarded kids that nobody but imbecile successors of physicians in private practice had applied for admission.
There had been NOT a single classmate who chose willingly against his/her common sense to go to the Do-Teihen(exclusively bottom-leveled medical school, currently also known as Gachi'Ura),
which would have cost outrageous money and its graduates are destined to be called Uraguchi morons who bought thier way into the Do-Teihen, by thier colleagues and even by thier own clients.
Although people won't call them names to their face,
certain 80-90% people of about my age have been yet scorning and sneering at Uraguchi graduates, speaking in the back of our mind,
" Uraguchi morons shall not behave like somebody."
We never speak out face to face in real life.
789:卵の名無しさん
19/01/24 18:59:58.86 oaTB7QH8.net
テロ工作機関が日本にあることが判明
福山友愛病院
日本人への薬物大量投与テロも発覚
また他は、
単純に、朝鮮殺戮殺人革命の為に、
日本人被害者側を、
カルト医師が医師免許を犯罪に使用し、
偽造カルテをでっち上げているにすぎない
でっち上げる理由は、カルト会社の犯罪にごねている、
だから、
警察から犯罪ライセンス与えられているので、
虚偽説明して騙して閉鎖病棟に閉じ込め、
携帯とりあげ、
監禁罪やり精神保健福祉法違反をやり偽造カルテ書きますね
それだけw
だって、警察がグルで共同組織犯罪をしているから、
犯罪を犯罪でないことにしてくれるんだもの
日本人被害者側を拉致していることがバレたしな
それをやらかしていたので、
福山友愛病院には電凸し、
インターネット上に書いていることを通告もしている
ケースワーカーのイソムラという女性曰わく、
カミサカ先生インターネットに書かれて辛いみたいで～・・・
日本国内にいるテロリストが辛くなることは、
日本人の喜びだっつーのww
790:卵の名無しさん
19/01/24 20:31:49.63 H5DL72sh.net
Rnorm<-function(n,m,s){
x=rnorm(n,m,s)
m + s*(x-mean(x))/sd(x)
}
N=1e3
x=replicate(5,Rnorm(N,50,10))
f <- function(){ # simulated test results of 3 or 5 subjects
y=numeric()
for(i in 1:5) y[i]=sample(x[,i],1)
y3=sum(y[1:3])
y5=sum(y[1:5])
c(y3,y5)
}
p2h <- function(x){ # points to hensa-chi
m=mean(x)
s=sd(x)
x=50+10*(x-m)/s
}
re=replicate(1e5,f())
y3=re[1,] # points of 3 subject test
y5=re[2,]
z3=p2h(y3) # hensa-ti of 3 subject test
z5=p2h(y5)
791:卵の名無しさん
19/01/24 20:32:02.08 H5DL72sh.net
hc <- function(h,d=0.5){ # hensa-chi conversion
idx=which(h-d<z5 & z5<h+d) # hensa-chi range in 5 subject test
mean(z3[idx]) # mean of its hensa-chi in 3 subject
}
hc=Vectorize(hc)
hh=40:80
plot(hh,hc(hh),bty='l',asp=1,xlab='5 subject',ylab='3 subject',
pch=19)
abline(a=0,b=1,col='gray')
792:卵の名無しさん
19/01/24 20:51:53.08 H5DL72sh.net
実行結果
URLﾘﾝｸ(i.imgur.com)
793:卵の名無しさん
19/01/27 05:17:22.65 Ba4uprm1.net
札幌ひばりが丘病院
麻薬取締法違反で書類送検
URLﾘﾝｸ(video.fc2.com)
794:卵の名無しさん
19/01/28 11:16:13.00 C5thEPa8.net
ゆっくり霊夢はFランク大学の就職課に就職したようです
URLﾘﾝｸ(www.youtube.com)
795:卵の名無しさん
19/01/28 17:29:00.31 AVoD+Cyq.net
働いたら負け
796:卵の名無しさん
19/01/28 21:53:43.21 +PCpJyrP.net
It is common knowledge among doctors and patients that Do-Teihen(exclusively bottom-leveled medical school) graduates mean morons who bought their way to Gachi'Ura(currently called by themselves)
According to the experience of entrance exam to medical school in the era of Showa, when the sense of discrimination against
privately-founded medical schools were more intense than it is now,
all such schools but for Keio had been so compared to some specialized institution for educable mentally retarded kids that nobody but imbecile successors of physicians in private practice had applied for admission.
There had been NOT a single classmate who chose willingly against his/her common sense to go to the Do-Teihen(exclusively bottom-leveled medical school, currently also known as Gachi'Ura),
which would have cost outrageous money and its graduates are destined to be called Uraguchi morons who bought thier way into the Do-Teihen, by thier colleagues and even by thier own clients.
Although people won't call them names to their face,
certain 80-90% people of about my age have been yet scorning and sneering at Uraguchi graduates, speaking in the back of our mind,
" Uraguchi morons shall not behave like somebody."
We never speak out face to face in real life.
797:卵の名無しさん
19/02/09 22:33:46.30 F1R4fVG9.net
f=function(p,q=0.99){
a=13*q
c=13*(1-q)
b=13000*p
d=13000*(1-p)
(a/(a+b)-c/(c+d)) / (a/(a+b))
}
p=q=seq(0.01,0.99,0.1)
z=outer(p,q,f)
contour(p,q,z)
g=Vectorize(f)
g((1:19)/20)
798:卵の名無しさん
19/02/10 12:37:21.35 9eXx3CuI.net
options(digits=3)
f=function(p,q=0.99){ # p:pyroli in population, q:pyroli in cancer
a=13*q
c=13*(1-q)
b=13000*p
d=13000*(1-p)
x=(a/(a+b)-c/(c+d)) / (a/(a+b)) # attributable to pyroli(+)
y=(c/(c+d)-a/(a+b)) / (c/(c+d)) # qtteibutabel to pyroli(-)
round(ifelse(x>0,x,y),3)
}
g=Vectorize(f)
g((1:19)/20)
p=seq(0.01,0.99,0.01) # pylori in population
q=seq(0.01,0.99,0.01) # pylori in cancer
z=outer(p,q,f)
image(p,q,z,bty='l',nlevels=20,xlab='pyroli in population',
ylab='pyroli in cancer')
contour(p,q,z,bty='l',nlevels=20,xlab='pyroli in population',
ylab='pyroli in cancer',add=T)
contour(p,q,z,bty='l',nlevels=20,xlab='pyroli in population',
ylab='pyroli in cancer')
source('tools.R')
Persp(p,q,z)
Persp3d(p,q,z)
799:卵の名無しさん
19/02/10 12:37:47.54 9eXx3CuI.net
'
がん 非がん
感染 a b
未感染 c d
attributable risk : x
x={a/(a+b)-c/(c+d)} / {a/(a+b)}
'
800:卵の名無しさん
19/02/10 19:40:18.09 9eXx3CuI.net
q1=0.9934
q2=0.9958
f <- function(p,q){
a=13*q
c=13*(1-q)
b=13000*p
d=13000*(1-p)
(a/(a+b)-c/(c+d)) / (a/(a+b))
}
x0=0.985
f1=Vectorize(function(p)f(p,q=q1))
f2=Vectorize(function(p)f(p,q=q2))
curve(f1(x))
abline(h=0.985)
uniroot(function(x)f1(x)-x0,c(0,0.9))$root
curve(f2(x))
abline(h=0.985)
uniroot(function(x)f2(x)-x0,c(0,0.9))$root
801:卵の名無しさん
19/02/11 12:56:16.04 rzg7XY57.net
一般論として医療統計って全然科学じゃないって、研究者の中ではバカにされてますよ。
802:卵の名無しさん
19/02/11 14:44:26.58 Pny0jm5M.net
統計で使われる数学って難し過ぎてついていけないよね。
AICとかわからないままにソフトで出させているだけ。
803:卵の名無しさん
19/02/13 12:53:38.65 hLnbSCPn.net
It is common knowledge among doctors and patients that Do-Teihen(exclusively bottom-leveled medical school) graduates mean morons who bought their way to Gachi'Ura(currently called by themselves)
According to the experience of entrance exam to medical school in the era of Showa, when the sense of discrimination against
privately-founded medical schools were more intense than it is now,
all such schools but for Keio had been so compared to some specialized institution for educable mentally retarded kids that nobody but imbecile successors of physicians in private practice had applied for admission.
There had been NOT a single classmate who chose willingly against his/her common sense to go to the Do-Teihen(exclusively bottom-leveled medical school, currently also known as Gachi'Ura),
which would have cost outrageous money and its graduates are destined to be called Uraguchi morons who bought thier way into the Do-Teihen, by thier colleagues and even by thier own clients.
Although people won't call them names to their face,
certain 80-90% people of about my age have been yet scorning and sneering at Uraguchi graduates, speaking in the back of our mind,
" Uraguchi morons shall not behave like somebody."
We never speak out face to face in real life.
804:卵の名無しさん
19/02/16 13:25:47.75 QaizKjLi.net
>>761
裏口馬鹿にもわかるように原文もつけなきゃ。
原文はこれな。
私は昭和の時代に大学受験したけど、昔は今よりも差別感が凄く、慶応以外の私立医は特殊民のための特殊学校というイメージで開業医のバカ息子以外は誰も受験しようとすらしなかった。
常識的に考えて、数千万という法外な金を払って、しかも同業者からも患者からもバカだの裏口だのと散々罵られるのをわかって好き好んで私立医に行く同級生は一人もいませんでした。
本人には面と向かっては言わないけれど、俺くらいの年代の人間は、おそらくは８－９割は私立卒を今でも「何偉そうなこと抜かしてるんだ、この裏口バカが」と心の底で軽蔑し、嘲笑しているよ。当の本人には面と向かっては絶対にそんなことは言わないけどね。
805:卵の名無しさん
19/02/17 22:18:47.40 rglxpU9N.net
It is common knowledge among doctors and patients that Do-Teihen(exclusively bottom-leveled medical school) graduates mean morons who bought their way to Gachi'Ura(currently called by themselves)
According to the experience of entrance exam to medical school in the era of Showa, when the sense of discrimination against
privately-founded medical schools were more intense than it is now,
all such schools but for Keio had been so compared to some specialized institution for educable mentally retarded kids that nobody but imbecile successors of physicians in private practice had applied for admission.
There had been NOT a single classmate who chose willingly against his/her common sense to go to the Do-Teihen(exclusively bottom-leveled medical school, currently also known as Gachi'Ura),
which would have cost outrageous money and its graduates are destined to be called Uraguchi morons who bought thier way into the Do-Teihen, by thier colleagues and even by thier own clients.
Although people won't call them names to their face,
certain 80-90% people of about my age have been yet scorning and sneering at Uraguchi graduates, speaking in the back of our mind,
" Uraguchi morons shall not behave like somebody."
We never speak out face to face in real life.
806:卵の名無しさん
19/02/19 19:34:53.54 5dzWj+iV.net
俺は医者ではない。
807:卵の名無しさん
19/02/19 22:25:59.08 j+jAYShv.net
医学統計は苦手だ
できなきゃ論文読めんけど
808:卵の名無しさん
19/02/20 09:47:48.92 P0P24ylG.net
>>765
>456の問題点を指摘できないなら不正統計に騙される素養があるね。
809:卵の名無しさん
19/02/27 17:12:27.40 L7SdQ9d/.net
まあ落ち着いて
萌え動画でもどうぞ
URLﾘﾝｸ(m.youtube.com)
810:卵の名無しさん
19/03/01 20:33:40.05 CpB3LdCN.net
「もし君の妻子や友人にいつまでも生きていてもらいたいと思うならば、君は愚かである。というのは、君は、君の力の及びうる範囲内にないものまでも君の力の及びうる範囲内にあることを欲し、よそのものまでも君のものであることを欲しているからである。」
811:卵の名無しさん
19/03/04 13:04:58.69 PhmNbZK3.net
>>767 これでも見て落ちつけ https://www.youtube.com/watch?v=mB3ZVIc4EF4
813:卵の名無しさん
19/03/14 16:45:10.23 JBQs0yPH.net
韓国で大麻が医療解禁した
国際条約では臨床試験や学術研究(輸出輸入など)が認められているんで一度訪ねてみては?
URLﾘﾝｸ(livedoor.blogimg.jp)
URLﾘﾝｸ(livedoor.blogimg.jp)
韓国で使われてる大塚製薬の大麻製剤は日本の医師会員向けに和訳も開始されている。
URLﾘﾝｸ(news.nicovideo.jp)
814:卵の名無しさん
19/03/18 15:05:39.58 pTprsItO.net
まあ落ち着いて
萌え動画でもどうぞ
URLﾘﾝｸ(www.youtube.com)
815:卵の名無しさん
19/03/27 22:03:49.30 Kh2JcAUc.net
統計楽しいね
816:卵の名無しさん
19/03/28 00:45:03.92 8PkLLiIS.net
URLﾘﾝｸ(www3.nhk.or.jp)
インフルエンザの新しい治療薬「ゾフルーザ」を投与されたＡ香港型のインフルエンザ患者30人を調べたところ、
70％余りに当たる22人から、この薬が効きにくい耐性ウイルスが検出されたことが国立感染症研究所の調査で分かりました。
調査件数は多くないものの、専門家は現在のような使用を続けると、耐性ウイルスが広がるおそれがあるとして使用基準を見直すべきだと指摘しています。
URLﾘﾝｸ(www3.nhk.or.jp)
問題：耐性発生率の95%信頼区間は？ 👀
Rock54: Caution(BBR-MD5:1341adc37120578f18dba9451e6c8c3b)
817:卵の名無しさん
19/03/28 00:49:17.10 8PkLLiIS.net
95%CI 0.541 0.877なので過半数に耐性出現と言っていいな。
818:卵の名無しさん
19/03/28 01:30:14.34 uQgP69pk.net
耐性出現率が過半数である確率はいくらか？
819:卵の名無しさん
19/06/25 16:00:08.86 ZxTc5GC6.net
>>775
あんた医科歯科大卒だろ？
学部1～２回生の時に線形代数と微分積分はやった？
820:卵の名無しさん
19/06/25 16:01:39.44 ZxTc5GC6.net
線形代数と微分積分を知らんと統計学の勉強は無理だね。
まあ、上を見れば切りがないけどね。
821:卵の名無しさん
19/06/28 10:53:20.90 30gP48kS.net
>>777
Rが使えれば統計扱えるよ。
今どき分布表で補間値計算なんてしないし。
非心ｔ分布もRが出してくれる。
URLﾘﾝｸ(www3.nhk.or.jp)
インフルエンザの新しい治療薬「ゾフルーザ」を投与されたＡ香港型のインフルエンザ患者30人を調べたところ、70％余りに当たる22人から、この薬が効きにくい耐性ウイルスが検出されたことが国立感染症研究所の調査で分かりました。
調査件数は多くないものの、専門家は現在のような使用を続けると、耐性ウイルスが広がるおそれがあるとして使用基準を見直すべきだと指摘しています。
耐性化率が50%以上である確率は
pbeta(0.5,1+22,1+8,lower=F)
[1] 0.9946631
822:卵の名無しさん
19/06/28 12:35:11.06 30gP48kS.net
こんな質問があったからRを使って答えておいた。
他の人も同じ数値出していたな。
655 名無しさん＠１周年 sage 2019/06/28(金) 12:05:50.17 ID:vuCIc6fU0
偏差値72.5ってIQにするとどれくらいなの？
URLﾘﾝｸ(asahi.2ch.net)
823:卵の名無しさん
19/06/28 20:16:59.55 30gP48kS.net
統計学のための数学
あんまり公式を使うことないなぁ。
定積分も数値積分でやっちゃうから
URLﾘﾝｸ(www.data-arts.jp) <
824:卵の名無しさん
19/10/02 15:20:15.02 nv8GYZZt.net
kk=(88+96)/2/100 #TP: kessei kando
kt=(89+100)/2/100 #TN: kessei tokuido
pLHk=kk/(1-kt) #TP/FP
nLHk=(1-kk)/kt #(1-TP)/(1-FP)=FN/TN
bk=(90+98)/2/100 #ben kando
bt=(87+100)/2/100 #ben tokuido
pLHb=bk/(1-bt)
nLHb=(1-bk)/bt
prep=0.5 # pre-possibility
preo=prep/(1-prep) # pre-odds
posto=preo*nLHk*pLHb # kessei(-)&ben(+)
(postp=posto/(1+posto))
#
posto=preo*nLHk*nLHb # kessei(-)&ben(-)
(postp=posto/(1+posto))
825:卵の名無しさん
19/10/02 17:26:13.36 nv8GYZZt.net
#
# 尿素呼気試験(B) 90～100 80～99
# 血清抗体(S) 88～96 89～100
# 尿中抗体(U) 89～97 77～95
# 便中抗原(F) 90～98 87～100
U=c(90/2+100/2,80/2+99/2)/100
S=c(88/2+96/2,89/2+100/2)/100
U=c(89/2+97/2,77/2+95/2)/100
F=c(90/2+98/2,87/2+100/2)/100
DOR <- function(T){
TP=T[1] ; FN=1-TP
TN=T[2] ; FP=1-TN
pLR=TP/FP
nLR=FN/TN
pLR/nLR
}
DOR(U)
DOR(S)
DOR(U)
DOR(F)
826:卵の名無しさん
19/10/02 17:34:50.26 nv8GYZZt.net
B=c(90/2+100/2,80/2+99/2)/100
S=c(88/2+96/2,89/2+100/2)/100
U=c(89/2+97/2,77/2+95/2)/100
F=c(90/2+98/2,87/2+100/2)/100
DOR <- function(T){
TP=T[1] ; FN=1-TP
TN=T[2] ; FP=1-TN
pLR=TP/FP
nLR=FN/TN
pLR/nLR
}
> DOR(B)
[1] 161.9524
> DOR(S)
[1] 197.5909
> DOR(U)
[1] 81.61224
> DOR(F)
[1] 225.359
>
827:卵の名無しさん
19/11/25 07:37:28 ZvHlCyJI.net
# ゴルゴ13は100発100中
# ゴルゴ14は10発10中
# ゴルゴ15は1発1中
# とする。
# 各々10000発撃ったとき各ゴルゴの命中数の期待値はいくらか？
G13 <- function(N,n,r){
pp=0:N
f <- function(x) choose(x,r)*choose(N-x,n-r)/choose(N,n)
sum(pp*f(pp)/sum(f(pp)))
}
G13(10000,100,100)
G13(10000,10,10)
G13(10000,1,1)
#.n発r中の狙撃手が.N発狙撃するときの命中数を返す
Go13 <- function(.N, .n, r, k=10^3){ # k:シミュレーション回数
f <-function(S,N=.N,n=.n){
y=c(rep(1,S),rep(0,N-S))
sum(sample(y,n))
}
xx=r:.N
SS=NULL
for(i in 1:k){
x=sapply(xx,f)
SS=c(SS,which(x==r)-1+r)
}
print(summary(SS))
invisible(SS)
}
828:卵の名無しさん
19/11/27 20:32:52 s7+UeQPF.net
a=c(1,10,100)
curve(dbeta(x,1+a[1],1),bty='l')
curve(dbeta(x,1+a[2],1),bty='l')
curve(dbeta(x,1+a[3],1),bty='l')
1000*qbeta(0.95,a+1,1,lower=F)
829:卵の名無しさん
19/11/27 20:34:40 s7+UeQPF.net
> a=c(1,10,100)
> curve(dbeta(x,1+a[1],1),bty='l')
> curve(dbeta(x,1+a[2],1),bty='l')
> curve(dbeta(x,1+a[3],1),bty='l')
> 10000*qbeta(0.95,a+1,1,lower=F)
[1] 2236.068 7615.958 9707.748
>
830:卵の名無しさん
19/11/30 09:56:41 bijshzSl.net
# 安倍総理招待枠(60-****)が招待状から
# 254x 2409 357xだという。
# URLﾘﾝｸ(youtu.be)
# 通し番号がふられているとして
# 357xを3570と低めに見積もる。
#>いや、それどころか、今年はそれ以上の数字も確認されている。
#選挙で安倍首相が遊説をおこなうとかなりの頻度で目撃されている熱烈な支持者であるM氏という人物がいるのだが、
#氏が菅官房長官や杉田水脈議員とのツーショット写真などとともにSNSにアップしていた
# 今年の「桜を見る会」の受付票には「60-4908」とナンバリングされているのだ。
#ｽﾚﾘﾝｸ(seijinewsplus板)
# 他�
831:ﾉ今年は254? 2409 357? という番号が国会で発言されている。 #参加者の18200の半数9100が官邸枠の可能性があると仮定して官邸枠招待者数の期待値と95%CIを算出せよ。 Nmin=4908 Nmax=9100 n=Nmin:Nmax m=4 pmf=choose(Nmax-1,m-1)/choose(n,m) #Pr(max=60|n) plot(n,pmf,ylab='probability') pdf=pmf/sum (pmf) plot(n,pdf,ylab='density',bty='l') sum(n*pdf) #E(n) plot(n,cumsum(pdf)) abline(h=0.95,lty=3) idx=min(which(cumsum(pdf)>0.95)) n[idx]
832:卵の名無しさん
19/11/30 10:28:03.80 bijshzSl.net
> Nmin=4908
> Nmax=9100
> n=Nmin:Nmax
> m=4
> pmf=choose(Nmax-1,m-1)/choose(n,m) #Pr(max=60|n)
> plot(n,pmf,ylab='probability')
> pdf=pmf/sum (pmf)
> plot(n,pdf,ylab='density',bty='l')
> sum(n*pdf) #E(n)
[1] 6191.377
> plot(n,cumsum(pdf))
> abline(h=0.95,lty=3)
> idx=min(which(cumsum(pdf)>0.95))
> n[idx]
[1] 8406
>
833:卵の名無しさん
19/12/09 08:25:56.61 zd3bijjh.net
# サイコロ
# 正６面体のサイコロがある．４面は青色、２面は赤色である．
# このサイコロを合計20回振るとき、最も起こりそうな順番はどれか？
# １．赤 青 赤 赤 赤
# ２．青 赤 青 赤 赤 赤
# ３．青 赤 赤 赤 赤 赤
sim <- function(){
a=sample(0:1,20, replace=TRUE, prob=c(4,2))
b=as.character(a)
c=paste(b,collapse="")
s1=paste(c(1,0,1,1,1),collapse="")
s2=paste(c(0,1,0,1,1,1),collapse="")
s3=paste(c(0,1,1,1,1,1),collapse="")
res=c(grepl(s1,c),grepl(s2,c), grepl(s3,c))
return(res)
}
k=1e5
re=replicate(k,sim())
mean(re[1,])
mean(re[2,])
mean(re[3,])
834:卵の名無しさん
19/12/10 08:45:53 k/8kaoYw.net
seqNp <- function(N=100,K=5,p=0.5){
if(p==0) return(0)
if(N==K) return(p^K)
q=1-p
a=numeric(N) # a(n)=P0(n)/p^n , P0(n)=a(n)*p^n
for(i in 1:K) a[i]=q/p^i # P0(i)=q
for(i in K:(N-1)){ # recursive formula
a[i+1]=0
for(j in 0:(K-1)){
a[i+1]=(a[i+1]+a[i-j])
}
a[i+1]=q/p*a[i+1]
}
P0=numeric(N)
for(i in 1:N) P0[i]=a[i]*p^i # P0(n)=a(n)*p^n
MP=matrix(rep(NA,N*K),ncol=K)
colnames(MP)=paste0('P',0:(K-1))
MP[,'P0']=P0
MP[1,'P1']=p
for(i in (K-2):K) MP[1,i]=0
for(k in 2:K){
for(i in 1:(N-1)) MP[i+1,k]=p*MP[i,k-1]
} # Pk(n+1)=p*P(k-1)(n)
ret=1-apply(MP,1,sum)
ret[N]
}
835:卵の名無しさん
19/12/10 08:46:19 k/8kaoYw.net
seqNpJ <- function(N,K,p) seqNp(N,K,p)-seqNp(N,K+1,p)
seqNpJ(100,5,0.5)
seq2pCI <- function(N,K,alpha=0.05,Print=T){
vp=Vectorize(function(p)seqNp(N,K,p)-seqNp(N,K+1,p))
if(Print){curve(vp(x),lwd=2,bty='l',xlab='Pr[head]',ylab=paste('Pr[max',K,'-head repetition]'))
abline(h=alpha,lty=3)}
peak=optimize(vp,c(0,1),maximum=TRUE)$maximum
mean=integrate(function(x)x*vp(x),0,1)$value/integrate(function(x)vp(x),0,1)$value
lwr=uniroot(function(x,u0=alpha) vp(x)-u0,c(0.01,peak))$root
upr=uniroot(function(x,u0=alpha) vp(x)-u0,c(peak,0.99))$root
c(lower=lwr,mean=mean,mode=peak,upper=upr)
}
seq2pCI(100,5,0.05,T)
vs=Vectorize(function(K)seq2pCI(N=100,K,alpha=0.05,Print=F))
x=2:23
y=vs(x)
head(y)
y=y*100
plot(x,y['mean',],bty='l',pch=19,ylim=c(0,100),
xlab="最大連続数",ylab="推定裏口入学者数")
points(2:23,y['mode',],bty='l')
segments(x,y['lower',],x,y['upper',])
legend('right',bty='n',pch=c(19,1),legend=c("期待値","最頻値"))
836:卵の名無しさん
19/12/10 08:46:54 k/8kaoYw.net
# pdfからcdfの逆関数を作ってhdiを表示させて逆関数を返す
# 0,1での演算を回避 ∫[1/nxxx,1-1/nxx]dxで計算
pdf2HDI <- function(pdf,xMIN=0,xMAX=1,cred=0.95,nxx=1001){
xx=seq(xMIN,xMAX,length=nxx)
xx=xx[-nxx]
xx=xx[-1]
xmin=xx[1]
xmax=xx[nxx-2]
AUC=integrate(pdf,xmin,xmax)$value
PDF=function(x)pdf(x)/AUC
cdf <- function(x) integrate(PDF,xmin,x)$value
ICDF <- function(x) uniroot(function(y) cdf(y)-x,c(xmin,xmax))$root
hdi=HDInterval::hdi(ICDF,credMass=cred)
c(hdi[1],hdi[2])
}
# N試行で最大K回連続成功→成功確率ｐの期待値、最頻値と95%HDI
# max K out of N-trial to probability & CI
mKoN2pCI <- function(N=100 , K=4 , conf.level=0.95){
pmf=Vectorize(function(p)seqNp(N,K,p)-seqNp(N,K+1,p))
mode=optimize(pmf,c(0,1),maximum=TRUE)$maximum
auc=integrate(pmf,0,1)$value
pdf=function(x) pmf(x)/auc
mean=integrate(function(x)x*pdf(x),0,1)$value
curve(pdf(x),lwd=3,bty='l',xlab='Pr[head]',ylab='density')
lu=pdf2HDI(pdf,cred=conf.level)
curve(pdf(x),lu[1],lu[2],type='h',col='lightblue',add=T)
c(lu[1],mean=mean,mode=mode,lu[2])
}
mKoN2pCI(100,4) # one minute plz
mKoN2pCI(100,5)
837:卵の名無しさん
19/12/10 08:47:01 k/8kaoYw.net
# N試行で最大K回連続成功→成功確率ｐの期待値、最頻値と95% Quantile
# max K out of N-trial to probability & CI_quantile
mKoN2pCIq <- function(N=100 , K=4 , alpha=0.05){
pmf=Vectorize(function(p)seqNp(N,K,p)-seqNp(N,K+1,p))
mode=optimize(pmf,c(0,1),maximum=TRUE)$maximum
auc=integrate(pmf,0,1)$value
pdf=function(x) pmf(x)/auc
curve(pdf(x),bty='l')
mean=integrate(function(x)x*pdf(x),0,1)$value
cdf=function(x) MASS::area(pdf,0,x)
vcdf=Vectorize(cdf)
lwr=uniroot(function(x)vcdf(x)-alpha/2,c(0,mode))$root
upr=uniroot(function(x)vcdf(x)-(1-alpha/2),c(mode,1))$root
c(lower=lwr,mean=mean,mode=mode,upper=upr)
}
mKoN2pCIq(100,4)
mKoN2pCIq(100,5)
838:卵の名無しさん
19/12/10 08:48:13 2lYn70em.net
## simulation
mKoN2pCIga<-function(N=100,K=4,Print=T){
pmf=Vectorize(function(p)seqNp(N,K,p)-seqNp(N,K+1,p))
auc=integrate(pmf,0,1)$value
pdf=function(x) pmf(x)/auc
(mode=optimize(pmf,c(0,1),maximum=TRUE)$maximum)
(mean=integrate(function(x)x*pdf(x),0,1)$value)
(var=integrate(function(x)(x-mean)^2*pdf(x),0,1)$value)
(sd=sqrt(var))
#正規分布で近似してみる
if(Print){
curve(pdf,bty='l',col='red')
curve(dnorm(x,mean,sd),add=T,col='blue')}
c(lower=qnorm(0.025,mean,sd),mode=mode,mean=mean,upper=qnorm(0.975,mean,sd))
}
mKoN2pCIga(100,5)
vmK=Vectorize(function(K) mKoN2pCIga(100,K))
(y=cbind(0,vmK(2:30)))
plot(1:30,y['mean',],pch=19)
points(1:30,y['mode',])