臨床統計もおもしろいですよ、その2at HOSP
臨床統計もおもしろいですよ、その2 - 暇つぶし2ch541:卵の名無しさん
18/12/16 11:42:14.17 /hMzfoFS.net
Nt <- 100
set.seed(123)
Ti <- sample(20:40, Nt, replace=TRUE)
Xt <- rlnorm(Nt, 100, 15)
Yt <- rbinom(Nt, size=Ti, prob=0.8)
20-40人からなるド底辺シリツ学生のグループを100組集め


542:る。 各グループを代表する知能指数が正規分布をしている。 裏口の確率は80%である。 その線形回帰は glm(Yt ~ Xt, family=poisson(link="log"), offset=log(Ti))



543:卵の名無しさん
18/12/16 14:16:40.36 /hMzfoFS.net
source('tmp.tools.R')
library(Mass)
set.seed(71)
size <- 10000
x <- -seq(1/size,1,by=1/size)
shape <- 10.7
mu <- exp(3.7*x+2)
y1 <- sapply(mu, function(m)rnegbin(1, mu=m, theta=shape))
pairsxyz(cbind(x,y1,mu))
model_glmnb <- glm.nb(y1 ~ x)
summary(model_glmnb)
plot(x,y1,pch=19,col=rgb(.01,.01,.01,.02))
points(x,fitted(model_glmnb,seq(0,1,len=size)),col=0)
lines(x,model_glmnb$fitted)
points(x,mu,col=rgb(0,0,1,0.01),cex=0.5)
#
lambda <- sapply(mu,function(x) rgamma(1,shape=shape,scale=x/shape))
y2 <- sapply(lambda,function(x) rpois(1,x))
hist(y1,f=F,breaks=50)
hist(y2,f=F,breaks=50, add=T,col=5)

544:卵の名無しさん
18/12/16 20:30:27.72 /hMzfoFS.net
# ネットワークノード1..j..N
# ネットワーク内の客の数1..k..K、
# ノードjでの平均客数L[k,j]、通過時間W[k,j]、到着率lambda[k,j]
# mu[j]:サービス率,r[i,j]ノードiからjへの移行確率
# theta[j] = Σtheta[i](i=1,N) * r[i,j] の解がtheta
rm(list=ls())
theta=c(0.6,0.4)
mu=c(0.6,0.4)
N=2
K=10
L=W=lambda=matrix(rep(NA,N*K),ncol=N)
k=1
for(j in 1:N){
W[k,j]=1/mu[j]*(1+0)
}
tmp=numeric(N)
for(j in 1:N){
tmp[j]=theta[j]*W[k,j]
}
den=sum(tmp)
lambda[k,1]=k/den
for(j in 2:N){
lambda[k,j]=theta[j]*lambda[k,1]
}
for(j in 1:N){
L[k,j]=lambda[k,j]*W[k,j]
}

545:卵の名無しさん
18/12/16 20:30:52.05 /hMzfoFS.net
#
for(k in 2:K){
for(j in 1:N){
W[k,j]=1/mu[j]*(1 + L[k-1,j])
}
tmp=numeric(N)
for(j in 1:N){
tmp[j]=theta[j]*W[k,j]
}
den=sum(tmp)
lambda[k,1]=k/den
for(j in 2:N){
lambda[k,j]=theta[j]*lambda[k,1]
}
for(j in 1:N){
L[k,j]=lambda[k,j]*W[k,j]
}
}
L
W
lambda

546:卵の名無しさん
18/12/17 04:38:09.36 GOYrgwsC.net
minimize 2*b+4*sin(θ)/(cos(θ)-1/2)*α where (θ-sin(θ)/2+(sin(θ)/(cos(θ)-1/2))^2(α-sin(α))/2+(1-b)*sin(θ)/2=pi/8 ,0<b<1,0<θ<pi/2

547:卵の名無しさん
18/12/17 04:45:32.12 GOYrgwsC.net
minimize (2*b+4rα) where (θ-sinθ)/2+r^2(α-sinα)/2+(1-b)sinθ/2=pi/8 ,0<b<1,0<θ<pi/2

URLリンク(www.wolframalpha.com)(2*b%2B4r%CE%B1)++where+(%CE%B8-sin%CE%B8)%2F2%2Br%5E2(%CE%B1-sin%CE%B1)%2F2%2B(1-b)sin%CE%B8%2F2%3Dpi%2F8+,0%3Cb%3C1,0%3C%CE%B8%3Cpi%2F2

Standard computation time exceeded...
Try again with Pro computation time

548:卵の名無しさん
18/12/17 05:14:49.85 GOYrgwsC.net
rangeB = map (/1000) [1..1000]
rangeTheta = map (\x -> x * pi/2/1000) [1..1000]
rangeAlpha = map (\x -> x * pi/1000) [1..1000]
rangeR = map (\x -> x * 1/1000) [1000..10000]
re = [2*b+4*r*α| b<-rangeB,θ<-rangeTheta, α<-rangeAlpha,r<-rangeR,(θ-sin(θ))/2+r^2*(α-sin(α))/2+(1-b)*sin(θ)/2==pi/8]
main = do
print $ minimum re

549:卵の名無しさん
18/12/17 05:39:04.50 GOYrgwsC.net
rangeB = map (/1000) [1..1000]
rangeTheta = map (\x -> x * pi/2/1000) [1..1000]
rangeAlpha = map (\x -> x * pi/1000) [1..1000]
re = [(2*b+4*sin(θ)/(cos(θ)-1/2)*α)| b<-rangeB,θ<-rangeTheta, α<-rangeAlpha,(θ-sin(θ))/2+(sin(θ)/(cos(θ)-1/2))^2*(α-sin(α))/2+(1-b)*sin(θ)/2==pi/8]
main = do
print $ minimum re

550:卵の名無しさん
18/12/17 05:47:22.52 GOYrgwsC.net
URLリンク(ja.wolframalpha.com)(2*b%2B4r%CE%B1)++where+(%CE%B8-sin%CE%B8)%2F2%2Br%5E2(%CE%B1-sin%CE%B1)%2F2%2B(1-b)sin%CE%B8%2F2%3Dpi%2F8+,0%3Cb%3C1,0%3C%CE%B8%3Cpi%2F2,0%3C%CE%B1%3C2*pi

551:卵の名無しさん
18/12/17 05:47:36.95 GOYrgwsC.net
minimize(2*b+4rα) where (θ-sinθ)/2+r^2(α-sinα)/2+(1-b)sinθ/2=pi/8 ,0<b<1,0<θ<pi/2,0<α<2*pi

552:卵の名無しさん
18/12/17 13:52:12.63 FuhGkKV8.net
向学心のないコミュ障事務員のスレ
事務員と見せかけてその正体は大学受験に失敗したニートだろうな
二十歳前くらいの大学受験の知識で一生わたっていけると勘違いしているらしい
いや勘違いというか自分にそう言い聞かせなければいけないほど
人生終わってんだろうな
ハゲワロス

553:卵の名無しさん
18/12/17 13:59:09.34 N+SAyf9z.net
>>518
こういう野糞投稿しかできないのはド底辺シリツかな?
ド底辺スレから部屋割計算できなくて逃げ出したアホと同一人物?

554:卵の名無しさん
18/12/17 14:01:55.57 N+SAyf9z.net
>>518
向学心を語るなら>499の計算してからな。
ド底辺頭脳クン。

555:卵の名無しさん
18/12/17 14:02:51.74 FuhGkKV8.net
人生についての向学心は大学受験なんていう狭い領域で机上の空論をこねくり回すことではない
だがいつまでも不合格の受験生は、その領域から抜け出せなくなるんだなW
向学心のない事務員にふさわしいスレを持ってきてやったぞ
【因果応報】医学部再受験【自業自得】
スレリンク(kouri板)

556:卵の名無しさん
18/12/17 14:07:07.62 N+SAyf9z.net
>>521
>499は受験に無関係な臨床問題だよ。
答えてみ!
はよ!はよ!
向学心を語るなら>499の計算してからな。
ド底辺頭脳クン。

557:卵の名無しさん
18/12/17 14:08:05.38 FuhGkKV8.net
このスレもお勧めだぞW
【悲報】ワイの医学部志望の妹、精神壊れる
URLリンク(vippers.jp)

558:卵の名無しさん
18/12/17 14:09:08.10 N+SAyf9z.net
>>523
>>518
向学心を語るなら>499の計算してからな。
ド底辺頭脳クン。

559:卵の名無しさん
18/12/17 14:14:33.64 N+SAyf9z.net
ド底辺頭脳でもこれくらいできるかな?
12^2 = 144 ........................ 21^2 = 441
13^2 = 169 ........................ 31^2 = 961
二桁の数でこうなるのは他にある?
向学心も頭脳もないから答えられないだろうな。

560:卵の名無しさん
18/12/17 14:26:10.34 FuhGkKV8.net
さあそろそろ今日の宿題の問題を出題しといてやろうか
事務員には無理だろうけどな
さて問題です
ラブアンドピースとアヘ顔ダブルピースの相似性について原稿用紙30枚以上で考察せよ

561:卵の名無しさん
18/12/17 14:40:27.65 qP6ufigZ.net
>>525
2行で答がでる。
ド底辺頭脳には無理。
[(x,y)|x<-[1..9],y<-[1..9],let z=(10*x+y)^2,x<y,let a=z `div` 100,let b=z `mod` 100 `div` 10, let c=z `mod` 10,100*c+10*b+a==(10*y+x)^2]
[(x,y)|x<-[0..9],y<-[0..9],x<y,let z=(10*x+y)^2,let a=z `div` 100, let a1=a `div`10,let a2=a `mod` 10,
let b=z `mod` 100 `div` 10, let c=z `mod` 10,1000*c+100*b+a2*10+a1==(10*y+x)^2]

562:卵の名無しさん
18/12/17 14:43:19.14 qP6ufigZ.net
>>526
数値計算とか統計処理の必要な問題すらも作れないのは統計数理の素養がない馬鹿だからだろ。
臨床医ならこれくらいの計算ができる教養がほしいね。できなきゃド底辺シリツ医並の頭脳。

インフルエンザの迅速キットは特異度は高いが感度は検査時期によって左右される。
ある診断キットが開発されたとする。
このキットは特異度は99%と良好であったが、
感度については確かな情報がない。
事前確率分布として一様分布を仮定する。
50人を無作為抽出してこの診断キットで診断したところ40人が陽性であった。
この母集団の有病率の期待値と95%CIはいくらか?
またこの診断キットの感度の期待値と95%CIはいくらか

563:卵の名無しさん
18/12/17 19:41:34.06 FuhGkKV8.net
やはり低学歴ニートには答えられないんだな
参加賞で1曲あげといてやろう
URLリンク(www.youtube.com)

564:卵の名無しさん
18/12/17 19:55:33.62 N+SAyf9z.net
>>529
>>518
向学心を語るなら>499の計算してからな。
ド底辺頭脳クン。

565:卵の名無しさん
18/12/17 20:34:35.71 qP6ufigZ.net
source('tools.R')
ab=Mv2ab(0.70,0.01)
pbeta(0.80,ab[1],ab[2])-pbeta(0.60,ab[1],ab[2])
pbeta(0.90,ab[1],ab[2])-pbeta(0.50,ab[1],ab[2])
x=rbeta(1e5,ab[1],ab[2])
o=x/(1-x)
TP=0.80
FN=1-TP
TN=0.99
(nLR=FN/TN)
po=o*nLR^2
pp=po/(1+po)
BEST::plotPost(pp,compVal=0.10)

566:卵の名無しさん
18/12/17 20:34:41.41 qP6ufigZ.net
source('tools.R')
ab=Mv2ab(0.70,0.01)
pbeta(0.80,ab[1],ab[2])-pbeta(0.60,ab[1],ab[2])
pbeta(0.90,ab[1],ab[2])-pbeta(0.50,ab[1],ab[2])
x=rbeta(1e5,ab[1],ab[2])
o=x/(1-x)
TP=0.80
FN=1-TP
TN=0.99
(nLR=FN/TN)
po=o*nLR^2
pp=po/(1+po)
BEST::plotPost(pp,compVal=0.10)

567:卵の名無しさん
18/12/17 20:36:09.78 qP6ufigZ.net
インフルエンザ迅速検査キットの感度が80%、偽陽性率1%とする。
臨床的にはインフルエンザを疑うが迅速検査陰性、
不要な薬は使いたくない、と患者の要望。
検査が早期すぎるための偽陰性かと考えて翌日再検したが、こちらも陰性。
「私がインフルエンザである可能性はどれくらいありますか?
可能性が10%以上なら仕事を休んで抗ウイルス薬希望します。」
という。
臨床所見から疑われる有病率が70%程度(最頻値70%標準偏差10%)としたときにこの患者が仕事を休むことになる確率はいくらか?

568:卵の名無しさん
18/12/18 08:26:52.13 3rdYoIBg.net
バグに気づくのに時間がかかった、数学板に上げる前に気づいてよかった。
PQR = function(b,s){ # OQ=b ∠POQ=s
r=sin(s)/(cos(s)-1/2)
a=atan(sqrt(1^2+b^2-2*1*b*cos(s))/r) # <- WRONG!!
(s-sin(s))/2+r^2*(a-sin(a))/2+(1-b)*sin(s)/2
# (pi/2-s)/2 + b*sin(s)/2 - r^2*(a/2-sin(a)/2)
}
bb=seq(0,1,len=101)
ss=seq(0,pi/2,len=101)
z=outer(bb,ss,PQR)
image(bb,ss,z,xlab='b(OQ)',ylab='s(rad)')
contour(bb,ss,z,bty='l',nlevels=20,add=T, lty=3)
contour(bb,ss,z,levels=pi/8,lwd=2, add=T)
Persp3d(bb,ss,z)
b2s <- function(b){
uniroot(function(x) PQR(b,x)-pi/8, c(0,pi/2))$root
}
b2s=Vectorize(b2s)
L = function(b){ # 2b + 4ra
s=b2s(b)
r=sin(s)/(cos(s)-1/2)
a=atan(sqrt(1^2+b^2-2*1*b*cos(s))/r)
2*b + 4*r*a
}
L=Vectorize(L)
curve(L(x),0,pi/6,bty='l')
L(0)

569:卵の名無しさん
18/12/18 09:56:44.79 iW+brusm.net
こんな素朴な質問したら理解できない長文レスが数学化卒からきた。
4円弧モデルが最小ってどうすれば言えるの?
URLリンク(i.imgur.com)
URLリンク(i.imgur.com)
URLリンク(i.imgur.com)
の順に小さくなるのは計算できたけど。

570:卵の名無しさん
18/12/18 15:48:23.96 3rdYoIBg.net
暇つぶしにサムスカの太鼓持ち講演を聞いた。
相関係数を算出に全例調査しているのにp値を出していて
統計処理を理解していないアホだとわかったのが面白かった。

571:卵の名無しさん
18/12/18 15:50:29.76 3rdYoIBg.net
>>536
まあ、聞いているやつが(ド底辺シリツ医大卒のようなアホが多数で)p<0.001とか出せばエビデンスレベルが高いと誤解するだろうとわかっての所作なら策士だとは思うが。

572:卵の名無しさん
18/12/18 15:55:06.24 3rdYoIBg.net
こんなのがあんのね
電解質Na,K専用測定器 Fingraph
俺は夜間休日は血ガス測定で代用してたけど。

573:卵の名無しさん
18/12/18 21:56:13.92 RRjIHldX.net
3S Policy Causes the ResultTruman of Panama document and 3S policy
We are keeping the citizens in a cage named "Freedom of Falsehood".
The way you just give them some luxury and convenience.Then release sports, screen, sex (3S).
Because the citizens are our livestock. For that reason, we must make it longevity. It makes it sick
(with chemicals etc.) and keeps it alive.This will keep us harvesting. This is also the authority
of the victorious country.The medal and the Nobel prize is not all Raoul declarationExamples of 3S policy
Hypocrites and disguised society and Panama document "He moved here here two years ago, I was afraid
of looking, but I thought it was a polite person to greet" Hello "when I saw it ... It was not like I was waking
this kind of incident ... ... "(Neighboring residents)On November 7, Shosuke Hotta, 42, a self-proclaimed office
worker in Sagamihara City, Kanagawa Prefecture, was arrested on suspicion of abusing her daughter Sakurai Ai.
This is the nature of the Earthlings

574:卵の名無しさん
18/12/18 23:30:22.38 rfcy7s34.net
本日も事務員放送の終了時間となりました
URLリンク(m.youtube.com)

575:卵の名無しさん



576:sage
>>540 相関係数を算出に全例調査しているのにp値を出だすのが 何でアホと言われるかわかる?



577:卵の名無しさん
18/12/19 13:50:09.99 CUSSzVDj.net
>>541
わからないからp値だしているのか
わからない>540のような馬鹿につけこんでいるのか、
どっちだろうな。
サイコロが2回続けて1の目出ました。
その確率は1/6*1/6=0.02777で0.05未満だから
サイコロは有意差をもって歪(いびつ)である、
即ち、1の目のでる確率=1/6という帰無仮説は棄却される
というのは正しいか?
これに即答できんアホはこのスレに書く基礎学力なし。

578:卵の名無しさん
18/12/19 16:34:10.87 3mCL+IGn.net
さて問題です
ガチユリとガチムチの違いを原稿用紙30枚以上で述べよ

579:卵の名無しさん
18/12/19 16:43:43.16 CUSSzVDj.net
>>543
>>540
相関係数を算出に全例調査しているのにp値を出だすのが
何でアホと言われるかわかる?
サイコロが2回続けて1の目出ました。
その確率は1/6*1/6=0.02777で0.05未満だから
サイコロは有意差をもって歪(いびつ)である、
即ち、1の目のでる確率=1/6という帰無仮説は棄却される
というのは正しいか?
オマエ、これに即答できんアホだろ。

580:卵の名無しさん
18/12/19 18:33:47.66 3mCL+IGn.net
問題に即答できんとは
さてはニートだな

581:卵の名無しさん
18/12/19 19:18:36.88 CUSSzVDj.net
>>545
統計スレで統計問題を書けないとはド底辺シリツ頭脳だな。

582:卵の名無しさん
18/12/19 20:12:36.56 3mCL+IGn.net
開業医でないのに開業スレに書きこむ
非医師のニート脳が何か言ってるなW

583:卵の名無しさん
18/12/19 21:12:41.91 CUSSzVDj.net
>>547
開業医スレに
こういうのを書いたら礼を言われましたよ。
URLリンク(egg.2ch.net)

584:卵の名無しさん
18/12/19 21:13:24.16 CUSSzVDj.net
>>547
サイコロが2回続けて1の目が出ました。
その確率は1/6*1/6=0.02777で0.05未満だから
サイコロは有意差をもって歪(いびつ)である、
即ち、1の目のでる確率=1/6という帰無仮説は棄却される
というのは正しいか?
オマエ、これに即答できんアホだろ。

585:卵の名無しさん
18/12/19 21:24:51.44 CUSSzVDj.net
>>547
ニート脳の学力すらないどアホがオマエな。
オマエ、これに即答できんアホだろ。
統計スレなんだから
これに答えてみ!はよ!はよ!
サイコロが2回続けて1の目が出ました。
その確率は1/6*1/6=0.02777で0.05未満だから
サイコロは有意差をもって歪(いびつ)である、
即ち、1の目のでる確率=1/6という帰無仮説は棄却される
というのは正しいか?

586:卵の名無しさん
18/12/19 21:25:04.74 3mCL+IGn.net
非医師ニートくんは
スレ1の説明も読めないようだ
スレリンク(hosp板:1番)

587:卵の名無しさん
18/12/19 21:28:31.22 CUSSzVDj.net
>>551
スルーできないアホがオマエじゃん。
これに即答できんアホだろ。
統計スレなんだから
これに答えてみ!はよ!はよ!
サイコロが2回続けて1の目が出ました。
その確率は1/6*1/6=0.02777で0.05未満だから
サイコロは有意差をもって歪(いびつ)である、
即ち、1の目のでる確率=1/6という帰無仮説は棄却される
というのは正しいか?

588:卵の名無しさん
18/12/19 21:49:54.61 CUSSzVDj.net
>>551
URLリンク(egg.2ch.net)
までの議論(>77が俺)は
非医者ニートとの会話だと思う?
相当頭が悪そうだな。
これに即答できんアホだろ。統計スレなんだから
これに答えてみ!はよ!はよ!
サイコロが2回続けて1の目が出ました。
その確率は1/6*1/6=0.02777で0.05未満だから
サイコロは有意差をもって歪(いびつ)である、
即ち、1の目のでる確率=1/6という帰無仮説は棄却される
というのは正しいか?

589:卵の名無しさん
18/12/19 23:57:49.27 3mCL+IGn.net
今日もニートの歌
URLリンク(m.youtube.com)

590:卵の名無しさん
18/12/20 04:19:04.11 b33W7V5i.net
>>554
URLリンク(egg.2ch.net)
の議論(>77が俺)は
非医者ニートとの会話だと思う?
相当頭が悪そうだな。
これに即答できんアホだろ。統計スレなんだから
これに答えてみ!はよ!はよ!
サイコロが2回続けて1の目が出ました。
その確率は1/6*1/6=0.02777で0.05未満だから
サイコロは有意差をもって歪(いびつ)である、
即ち、1の目のでる確率=1/6という帰無仮説は棄却される
というのは正しいか?

591:卵の名無しさん
18/12/20 08:39:31.68 YfbLPysN.net
>>547
開業医スレに
こういうのを書いたら礼を言われましたよ。
URLリンク(egg.2ch.net)
非医師ニートに教えてもらって例を言ってるのも非医師ニートなのか。二人から褒めるレスが返って来たんだが。

統計スレなんだから!これに答えてみ!はよ!はよ!
サイコロが2回続けて1の目が出ました。
その確率は1/6*1/6=0.02777で0.05未満だから
サイコロは有意差をもって歪(いびつ)である、
即ち、1の目のでる確率=1/6という帰無仮説は棄却される
というのは正しいか?

592:卵の名無しさん
18/12/20 19:49:02.89 LEZzzwW7.net
事務員の日常
URLリンク(www.nicovideo.jp)

593:卵の名無しさん
18/12/20 20:49:22.03 s4OL7Lza.net
>>557
URLリンク(egg.2ch.net)
の議論(>77が俺)は
非医者ニートとの会話だと思う?
相当頭が悪そうだな。
これに即答できんアホだろ。統計スレなんだから
これに答えてみ!はよ!はよ!
サイコロが2回続けて1の目が出ました。
その確率は1/6*1/6=0.02777で0.05未満だから
サイコロは有意差をもって歪(いびつ)である、
即ち、1の目のでる確率=1/6という帰無仮説は棄却される
というのは正しいか?

594:卵の名無しさん
18/12/20 22:04:13.70 aEbz7irJ.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 pr


595:actice 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.



596:卵の名無しさん
18/12/21 07:36:37.26 RO6zf9jd.net
リンゲル液を凍らせた。
室温に放置して固相と液相が共存している時の液相の浸透圧は
元のリンゲル液より高いか低いかを即答できないのが
ド底辺シリツ医な、オマエも即答できんのじゃね。
ある学会でリンゲル液と同じと答えたアホがいたな。

597:卵の名無しさん
18/12/21 08:36:03.18 9yiQ3K91.net
>>559
ド底辺シリツ医にもわかるように原文をつけておけよ。
>>
私は昭和の時代に大学受験したけど、昔は今よりも差別感が凄く、特殊民のための特殊学校というイメージで開業医のバカ息子以外は誰も受験しようとすらしなかった。
常識的に考えて、数千万という法外な金を払って、しかも同業者からも患者からもバカだの裏口だのと散々罵られるのをわかって好き好んでド底辺医に行く同級生は一人もいませんでした。
本人には面と向かっては言わないけれど、俺くらいの年代の人間は、おそらくは8-9割はド底辺医卒を今でも「何偉そうなこと抜かしてるんだ、この裏口バカが」と心の底で軽蔑し、嘲笑しているよ。
当の本人には面と向かっては絶対にそんなことは言わないけどね。
<<

598:卵の名無しさん
18/12/21 12:14:03.53 oQHNIYqE.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.

599:卵の名無しさん
18/12/21 17:03:49.74 JdSFbP47.net
 「イグっ&#9829; イグううぅーー&#9829; 鋼太郎、池袋にイっちゃうのおおぉぉ!」

600:卵の名無しさん
18/12/21 19:52:00.45 b22UegL5.net
age

601:卵の名無しさん
18/12/22 14:08:40.70 Ihwj4uih.net
URLリンク(www.nicovideo.jp)

602:卵の名無しさん
18/12/23 12:00:45.56 9jfRvGpH.net
f=c(179,28,19,13,9,7,6,5,5,4,4,4,5,7,9,10,11,12,14,16,18,19,19,20,21,22,23,24,25,27,28,30,
31,34,37,39,41,43,46,51,57,63,70,77,83,91,99,109,119,130,142,155,167,178,190,202,216,
233,249,265,282,302,326,354,387,420,457,502,549,598,648,700,761,836,927,1026,1142,1284,
1455,1651,1862,2089,2341,2625,2934,3264,3598,3923,4233,4508,4740,4893,4973,5007,4999,
4729,4314,3797,3222,2634,2071,1566,1135,788,523,761)
m=c(191,31,21,13,10,8,8,8,7,7,7,8,9,11,14,17,21,26,32,37,42,46,49,50,50,50,49,49,50,52,55,
58,60,63,65,67,71,76,82,89,97,104,112,122,134,148,165,183,203,224,246,268,294,324,357,
391,425,461,502,549,601,659,722,792,872,958,1052,1147,1239,1331,1433,1546,1663,1783,
1905,2026,2167,2333,2532,2750,2973,3195,3414,3630,3827,4000,4133,4200,4194,4104,3916,
3681,3388,3046,2669,2272,1875,1494,1145,841,589,392,246,144,79,71)
LE <-function(ndx,Y,N0=10^5){ # life expectancy
n=length(ndx)
lx=numeric(n)
lx[1]=N0
for(i in 1:(n-1))
lx[i+1] <- lx[i] - ndx[i]
nqx=ndx/lx
nLx=numeric(n)
for(i in 1:n)
nLx[i] <- mean(c(lx[i],lx[i+1]))
nLx[n]=0
Tx=rev(cumsum(rev(nLx)))
le=Tx/lx
return(round(le[Y+1],1))
}

603:卵の名無しさん
18/12/23 12:02:35.17 9jfRvGpH.net
Y=1:105
sapply(Y,function(x) LE(f,x))
sapply(Y,function(x) LE(m,x))
Y[sapply(Y,function(x) LE(m,x))<5]
sapply(Y,function(x) LE(f,x))
sapply(Y,function(x) LE(m,x))
Y=1:105
Y[sapply(Y,function(x) LE(f,x))<5]
Y[sapply(Y,function(x) LE(m,x))<5]

604:卵の名無しさん
18/12/23 18:51:20.54 u874lVPY.net
メリークリスマス
URLリンク(m.youtube.com)

605:卵の名無しさん
18/12/24 13:54:15.39 muWLHtMs.net
平均余命(縦線)と「平均寿命-現年齢」(点線の斜線)をグラフにしてみた。
URLリンク(i.imgur.com)

606:卵の名無しさん
18/12/24 19:08:03.81 L4b+RVwf.net
URLリンク(m.youtube.com)

607:卵の名無しさん
18/12/24 20:36:41.73 6ZVUEdxI.net
他スレで喚いて荒らしていても
算数クイズくんが
底辺から抜け出すことは
困難を極めます

608:卵の名無しさん
18/12/24 22:12:45.55 lZLAiNKc.net
開業医スレに
こういうのを書いたら礼を言われましたよ。
URLリンク(egg.2ch.net)

609:卵の名無しさん
18/12/24 22:21:17.70 lZLAiNKc.net
>>571
これ理科のクイズ。
リンゲル液を凍らせた。
室温に放置して固相と液相が共存している時の液相の浸透圧は
元のリンゲル液より高いか低いかを即答できないのが
ド底辺シリツ医な、オマエも即答できんのじゃね。
ある学会でリンゲル液と同じと答えたアホがいたな。

610:卵の名無しさん
18/12/25 09:25:20.23 Lx1BgrzL.net
平均余命の信頼区間を計算しようと思ったが
各年齢での人口データがないと無理だな。

611:卵の名無しさん
18/12/26 00:18:55.39 rHMA/ZVK.net
さて国語の問題です
信じて送り出したフタナリ事務員が農家の叔父さんの・・・・・
続きを原稿用紙30枚以上で作成せよ

612:卵の名無しさん
18/12/26 20:53:45.05 0RdixcRl.net
clinical sense の症例が配信されてきた。
Devic病と思ったら経過中にレルミット兆候が出てMSに診断変更。
インターフェロンβ→ナタリズマブの前にステロイドパルスがないのが解せなかったな。

613:卵の名無しさん
18/12/26 23:30:42.66 wApOj4Hd.net
ところでMS法人って
みんな作ってる?

614:卵の名無しさん
18/12/27 14:17:15.36 7ElT0ESI.net
Welch.test=function(n1,n2,m1,m2,sd1,sd2){
T=(m1-m2)/sqrt(sd1^2/n1+sd2^2/n2)
df=(sd1^2/n1+sd2^2/n2)^2 / (sd1^4/n1^2/(n1-1)+sd2^4/n2^2/(n2-1))
p.value=2*pt(abs(T),df,lower.tail = FALSE)
return(c(p.value=p.value,t.statistic=T))
}
# white zoneの視認性と腺管密度の関係
# white zoneが鮮明な癌 窩間部の長さ136±45 (p<0.05)
# white zoneが不鮮明な癌 窩間部の長さ77±22 (p<0.05)
m1=136
sd1=45
m2=77
sd2=22
n2p <- function(n){
Welch.test(n,n,m1,sd1,m2,sd2)[1]
}
n2p=Vectorize(n2p)
plot(2:15,sapply(2:15,n2p),bty='l',pch=19)
abline(h=0.05)

615:卵の名無しさん
18/12/27 14:47:38.86 7ElT0ESI.net
# se=sd/√n , sd=se*√n
n2p_se <- function(n){
sd1=sd1*sqrt(n)
sd2=sd2*sqrt(n)
Welch.test(n,n,m1,sd1,m2,sd2)[1]
}
n2p_se=Vectorize(n2p_se)
n2p_se(5)
plot(n,sapply(n,n2p_se),bty='l',pch=19,ylim=c(0,1))
abline(h=c(0.01,0.05),col=8)
n2p_se(18:22)

616:卵の名無しさん
18/12/27 15:17:44.83 7ElT0ESI.net
# sim
# white zoneが鮮明な癌 窩間部の長さ136±45 (p<0.05)
# white zoneが不鮮明な癌 窩間部の長さ77±22 (p<0.05)
n=7
y1=scale(rnorm(n))*45 + 136
y2=scale(rnorm(n))*22 + 77
library(BEST)
BESTout=BESTmcmc(y1,y2)
# Look at the result:
BESTout
summary(BESTout)
plot(BESTout)
plot(BESTout, "sd")
plotPostPred(BESTout)
plotAll(BESTout, credMass=0.8, ROPEm=c(-0.1,0.1),
ROPEeff=c(-0.2,0.2), compValm=0.5)
plotAll(BESTout, credMass=0.8, ROPEm=c(-0.1,0.1),
ROPEeff=c(-0.2,0.2), compValm=0.5, showCurve=TRUE)
summary(BESTout, credMass=0.8, ROPEm=c(-0.1,0.1), ROPEsd=c(-0.15,0.15),
ROPEeff=c(-0.2,0.2))
pairs(BESTout)

617:卵の名無しさん
18/12/27 15:20:39.76 7ElT0ESI.net
white zoneの視認性と腺窩の深さの関係
white zoneが鮮明な癌 腺窩の深さ 180±47 (p<0.001)
white zoneが不鮮明な癌 腺窩の深さ 81±20 (p<0.001)

618:卵の名無しさん
18/12/27 15:25:25.68 7ElT0ESI.net
>>581
m1=180
sd1=47
m2=81
sd2=20
n2p <- function(n){
Welch.test(n,n,m1,sd1,m2,sd2)[1]
}
n2p=Vectorize(n2p)
n=2:10
plot(n,sapply(n,n2p),bty='l',pch=19)
abline(h=c(0.001,0.01),col=8)
n2p(5:10)

619:卵の名無しさん
18/12/27 15:30:09.78 7ElT0ESI.net
m1=180
sd1=47
m2=81
sd2=20
n2p <- function(n){
Welch.test(n,n,m1,sd1,m2,sd2)[1]
}
n2p=Vectorize(n2p)
n=2:10

# se=sd/√n , sd=se*√n
n2p_se <- function(n){
sd1=sd1*sqrt(n)
sd2=sd2*sqrt(n)
Welch.test(n,n,m1,sd1,m2,sd2)[1]
}
n2p_se=Vectorize(n2p_se)
n=25:35
plot(n,sapply(n,n2p_se),bty='l',pch=19,ylim=c(0,0.05))
abline(h=c(0.001,0.01),col=8)
n2p_se(27:32)

620:卵の名無しさん
18/12/27 15:30:31.81 7ElT0ESI.net
サンプルサイズが推定できて、気分が( ・∀・)イイ!!

621:卵の名無しさん
18/12/27 21:28:20.91 ZlVEHwaK.net
アイドル映像でも見て もちつけ
URLリンク(m.youtube.com)

622:卵の名無しさん
18/12/27 22:51:46.90 gYAr3TuJ.net
qqct.jp
URLリンク(bvdeur67.secure.ne.jp)

623:卵の名無しさん
18/12/28 00:56:18.13 Oc+hTpm3.net
日本にも、
ついに朝鮮殺戮殺人カルト宗教
殺人学会の、
テロ工作拠点が明らかになった!
テロ工作拠点
福山友愛病院
なんと日本警察に朝鮮殺戮殺人学会が侵入し、
カルトによる被害者側を警察が拉致して連れていくテロ工作拠点の一角!!
警察と完全犯罪達成の為に、共同組織犯罪をやっていて、
与えられていたのは、
犯罪ライセンス!!!!!!

薬物大量投与テロ発生!!
なんと故意に八倍!!!
カミサカというテロ医師による、
偽造カルテ作成発覚!!
なんとその為に監禁罪をやっていた!!!

624:卵の名無しさん
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
4分割ですらの正答出せなかった。
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ピストル(当時のお金の単位)の金を賭けて勝負したとする。
 そしてどちらかが先に3点を得たものを勝ちとし、勝った方がかけ金の総額64ピストルをもら
 えるとする。ところが甲が2点、乙が1点を得たとき、勝負が中止になってしまった。
 このとき、二人のかけ金の総額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
ド底辺シリツ医にもわかるように原文をつけておけよ。
>>
私は昭和の時代に大学受験したけど、昔は今よりも差別感が凄く、特殊民のための特殊学校というイメージで開業医のバカ息子以外は誰も受験しようとすらしなかった。
常識的に考えて、数千万という法外な金を払って、しかも同業者からも患者からもバカだの裏口だのと散々罵られるのをわかって好き好んでド底辺医に行く同級生は一人もいませんでした。
本人には面と向かっては言わないけれど、俺くらいの年代の人間は、おそらくは8-9割はド底辺医卒を今でも「何偉そうなこと抜かしてるんだ、この裏口バカが」と心の底で軽蔑し、嘲笑しているよ。
当の本人には面と向かっては絶対にそんなことは言わないけどね。
<<>>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が標準正規分布に従う時、次の値を求めよ。
1:P{Z>u1}=0.05を満たすu1の値 → u1=1.645
2:P{Z<u2}=0.005を満たすu2の値 → u2=-2.576
3: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)に従う時、次の確率を求めよ。
1:P{X>20} → 0.0228
2:P{X<5} → 0.1587
3: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
ド底辺シリツ医にもわかるように原文をつけておけよ。
>>
私は昭和の時代に大学受験したけど、昔は今よりも差別感が凄く、特殊民のための特殊学校というイメージで開業医のバカ息子以外は誰も受験しようとすらしなかった。
常識的に考えて、数千万という法外な金を払って、しかも同業者からも患者からもバカだの裏口だのと散々罵られるのをわかって好き好んでド底辺医に行く同級生は一人もいませんでした。
本人には面と向かっては言わないけれど、俺くらいの年代の人間は、おそらくは8-9割はド底辺医卒を今でも「何偉そうなこと抜かしてるんだ、この裏口バカが」と心の底で軽蔑し、嘲笑しているよ。
当の本人には面と向かっては絶対にそんなことは言わないけどね。

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
戸田A哲君もアメリカで経済学助狂をやっているらしい。
バブル崩壊で安定を求めて理3に入ってしまったらしい。
無駄に齢とっちまったな。

651:卵の名無しさん
19/01/07 20:32:18.68 31hT7o17.net
>>610
正月休みで1週間ぶりに内視鏡したが違和感があるなぁ。
平日は週3やってる。

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
正月休みで1週間内視鏡から離れたら却って体調が悪いなぁ。
金のために不愉快な仕事をしないというのが最大幸福だな。
おかげで自宅で安眠できる。

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
では模範解答をお願いしますね。
原稿用紙30枚以上必要らしいけど。
はよ、はよ!

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
算数クイズができると論文の検証ができるわけね。
内視鏡の洗浄待機時間に 胃と腸 の所見用語集2017を読んでいたらこんな記載があった。(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なのだろう。
±のあとの数字は標準偏差なのか標準誤差なのか判然としない。
検定はt検定(等分散を前提とはできないのでWelch法を使用したとする)で
サンプルサイズは等しいとして±の後の数字は標準偏差あるとすると
> n2p(5:9)
p.value p.value p.value p.value p.value
0.055446074 0.032947710 0.019850610 0.012077956 0.007403284
サンプルサイズは僅か6~8と推定される。

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
どアホにもわかる計算の例示として税込みガソリン価格を上げただけということすらわからんらしいな。
こういう問題に算数が必要と書いてもオマエの頭じゃ正解だせないだろ?
診療所から病院に患者を救急搬送する。
病院から医師搭乗の救急車が診療所に向かっており10時到着予定と連絡が入った。
患者の病態が悪化したら、診療所の普通車で病院に向かい救急車と出会ったら
救急車に患者を移して搬送し病院到着を急ぐという計画を立てた。
普通車から救急車への患者の乗り換えで10分余分に時間がかかる。
道路事情から救急車は病院から診療所への道は
平均時速60kmで、逆方向は平均時速45kmで定速走行する。


691:診療所の普通車は信号待ちもあり平均時速30kmで定速走行する。 何時以降の病態悪化は診療所の車を使わずに救急車の到着を待つ方が病院に早く着くか?



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
どアホにもわかる計算の例示として税込みガソリン価格を上げただけということすらわからんらしいな。
こういう問題に算数が必要と書いてもオマエの頭じゃ正解だせないだろ?
診療所から病院に患者を救急搬送する。
病院から医師搭乗の救急車が診療所に向かっており10時到着予定と連絡が入った。
患者の病態が悪化したら、診療所の普通車で病院に向かい救急車と出会ったら
救急車に患者を移して搬送し病院到着を急ぐという計画を立てた。
普通車から救急車への患者の乗り換えで10分余分に時間がかかる。
道路事情から救急車は病院から診療所への道は
平均時速60kmで、逆方向は平均時速45kmで定速走行する。診療所の普通車は信号待ちもあり平均時速30kmで定速走行する。
何時以降の病態悪化は診療所の車を使わずに救急車の到着を待つ方が病院に早く着くか?

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
どアホにもわかる計算の例示として税込みガソリン価格を上げただけということすらわからんらしいな。
こういう問題に算数が必要と書いてもオマエの頭じゃ正解だせないだろ?
診療所から病院に患者を救急搬送する。
病院から医師搭乗の救急車が診療所に向かっており10時到着予定と連絡が入った。
患者の病態が悪化したら、診療所の普通車で病院に向かい救急車と出会ったら
救急車に患者を移して搬送し病院到着を急ぐという計画を立てた。
普通車から救急車への患者の乗り換えで10分余分に時間がかかる。
道路事情から救急車は病院から診療所への道は
平均時速60kmで、逆方向は平均時速45kmで定速走行する。診療所の普通車は信号待ちもあり平均時速30kmで定速走行する。
何時以降の病態悪化は診療所の車を使わずに救急車の到着を待つ方が病院に早く着くか?

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
どアホにもわかる計算の例示として税込みガソリン価格を上げただけということすらわからんらしいな。
こういう問題に算数が必要と書いてもオマエの頭じゃ正解だせないだろ?
診療所から病院に患者を救急搬送する。
病院から医師搭乗の救急車が診療所に向かっており10時到着予定と連絡が入った。
患者の病態が悪化したら、診療所の普通車で病院に向かい救急車と出会ったら
救急車に患者を移して搬送し病院到着を急ぐという計画を立てた。
普通車から救急車への患者の乗り換えで10分余分に時間がかかる。
道路事情から救急車は病院から診療所への道は
平均時速60kmで、逆方向は平均時速45kmで定速走行する。診療所の普通車は信号待ちもあり平均時速30kmで定速走行する。
何時以降の病態悪化は診療所の車を使わずに救急車の到着を待つ方が病院に早く着くか?

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
どアホにもわかる計算の例示として税込みガソリン価格を上げただけということすらわからんらしいな。
こういう問題に算数が必要と書いてもオマエの頭じゃ正解だせないだろ?
診療所から病院に患者を救急搬送する。
病院から医師搭乗の救急車が診療所に向かっており10時到着予定と連絡が入った。
患者の病態が悪化したら、診療所の普通車で病院に向かい救急車と出会ったら
救急車に患者を移して搬送し病院到着を急ぐという計画を立てた。
普通車から救急車への患者の乗り換えで10分余分に時間がかかる。
道路事情から救急車は病院から診療所への道は
平均時速60kmで、逆方向は平均時速45kmで定速走行する。診療所の普通車は信号待ちもあり平均時速30kmで定速走行する。
何時以降の病態悪化は診療所の車を使わずに救急車の到着を待つ方が病院に早く着くか?

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
内視鏡の洗浄待機時間に 胃と腸 の所見用語集2017を読んでいたらこんな記載があった。(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なのだろう。
±のあとの数字は標準偏差なのか標準誤差なのか判然としない。
検定はt検定(等分散を前提とはできないのでWelch法を使用したとする)で
サンプルサイズは等しいとして±の後の数字は標準偏差あるとすると
> n2p(5:9)
p.value p.value p.value p.value p.value
0.055446074 0.032947710 0.019850610 0.012077956 0.007403284
サンプルサイズは僅か6~8と推定される。
n=7
p.value=0.019850610
としてFalse Positive Riskを計算してみる(非心t分布でのコードは既出)
# 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)
あるんだけど、
非心t分布にはなんにも言及がないんだよなぁ。
非心分布の解説本ってこれを凌駕するのはないな。
サンプルサイズの決め方 (統計ライブラリー)
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
虚飾値がmの値によらないのは興味深いな。

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
カネが欲しいから医学部にいくんや
文句あっか


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