[過去ﾛｸﾞ] 臨床統計もおもしろいですよ、その2 (1002ﾚｽ)
上下前次1-新
このｽﾚｯﾄﾞは過去ﾛｸﾞ倉庫に格納されています｡
次ｽﾚ検索 歴削→次ｽﾚ 栞削→次ｽﾚ 過去ﾛｸﾞﾒﾆｭｰ
912: 、 2020/02/12(水)20:50 ID:dWksx8Sa(1) AAS
こんな凄いスレッドあったんですね、読み返してみます、
913(1): 2020/02/12(水)21:02 ID:iCMxhp27(1) AAS
ここは国試浪人の事務員のスレ
914: 2020/02/12(水)21:21 ID:6n8f5vvH(5/6) AAS
>>913
手コキ問題の答はまだか？
915: 2020/02/12(水)22:17 ID:6n8f5vvH(6/6) AAS
>>911
まず、シミュレーションプログラムを作る。
p=1-c(10,16,34,9,10,26)/125 #進級確率
q=1-p
sim <- function(){
g=1 # 学年,7:卒業
y=0 # 在学年数
while(g<7 & y<12){
g=g+sample(1:0,1,prob=c(p[g],q[g]))
y=y+1
省5
916: 2020/02/13(木)13:57 ID:DpKqMtsQ(1/3) AAS
p=1-c(10,16,34,9,10,26)/125 #進級確率
q=1-p
sim <- function(){
g=1 # 学年,7:卒業
y=0 # 在学年数
while(g<7 & y<12){ # 卒業前、在学１２年未満なら
g=g+sample(1:0,1,prob=c(p[g],q[g])) # 進級判定して
y=y+1 # 在学年数を増やす
}
c(g,y) # 学年と在学年数を返す
省7
917: 2020/02/13(木)13:58 ID:DpKqMtsQ(2/3) AAS
> nrow(R12[R12[,1]<7,])/k # 退学確率
[1] 0.001039
918: 2020/02/13(木)21:06 ID:DpKqMtsQ(3/3) AAS
p=1-c(10,16,34,9,10,26)/125 # 進級確率
q=1-p # 留年確率
p6=prod(p) # 6年で卒業確率
f <- function(x) prod(q[x])
pf=c(6,p6)
for(i in 7:12){
cmb=combinations(6,i-6,rep=T)
pf=rbind(pf,c(i,p6*sum(apply(cmb,1,f))))
}
pf
省2
919(1): 2020/02/14(金)01:05 ID:9rZr8ZNE(1/4) AAS
library(gtools)
p=1-c(10,16,34,9,10,26)/125 # 進級確率
q=1-p # 留年確率
p6=prod(p) # 6年で卒業確率
f <- function(x) prod(q[x]) # x:留年学年の配列→その組合せでの確率
pf=c(6,p6) # 在学年数,その確率
for(i in 7:12){#在学年数i 7~12において
cmb=combinations(6,i-6,rep=T) # 1:6から留年学年の組合せi-6個を重複を許して選ぶ
pf=rbind(pf,c(i,p6*sum(apply(cmb,1,f)))) #その組合せでの確率の総和に6年分の進級確率を乗ずる
}
省4
920: 2020/02/14(金)01:13 ID:9rZr8ZNE(2/4) AAS
理論値
> (yg=sum(pf[,1]*pf[,2])/P) # 卒業までの平均在学年数 Σ x*pdf
[1] 7.027985
シミュレーション値
> mean(R7[,2]) # 卒業生の在学年数
[1] 7.027643
合致しているから、どちらも正しそう。
921(1): 2020/02/14(金)18:43 ID:x9TXrlfR(1) AAS
思えば東京に住んでる奴はそれだけで経済的に優遇されてるな
田舎に住んでる貧乏人は学力があっても地元国立一択だったり
学力があっても就職コースだったりするからなぁ
922: 2020/02/14(金)20:30 ID:9rZr8ZNE(3/4) AAS
>>921
東京で地元の国立に進学するのは容易ではないと思う。
923: 2020/02/14(金)20:33 ID:9rZr8ZNE(4/4) AAS
>>919
この問題は数学板に投稿してもプログラムの練習にしかならんな。
924(1): 2020/02/14(金)22:03 ID:wJciaHaO(1) AAS
難し、ネットで医療統計のページ見つけたが、
エクセルシートが出てくるだけだし、χ２検定だけでも
俺みたいなのに判るように解説してほしい、
MRが言う当社の薬剤は他のよりも優れていない確率は5％
って言うのはホントなんだろうか？違うよね
925: 2020/02/15(土)07:37 ID:6ICubIld(1) AAS
>>924
Χ２乗検定が知りたいのか
p値で判定することの欺瞞性を知りたいのか
非劣性検定の問題点を知りたいのか
どれ？
926: 2020/02/15(土)15:24 ID:BqM3deVD(1) AAS
具体的にはグーフィス（持田）の製品紹介にある
グーフィスはプラセボ群に比して有意に大きな値をしめしました（ｐ＜0.00０１共分散解析）
次行にはほぼ同じ文面で同じ数字でFisherの正確検定
更にしたには Wilconxonの順位和検定とあります
なんのことやら、
927: 2020/02/15(土)16:05 ID:9KqLKmld(1/2) AAS
分散分析は郡内分散と群間分散の比がF分布に従うことをつかった検定。
（この比をF ratioというのが色っぽいｗ）
統計の教科書読めばFisherもWilcoxonも載っているから独習してください。
頻度主義統計を理解する前に、
p値で判定することの欺瞞性や非劣性検定の問題点を論じても混乱するだけだと思います。
928: 2020/02/15(土)19:46 ID:9KqLKmld(2/2) AAS
新型コロナウイルスが新たに６７人で陽性であった。
検査可能数を100~1000と仮定したときの検査数の期待値と９５％信頼区間を求めよ
"
f <- function(N,conf=0.95,print=FALSE){
# N=100 ; conf=0.95
m=1
r=67
n=r:N
pmf=choose(r-1,m-1)/choose(n,m) #Pr(max=60|n)
# plot(n,pmf,ylab='probability',bty='l')
省13
929: 2020/02/16(日)08:15 ID:dI8Llzzn(1) AAS
コロナ患者はおことわり
当たり前だよなぁ
普通の病院に来られても検査はできん
医師なら当然知ってることだ
930: 2020/02/18(火)12:09 ID:OUmetWs8(1/2) AAS
製薬会社の社員って基本的に自社製品の宣伝しかしないよね
金儲けのため病院訪問してんだから当たり前かも知らんが
コンプ薬屋クンとかサラリーマン時代はどんな営業してたのかねぇ？
931: 2020/02/18(火)19:48 ID:OUmetWs8(2/2) AAS
人生に失敗してる奴ほど反日野党に肩入れするよねW
上下前次1-新書関写板覧索設栞歴
あと 71 ﾚｽあります
ｽﾚ情報 赤ﾚｽ抽出 画像ﾚｽ抽出 歴の未読ｽﾚ AAｻﾑﾈｲﾙ
ぬこの手 ぬこTOP 0.010s