底辺私立医大を卒業した医者って頭悪いよね？ Part16

タイトルまんまです

すなわち、タイトルの正しさを検証するスレです。

※前スレ
底辺私立医大を卒業した頭悪いよね？
http://egg.2ch.net/test/read.cgi/hosp/1472273747/
底辺私立医大を卒業した医者って頭悪いよね？ Part2
http://egg.2ch.net/test/read.cgi/hosp/1520807375/
底辺私立医大を卒業した医者って頭悪いよね？ Part3
http://egg.2ch.net/test/read.cgi/hosp/1522967067/
底辺私立医大を卒業した医者って頭悪いよね？ Part4
http://egg.2ch.net/test/read.cgi/hosp/1526677794/
底辺私立医大を卒業した医者って頭悪いよね？ Part5
http://egg.2ch.net/test/read.cgi/hosp/1529666055/
底辺私立医大を卒業した医者って頭悪いよね？ Part6
http://egg.5ch.net/test/read.cgi/hosp/1531305684/
底辺私立医大を卒業した医者って頭悪いよね？ Part7
http://egg.2ch.net/test/read.cgi/hosp/1532472024/
底辺私立医大を卒業した医者って頭悪いよね？ Part8
http://egg.2ch.net/test/read.cgi/hosp/1533510399/
底辺私立医大を卒業した医者って頭悪いよね？ Part9
http://egg.2ch.net/test/read.cgi/hosp/1535547275/
底辺私立医大を卒業した医者って頭悪いよね？ Part10
https://egg.2ch.net/test/read.cgi/hosp/1537867846/
底辺私立医大を卒業した医者って頭悪いよね？ Part11
https://egg.5ch.net/test/read.cgi/hosp/1544320163/
底辺私立医大を卒業した医者って頭悪いよね？ Part12
https://egg.5ch.net/test/read.cgi/hosp/1550969602/
底辺私立医大を卒業した医者って頭悪いよね？ Part13
https://egg.5ch.net/test/read.cgi/hosp/1575242106/
底辺私立医大を卒業した医者って頭悪いよね？ Part14
https://egg.5ch.net/test/read.cgi/hosp/1578616997/
底辺私立医大を卒業した医者って頭悪いよね？ Part15
https://egg.5ch.net/test/read.cgi/hosp/1584050953/

ド底辺シリツ医大卒ってこれだよな

底辺私立医大卒が国立大学医学部卒を妬むスレ [転載禁止](c)2ch.net
http://potato.5ch.net/test/read.cgi/hosp/1424846816/695

695 名前：卵の名無しさん[] 投稿日：2016/09/25(日) 22:04:25.53 ID:FzuOCpqV
匿名先生
精神科
▲
1
評価

私立出でも、素晴らしい医師がいるし、国立でも、どうしようもないのもいる。
そんなのは、あくまで個人の資質差と思ってました。

が、

今いる病院（医師の大部分が、いわゆる底辺私立）に来て、その考えは変わりました。

とある、気に食わない医師の白衣やサンダルを隠す、机の上にゴミを置く、など、小学生じみた嫌がらせをしている医師というのを、はじめて見ました。

国立出の変人も多々みてきましたが、こんな、小学生のイタズラみたいなことをする人は、いませんでした。

やはり、底辺私立は、それなりなんだと思いました。

2016年09月02日 00時21分09秒

人種や民族など本人が選択できないことを根拠とする区別は差別の批判を免れないが、
ド底辺シリツ医大進学は本人の選択だよなぁ。

これな！

不朽の名投稿
>>
私は昭和の時代に大学受験したけど、昔は今よりも差別感が凄く、特殊民のための特殊学校というイメージで開業医のバカ息子以外は誰も受験しようとすらしなかった。
常識的に考えて、数千万という法外な金を払って、しかも同業者からも患者からもバカだの裏口だのと散々罵られるのをわかって好き好んでド底辺医に行く同級生は一人もいませんでした。
本人には面と向かっては言わないけれど、俺くらいの年代の人間は、おそらくは８−９割はド底辺医卒を今でも「何偉そうなこと抜かしてるんだ、この裏口バカが」と心の底で軽蔑し、嘲笑しているよ。
当の本人には面と向かっては絶対にそんなことは言わないけどね。
<<

東京医大の事件は裏口入学が現在進行形であること如実にしめした事件だよね。
シリツ医の使命は裏口入学撲滅国民運動の先頭に立つことだよ。
裏口入学の学生を除籍処分にしないかぎり、信頼の回復はないね。つまり、いつまで経ってもシリツ医大卒＝裏口バカと汚名は拭えない。シリツ出身者こそ、裏口入学に厳しい処分せよを訴えるべき。
裏口入学医師の免許剥奪を！の国民運動の先頭に立てばよいぞ。
僕も裏口入学とか、言ってたら信頼の回復はない。

シリツ医のこれが現実

１次方程式もできないド底辺特殊シリツ医大卒の記録
http://imagizer.imageshack.com/img923/2715/RosCsf.jpg


何度読んでも馬鹿すぎる。
男女別の割合と全体での割合から男女比が計算できるとも思わないとは。
なんでこんなのが大学に入れるわけよ？
裏口入学以外に説明がつく？
中学生でも解ける一次方程式の問題だろ。

それすらできない馬鹿が自信を持って発言。
>患者の男女比が必要なのもわからないのか？
だとさ。
http://imagizer.imageshack.com/img923/9687/zNivZW.jpg


0.2575549
と答を書いてやったら
>単位も書かずに答えだとか…
ド底辺シリツ医大では確率に単位があるらしいぞw
何でこんな馬鹿が大学に入れるわけ？
裏口入学以外に説明がつく？
http://imagizer.imageshack.com/img923/5683/vMw9nv.jpg

算数の掛け算すら怪しいのがシリツ医がこういう事故を起こす。
氷山の一角だと俺は思っている。

http://i.imgur.com/ArPaux9.png

底辺シリツ裏口卒業のガイジ


https://video.twimg.com/ext_tw_video/1222876611202039813/pu/vid/960x540/6bltLho4FsG07Wk2.mp4
山添「つまり趣旨は伝わってるんですよ、総理、いかがですか」
安倍
「ですが、しゅし、趣旨については伝わったという話をいたしましたが、こっ…こっ…ﾁｯ…
このですね、この、ｲﾜﾊﾞｯ、ｻｧﾝｷﾞﾝの……………が、書いた、いるですね、この　ｫｰ
ｷﾞｮｳｾｲｷｶﾝﾉﾎﾕｳｽﾙｼﾞｮｳﾎｳﾉｺｳｶｲﾆｶﾝｽﾙﾎｳﾘﾂﾆﾓﾄﾁﾞｭｲﾃ、名簿全体を、公開されることもありますと
いうこと、では、ｷきていない、ということは申し上げた、わけで、あります
で、そもそもですね！えー、この…この、最初の、ないか、えー、れんらく、うー
事務連絡で、かいたあることもですね、えー、ｶｲ請求の対象！とされることも、あります！
とこうかいて、ｳｪ、ある、わけぇ、で、ｱﾘﾏｼﾃ、対象とされることもありますということとですね、
ｴ、ｶｲｼﾞ、ｴ、ｱﾉﾑ、名簿ｼﾞｪﾝﾀｲを公開される、ということとは、これ、ｴｰ、…f…fｾ、請求の対象ﾄｾ！
対象ですから！対象とされたってこと、でありますから、そこが、ﾏｧ違うということだとｵﾓﾏｽｶﾞ
いずれぇいたしましてもですね、ｲｽﾞﾚｪいたしましても、いわば、ｱ､そうい、う、方、をですね、
ｴｰ、選んでくださいと、いうことについては、ｿﾚｧ、当然、ｴｷﾃｲﾙということは
先ほども申し上げた趣旨、にﾁｭｲﾃﾊ、申し上げている通りであります」
山添「ちょっとはっきり趣旨が分かりませんけども」

N=950
y=c(rep(NA,N-1),30) # c(NA,NA,....,30)
thresh=cbind(rep(15,N),rep(25,N),rep(30,N)) # 15,25,30
ybin=c(rep(1,N-1),3) # 0:[,15), 1:[15,25), 2:[25,30), 3:[30,]
data=list(N=N,y=y,ybin=ybin,thresh=thresh)
inits=list(theta=0.5)
model.string='
model{
for (i in 1:N){
ybin[i] ~ dinterval( y[i] , thresh[i, ] )
y[i] ~ dbin(theta,50)
}
theta ~ dbeta(1,1)I(0.25,1)
}
'

https://image.slidesharecdn.com/stanintermediate-160812024702/95/stan-41-638.jpg?cb=1470970102

>>8

https://kosugitti.github.io/JSSP2018Summer/jssp05.html#%E6%89%93%E3%81%A1%E5%88%87%E3%82%8A%E3%83%87%E3%83%BC%E3%82%BF

fakekit <- function(
prev, # prevalence
sen=0, # sensitivity
spc=1){ # specificity

data.frame(p = prev*sen + (1-prev)*(1-spc),
PPV=sen*prev/(sen*prev+(1-prev)*(1-spc)),
NPV=(1-prev)*spc/((1-prev)*spc+prev*(1-sen)),
precision=(prev*sen+(1-prev)*spc)/
((prev*sen+(1-prev)*spc+(1-prev)*(1-spc)+(prev*(1-sen)))),
pLR=sen/(1-spc),
nLR=(1-sen)/spc,
DOR=pLR/nLR)}

rm(list=ls())
source('tools.R')

Prev2PV <- function(
prev=0.003,
SEN=0.5, # 感度最頻値 beta(13.6991,9.4661)でmode 0.6, sd=0.1
SD1=0.2,# 感度標準偏差
SPC=0.95, # 特異度最頻値 beta(36.172,4.908)でmode 0.9 sd=0.05
SD2=0.05, # 特異度標準偏差
MODE=FALSE, # 最頻値表示,期待値表示
N.ITER=NULL, # JAGSのサンプリング回数
THIN=NULL, # THIN毎にデータ採用
print=TRUE){ # グラフを描くか?
library(rjags)
library(BEST)
sn=Mv2ab(SEN,SD1^2) # Mode variance -> Beta(a,b)のパラメータに
sp=Mv2ab(SPC,SD2^2)
# n : 日々の検査人数配列 x:日々の陽性人数配列
modelstring='
model
{
p = prev*sen + (1-prev)*(1-spc) ;
PPV=sen*p/(sen*p+(1-p)*(1-spc));
NPV=(1-p)*spc/((1-p)*spc+p*(1-sen));
precision=(p*sen+(1-p)*spc)/((p*sen+(1-p)*spc+(1-p)*(1-spc)+(p*(1-sen))));
pLR=sen/(1-spc);
nLR=(1-sen)/spc;
DOR=pLR/nLR;
sen ~ dbeta(sn[1],sn[2])
spc ~ dbeta(sp[1],sp[2])
}
'

writeLines(modelstring,'TEMPmodelj4.txt')
dataList=list(prev=prev,sn=sn,sp=sp)
jagsModel = jags.model( file="TEMPmodelj4.txt" ,data=dataList, n.chains=4,quiet=TRUE)
update(jagsModel)
codaSamples=coda.samples(jagsModel,
variable.names=c("PPV","NPV","precision","pLR","nLR","DOR"),
n.iter=ifelse(is.null(N.ITER),1e5,N.ITER),
thin=ifelse(is.null(THIN),1,THIN))
js=as.data.frame(as.matrix(codaSamples))
par(mfrow=c(3,2))
BEST::plotPost(js$PPV,xlab='PPV',col='pink')
BEST::plotPost(js$NPV,xlab='NPV')
BEST::plotPost(js$precision,xlab='precision')
BEST::plotPost(js$pLR,xlab='pLR')
BEST::plotPost(js$nLR,xlab='nLR')
BEST::plotPost(js$DOR,xlab='DOR')
invisible(js)
}

他の病気で入院してきた患者が自覚症状なしでコロナ持ってくるんだぜ
さらにそれで感染した同僚が無症候でまわりにばら撒くんだぜ

いくらプロでもこんなん防げねーよ全患者・全職員に対し常時完全防備なんてできねーし
なんならいっそ一般病棟よりコロナ病棟入ったほうが全員が最初からわかって対策できるぶん感染リスク少ねーと思うわ

rm(list=ls())

library(rjags)
# model string
modelstring = "
model {
for ( i in 1:N ) {
ybin[i] ~ dinterval( y[i] , thresh[i,] )
y[i] ~ dnorm( mu , tau )
}
tau ~ dt(0,pow(2.5,-2),1)T(0,)
# tau ~ dgamma(.001,.001)
sigma = sqrt(1/tau)
mu ~ dnorm(0,1E-6)
}
"

# Generate some random data Y:
N=500 ; MU=100 ; SD=15
set.seed(19791017) ; Y=MU + SD*scale(rnorm(N))
# Generate cencored data y:
# threshold 0 1 2 3
# censored 94-100 106-
fn <- function(
t1=94,
t2=100,
t3=106){
y=Y=as.numeric(Y)
thresh=matrix(rep(c(t1,t2,t3),N),ncol=3,b=T)
ybin=numeric(N)
idx0=which(Y<t1)
ybin[idx0]=0
idx1=which(t1<=Y & Y<=t2)
ybin[idx1]=1
idx2=which(t2<Y & t3<Y)
ybin[idx2]=2
idx3=which(t3<=Y)
ybin[idx3]=3
y[c(idx1,idx3)]=NA
list(y=y,ybin=ybin,thresh=thresh)
}
re=fn()
data=list(y=re$y,ybin=re$ybin,thresh=re$thresh,N=N)
writeLines(modelstring,con="model.txt")
jagsModel=jags.model('model.txt',data=data,n.chains = 1)
codaSamples=coda.samples(jagsModel,n.iter=1e4,thin=1,
variable.names = c("mu","sigma") )
js=as.data.frame(as.matrix(codaSamples))

summary(js)
layout(matrix(1:2,2))
BEST::plotPost(js$mu,xlab='μ = 100',showMode = T)
BEST::plotPost(js$sigma,xlab='σ = 15',showMode = T)

Natureのこの論文はエクセルデータとRのコードがついていて自分で再現できるので入力の手間が省ける。
https://www.nature.com/articles/s41591-020-0869-5
感染させる確率分布をガンマ分布を平行移動させた分布として想定して、データから最尤法でパラメータ算出している。

そのプログラムをみていくと
感染させる確率分布を
#--- infectiousness, gamma distribution ---
# gpar[1:2]: hyper-parameters (gamma)
# x : infection time of infectee w.r.t onset time of infector
f.Xc = function(x, gpar) { dgamma(x, gpar[1], gpar[2]) }
ガンマ分布にしているけど
一人が一人を感染させるモデルだから
形状パラメータgpar[1]は1で固定、つまり、指数分布だと思うんだがどうだろ？待ち時間の分布と同じ考え。

w.r.t = with reference to

ひとりが別の一人に一回感染させるというモデルなので待ち時間の分布の指数分布でいいと思う。
ガンマ分布の形状パラメータ＝１とおけば指数分布になる。
プログラムを書き直してグラフにすると
https://i.imgur.com/paErvfa.png
発症前に感染させている確率は47％と原著とあまりかわらないが、そのピークは発症1.6日前という結果になった。

算数の掛け算すら怪しいのがシリツ医がこういう事故を起こすわけなんだな。

http://i.imgur.com/ArPaux9.png

私立医大を増やすとこうなる。

算数の掛け算すら怪しいのがシリツ医がこういう事故を起こすわけなんだな。

://i.imgur.com/ArPaux9.png

5卵の名無しさん2020/03/13(金) 07:13:02.02ID:kiP5MRUD
前スレからのド底辺シリツ医への宿題

若い女性研修医（嘘つきかどうかは不明）から
「あなたのいうことが正しければ手コキかフェラをしてあげる、そうでなければセンズリを命じる」と言われた。
フェラをしてもらうには何と言えばいいか？
答だせないってアホなの？　
フェラよりセンズリが好きなの？

997卵の名無しさん2020/04/27(月) 19:05:04.10ID:ek/DjIvE
>>5
答えてやろうｗ
正直に「俺はド底辺私立医大卒の医者です」と言うｗ

前スレの続きだがこれで答えは正解？

阪大卒業しました（恥ずかしいっす）

InstalledPackages=installed.packages()[,'Package']
install.packages(InstalledPackages)

>>21

アホ晒して楽しい？

" https://www.medrxiv.org/content/10.1101/2020.04.16.20067835v1.full.pdf
saliva+ saliva-
swab+ 26 3
swab- 8 1
"
prop.test(c(9,4),c(38,38))
poisson.test(c(9,4),c(38,38),alt='g')

https://www.researchgate.net/publication/340869568_Saliva_Sample_as_a_Non-Invasive_Specimen_for_the_Diagnosis_of_Coronavirus_Disease-2019_COVID-19_a_Cross-Sectional_Study/link/5ea19006458515ec3aff8f36/download

https://www.researchgate.net/publication/340869568_Saliva_Sample_as_a_Non-Invasive_Specimen_for_the_Diagnosis_of_Coronavirus_Disease-2019_COVID-19_a_Cross-Sectional_Study/link/5ea19006458515ec3aff8f36/download

https://i.imgur.com/3GWmZz3.png

を見つけたのでKappa値や信頼区間を出して遊んでみよう〜っと。

>>27
仮想データで遊ぶより、面白いな。

swabと唾液でまあ、結果が合致している。

> print(fit.kappa,pars=c('kappa','PABAK'))
Inference for Stan model: kappa.
4 chains, each with iter=2000; warmup=1000; thin=1;
post-warmup draws per chain=1000, total post-warmup draws=4000.

mean se_mean sd 2.5% 25% 50% 75% 97.5% n_eff Rhat
kappa 0.81 0 0.07 0.65 0.77 0.81 0.86 0.92 3540 1
PABAK 0.93 0 0.02 0.88 0.92 0.93 0.95 0.97 3342 1

JAGSでMCMCしてもほぼ同じ結果。
https://i.imgur.com/XKT53Fy.png

キスも……慶応病院研修医「集団感染」
40人懇親会で何が起きていたか？？？

文春オンライン
https://bunshun.jp/articles/-/37104?page=1

うりゅう爺さん今日も絶好調ｗｗｗｗ

きちんと勉強机の前に座って勉強しないとまた浪人するよｗｗｗｗ

で25浪くらいかｗｗｗｗ

そろそろ就職して年老いた親を安心させないとｗｗｗｗ

医科歯科大卒を自称してもウソはバレバレだからｗｗｗｗ

スレタイ読んでここに投稿するシリツ医ってアホだと思うぜ。

内視鏡バイト先が緊急事態宣言解除まで休診だけど６〜１０割給与がでるのでなんだか悪いな。
バイト医を解雇する病院とかでてきているらしいなぁ。

国立卒は謙虚な人がおおいなぁ。

新型コロナウイルス対策 in 医者板 Part.8
https://egg.5ch.net/test/read.cgi/hosp/1587582729/771

771 名前：卵の名無しさん[] 投稿日：2020/04/28(火) 15:49:26.22 ID:kb7Uq7l8
>>766
＞こういう計算ができる歯医者は国立卒だろうな。

国立の阪大卒の歯医者ですが全く分かりません。（石神井の学院出身）
歯医者の知性の低さを侮らないで下さい。
医学部に行けなかったから歯学部に行ったのでは無く
親の病院を継ぎたく無いから歯学部に入ったリヴァタリアンなのですから
お医者さんのような社会的動機的知性を持っている歯医者なんていませんよ。（ダブルライセンスの先生は別として）
やる気ない人生を楽しむ為に歯医者になったんですから。でもお医者様に関してはリスペクトしてます。

Suppose that an urn contains n1 white balls and n2 black balls. A total of m1 balls are
drawn from the urn without replacement. Then y1, the number of white balls drawn from
the urn, has a hypergeometric distribution
y1 ~ dhyper(n1,n2,m1,1)

The non-central hypergeometric distribution with odds ratio parameter ψ≠ 1 arises when
m1 is the sum of two binomial distributions
m1 ~ sum(y1, y2)
y1 ~ dbin(p1, n1)
y2 ~ dbin(p2, n2)
If m1 is observed, the conditional distribution of y1 has a non-central hypergeometric distribution
with odds ratio parameter
ψ =　p1(1 - p2) / (1 - p1)p2
There are two important things to note. Firstly, the value of m1 is not fixed a priori but is a
random variable. Secondly, the hypergeometric distribution does not depend on the absolute
values of the underlying binomial probability parameters p1; p2 but only their odds ratio .

"
ド底辺シリツ医大にある年の裏口受験票があることが判明した。
枚数は１０００枚以下であることは判明しているが、その数を推定したい。
裏口入学調査員が無作為に１０枚選んで番号を記録して元に戻した。
別の裏口入学調査員が無作為に２０枚選んで番号を記録した。
二人の調査員の記録した番号を照合すると３枚の番号が一致していた。
この情報からこの年の裏口受験証の数の最頻値、期待値とその９５％信頼区間を求めよ。
"

〇〇大学入試問題


国語　次の漢字の読み方を答えよ
血管　細菌　病気

数学
3 +5=
7+2=
6-4=

英語　アルファベットをAからZまで大文字で全て書け


社会
江戸幕府を開いた人は次のうち誰？
A.徳川家康
B.豊臣秀吉
C.織田信長

理科
時速50Kmの自家用車が15Kmを走るのに要する時間を答えなさい。

JAGSもプログラムのコメントに日本語を使うとエラーがでるな。

シリツ医に計算させて医療に従事させるみたいなもんだな。

算数の掛け算すら怪しいのがシリツ医がこういう事故を起こす。

http://i.imgur.com/ArPaux9.png

だってマスク要りません
検査やりません
外国人どんどん来てだもん
初期行動がこれだから、いくら問題になろうとどっかでコソコソ動いて本筋に戻したいのよ準備してたから

ほとぼり冷めたと思ったのか感染研辺り奇妙な事言い出し始めた

俺と同じような計算している人いるんだなぁ。

http://statmodeling.hatenablog.com/entry/covid19-estimate-total-number-of-positives-in-japan

COVID-19 日本国内の潜在的な陽性者数を推定する試み

比較物、用いる価値観によるとしか
実際の医療はパターンばっかりの中に器用さが求められたり、実直さが必要だったりいろんな才能が必要だと思うよ
ゲームが上手いと手術が上手い事が多いかもしれないので私立の方がいいかもしれん
第一出身大学は聞かないとわからない
いちいち入試問題を解かすわけじゃ無いのでね
私立は入試のボーダーは国立より低いだろうけど低いなりに競争はある
私立にも入れない人か、学歴コンプが気にするんだろう

# http://statmodeling.hatenablog.com/entry/covid19-estimate-total-number-of-positives-in-japan
rm(list=ls())
library(rjags)
modelString='
model {
for (a in 1:A) {
n.pos[a] ~ dbinom(p.pos, N[a])
n.sym[a] ~ dbinom(q[a], n.pos[a])
n.asy[a] <- n.pos[a] - n.sym[a]
Y.sym[a] ~ dbinom(p.obs.sym, n.sym[a])
Y.asy[a] ~ dbinom(p.obs.asy, n.asy[a])
Y.sym.DP[a] ~ dbinom(q[a], N.pos.DP[a])
logit(q[a]) <- q.x[a]
}
n.pos.sum <- sum(n.pos)
# positive rate = prevalence*sensitivity + (1 - prevalence)*(1 - specificity)
# prevalence = (specificity + positive rate - 1) / (specificity + sensitivity - 1)
prevalence <- (spc + p.pos -1)/(spc + sen -1)
sen ~ dbeta(sn[1],sn[2])
spc ~ dbeta(sp[1],sp[2])
p.pos ~ dbeta(0.5, 1)
p.obs.sym ~ dbeta(4, 2)
p.obs.asy ~ dbeta(2, 4)
q.x[1] ~ dnorm(0, 1e-4)
q.x[2] ~ dnorm(0, 1e-4)
for (a in 3:A) {
q.x[a] ~ dnorm(2*q.x[a-1] - q.x[a-2], tau)
}
tau ~ dt(0,pow(2.5,-2),1)T(0,)
}
'

# https://www.stat.go.jp/data/jinsui/pdf/202004.pdf

SEN=0.5 ; SD1=0.2
SPC=0.9 ; SD2=0.05
sn=Mv2ab(SEN,SD1^2) # Mode variance -> Beta(a,b)のパラメータに
sp=Mv2ab(SPC,SD2^2)
pop=c(472,506,534,577,639,627,667,747,855,982,861,778,745,843,899,716,537,369,181,52)
A=10
N=(pop[(1:10)%%2==1] + pop[(1:10)%%2==0])*1e4
plot(N,type='h',bty='l')
Y.sym.DP=c(0,2,25,27,19,28,76,95,27,2)
N.pos.DP=c(1,5,28,34,27,59,177,234,52,2)
Y.sym=c(19,20,256,212,296,334,280,212,114,4)
Y.asy=c(6,12,30,23,32,20,21,19,17,6)
plot(0:(A-1)*10,Y.sym.DP/N.pos.DP,type='h',bty='l',xlab='age')
lines((0:9*10)+1.5,Y.sym/(Y.sym + Y.asy),type='h',col='red')
dataList=list(N=N,Y.sym=Y.sym,Y.asy=Y.asy,
Y.sym.DP=Y.sym.DP,N.pos.DP=N.pos.DP,A=A,
sn=sn,sp=sp)
writeLines(modelString,'temp.txt')
jagsModel=jags.model('temp.txt',data=dataList, n.chains=4,n.adapt = 50000)
update(jagsModel,50000)
codaSamples=coda.samples(jagsModel, n.iter=50e4,thin=100,
variable.names = c('n.pos.sum', 'n.pos', 'p.pos', 'p.obs.sym', 'p.obs.asy', 'q'))
mat=as.matrix(codaSamples)
colnames(mat)
js=as.data.frame(mat)
attach(js)
summary(p.pos)
summary(n.pos.sum)

>>41
手先の器用な奴はシリツにもいるが（上皇の執刀医）、
頭が器用なシリツ医はまず、いない。特待生ならいるかもしれんが。
このスレのシリツ医は知恵遅れなのは明らか。

これには同意。

なんで、不正って分かってて退学にならないの？
大分の公務員試験は、本人は全く知らずに、親が頼んでたのが分かっただけで、勤務してたのに、解雇になったのに。
ありえないよ。

東京医大の事件は裏口入学が現在進行形であること如実にしめした事件だよね。
シリツ医の使命は裏口入学撲滅国民運動の先頭に立つことだよ。
裏口入学の学生を除籍処分にしないかぎり、信頼の回復はないね。つまり、いつまで経ってもシリツ医大卒＝裏口バカと汚名は拭えない。シリツ出身者こそ、裏口入学に厳しい処分せよを訴えるべき。
裏口入学医師の免許剥奪を！の国民運動の先頭に立てばよいぞ。
僕も裏口入学とか、言ってたら信頼の回復はない。

僕には友達がいません　誰か友達になってくれないかな〜

by 浩之

感度（30-70％）特異度（90-100%）の確率分布をいれて、
http://statmodeling.hatenablog.com/entry/covid19-estimate-total-number-of-positives-in-japan
の設定を踏襲して、データを更新して
COVID-19 日本国内の陽性率を推定してみると

推定陽性率は
> summary(p.pos) ; HDInterval::hdi(p.pos)[1:2]
Min. 1st Qu. Median Mean 3rd Qu. Max.
0.000031 0.000067 0.004681 0.005229 0.009325 0.016658
lower upper
0.00003082 0.01408833

推定有病率は
> summary(prevalence) ; HDInterval::hdi(prevalence)[1:2]
Min. 1st Qu. Median Mean 3rd Qu. Max.
0.00000 0.00166 0.00529 0.01040 0.01129 0.16078
lower upper
0.0000005693 0.0305640211

ド底辺裏口シリツ医をコロナ病院に徴用すべきだな。便所掃除とかゴミ捨てとかに使えば、少しは裏口シリツ医でも役に立つと思われるかもなぁ。
算数の掛け算すら怪しいから、治療に当たらせては駄目。医療ミスをコロナのせいにしてしまうから。

http://i.imgur.com/ArPaux9.png

さっさと解除しろよ

「緊急事態宣言で失う命は救う命より多い」
http://ikedanobuo.livedoor.biz/archives/52040812.html

https://twitter.com/HY29522792/status/1255422964734070785
https://twitter.com/5chan_nel (5ch newer account)

https://twitter.com/big3dead/status/1255423465202552833
https://twitter.com/5chan_nel (5ch newer account)

読め
https://twitter.com/ikedanob/status/1254976443072737281
https://twitter.com/5chan_nel (5ch newer account)

医師会碌なもんじゃねえな
https://twitter.com/ikedanob/status/1254973021074288640
https://twitter.com/5chan_nel (5ch newer account)

https://twitter.com/SubeSube_Ginko/status/1255144973135048707
https://twitter.com/5chan_nel (5ch newer account)

反日のクズはチョンだけではないようだ
https://twitter.com/ikedanob/status/1255144208433770502
https://twitter.com/5chan_nel (5ch newer account)

https://twitter.com/ikedanob/status/1255041050508288000
https://twitter.com/5chan_nel (5ch newer account)

https://twitter.com/MalusPumila4/status/1255046450817294338
https://twitter.com/5chan_nel (5ch newer account)

適確です
https://twitter.com/suzukimasayuki/status/1255053567578173536
https://twitter.com/5chan_nel (5ch newer account)

リアルでは誰も相手をしてくれません　匿名のネットでは誰かが相手をしてくれるので嬉しいです

by 浩之

https://twitter.com/ShinodaHideaki/status/1255056644339187712
https://twitter.com/5chan_nel (5ch newer account)

>>60
それって某医科歯科爺さんのことですか？ｗ

うりゅうよ、これを読めｗ
https://twitter.com/ikedanob/status/1254964694961381376
https://twitter.com/5chan_nel (5ch newer account)

もっともうりゅうは医者じゃないけどねｗｗｗｗｗｗ

https://twitter.com/Tokyo92907342/status/1254966606075006976
https://twitter.com/5chan_nel (5ch newer account)

僕は福井県筑後市瓜生町に住んでいます　リアルでの友達が欲しいです

並列コードのバグ取りは大変だったな・・・

友達が一人もいないので自分と会話をするしかありません　寂しいです

by 浩之

アホ私立医大に裏口入学した奴は死刑でいいと言っておく

んなのもん正々堂々と試験受けて合格すればいいだろうに

僕の不幸の原因は全て僕のせいです　認めるのは辛いですが、、、

by浩之

誰もうりゅうの不幸の原因なんて興味あらへんがなｗ

「私大医学部は金持ちのバカばかり」伝説はもはや過去の話
「私大医学部は金持ちのバカばかり」伝説はもはや過去の話
「私大医学部は金持ちのバカばかり」伝説はもはや過去の話

https://ent.smt.docomo.ne.jp/article/527223

んなこたーないｗ
アホばっかやでｗ

処女にフェラを仕込むのはさほど難しくないが（N＝４）、アホにアホと認識させるのはむつかしい

自分の不幸をどうしても他人のせいにしてしまうほど　僕は不幸です

by 浩之

自分の不幸をどうしても他人のせいにしてしまうほど　僕は不幸です

by 浩之

以前はよく「社会勉強」に出かけたが今は無理やねｗ
嫁はんが厳しいからのう

恥ずかしながら　僕はモテないので今まで女性と付き合ったことがありません　女性の前でどう振る舞えばいいのかわからないのです　あ、年齢は今50歳台後半です

by 浩之

>>74
臨床をやっている女医は不細工か不倫しているかのどちらかである。
(エビデンスレベルV、個人の経験 n =21)

>>77
当直がなくなると、外泊の口実がなくなるからねｗ

恥ずかしながら　僕はモテないので今まで女性と付き合ったことがありません　女性の前でどう振る舞えばいいのかわからないのです　あ、年齢は今50歳台後半です

by 浩之

>>79
女医と看護師は地雷原ｗ
ドエライ目に会うぞｗ

>>80
学会が口実としては最高ｗ
まずバレないｗ
事前調査をしないと「社会勉強」で地雷を踏むｗ

僕の幸せってなんだろう？僕は今まで何をしてきたのだろうか？失敗してしまった人生

by浩之

人種や民族など本人が選択できないことを根拠とする区別は差別の批判を免れないが、
ド底辺シリツ医大進学は本人の選択だよなぁ。

これな！

不朽の名投稿
>>
私は昭和の時代に大学受験したけど、昔は今よりも差別感が凄く、特殊民のための特殊学校というイメージで開業医のバカ息子以外は誰も受験しようとすらしなかった。
常識的に考えて、数千万という法外な金を払って、しかも同業者からも患者からもバカだの裏口だのと散々罵られるのをわかって好き好んでド底辺医に行く同級生は一人もいませんでした。
本人には面と向かっては言わないけれど、俺くらいの年代の人間は、おそらくは８−９割はド底辺医卒を今でも「何偉そうなこと抜かしてるんだ、この裏口バカが」と心の底で軽蔑し、嘲笑しているよ。
当の本人には面と向かっては絶対にそんなことは言わないけどね。
<<

東京医大の事件は裏口入学が現在進行形であること如実にしめした事件だよね。
シリツ医の使命は裏口入学撲滅国民運動の先頭に立つことだよ。
裏口入学の学生を除籍処分にしないかぎり、信頼の回復はないね。つまり、いつまで経ってもシリツ医大卒＝裏口バカと汚名は拭えない。シリツ出身者こそ、裏口入学に厳しい処分せよを訴えるべき。
裏口入学医師の免許剥奪を！の国民運動の先頭に立てばよいぞ。
僕も裏口入学とか、言ってたら信頼の回復はない。

>>82
学生時代には女子医大の学生には息子が大変お世話になりました。

美人な女医ってたまにいるけどね、いるんだけどさ・・・
トンデモない性癖の持ち主だったな（N＝２）
具体的には書かれへんｗ

自分の不幸をどうしても他人のせいにしてしまうほど　僕は不幸です

by 浩之

>>85
俺の息子は大変に鍛えられましたｗ

人種や民族など本人が選択できないことを根拠とする区別は差別の批判を免れないが、
ド底辺シリツ医大進学は本人の選択だよなぁ。

ド底辺シリツ医大にしか入れない子供だとわかったとき両親はどれだけ悲しんだことだろう。

症例提示

私立医学部医学科ｋａｒｔｅ213 [無断転載禁止]&#169;2ch.net
https://tamae.5ch.net/test/read.cgi/kouri/1454692118/319

319 名前：大学への名無しさん[] 投稿日：2016/02/11(木) 22:52:02.21 ID:50TBuVXh0
恥ずかしながら実はうちの親がそうなんです
私立洗顔でした
本人は隠したがっていたけど祖父が公言してたからね
祖父は息子（うちの父親）の学力など知らずに進路指導で当たり前に自分の出身医大（国立）を希望させ父親の担任に無理だと断言されこれが一生で最大の屈辱だったそうです
祖父が何十年もたってるのに死ぬまで祖母をその事で叱りつけてました
お前がしっかり勉強させないからこうなったっとね
自分は教育に失敗したが祖父の口癖
これが私大医学部ＯＢの悲しい現実
これが現実です

無理をして多浪してまで医学部に入学したがまさか無職でニートになるとは思わなかった　背伸びしちゃったな　今は医者ではなく患者だもんな

by 浩之

医者になりたきゃフツーに受験してフツーに合格すればいいんですよ
地方国立ぐらいで十分だからさ
なにも大金費やして「特別受験」する必要はない
log（対数）を10グラムと読んでしまうアホなんか使い物になるわけあらへん

僕は福井県筑後市瓜生町に住んでいます　リアルでの友達が欲しいです

by 浩之

>>1
入りたくても入れなかったんだね
真っ赤だよ

誰かリアルで友達になってくれませんか？

by 浩之

もうすぐ還暦が来る　余生は5ちゃんねるで書き込んで遊ぶつもりです　皆さん僕とネットで友達になってくださいね　よろしくお願いします

by 浩之

お金は精神障害者手帳でなんとかなるかな？チョット不安

by 浩之

これもシリツの悪寒

誤って10倍量の薬を処方 子供が救急搬送される　院外の薬局が｢投与量が多いのでは｣と問い合わせるも、間違いに気づかず・尼崎の医者 [水星虫★]
https://asahi.5ch.net/test/read.cgi/newsplus/1588064644/

>>91
数字以前に、手コキ問題にも答が出せないのがド底辺シリツ医。
同僚も上司も裏口だらけなのか、本人が馬鹿すぎて相手にされないんだろう。

実情はかかりつけ医が保健所に電話すると病院は紹介される。(保健所は責任回避)
病院はまずCT撮る。肺炎あればPCR検査。なければ帰す。(病院は通達通りで責任回避)
紹介した開業医にはコロナ肺炎あれば結果が保健所から通知あり。
肺炎なけれぱ紹介のガイドラインを書いてある通知が初めの1回目から全部再送られてくる。
よく読めよ、と。

だから誰も悪くないんだよ。


これに似ているなぁ。

ド底辺シリツ医大が悪いのではない、本人の頭が悪いんだ。

うりゅうひろゆき先輩、

国家試験のお勉強、

すすんでますか？

何年も勉強してんだから、

そろそろ.ね…

応援してますよ


それはともかく、

全部、

自演ってことは

ないです…

よ、ね、？

なんか、昔からのこと考えたら

やりそうなんで

>>101
誰だよ、それ？

俺、医科歯科卒だから、神経解剖学の故・萬年 甫教授の飼い猫の名前も覚えているぞ。
絵画も趣味で私の画風はマンネリズムと言って学生を笑わせていたのもいい思い出した。

んで、あんたどこ卒？

え？何書いてるの？
重症じゃないのか？、うりゅう先輩！

じゃ
うりゅうさん、テストだよ

100から7ひいたらいくつかな？

それと、
思いつく野菜を5つ、書いてみようか

内視鏡検査休止（６〜１０割の賃金補償付き）で纏まった時間がとれるで、ベイズの階層モデルの学習中。
コロナ関連で公表された数値を題材にして遊べて( ・∀・)ｲｲ!!
仮想データではどうも味気ないから。

前々スレからのド底辺シリツ医への宿題

若い女性研修医（嘘つきかどうかは不明）から
「あなたのいうことが正しければ手コキかフェラをしてあげる、そうでなければセンズリを命じる」と言われた。
フェラをしてもらうには何と言えばいいか？
答だせないってアホなの？　
フェラよりセンズリが好きなの？

>>103
んで、あんたどこ卒？

もうすぐ還暦が来る　余生は5ちゃんねるで書き込んで遊ぶつもりです　皆さん僕とネットで友達になってくださいね　よろしくお願いします

by 浩之

悪口でもいいですから相手してください　寂しくて寂しくて辛いんです

by 浩之

それより

うりゅう先輩

100から7引くと

いくつかな？

93です　すみません

by 浩之

特に私立大学医卒の人が順調に医師になっていくのが悔しい　私立医卒の人に負けてしまった　悔しくてたまりません

by 浩之

よしよし、いいぞ！うりゅう先輩！

そこから

また

7を引くと？

いくつかな？

86？ちょっと苦しくなってきました　すみません

by 浩之

特に私立大学医卒の人が順調に医師になっていくのが悔しい　私立医卒の人に負けてしまった　悔しくてたまりません

by 浩之

>>

検査結果では、一般市民の百四十七人の４・８％にあたる七人が陽性、医療従事者五十五人のうち９・１％の五人が陽性だった。
市民・医療従事者を合計した二百二人全体では５・９％の十二人（男女とも六人）が陽性だった。

<<

このデータから、一般市民と医療従事者で陽性率に有意差があるかを検定してみよう。

ド底辺シリツ医には無理だろね。

こういう実測データの統計処理は仮想データと違って役立つね。

ド底辺シリツ医が役立たないのは自明だが。

みんなが順調に医師になっていくのがたまらなく妬ましいし悔しくてたまりません
わかってください

by 浩之

うりゅう先輩、

がんばれ

そこから

また

7を引くと…

いくつかな？

難しいの？

時間かけたらダメですよ

わからなくなってきましたよ　電卓使っていいですか？　すみません

by 浩之

>>116
誰だよ、それ？
んで、あんたどこ卒？

俺、医科歯科。入学は二期校時代。

内視鏡バイト先は休診だけど、有給休暇分は１０割以後は６割の給与が補償されて( ・∀・)ｲｲ!!

内視鏡バイト先は休診だけど、有給休暇分は１０割以後は６割の給与が補償されて( ・∀・)ｲｲ!!

ド底辺シリツ医大の面接試験に手コキ問題を出しても正答できない馬鹿だらけなので○×式で試験をすることになった。
15人を面接して40問の問題の正答数が21,17,21,18,22,31,31,34,34,35,35,36,39,36,35
であった。当てずっぽうの正解確率が0.5なので正解率の期待値が0.5以上である受験生を選択したい。
何点以上の受験生を選択すればよいか？

Score = c(21,17,21,18,22,31,31,34,34,35,35,36,39,36,35)
N = length(Score) # number of people : 15
Q = 40 # number of questions : 40
dataList=list(Score=Score,N=N,Q=Q)

modelString='
model{
for(i in 1:N){
z[i] ~ dbern(0.5) # 0:nerd 1:bright
}
nerd <- 0.5
bright ~ dunif(0.5,1) # dbeta(1,1)T(0.5,1)

for(i in 1:N){
theta[i] <- ifelse(z[i]==0,nerd,bright)
Score[i] ~ dbin(theta[i],Q)
}
}
'
writeLines(modelString,'temp.txt')
jagsModel=jags.model('temp.txt',data=dataList,n.chains = 4)
codaSamples=coda.samples(jagsModel,n.iter=1e6,thin=10,
variable.names = c('bright','theta'))
js=as.data.frame(as.matrix(codaSamples))

とりえだった学力でも私立大学卒の人に負けてしまった　もう一度学生からやり直したい

by 浩之

友達が一人もいないので自分と会話をするしかありません　寂しいです

by 浩之

"検査結果では、一般市民の百四十七人の４・８％にあたる七人が陽性、
医療従事者五十五人のうち９・１％の五人が陽性だった。
https://www.tokyo-np.co.jp/s/article/2020043090070748.html

"

"検査結果では、一般市民の百四十七人の４・８％にあたる七人が陽性、
医療従事者五十五人のうち９・１％の五人が陽性だった。
https://www.tokyo-np.co.jp/s/article/2020043090070748.html

"
https://i.imgur.com/kP3EYrt.jpg

リスク比が１を跨ぐから有意差はないな。

>>122
視覚化した方がわかりやすいな。
内視鏡バイトが休診で勉強が捗る。
まあ、寝当直病院の方が飲酒しないからいいんだが、
当直室のPCは遅くてstanでのコンパイルには不向き。

金がないこと　女性にモテないこと　精神病　友達がいないこと　仕事が見つからないこと　で苦しんでいます　誰か助けてください

by 浩之

そんなことより

うりゅう先輩

思いつく野菜を5つ言え

あとから

また

質問しますよ

さあ

はい　お答えします

にんじん　キュウリ　じゃがいも　白菜　タマネギ　です

あっていますか？

by 浩之

# data
pop=c(472,506,534,577,639,627,667,747,855,982,861,778,745,843,899,716,537,369,181,52)
N=(pop[(1:10)%%2==1] + pop[(1:10)%%2==0])*1e4
A=10
plot(N,type='h',bty='l')
Y.sym.DP=c(0,2,25,27,19,28,76,95,27,2)
N.pos.DP=c(1,5,28,34,27,59,177,234,52,2)
Y.sym=c(19,20,256,212,296,334,280,212,114,4)
Y.asy=c(6,12,30,23,32,20,21,19,17,6)


COVID-19 日本国内の潜在的な陽性者数を推定する試み
http://statmodeling.hatenablog.com/entry/covid19-estimate-total-number-of-positives-in-japan
のデータを更新して実行してみた。

model {
for (i in 1:A) { # 1:age 0-10 2:age 10-19
n.pos[i] ~ dbinom(p.pos, N[i])
n.sym[i] ~ dbinom(q[i], n.pos[i]) # q[i] : Pr[symptomatic|PCR+]
n.asy[i] <- n.pos[i] - n.sym[i]
Y.sym[i] ~ dbinom(p.obs.sym, n.sym[i])
Y.asy[i] ~ dbinom(p.obs.asy, n.asy[i])
Y.sym.DP[i] ~ dbinom(q[i], N.pos.DP[i])
logit(q[i]) <- q.x[i]
# q[i] <- ilogit(q.x[i])
}
n.pos.sum <- sum(n.pos)
p.pos ~ dbeta(0.5,0.5) # original : dbeta(0.5, 1)
p.obs.sym ~ dbeta(4, 2) # mean 0.667 sd 0.1781742
p.obs.asy ~ dbeta(2, 4) # mean 0.333 sd 0.1781742
q.x[1] ~ dnorm(0, 1e-4) # N(m=0,sd=100)
q.x[2] ~ dnorm(0, 1e-4)
for (j in 3:A) {
q.x[j] ~ dnorm(2*q.x[j-1] - q.x[j-2], tau)
}
tau ~ dt(0,pow(2.5,-2),1)T(0,)
# original : tau = pow(sigma,-2) ; sigma ~ dunif(0,5)
}

> summary(n.pos.sum) ; HDInterval::hdi(n.pos.sum)[1:2]
Min. 1st Qu. Median Mean 3rd Qu. Max.
3839 5040 5788 16264 7200 392979
lower upper
3839 79969

現時点での推定陽性者数の期待値は16000人という結果が得られた。

一般人はド底辺シリツ医大生をご存知ないのだろうか？

【速報】#安倍晋三内閣総理大臣 「だって…こ、これに書いて、これに、これに、これに書いてないじゃないですか！？」 ★6 [１号★]
https://asahi.5ch.net/test/read.cgi/newsplus/1588244172/988

988 名前：不要不急の名無しさん[] 投稿日：2020/04/30(木) 21:47:38.22 ID:tHduSa9z0
安倍は正直そこらの大学生より遥かにアホなんだから安倍に質問するとか無駄なことをやめろよ
まともな奴が官僚に直接質問する場にしろ

1 名前:赤紙が来た :2020/04/26(日) 14:06:53.09 ID:sCYq+H5P
【訓示】
現在、激烈極まる対新型コロナ戦の真っただ中に在る事は諸君らも知る通りである！
しかし、国民の健康を防衛するための医師の絶対数が不足しておる！
特に、発熱センターで検体採取など雑務をこなし、最前線で働く医師が不足しておる！

そこで、この度わが国民医療防衛隊では、最前線に赴く医師の志願を広く求める事となった！
皆存分に「発熱センターへの出務を希望する」に〇を付け提出して頂きたい！

我々医師会上層部と国民は諸君ら開業医にも期待しておる！
余談だが、健康の保証はできない、感染時の補償も不明、その後の保険についても不透明である！

異常！

【新型コロナ】「赤紙」報告本部【召集令状】
https://egg.5ch.net/test/read.cgi/hosp/1587877613/

>>137
ド底辺シリツ医を徴用すべき。

但し、頭脳労働に従事させてはならない。喀痰吸引とか体位変換要員でいい。


算数の掛け算すら怪しいのがシリツ医がこういう事故を起こす。


http://i.imgur.com/ArPaux9.png

上級国民の中から安倍よりアホを探した結果だという。
ガイジじゃね？
https://i.imgur.com/IkE6Hnt.gif

ド底辺シリツ医を徴用すべき。

但し、頭脳労働に従事させてはならない。
喀痰吸引とか体位変換要員でいい。

算数の掛け算すら怪しいシリツ医がこういう事故を起こす。


://i.imgur.com/ArPaux9.png

専門医自慢のウチの私立卒医の会話。

「全盲なのに眼振が出てる患者がいるんだけど何故なんだろう？」

うりゅうひろゆき先輩！

寝ないで何してんですか？！

俳諧したら、

うりゅうウイルスを

撒き散らかしますよ！

うりゅうひろゆき先輩、

もう

還暦すぎてんだから！


そんだけ、元気なら、


じゃあ、そんな事より

画用紙に、クレヨンで

時計の絵を描けるかな？

針は

10時10分

描いてみようか？

恥ずかしながら　僕はモテないので今まで女性と付き合ったことがありません　女性の前でどう振る舞えばいいのかわからないのです　あ、年齢は今50歳台後半です 友達が一人もいないので自分と会話をするしかありません　寂しいです

by 浩之

うりゅう先輩、

そんなことどうでもいいよ

画用紙用意しなさいよ

時計のお絵かきしなさい

みんなが順調に医師になっていくのがたまらなく妬ましいし悔しくてたまりません
わかってください

by 浩之

A B C D E F G H I J K L M N O P Q R S T
P1 1 1 1 1 0 0 1 1 0 1 0 0 NA 0 0 1 0 1 0 0
P2 0 1 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
P3 0 0 1 0 0 0 1 1 0 0 0 0 1 0 0 0 0 0 0 0
P4 0 0 0 0 0 0 1 0 1 1 0 0 0 0 0 0 0 0 0 0
P5 1 0 1 1 0 1 1 1 0 1 0 0 1 0 0 0 0 1 0 0
P6 1 1 0 1 0 0 0 1 0 1 0 1 1 0 0 1 0 1 0 0
P7 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0
P8 0 0 0 0 NA 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
P9 0 1 1 0 0 0 0 1 0 1 0 0 1 0 0 0 0 1 0 1
P10 1 0 0 0 0 0 1 0 0 1 0 0 1 0 0 0 0 NA 0 0

の欠測データが１である確率とその95%信頼区間を求めよ。

model{
for(i in 1:np){
for(j in 1:nq){
k[i,j] ~ dbern(theta[i,j])
theta[i,j]=p[i]*q[j]
}
}
# prior
for(i in 1:np){
p[i] ~ dbeta(0.5,0.5)
}
for(j in 1:nq){
q[j] ~ dbeta(0.5,0.5)
}
# posterior
for(i in 1:np){
for(j in 1:nq){
predK[i,j] <- p[i]*q[j]
}
}
}

周辺尤度を出して
https://i.imgur.com/kLWdBcB.png
https://i.imgur.com/taLUU3H.png
これから求めた結果はこうなるな。

https://i.imgur.com/BRAI9A6.png

点推定での解説
The inferred probability for M-1 being
correct is approximately 0.74. For item-person combination E-8 the reverse holds.
The item is difficult and the participant is inattentive. Consequently, the inferred
probability for E-8 being correct is approximately 0.01. Finally, combination R-10
is middle-of-the-road on both dimensions, and the inferred probability for it being
correct is 0.41.
These inferred probabilities are directly related to the knowledge about
each participant and item. In fact, if you multiply the estimated p’s and q’s
you can recover the inferred probabilities. For instance, p[1] is approximately
0.88, and q[13] (the M) is approximately 0.84. The multiplication of these
probabilities yields .74. The same calculations may be performed for the other
missing item-participant calculations.

次にこれを実装

The definition of the accuracy for a person on a question in terms
of the product ij = piqj is very simple to understand, but other models
of the interaction between person ability and question difficulty are used in
psychometric models. For example, the Rasch model (e.g., Andrich, 1988)
uses ij = exp (pi ? qj) / (1 + exp (pi ? qj)). Change the graphical model to
implement the Rasch model.

>>150
モデルを変えると随分と分布が変化するなぁ。
https://i.imgur.com/EajssHe.png

まあ、ド底辺シリツ医はどこでもドアホのままだが。

ニュー速＋などで、医師をコロナ診療に徴用せよという一般人からの意見が多いので
俺は、こういうレスをつけている。

裏口容疑者のド底辺シリツ医を徴用すべき。

但し、頭脳労働に従事させてはならない。

喀痰吸引とか体位変換要員でいい。

算数の掛け算すら怪しいシリツ医がこういう事故を起こす。


://i.imgur.com/ArPaux9.png

特に私立大学医卒の人が順調に医師になっていくのが悔しい　私立医卒の人に負けてしまった　悔しくてたまりません

by 浩之

医師8(P1-P8)人に8問の問題(A-H)を解いてもらった。
問題は医学領域の問題と裏口関連の問題である。
国立医と裏口シリツ医がいるとする。
国立医は医学領域問題への正解率は高いが裏口関連の問題には正解率が低い。
裏口シリツ医は医学領域の問題の正解率は低いが裏口関連の問題には正解率が高い。

以下が成績である。P1は裏口シリツ医でAは裏口関連問題であることが判明している。

A B C D E F G H
P1 1 0 0 1 1 0 0 1
P2 1 0 0 1 1 0 0 1
P3 0 1 1 0 0 1 0 0
P4 0 1 1 0 0 1 1 0
P5 1 0 0 1 1 0 0 1
P6 0 0 0 1 1 0 0 1
P7 0 1 0 0 0 1 1 0
P8 0 1 1 1 0 1 1 0

P1の他に裏口シリツ医は誰か、Aの他に裏口関連問題はどれかを答えよ。

いよっ　自称医科歯科卒の偽医者さん　まだお仕事は見つかりませんか？

僕は福井県筑後市瓜生町に住んでいます　リアルでの友達が欲しいです

by 浩之

>>101
誰だよ、それ？

俺、医科歯科卒だから、神経解剖学の故・萬年 甫教授の飼い猫の名前も覚えているぞ。
絵画も趣味で私の画風はマンネリズムと言って学生を笑わせていたのもいい思い出した。

んで、あんたどこ卒？

俺、医科歯科卒だから、神経解剖学の故・萬年 甫教授の飼い猫の名前も覚えているぞ。
絵画も趣味で私の画風はマンネリズムと言って学生を笑わせていたのもいい思い出した。


w

もうすぐ還暦が来る　余生は5ちゃんねるで書き込んで遊ぶつもりです　皆さん僕とネットで友達になってくださいね　よろしくお願いします　恥ずかしいですが国試合格できなくて　 医師にはなれなかったですけど

by 浩之

うりゅうひろゆき先輩、

またおかしなこと書き出したな

いつも同じ、ネコを解剖したが

教授がどうとか、


てかどうでもいいよ、そんな話

そんなことより、

うりゅうさん、

昨日、聞いた野菜、

なにを言ったか、思い出してみようよ

順番にね

焦んないで、

ゆっくりでもいいからね

自分で言ったんだよ？

できるよね？

すみません　

じゃがいも　にんじん　キャベツ　タマネギ　白菜

ではなかったでしょうか？チョット違うような気もしますが、、

by 浩之

しかし、医科歯科卒と教えたら発狂するみたいだな。

何人（なにじん）ですかと聞かれたら、俺は普通に日本人です　と答えられる。
べつに自慢でもなんでもなくて、日常会話のレベル。

臨床医やってる時点で負け犬のような気が・・・・・

Alzheimer’s recall test cheating

Ifthe focus is purelyon prediction, other statistical approaches, including especially ones that combinethe best aspects of generative and discriminative modeling, may be superior (e.g.,Lasserre, Bishop, & Minka, 2006).
The advantage of the generative model is in providing details about the underly-ing processes assumed to produce the data, particularly by quantifying uncertainty.

>>162
接客業は賎業。
とりわけド底辺シリツ医はナマポからも蔑まれる。

アルツハイマー患者の再生テストのいかさま　というモデルを学習中。

model{
phi ~ dbeta(5,5) # non-informative prior for grouping
for(i in 1:np){
grp[i] ~ dbern(phi) # 0:bona fide 1:cheating
score[i] ~ dbin(theta[i],nq)
thetatmp[i] ~ dbeta(alpha[grp[i]+1],beta[grp[i]+1])
theta[i] <- max(0.01, min(0.99,thetatmp[i])) # theta [0.01-0.99]
}
alpha[1] <- mubon*lambdabon # a/(a+b) * (a+b) == a
beta[1] <- (1-mubon)*lambdabon # (1-a/(a+b)) * (a+b) == b

logit(muche) <- logit(mubon) + mudiff # cheating uses more words
alpha[2] <- muche*lambdache # muche : mean of cheating
beta[2] <- (1-muche)*lambdache # lambdache : lambda for cheating

mubon ~ dbeta(1,1)
mudiff ~ dnorm(0,0.5)T(0,)
lambdabon ~ dunif(5,40)
lambdache ~ dunif(5,40)

for(i in 1:np){
pct[i] <- ifelse(grp[i]==truth[i],1,0)
}
pc <- sum(pct)
}

ここ最近
明らかにうりゅうひろゆき先輩以外の奴が入ってきて
先輩がやり込められているな
誰かも書いていたが、
やはりこのスレの前半は
全部同じオッサン、つまり、うりゅうひろゆき先輩が
全部自演してたことが判明した

par(cex.main = 1.5, mar = c(5, 6, 4, 5) + 0.1, mgp = c(3.5, 1, 0), cex.lab = 1.5,
font.lab = 2, cex.axis = 1.3, bty = "n", las=1)

plot(score,summary(samples)$summary[paste("grp[", 1:np, "]", sep=""), 'mean'], ylim=c(0,1),
xlim=c(0,nq), lwd=2, pch=4, xlab= "Number of Items Recalled Correctly",
ylab="Cheater Classification")

summaryで表示されたmeanを使ってグラフ化かぁ。

ms=rstan::extract(samples)
mean.grp=apply(ms$grp,1,mean)
plot(score,mean.grp)
とやっていたけど、1行で済むなぁ。可読性は悪いが。

他人のコードを読むのは勉強になる。
内視鏡休診でまとまった時間がとれるので統計モデルの勉強中.
寝当直勤務の方が飲酒しないから捗るんだがメモリ4G corei3なので処理の待ち時間で思考が中断されるのが難点。

ド底辺シリツ医ってコロナ待機で纏まった時間がとれるのに勉強に充てないんだなぁ。
やっぱりセンズリ専従医みたいだね。

宿題しないと安倍晋三みたいなアホになるぞいないの？
フェラよりセンズリが好きなの？

前々スレからのド底辺シリツ医への宿題

若い女性研修医（嘘つきかどうかは不明）から
「あなたのいうことが正しければ手コキかフェラをしてあげる、そうでなければセンズリを命じる」と言われた。
フェラをしてもらうには何と言えばいいか？


答だせないってアホなの？　
フェラよりセンズリが好きなの？
助言してくれる友人も上司もいないの？
みんな裏口シリツ医？

んで、あんたどこ卒？

うりゅうひろゆき先輩の

ますます重症だな

治癒効果が見られんなぁ

またおなはしもどすけど

質問しますよ

１００から７ひいたら、

いくつかな？

rm(list=ls())

library(rstan)
Sys.setenv(LOCAL_CPPFLAGS = '-march=corei7 -mtune=corei7')
rstan_options(auto_write = TRUE)
options(mc.cores = parallel::detectCores())

options(scipen = 5)

model = "
// Cheating Latent Mixture Model
data {
int<lower=1> nq; // number of questions
int<lower=1> np; // number of people
int<lower=1,upper=nq> score[np];
int<lower=0,upper=1> truth[np]; // 1:cheating 0:bona fide
}
parameters {
real<lower=0,upper=1> phi; // grouping probability ~ bernoulli[phi]
real<lower=0,upper=1> mubon; // mean of bona fide
real<lower=0> mudiff; // == logit(muche) - logit(mubon)
real<lower=5,upper=50> lambdabon; // shape1+shape1 in beta distuributui
real<lower=5,upper=50> lambdache; // bon : bona fide, che : cheating
vector<lower=0,upper=1>[np] theta; // ~ dbeta(nq,theta)
}

transformed parameters {
vector[2] loglik[np]; // log-likelihood
vector<lower=0>[2] alpha; // shape1 parameter
vector<lower=0>[2] beta; // shape2 parameter
real<lower=0,upper=1> muche; // mean of cheating score

// Additivity on Logit Scale
//muche = inv_logit(logit(mubon) + mudiff); // mudiff = logit(muche) - logit(mubon)
muche = exp(logit(mubon) + mudiff)/(1 + exp(logit(mubon) + mudiff));

// Transformation to Group Mean and Precision
alpha[1] = mubon * lambdabon;
beta[1] = lambdabon * (1 - mubon);
alpha[2] = muche * lambdache;
beta[2] = lambdache * (1 - muche);

// Data are Binomial with Rate Given by
// Each Person's Group Assignment
for (i in 1:np) { // log1m(x)==log(1-x)
loglik[i,1] = log(1-phi) + beta_lpdf(theta[i] | alpha[1], beta[1]);
loglik[i,2] = log(phi) + beta_lpdf(theta[i] | alpha[2], beta[2]);
}
}

model {
// Priors
mubon ~ beta(1, 1); // can be removed
mudiff ~ normal(0, 1 / sqrt(.5))T[0,]; // Constrained to be Positive
// Relatively Uninformative Prior on Base Rate
phi ~ beta(5, 5);

for (i in 1:np)
target += log_sum_exp(loglik[i]);
// increment_log_prob(log_sum_exp(loglik[i]));

score ~ binomial(nq, theta);
}
generated quantities {
int<lower=0,upper=1> grp[np];
real pc;
vector[np] pct;

for (i in 1:np) {
vector[2] prob;
prob = softmax(loglik[i]);
// Each Person Belongs to One of Two Latent Groups
grp[i] = bernoulli_rng(prob[2]);
// Correct Count
pct[i] = (grp[i]==truth[i])? 1:0 ;
}
pc = sum(pct);

}"

川崎卒はダメ

stanは変数宣言も面倒、まあ、これはCで慣れているけど
離散量のMCMCに対数尤度を書くのは面倒。
,　が｜だったり、_rngを付加とか、
そのうち　deprecated とかのメッセージで新形式に更新しろと、手間がかかる。
整数を推測するにはJAGSの方がいい。
まあ、WINBUGの文法を踏襲しているから、

>>173
川崎は学生の不祥事は多いが医者の不祥事は少ない印象。
このスレのド底辺シリツ医は聖魔か藤田だと思う。

>>169
もっと、賢げな問題にしたらどう？
こういうのとか、７の１００乗の末尾１０桁の数字を答えよ。

センズリ専従医じゃ無理なのか？
助言してくれる友達もいないの？？

んで、あんたどこ卒？

ここで無料で読める分は読破してしまった。
https://bayesmodels.com/
緊急事態宣言は延長されるみたいだし、Kindle版があるみたいだから買うかな。
本文は斜め読みして、プログラムコードを読む方が速いし、電子図書の方が入力が楽。
まあ、サポートの終了したWINBUGのコードだから、JAGSやstanに移植するのがプログラムの勉強にもなる。
まとまった時間が取れそうだから丁度いいな。

内視鏡をしばらく触ってないと間隔が鈍りそう。
ド底辺シリツ医は毎日センズリらしいが。

The marginal likelihood is computed by averaging the likelihood across the parameter
space, with prior probabilities acting as averaging weights. Thus, in order
to determine how well a model predicted the data, we need to take into account all
predictions that the model made and weight these by their prior probability. We
can restate this mathematically by saying that the marginal likelihood is obtained
by averaging out the model parameters in accordance with the law of total probability.

p (D | M1) =∫p (D | θ,M1) p (θ | M1) dθ.

うりゅうひろゆき先輩

次のお仕事

画用紙を用意しなさい

お絵かきの

時間ですよ

時計の絵を描けるかな

そして

針は

10時10分

かけるかな？

緊急事態宣言が延長になったら、内視鏡バイトも休診だけど給与は６〜１０割は給付される。
これまでの実績がこういうとこで活きてくるなぁ。
ここのド底辺シリツ医は解雇されんじゃないかな？

>>175
提供卒は意外にまとも

こういう投稿があったから、

＞地方住まいの開業医たちを 応援に来させたいのかもしれんね

このレスを貼ってきた。

ド底辺シリツ医を徴用すべき。
但し、頭脳労働に従事させてはならない。喀痰吸引とか体位変換要員でいい。
算数の掛け算すら怪しいのがシリツ医がこういう事故を起こす。

http://i.imgur.com/ArPaux9.png

>>182
このスレのド底辺シリツ医は聖魔か藤田じゃないかな？
あんまり、ド底辺シリツ医と仕事をする機会がなかったから判断材料に乏しい。
まあ、理一に受かるような頭脳をもったシリツ医は思いつかん。

医師のふりぐらいさせてくださいよ　医師になるのが夢だったんですから　お願いします

by 浩之

副院長の息子が医科歯科の後輩というのもプラスに作用したなぁ。

英語を書いたり　統計を書いたりするのは　少しでも医師らしく　また　頭が良い雰囲気を出すためなんです　実際の僕は　国試に合格できなくて医師になれなかった馬鹿ものです　ほんとうにごめんなさい　本当の僕は無職でニートの社会的底辺です

by 浩之

しかし、ド底辺シリツ医って医科歯科卒と教えたら発狂するみたいだな。

何人（なにじん）ですかと聞かれたら、俺は普通に日本人です　と答えられるだろ。
べつに自慢でもなんでもなくて、日常会話のレベル。

しかし、ド底辺シリツ医って医科歯科卒と教えたら発狂するみたいだな。





w

恐れていたことが現実になったな。

【埼玉】骨折で手術の女性死亡、院内感染か　クラスター発生とみられる病院に入院 [ばーど★]
https://asahi.5ch.net/test/read.cgi/newsplus/1588378595/1

COVID-19時代に必要な三原則

（１）女をみたら妊娠と思え。

（２）患者をみたらコロナと思え。

（３）シリツ医をみたら裏口医と思え。

裏口シリツ医をCOVID-19の診療にあたらせればいい。

但し、頭脳労働に従事させてはならない。喀痰吸引とか体位変換要員でいい。

算数の掛け算すら怪しいのがシリツ医がこういう事故を起こす。

http://i.imgur.com/ArPaux9.png

馬鹿を医師から淘汰するのに丁度いいと思う。

医師になれた人が羨ましい　僕ももう少し頭がよければなー　叶わぬ夢だった

by 浩之

俺、医科歯科卒だから、神経解剖学の故・萬年 甫教授の飼い猫の名前も覚えているぞ。
絵画も趣味で私の画風はマンネリズムと言って学生を笑わせていたのもいい思い出した。







w

>>194
萬年先生は絵が上手だったね。
医科歯科の卒業生なら知っているよ。
毎年、同じ冗談を言っていたから先輩から伝わっていて、
さあ出るぞと学生も身構えて笑っていたな。
でたー、て感じで。

俺、医科歯科卒だから、神経解剖学の故・萬年 甫教授の飼い猫の名前も覚えているぞ。
絵画も趣味で私の画風はマンネリズムと言って学生を笑わせていたのもいい思い出した。







w

リアルでは誰も相手をしてくれません　匿名のネットでは誰かが相手をしてくれるので嬉しいです

by 浩之

うりゅうひろゆき

先輩

症状酷くなってないか？

妙なこと書いていないで

時計のお絵かき

しなさいよ

それで、

画用紙、

買ってきた？

いえ　まだ買ってきていません　すみません　金がないんです

by 浩之

■ 神奈川県 ３人死亡 新たに15人感染 計1038人に 新型コロナ

5月1日 18時48分

横浜市は５月１日、新型コロナウイルスへの集団感染が起きている医療機関で感染したとみられる70代と80代の男女合わせて３人が死亡したと発表しました。
また、神奈川県内では合わせて男女15人がウイルスに新たに感染したことが確認されたと発表されました。
横浜市によりますと、死亡したのは80代の男性２人と70代の女性の合わせて３人で、４月28日から29日にかけて、新型コロナウイルスによる肺炎で死亡したということです。

３人は、いずれも旭区にある「聖マリアンナ医科大学 横浜市西部病院」に、別の病気で入院していましたが、病院内で感染したとみられるということです。

また横浜市は、この病院に勤める女性の看護師３人や小学生の女児を含む合わせて男女７人の感染が、新たに確認されたと発表しました。

さらに神奈川県は、小田原市や湯河原町などに住む30代から90代以上の男女４人の感染が確認されたと発表し、川崎市は入院患者の感染が相次いでいる川崎協同病院の入院患者を含む男女３人の新たな感染が確認されたと発表しました。

そして、茅ヶ崎市も市内に住む50代の小売店従業員の女性１人の感染確認を発表しました。

５月１日、神奈川県内ではこれまでに合わせて15人の感染の確認が発表されたことになります。

ＮＨＫのまとめでは、神奈川県内で発表された感染者は、これで1038人となり、このうち37人が死亡しています。

>>200
ド底辺シリツ医は玉よけにも使えんのか？
やっぱりセンズリ専従じゃないかよ。

Holmes and Watson are involved in a rather simple game of darts,
in which, with each dart, the player tries to score as many points as possible.
The maximum score per dart is 60, and the minimum score is 0 (when the
dart lands outside the board). After 5 darts, Holmes scored {38, 10, 0, 0,
0} and Watson scored {20,20,20,18,16}.
How do you determine who is the
better player? Consider another game of darts, but now one of the players
gets to throw 50 times instead of 5. Explain how this scenario shows the
importance of averaging instead of maximizing in order to penalize complexity.

こちら、埼玉医大ですが何か？
w

>>
東京女子医大病院は４月６日付で、外勤中止を求める文書を医師に配布した。
大学病院で感染者を受け入れる日大は４月１０日、外勤中止を派遣先に告知。
順天堂大や帝京大は「状況に応じ（中止を）検討する」とした。
<<

さっさと、外勤禁止をしないと院内感染が拡大する一途だと思う。
発症２日前に感染させている確率が一番高いという統計解析がでているのだから。

ド底辺シリツ医は
裏口シリツ医をCOVID-19の診療にあたらせればいい。
但し、頭脳労働に従事させてはならない。喀痰吸引とか体位変換要員でいい。
算数の掛け算すら怪しいのがシリツ医がこういう事故を起こす。
://i.imgur.com/ArPaux9.png
馬鹿を医師から淘汰するのに丁度いいと思う。

>>203
埼玉医大卒なら、手コキ問題を解けるはずもないだろ。
即答して俺の認識を変えさせてみてくれ。

うりゅうひろゆき先輩

時計のお絵かきだよ

画用紙用意しなさい

10時10分の絵を描く事

ココは

うりゅうひろゆきの

負のオーラしか

ないな

コロナ騒動の

相乗効果で

なんか

きもち

悪い

笑

インバウンド、爆買い、民泊、アベノミクスの悪政の逆回転が始まっただけ。日本は中華支配から逃れることができる。まさにコロナは神風。

休みだ〜

萬年先生は絵が上手だったね。
医科歯科の卒業生なら知っているよ。
毎年、同じ冗談を言っていたから先輩から伝わっていて、
さあ出るぞと学生も身構えて笑っていたな。
でたー、て感じで。







w w w

ちょこっと　恥ずかしいね　


w

ココ覗いてたら、

うりゅうひろゆき先輩の

国試浪人が

うつりそうで

ホンマに

ヤダな

せっかくの

ツキが

どんどん

なくなりそう

です

わ

マジです

よ

いや、

ほんと

すみません　すみません　リアルでは誰も相手をしてくれません　匿名のネットでは誰かが相手をしてくれるので嬉しいです

by 浩之

萬年先生は絵が上手だったね。
医科歯科の卒業生なら知っているよ。
毎年、同じ冗談を言っていたから先輩から伝わっていて、
さあ出るぞと学生も身構えて笑っていたな。
でたー、て感じで。

強烈ですね　ちょこっと　寒気がします


by 浩之

しばらく内視鏡を握ってないなぁ。
ここのド底辺シリツ医って首になったんだろうな。

厚労省・専門家会議がPCR検査絞るように指示を出す
↓
ほとんどの自治体が指示を守って検査絞る
↓
一部の自治体の首長が検査しっかりするように指示、その自治体だけ検査しっかりする
↓
検査しっかりするように指示する首長が増えてくる

自民系保守政治家の場合、政府の指示守って検査数増やさないことが多い気がする
ただし和歌山県とかは自民系知事なのに検査しっかりしてる

地方で車社会なら、検査絞っても勝手に終息するが、
地方でも電車があって繁華街が発展してるところは、
勝手に終息せずに、検査絞ったところで感染拡大していく

国難から日本を守ったの北条時宗と菅直人。
国難を春節中国人ウェルカムで招いたのが安倍。
原発の冷却電源喪失は起こらないという神話を政府の公式見解にしたのも安倍。

熊大医学部卒の高木化け物不細工日本で一番ブス化け物不細工ドブス超絶不細工女医よりマシ

ここで盛んに書き込まれてる、うりゅうひろゆきとかいう人、有名人なの？

>>205
やはり、ド底辺シリツ医には手コキ問題は答えられず。
ここに書いているド底辺シリツ医ってすでに首になってんじゃないかなぁ？
いくらなんでもカンフェレンス参加者に答を教えてもらえるだろうと思う。シリツ医大でも教員は国立卒のことがおおいし。
帝京は東大卒が多いね。

>>219
ド底辺シリツ医の脳内仮想人物じゃないの？
なんのソースも出してないし。

啓発活動建設的提案を投稿してきた。

これ！
https://asahi.5ch.net/test/read.cgi/newsplus/1588408955/867

(2) 要請を受けていなくても、休業手当について60％を超えて支給する場合には、その部分に係る助成率を100％にする
こととしました。
この制度のおかげで今月は不労所得になりそうだが、内視鏡できないのはつまらんな。

アビガン希望者に投与できるようにするだけで経済は回り、自殺者も待機死者もなくなる。

それをしないのは改憲、戦争へ持っていきたいから。アメリカ・ディープステートは米中戦争を望んでいる。あべはディーステートの子分。
改憲して自衛隊をアメリカに売り渡すつもり。

種苗法も改悪して、日本人から自給自足も奪おうとしている。
あべに賛同してるのは先祖が半島系や戦後密入国してきた者。
マスコミ、政界、経済界に多く進出している。

もうすぐ還暦が来る　余生は5ちゃんねるで書き込んで遊ぶつもりです　皆さん僕とネットで友達になってくださいね　よろしくお願いします
by 浩之

医師のふりぐらいさせてくださいよ　医師になるのが夢だったんですから

by 浩之

みんなが順調に医師になっていくのがたまらなく妬ましいし悔しくてたまりません
わかってください

by 浩之

>>226
ド底辺シリツ医は首になって暇なんだな。

日当直中だけど、新患は診なくてすむので楽。

同一法人の老人ホームの入所者が帯状疱疹でやってきたくらい。

アメナメビルが院内にあってよかった。

ここに書いているド底辺シリツ医ってすでに首になってんじゃないかなぁ？
いくらなんでもカンフェレンス参加者に答を教えてもらえるだろうと思う。

俺は、今月は不労所得だけで食えるが、内視鏡を握れないのはつまらんな。

医師のふりぐらいさせてくださいよ　医師になるのが夢だったんですから

by 浩之

うりゅうひろゆきさん、

ますます重症ですってね

お絵かきもしてないし

ダメですね

ほんと

じゃ、こないだ

思いつく野菜

5つ挙げなさい

って言ったこと

覚えてますか？

同じの思い出して

もういちど挙げなさい

同じのですよ
時計のお絵かきは

0点ってことで

しかたないからな

しかしこりゃ困ったなあ

はい　すみません

じゃがいも　キャベツ　キュウリ　にんじん　イチゴ

です

違うような気もしますが

ご教授お願いします

by 浩之

ド底辺シリツ医は
裏口シリツ医をCOVID-19の診療にあたらせればいい。
但し、頭脳労働に従事させてはならない。喀痰吸引とか体位変換要員でいい。
算数の掛け算すら怪しいのがシリツ医がこういう事故を起こす。
://i.imgur.com/ArPaux9.png
馬鹿を医師から淘汰するのに丁度いいと思う。

人種や民族など本人が選択できないことを根拠とする区別は差別の批判を免れないが、
ド底辺シリツ医大進学は本人の選択だよなぁ。

これな！

不朽の名投稿
>>
私は昭和の時代に大学受験したけど、昔は今よりも差別感が凄く、特殊民のための特殊学校というイメージで開業医のバカ息子以外は誰も受験しようとすらしなかった。
常識的に考えて、数千万という法外な金を払って、しかも同業者からも患者からもバカだの裏口だのと散々罵られるのをわかって好き好んでド底辺医に行く同級生は一人もいませんでした。
本人には面と向かっては言わないけれど、俺くらいの年代の人間は、おそらくは８−９割はド底辺医卒を今でも「何偉そうなこと抜かしてるんだ、この裏口バカが」と心の底で軽蔑し、嘲笑しているよ。
当の本人には面と向かっては絶対にそんなことは言わないけどね。
<<

東京医大の事件は裏口入学が現在進行形であること如実にしめした事件だよね。
シリツ医の使命は裏口入学撲滅国民運動の先頭に立つことだよ。
裏口入学の学生を除籍処分にしないかぎり、信頼の回復はないね。つまり、いつまで経ってもシリツ医大卒＝裏口バカと汚名は拭えない。シリツ出身者こそ、裏口入学に厳しい処分せよを訴えるべき。
裏口入学医師の免許剥奪を！の国民運動の先頭に立てばよいぞ。
僕も裏口入学とか、言ってたら信頼の回復はない。

まあ、こっちでもいいような。
ド底辺シリツ医は
裏口シリツ医をCOVID-19の診療にあたらせればいい。
但し、頭脳労働に従事させてはならない。喀痰吸引とか体位変換要員でいい。
算数の掛け算すら怪しいのがシリツ医がこういう事故を起こす。
://i.imgur.com/ArPaux9.png
馬鹿を医師から淘汰するのに丁度いいと思う。

特に私立大学医卒の人が順調に医師になっていくのが悔しい　私立医卒の人に負けてしまった　悔しくてたまりません

by 浩之

うりゅうひろゆき先輩

野菜は？

はい　お答えします

野菜は

じゃがいも　キュウリ　にんじん　キャベツ　みかん

です

僕　変ですか？　

by 浩之

金がないこと　女性にモテないこと　精神病　友達がいないこと　仕事が見つからないこと　で苦しんでいます　誰か助けてください

by 浩之

恥ずかしながら　十万円貰えたら　性風俗で遊ぼうと思ってます　楽しみです

by浩之

お金がない！

by 浩之

タクシーチケット支給の職場だが却ってハイリスクかもしれん。

乗客から感染か　熊本のタクシー運転手死亡 [どどん★]
https://asahi.5ch.net/test/read.cgi/newsplus/1588492505/

そんなの、

うりゅうひろゆきさんには関係無いよ

お医者さんのお話でしょ

それよりも

画用紙は？

お絵かきできたの？

はい　すみません

お金がなくて　画用紙が買えませんでした

早く　仕事を見つけて

お金が欲しいです

なにか仕事をください

by 浩之

私立大学医卒で順調に医師になってゆく人が　羨ましくて　悔しくて　たまりません

僕は国試に合格できなくて　今は無職でニートの　社会的底辺です　どうか優しくし

てください　お願いします

by 浩之

SEIRモデルに再感染(免疫を失ってR からSへの変遷）があるとしてグラフ化してみた。

1万人に1人の感染者と9人の「保菌者」がいるとして、

再感染率０％のとき
https://i.imgur.com/XWzQN7I.png

再感染率０．１％（１０００人に１人の割合で免疫を失う）場合
https://i.imgur.com/EdxIhMT.png

再感染率０．１％で外出自粛等で接触が８割減になった場合
https://i.imgur.com/sDppUfY.png

第二波が来るのもみてとれる。

最近の報道

インドは、中国企業から購入した中共ウイルス（新型コロナウイルス）の迅速スクリーニング検査キットの精度がわずか5％だとして、約50万個の注文をキャンセルした。

....
神戸市中央区にある市立医療センター中央市民病院の医師などのグループは、ことし３月末から先月７日にかけて、新型コロナウイルス以外の理由で外来を受診した患者から無作為に1000人を選び、血液中に新型コロナウイルスに感染したあとにできる「抗体」があるか調べました。
グループによりますとその結果、3.3％にあたる33人から抗体が検出されたということです。

以上を知ったあるド底辺シリツベンチャー企業が
必ず陽性がでる試薬33個と必ず陰性がでる試薬967個を混ぜた1000試薬をセットにしてインドに売り込んだ。

問題
(1)このセットの感度　特異度　陽性的中率　陰性的中率を求めよ
(2)各試薬は1回しか検査できないとして
これがイカサマキットであることを証明する手段はあるか？

英語を書いたり　統計を書いたりするのは　少しでも医師らしく　また　頭が良い雰囲気を出すためなんです　実際の僕は　国試に合格できなくて医師になれなかった馬鹿ものです　ほんとうにごめんなさい　本当の僕は無職でニートの社会的底辺です

by 浩之

rm(list=ls())
library(rjags)
modelString='
model{
marginal_likelyhood = N*p^r*(1-p)^(n-r)
p ~ dbeta(a,b)
}
'
n=10 ; r=9
N=choose(n,r)
a=1 ; b=1
data=list(n=n,r=r,N=N,a=a,b=a)
writeLines(modelString,'tmp.txt')
jagsModel=jags.model('tmp.txt', data=data,n.adapt=1e5)
update(jagsModel)
codaSamples=coda.samples(jagsModel, n.iter=1e7, thin=1,
variable=c('marginal_likelyhood'))
js=as.data.frame(as.matrix(codaSamples))
mean(js$marginal_likelyhood)
BEST::plotPost(js$marginal_likelyhood)
(PDM1=choose(10,9)*(0.5^10)) # P(D|M1) M1:p=0.5
mean(js$marginal_likelyhood)/(choose(10,9)*(0.5^10)) # P(D|M2)

(PDM2=integrate(function(x,n=10,r=9,a=1,b=1) # P(D|M2) M2:p ~ beta
choose(n,r)*x^r*(1-x)^(n-r)*dbeta(x,a,b),0,1)$value)
PDM1/PDM2
PDM2/PDM1

二項分布の共役分布がβ分布であることを体感

likelihood <- function(x,n=10,r=9) choose(n,r)*x^r*(1-x)^(n-r)
prior <- function(x,a=1,b=1) dbeta(x,a,b)
posterior <- function(x,n,r,a,b) likelihood(x,n,r)*prior(x,a,b)/
integrate(function(x,n=10,r=9,a=1,b=1)
choose(n,r)*x^r*(1-x)^(n-r)*dbeta(x,a,b),0,1)$value
curve(posterior(x,10,9,1,1))
curve(dbeta(x,1+9,1+10-9),type='h',add=T)

>>247
この（2）は確実な方法は思いつかんなぁ。数学板でのレスを待とう。

民間病院での心掛け　：　よそいけや　ふいにとびこむ　しんころな
リハビリテーション病院で集団感染しているから他人事じゃないね。

今時、集まってカンファレンスをやっているって感染対策を考えたらありえんな。
昼休みの薬屋の説明会も全部中止になったし。

ド底辺シリツ医を裏口シリツ医をCOVID-19の診療にあたらせればいい。
但し、頭脳労働に従事させてはならない。喀痰吸引とか体位変換要員でいい。
算数の掛け算すら怪しいのがシリツ医がこういう事故を起こす。
://i.imgur.com/ArPaux9.png
馬鹿を医師から淘汰するのに丁度いいと思う。

このスレのド底辺シリツ医はすでに首になったのかな。
この時期にカンファレンスとか書いていたなぁ。

ネトウヨは安倍は改憲派の右派思想の政治家だからマンセーしていて
菅直人は左翼思想の政治家だからバッシングしているんだろうけどさ

もしもそういう右左の思想をまったく考慮しないで
福一のときの安倍のやったことを考えると
安倍って最低な人間性なんだよ

一生懸命に国民のために何かできないか、東電から原発の報告がまったく途絶えたので
これは自分で情報を見に行くしかないとヘリで視察を決意
当時、枝野にそんなことすれば野党の自民党に攻撃の口実を与えるだけになるから
やめてくださいと止められたのに
菅直人は俺や東電幹部は国民のために犠牲になってもいいんだと情報収集のために福一視察を強行する
それをもって安倍がメルマガで、菅直人の視察が水素爆発の原因になったとデマを発信する

情報が途絶えた原発事故の状況を知るために、つまり国民の命を守るために
自分が被ばくしても、犠牲になってもいいから情報収集しに視察した菅直人を
なんの根拠もないデマで貶めて、いまだに叩かれる原因をつくったのが安倍晋三


思想に対する偏見を捨てて、この二人を見比べると
人間性ではどうみても

菅直人＞＞＞＞＞＞＞＞＞＞安倍晋三（ただのクズ）

となるんだよな

こういう投稿があったので

675 名前：不要不急の名無しさん[sage] 投稿日：2020/05/04(月) 08:28:07.66 ID:969To04l0
>>24
自粛で干上がった歯科医師がやるように仕向けてる
ホテルも干上がるまではほとんど非協力的だったから同じ作戦だろ

こういう建設的提言をレスしておいた。


678 名前：不要不急の名無しさん[sage] 投稿日：2020/05/04(月) 08:33:13.75 ID:H6kiOvdG0
>>675
ド底辺シリツ医を裏口シリツ医をCOVID-19の診療にあたらせればいい。
但し、頭脳労働に従事させてはならない。喀痰吸引とか体位変換要員でいい。
算数の掛け算すら怪しいのがシリツ医がこういう事故を起こす。
://i.imgur.com/ArPaux9.png
馬鹿を医師から淘汰するのに丁度いいと思う。

裏口シリツ医をCOVID-19の診療にあたらせればいい。
但し、頭脳労働に従事させてはならない。
喀痰吸引とか体位変換要員でいい。
算数の掛け算すら怪しいのがシリツ医がこういう事故を起こす。
://i.imgur.com/ArPaux9.png
馬鹿を医師から淘汰するのに丁度いいと思う。

シリツ医が正規入学か調査したとする。2人しらべたら0であった。
シリツ医の正規入学確率の事前分布として確率は0.5に設定(M1)するのと、一様分布にする(M2)のではどちらが優れたモデルか？

周辺尤度(エビデンス)を比較する。
Model1の周辺尤度はn=2,r=0,p=0.5として

PDM1=choose(n,r)*p^r*(1-p)^(n-r)

Model2を一般化すると

model{
marginal_likelyhood = N*p^r*(1-p)^(n-r)
p ~ dbeta(a,b)
}
PDM2=integrate(function(x) # P(D|M2) M2:p ~ beta
choose(n,r)*x^r*(1-x)^(n-r)*dbeta(x,a,b),0,1)$value
つまり、
∫[0,1] nCr*p^r*(1-p)^(n-r)*dp
dbeta(1,1)は一様分布でpdf=1

PDM2/PDM1=1.33 Bayes Factor

やはり、裏口入学確率を５分５分と考える前提が誤りと判明ｗ

英語を書いたり　統計を書いたりするのは　少しでも医師らしく　また　頭が良い雰囲気を出すためなんです　実際の僕は　国試に合格できなくて医師になれなかった馬鹿ものです　ほんとうにごめんなさい　本当の僕は無職でニートの社会的底辺です

by 浩之

数式は斜め読みにしてプログラムコードを読んでいたけど、
コロナのおかげでまとまった時間がとれるので数式を理解しつつ自分でプログラムを書いていくと勉強になるなぁ。

コロナ診療をやっている大学病院などは職員の外勤を禁止してほしいね。
うちは回復期リハビリテーション病棟があるからここで感染すると大変。
DNARじゃなくて自宅へ退院が前提になっている患者群だからここを死守しなくちゃならん。

コロナ診療の人手不足は裏口シリツ医を強制配置すればいいと思う。
但し、頭脳労働に従事させてはならない。
喀痰吸引とか体位変換要員でいい。
算数の掛け算すら怪しいのがシリツ医がこういう事故を起こす。
://i.imgur.com/ArPaux9.png
馬鹿を医師から淘汰するのに丁度いいと思う。

みんなが順調に医師になっていくのがたまらなく妬ましいし悔しくてたまりません
わかってください

by 浩之

このスレの趣旨に合わせた問題。

ド底辺シリツ医の裏口入学確率を推定するのに
モデルM1は頭が悪いよねスレへの投稿回数から、
モデルM2はそれに加えてセンズリ回数から、算出するとする。
経験上M1のM2に対するベイズファクター(エビデンス強度の比率の指標)は1.5である。
裏口確率をモデル1は90%モデル2は80%と推定した。
このド底辺シリツ医の裏口入学確率の期待値を求めよ。

建設的提言を投稿してきました。

https://asahi.5ch.net/test/read.cgi/newsplus/1588544942/305

建設的提言を追加投稿してきました。

エレベータがひとつしかない個人病院にコロナ疑いの患者を診て院内感染を防げと言われても無理。
下手したら倒産。集中と集約が必要。

コロナ診療の人手不足は裏口シリツ医を強制配置すればいいと思う。
但し、頭脳労働に従事させてはならない。
喀痰吸引とか体位変換要員でいい。
算数の掛け算すら怪しいのがシリツ医がこういう事故を起こす。
://i.imgur.com/ArPaux9.png
馬鹿を医師から淘汰するのに丁度いいと思う。

馬鹿を医師から淘汰するのに丁度いいと思う。

のエビデンスとしてこれも追加した方がいいな。


シリツ医のこれが現実

１次方程式もできないド底辺特殊シリツ医大卒の記録
http://imagizer.imageshack.com/img923/2715/RosCsf.jpg

コロナ診療の人手不足は裏口シリツ医を強制配置すればいいと思う。
但し、頭脳労働に従事させてはならない。
喀痰吸引とか体位変換要員でいい。
算数の掛け算すら怪しいのがシリツ医がこういう事故を起こす。
://i.imgur.com/ArPaux9.png

シリツ医のこれが現実
１次方程式もできないド底辺シリツ医大卒の記録
://imagizer.imageshack.com/img923/2715/RosCsf.jpg

馬鹿を医師から淘汰するのに丁度いいと思う。

医師会から野外PCR検査検体採取する医師を募集していたなぁ。
唾液検体でもKappa値は十分なのに。
こういう仕事にはド底辺シリツ医でも徴用してほしいね。

頭脳労働に従事させてはならない。

喀痰吸引とか体位変換要員でいい。

算数の掛け算すら怪しいのがシリツ医がこういう事故を起こす。

://i.imgur.com/ArPaux9.png

シリツ医のこれが現実
１次方程式もできないド底辺シリツ医大卒の記録

://imagizer.imageshack.com/img923/2715/RosCsf.jpg

馬鹿を医師から淘汰するのに丁度いいと思う。

やっぱり、リハビリテーションって密着するから危ないな

>>

その病院は、もとは消化器・循環器・糖尿病等の専門医のいる病院だったが、
去年の暮れに山口県のリハビリ専業医療法人グループに譲渡され、年末で外来終了。
それまでいた医者もナースも用済みサヨナラ、４月１日付で理事長・院長の首もスゲ変えて
文字通り「生まれ変わった」病院風リハビリ施設。

リストラしまくりホヤホヤのリハ専病院がコロナで人手不足になったって、医者もナースも助けちゃくれんわな

発覚時のニュースで「医療の現場はそれほどに切迫して…」的なアホコメントがそこかしこで語られてて
あまりの取材力のなさに愕然としたわｗ

American astronomer
Carl Sagan (1934?1996) coined the exact phrase “extraordinary
claims require extraordinary evidence.” The maxim is incorporated in the
Bayesian computation where prior odds are combined with the Bayes factor
to yield posterior odds.

意訳：ド底辺シリツ医が裏口入学でないことを示すのは困難。

In contrast to p-value hypothesis testing, Bayesian statistics assigns
no special status to the null hypothesis and this means that Bayes factors can
be used to quantify evidence for the null hypothesis just as for any other hypothesis.
A second advantage of Bayes factors is that they allow one to monitor the evidence
as the data come in.

p hackingへの批判だなぁ。ド底辺シリツ医の知能指数も有意差なしという統計処理が頻度主義統計では可能なんだよな。
boot strap使ってp値の分布とか出せば簡単に反駁できるけど、ドヤ顔でｐ値を出されたら黙るアホが大杉。

It is entirely appropriate to collect data until a point has been proven or disproven,
or until the data collector runs out of time, money, or patience” (Edwards et al.,
1963, p. 193). This means that researchers are free to continue data collection in
case the evidence (i.e., the Bayes factor) is inconclusive. Likewise, they are free to
terminate data collection as soon as the interim evidence is sufficiently compelling.


This freedom of “optional stopping” is denied to the researchers who use p-values
to test their hypotheses (Wagenmakers, 2007).

p(D|θ,M1) p(θ|M1)
p(D|M1) = ------------------
p(θ|D,M1)

裏口私立医は頭の他に手先も不器用だと誤解している人がいるようだ。

【経済】一般病院、患者激減　新型コロナで減収も補償対象外　地域医療の崩壊懸念　｢今後1〜2カ月、この状況が続けば経営が厳しい｣　★2 [特選八丁味噌石狩鍋★]
https://asahi.5ch.net/test/read.cgi/newsplus/1588552251/581

581 名前：不要不急の名無しさん[] 投稿日：2020/05/04(月) 13:39:36.57 ID:iRJqsoMZ0
>>576
裏口私立じゃぁ、感染対策の手技、心配(ごめん
さすがにジョーク)

擁護する投稿をしてきた。

https://asahi.5ch.net/test/read.cgi/newsplus/1588552251/588

大学入試の偏差値なんか知性と関係ないよ

安部首相の取り巻きの偏差値秀才を見れば分かるだろうにｗ

いつまで受験偏差値に拘ってんだよｗ

だから浪人２５年なんだよｗ

おい、うりゅう爺さん聞いてるか？

いつまでも医科歯科卒って詐欺ってたらダメだよ

現実と戦わなきゃｗ

>>268
自分の頭で考えない調べないバカは死ねばいいんだよ

どうでもいいが藤田と川崎は何であんなにバカなんだ！
勘弁してくれよ！！！

建設的提案をしてきました。
https://asahi.5ch.net/test/read.cgi/newsplus/1588552251/756

うりゅうひろゆき先輩

画用紙用意して、

時計のお絵かき

できないのなら、

ちがう事して

リハビリしますね



あんたが１３年間も

学生として過ごした町で

ダイヤルQ2で

女の子達からの呼び出しを

真に受けて、ヤンキー中学生からカツアゲされた

思い出の港で

今まさにクラスターだ

一応医学部卒業してんだから

行って来い

Winter <- c(-0.05,0.41,0.17,-0.13,0.00,-0.05,0.00,0.17,0.29,0.04,0.21,0.08,0.37,0.17,0.08,-0.04,-0.04,0.04,-0.13,-0.12,0.04,0.21,0.17,
0.17,0.17,0.33,0.04,0.04,0.04,0.00,0.21,0.13,0.25,-0.05,0.29,0.42,-0.05,0.12,0.04,0.25,0.12)

Summer <- c(0.00,0.38,-0.12,0.12,0.25,0.12,0.13,0.37,0.00,0.50,0.00,0.00,-0.13,-0.37,-0.25,-0.12,0.50,0.25,0.13,0.25,0.25,0.38,0.25,0.12,
0.00,0.00,0.00,0.00,0.25,0.13,-0.25,-0.38,-0.13,-0.25,0.00,0.00,-0.12,0.25,0.00,0.50,0.00)
t.test(Winter,Summer)

> t.test(Winter,Summer)

Welch Two Sample t-test

data: Winter and Summer
t = 0.81514, df = 68.386, p-value = 0.4178
alternative hypothesis: true difference in means is not equal to 0
95 percent confidence interval:
-0.04943605 0.11772874
sample estimates:
mean of x mean of y
0.10756098 0.07341463

nx=length(x)
modelString='
model{
for(i in 1:nx){
x[i] ~ dnorm(mu,pow(sigma,-2))
}
mu = delta*sigma
sigma ~ dt(0,1,1)T(0,)
delta ~ dt(0,1,1)
# delta ~ dt(0,1,1)T(,0)
}
'
writeLines(modelString,'tmp.txt')
dataList=list(x=x,nx=nx)
jagsModel=jags.model('tmp.txt', data=dataList, n.chains=4)
update(jagsModel)
codaSamples=coda.samples(jagsModel,n.iter=1e4,thin=1, var=c('delta'))
coda2summary(codaSamples)
plot(codaSamples)
js=as.data.frame(as.matrix(codaSamples))
plotPost2(js$delta,xlab=bquote(delta),compVal = 0, ROPE=c(-0.1,0.1))
delta=density(js$delta,n=1e5)
plot(delta,xlab=bquote(delta),main='',xlim=c(-3,3),lty=2)
curve(dt(x,1),add=T)
abline(v=0,lty=3,col='gray')
prior=dt(0,1)
x=delta$x[which(round(delta$x,5)==0)]
y=delta$y[which(round(delta$x,5)==0)]
posterior=y
posterior/prior # BayesFactor by Savage-Dickey density ratio
points(x,y,pch=19)　; points(0,dt(0,1),pch=19)

Bayes Factor （計算はMCMC→ Savage-Dickey Density Ratio）で帰無仮説H0 mu=0　採択の根拠を定量化できる。

老人でも血中濃度は同等だから安全とか、ド底辺シリツ医でも知能指数は劣らない　という統計捏造に　
サンプルサイズを小さくするのは見え透いているが、

Bayes Factorで帰無仮説を定量的に採択という捏造は巧妙にできるなぁ。
複雑な階層モデルにするほど反証が難しくなる。

自分で統計処理できないド底辺シリツ医なら簡単に騙せる。

棄却前提の帰無仮説に採択の根拠を与えるBayesFactorにKruscheが反対する理由がよくわかる。

こういうツールがあるといいなと思うと、大抵、Rにパッケージがあるな。
MCMCサンプルから密度を計算

library(polspline)
fit=logspline(js$delta)
plot(fit,xlim=c(-3,3),lty=2)
curve(dt(x,1),add=T)
(d0=dlogspline(0,fit))
points(0,d0,pch=19)
points(0,dt(0,1),pch=19)
(BF01=d0/dt(0,1))

ベイズ係数によるPCR検査感度モデル比較

rm(list=ls())
par(bty='l')
library(rjags)
source('tools.R')
library(polspline)
# https://toyokeizai.net/sp/visual/tko/covid19/
P=14895/153581 # 2020/05/04
n=c(153581)
r=c(14895)
N=length(n)

sn1=Mv2ab(0.5,0.2^2) # 30-70%
sn2=Mv2ab(0.8,0.1^2) # 70-90%
sp=Mv2ab(0.95,0.25^2) # 92.5-97.5%
data1=list(n=n,N=N,r=r,sn=sn1,sp=sp)
data2=list(n=n,N=N,r=r,sn=sn2,sp=sp)

modelString='
model{
for(i in 1:N){
r[i] ~ dbin(p,n[i])
}
p = prev*sen + (1-prev)*(1-spc)
sen ~ dbeta(sn[1],sn[2])
spc ~ dbeta(sp[1],sp[2])
prev ~ dbeta(1,1)
}
'

writeLines(modelString,'tmp.txt')
n.iter=1e4
thin=1
jagsModel1=jags.model('tmp.txt',data1,n.chains=4,n.adapt = n.iter/5)
jagsModel2=jags.model('tmp.txt',data2,n.chains=4,n.adapt = n.iter/5)
update(jagsModel1)
update(jagsModel2)
var=c('p','sen','spc','prev')
codaSamples1=coda.samples(jagsModel1,n.iter=n.iter,thin=thin,var=var)
codaSamples2=coda.samples(jagsModel2,n.iter=n.iter,thin=thin,var=var)
coda2summary(codaSamples1)
coda2summary(codaSamples2)

js1=as.data.frame(as.matrix(codaSamples1))
js2=as.data.frame(as.matrix(codaSamples2))

fit1=logspline(js1$p)
fit2=logspline(js2$p)
(d1=dlogspline(0,fit1)) # 30-70%
(d2=dlogspline(0,fit2)) # 70-90%
plot(density(js1$p),xlim=c(0.096,0.098),xlab=bquote(p),main='')
lines(density(js2$p),col=2)
abline(v=P,lty=3)
d1/d2

バグ修正

fit1=logspline(js1$p)
fit2=logspline(js2$p)
plot(fit1,xlim=c(0.0965,0.0975),ylim=c(400,550))
par(new=T)
plot(fit2,col=2,xlim=c(0.0965,0.0975),ylim=c(400,550),axes=F)
abline(v=P,lty=3)
(d1=dlogspline(P,fit1)) # 30-70%
(d2=dlogspline(P,fit2)) # 70-90%
d1/d2 # Savage-Dickey densiti ratio = BF12

https://toyokeizai.net/sp/visual/tko/covid19/
のデータP=14895/153581 　（#　2020/05/04の時点)
を使って
PCR検査の感度30-70％のモデルM1と感度70-90％のモデルM2のどちらが信憑性があるか、ベイズファクターで計算してみる。
M1は感度が最頻値60%標準偏差10%、M2は最頻値80%標準偏差10%のベータ分布に設定
特異度はいずれも最頻値95%標準偏差2.5%に設定し、有病率は一様分布を仮定
陽性数は、陽性率(P)＝真陽性率＋偽陽性率=有病率=有病率*感度 + (1-有病率)*(1-特異度)の二項分布に従うとする。
事後確率分布は
https://i.imgur.com/81bK4KE.png
陽性率P=14895/153581=0.09698465での事後確率分布の密度比(Savage-Dickey density ratio)でベイズファクターを出すと
> d1/d2 # Savage-Dickey densiti ratio = BF12
[1] 1.007722

まあ、ちょっぴり、感度30-70％のモデルの方がいいかも、という結果。

離散量データから確率密度計算をやってくれるlogsplineは便利で(　・∀・)ｲｲ!!

https://toyokeizai.net/sp/visual/tko/covid19/
のデータP=14895/153581 　（#　2020/05/04の時点)
を使って
PCR検査の感度30-70％のモデルM1と感度70-90％のモデルM2のどちらが信憑性があるか、ベイズファクターで計算してみる。
M1は感度が最頻値60%標準偏差10%、M2は最頻値80%標準偏差10%のベータ分布に設定
特異度はいずれも最頻値95%標準偏差2.5%に設定し、有病率は一様分布を仮定
陽性数は、陽性率(P)＝真陽性率＋偽陽性率=有病率*感度 + (1-有病率)*(1-特異度)の二項分布に従うとする。
事後確率分布は
://i.imgur.com/81bK4KE.png
陽性率P=14895/153581=0.09698465での事後確率分布の密度比(Savage-Dickey density ratio)でベイズファクターを出すと
> d1/d2 # Savage-Dickey densiti ratio = BF12
[1] 1.007722

まあ、ちょっぴり、感度30-70％のモデルの方がいいかも、という結果。

離散量データから確率密度計算をやってくれるlogsplineは便利で(　・∀・)ｲｲ!!

こういう警告とともに旧版での計算に移行しているみたいなので、最初からoldlogsplineを使った方が時間の短縮になるな。


> fit1=logspline(js1$p)
Warning messages:
1: In logspline(js1$p) : too much data close together
2: In logspline(js1$p) : re-ran with oldlogspline
> fit2=logspline(js2$p)
Warning messages:
1: In logspline(js2$p) : too much data close together
2: In logspline(js2$p) : re-ran with oldlogspline

>>290
旧版での計算に修正

fit1=oldlogspline(js1$p)
fit2=oldlogspline(js2$p)
plot(fit1,xlim=c(0.0965,0.0975),ylim=c(400,550),
xlab='陽性率')
legend('topright',bty='n',legend=c('感度:30-70%','感度:70-90%'),lty=1:2,col=1:2,)
par(new=T)
plot(fit2,col=2,xlim=c(0.0965,0.0975),ylim=c(400,550),axes=F,lty=2)
abline(v=P,lty=3)
(d1=doldlogspline(P,fit1)) # 30-70%
(d2=doldlogspline(P,fit2)) # 70-90%
d1/d2 # Savage-Dickey densiti ratio = BF12

うりゅうひろゆき先輩

ますます重症みたいだが？

先輩が、エア医師遊びしても

誰も困らんから

別にええけど

画用紙買ってきましたか？

時計のお絵かきしよう

エア医師遊びは

時計のお絵かき出来てからしようか？

10時10分の…

いや、まず、画用紙買っきてもらいなさい

担当の

介護のおばさん、

呼んでごらん

k=1e5
x=runif(k)-runif(k)
var(x)

頻度差の分布を出すのに事前分布を一様分布
θ1 ~ dunif(0,1)
θ2 ~ dunif(0,1)
にしたときの
θ1 - θ2の分布って三角分布になるんだな。

体感するには、

library(triangle)
k=1e7
x=runif(k)-runif(k)
hist(x,freq=F,breaks=100,ann=F)
var(x)
mean(x)
curve(dtriangle(x,-1,1),add=T)
curve(dnorm(x,0,sqrt(1/6)),add=T)

ニュー速＋で医師叩きの投稿が続くので建設的提案をレスしてきました。

うりゅうひろゆき先輩

いつものリハビリですか？

エア医師遊びできてよかったですね

じゃこんどこそ

画用紙買ってきてもらいなさい

時計のお絵かきしようかな

出来ますか？

Bayes　Factorによる二項比率の比較を関数化

# s1,s2,n1,n2 -> BF10
BF01 <- function(s1,s2,n1,n2){
library(gmp)
bf=chooseZ(n1,s1)*chooseZ(n2,s2)/chooseZ(n1+n2,s1+s2)*(n1+1)*(n2+1)/(n1+n2+1)
BF01=as.numeric(bf)
return(c('in favor of null-hyposisis : H0'=BF01))
}
(bf01=BF01(s1=424,s2=5416,n1=777,n2=9072))
c('in favor of alternative hyposis : H1' = 1/as.numeric(bf01))

コンドーム装着率のBayes Factorによる検定

(bf01=BF01(s1=424,s2=5416,n1=777,n2=9072))

私は九大出身の臨床医だが
国立大学卒の医者でもできる奴もいれば
できない奴もいる
私立大学卒でもできる奴もいれば
できない奴もいるという
印象だ
要は医師になってからが
勝負だろう
強く言っておく

それから
国家試験浪人は
1年までだ
それ以上は
軽蔑される

俺の頃は秋にも国試があったから半年で通っていたな。

>>303
じゃあ、ド底辺シリツ医にかわって手コキ問題に答えてみ。
数学板に類題の答を書いたから検索すればみつかるかも。

うりゅうひろゆき先輩

画用紙買ってきましたか？

時計のお絵かき出来てからだよ

そんなの

日本には感染症対策の指揮官が居ないのが残念。
欧米もそんな感じ。
でもサーズとマーズを経験した香港・台湾・韓国には居た。

国民を分断して収奪する安倍らの集団の罠に嵌る椰子大杉。

熱がある人は電話してとか書いてあるのに直接診察室にはいってきたり、別の病気で来た人がコロナだったり、濃厚接触者であることを隠して病院来て診察受けたりして、院内感染させて外来休診させている。
診なくて外来受けてないんじゃなく、コロナが無防備のところに入ってきて、休診させられる。

うりゅうひろゆき先輩には

なんの関係もない話だけど

いつものエア医師遊びを

したいのなら

参加してもいいけど

別に

誰も困らんから

自由

でも、ね、

今日は

何月何日ですか？

答えられてからだよ

まずはね

>>303
でも私立卒のできないヤツは本当に出来ないよｗ
特に60歳以上の医者。

>>312
若いのも同じ。

算数の掛け算すら怪しいのがシリツ医がこういう事故を起こす。
http://i.imgur.com/ArPaux9.png

シリツ医のこれが現実
１次方程式もできないド底辺シリツ医大卒の記録
http://imagizer.imageshack.com/img923/2715/RosCsf.jpg

自称九州大学卒というのに
手コキ問題に答えられないアホの存在も確認されたんだが、
九州大学ってその程度の頭でも入れるのか？
それとも自称なのか？、

安倍に関わった奴は東大卒でも馬鹿になっていくな。

お仕事はなにをされてるんですか？

うりゅうひろゆき先輩

ぜんぜん治療が

できてないぞ

エア医師遊びしている場合か？

ヘルパーさんに叱られるよ

画用紙用意は？

でも良いよな〜、裏口底辺私立卒の医者は。
だって親がお金持ちじゃなかったら本来は医師じゃ無かった連中なんだぜ。
その上、受験科目も3科目と国立に比べて極端に少ないし。

国立大学2次試験不合格となって浪人が決まったヤツや国立医学部を諦めて
旧帝理工学部や早慶の理工系に鞍替えした学生たちが、底辺私立医大に合格して
意気揚々としてる連中を遠くから悔しそうな表情で、一言ボソッと「バカのくせに、、、」
と言ったことを昨日のことのように覚えてるわ。

本当に睨みつけるような表情してたな。

今の受験制度、やっぱり間違ってると思うわ。

お仕事はされてないんですか？