高校数学の質問スレ（医者・東大卒専用） Part438 (991ﾚｽ)
上下前次1-新
抽出解除 ﾚｽ栞
ﾘﾛｰﾄﾞ規制です｡10分ほどで解除するので､他のﾌﾞﾗｳｻﾞへ避難してください｡
44: 2024/08/15(木)10:18:26.54 ID:5ggTekS4(1/2) AAS
bWin <- function(x){# 1 beats 3, 3 beats 2, 2 beats 1
if(length(unique(x))!=2 ) return(0) # no winner
u=sort(unique(x))
if(all(u==c(1,2))) return(sum(x==2)) # how many winners who won by 2
if(all(u==c(2,3))) return(sum(x==3)) # how many winners who won by 3
if(all(u==c(1,3))) return(sum(x==1)) # how many winners who won by 1
}
janken.till.winner.sim <- function(n){ # janken till someone wins
if(n<2) return(NA)
x=sample(1:3,n,replace = TRUE) # janken by n people
count=1 # counter of janken game
while(bWin(x)==0){ # if draw try again
x=sample(1:3,n,replace = TRUE) # janken again till someone win
count=count+1
}
c(bWin(x),count) # return how many winners and counter
}
janken.till.chamipion.sim <- \(n){
re=janken.till.winner.sim(n)
count=re[2]
winner=re[1]
while(winner>1){
re=janken.till.winner.sim(n)
count=count+re[2]
winner=re[1]
}
count
}
res=replicate(1e5,janken.till.chamipion.sim(9))
summary(res)
hist(res,col=2,breaks='scott',freq=FALSE,xlab='count',main='')
BEST::plotPost(res,showMode = TRUE,breaks='scott')
244: 2024/10/26(土)21:19:37.54 ID:dZVTecny(1/3) AAS
n100=100
n10=10
n9=9
# P[9|m]
p9_m=\(m) choose(m,n9)*choose(n100-m,n10-n9)/choose(100,10)
p9_m=Vectorize(p9_m)
P9m=p9_m(1:100)
plot(1:100,P9m)
# P[m|9]=P[9|m]P[m]/P[9]
# P[9]=P[9|m]P[m]+P[9|!m]P[!m]
# P[m]=1/101
# P[!m]=100/101
# P[9|!m] = P9m[-m]
pm_9=\(m) p9_m(m)/101 / (p9_m(m)/101+sum(P9m[-m])*100/101)
pm_9=Vectorize(pm_9)
auc=sum(pm_9(1:100))
x=1:100
sum(x*pm_9(x)/auc)
plot(x,pm_9(x))
315: 2024/11/17(日)20:37:36.54 ID:WmoqT3ZI(2/2) AAS
数直線上に動点Pがあり、
時刻0には原点に居て、速度が位置の関数としてv(x)= ax^2+bx +c で与えられているときの
時刻tでのPの位置を求めよ。
343: 2024/11/29(金)16:54:45.54 ID:a8qsJeEV(1) AAS
画像ﾘﾝｸ
402: 2024/12/08(日)22:22:28.54 ID:Yt2vU9nT(2/2) AAS
>>401
検査しながら説明するとかほざいてたやつ？
麻酔はアンタがするんだろ？
検査中は質問とかできないしそんなのできる訳ないだろアホかよ
504: 01/03(金)22:47:45.54 ID:KUn71WA0(2/2) AAS
画像ﾘﾝｸ
694: 02/06(木)20:30:35.54 ID:r+g7n4SK(4/5) AAS
>>692
アンタはスレタイすらろくに読めないアホだろうが
705(5): 02/08(土)08:09:53.54 ID:JbtcoP/Y(1) AAS
裏口シリツ医にはとうてい無理な臨床応用問題
2025年版の当直医マニュアルが届いたので読んでいたら、
インフルエンザ感染咽頭の所謂イクラサインをAIに判定させる装置（商品名nodoca）に関して次のような記載があった。
>>
内視鏡用テレスコープを用いた咽頭画像等解析は発症12時間以内で感度、特異度ともに高いが
発症24時間を過ぎると抗原検査の方が感度が高い。
<<
どの程度の差があるのかと検索したらnodocaの添付文書に次のような図が掲載されていた。
画像ﾘﾝｸ
このグラフから数値を読み取って、好みの統計処理ソフトを用いて
nodocaと抗原検査の感度の比の分布を算出し、その比の95%CIが1を跨ぐかどうかで下記の主張が有意かどうかを判断せよ。
(1) 内視鏡用テレスコープを用いた咽頭画像等解析は発症12時間以内で感度が高い。
(2) 発症24時間を過ぎると抗原検査の方が感度が高い。
定性的な記述を定量的に裏付けられるかという問題。
感度の比の分布が作成できれば、臨床的に有意義な差を設定（例えば感度が1.01倍は意味がないが、1.5倍なら有用とか）してその確率も算出できる。
763: 03/25(火)12:07:55.54 ID:JvJVqzfB(1/6) AAS
"
赤玉4個、黒玉3個、白玉2個、青玉1個を空箱なしで3つの箱にいれる。
箱を区別しないとき、入れ方は何通りあるか？
"
rm(list=ls())
# vector から remove を除いた残りを返す
fn_rest <- function(vector,remove) Reduce(function(x, y) x[-match(y, x)], remove, init = vector)
fn_rest(c(1,2,2,3,3,3),c(2,3,3))
library(RcppAlgos)
v=rep(1:4,4:1)
pa=partitionsGeneral(10,3,repetition=TRUE) ; pa
pa[1,] # 1 1 8
c1=unique(comboGeneral(v,2)) ; c1
matrix(c("赤","黒","白","青")[c1],ncol=2) |> noquote()
re1=nrow(c1) ; re1
pa[2,] # 1 2 7
c2=unique(permuteGeneral(v,3)) ; c2
ans2=c2[apply(c2,1,\(x) x[2]<=x[3]),] ; ans2
matrix(c("赤","黒","白","青")[ans2],ncol=3) |> noquote()
re2=nrow(ans2) ; re2
pa[3,] # 1,3,6
c3=unique(permuteGeneral(v,4)) ; c3
ans3=unique(t(apply(c3,1,\(x) c(x[1],sort(x[2:4]))))) ; ans3
matrix(c("赤","黒","白","青")[ans3],ncol=4) |> noquote()
re3=nrow(ans3) ; re3
pa[4,] # 1 4,5
c4=unique(permuteGeneral(v,5)) ; c4
ans4=unique(t(apply(c4,1,\(x) c(x[1],sort(x[2:5]))))) ; ans4
matrix(c("赤","黒","白","青")[ans4],ncol=5) |> noquote()
re4=nrow(ans4) ; re4
pa[5,] # 2 2 6
c4=unique(comboGeneral(v,4)) ; c4
ans5=NULL
for(i in 1:nrow(c4)){
c22=c4[i,] ; c22
c2=unique(comboGeneral(c22,2)) ; c2
re22=cbind(c2,t(apply(c2,1, \(x) fn_rest(c22,x)))) ; re22
ans22=re22[apply(re22,1,\(x) sum(c(10,1)*x[1:2]) <= sum(c(10,1)*x[3:4])),] ; ans22
ans5=rbind(ans5,ans22)
}
matrix(c("赤","黒","白","青")[ans5],ncol=4) |> noquote()
re5=nrow(ans5) ; re5
pa[6,] # 2 3 5
c6=unique(permuteGeneral(v,5)) ; c6
ans6=unique(t(apply(c6,1,\(x) c(sort(x[1:2]),sort(x[3:5]))))) ; ans6
matrix(c("赤","黒","白","青")[ans6],ncol=5) |> noquote()
re6=nrow(ans6) ; re6
810: 04/30(水)08:07:51.54 ID:wedVH8wl(7/10) AAS
# Find data where the overall Fisher's exact test is not significant,
# but at least one pairwise Fisher's exact test (with Bonferroni correction) is significant.
res_no_overall_sig_pairwise_sig <- NULL
while (is.null(res_no_overall_sig_pairwise_sig) || res_no_overall_sig_pairwise_sig$fisher_pg > alpha || res_no_overall_sig_pairwise_sig$min_pairwise_p_bonf > alpha) {
res_no_overall_sig_pairwise_sig <- sim_fisher()
}
cat("Data where overall Fisher's test is not significant, but pairwise is:\n")
print(res_no_overall_sig_pairwise_sig)
cat("\n")
# Find data where the overall Fisher's exact test is significant,
# but all pairwise Fisher's exact tests (with Bonferroni correction) are not significant.
res_overall_sig_no_pairwise_sig <- NULL
while (is.null(res_overall_sig_no_pairwise_sig) || res_overall_sig_no_pairwise_sig$fisher_pg < alpha || res_overall_sig_no_pairwise_sig$min_pairwise_p_bonf < alpha) {
res_overall_sig_no_pairwise_sig <- sim_fisher()
}
cat("Data where overall Fisher's test is significant, but pairwise is not:\n")
print(res_overall_sig_no_pairwise_sig)
options(warn = 0)
873: 06/12(木)18:20:13.54 ID:QiRqli9X(1) AAS
>>871
バグだと思ったのにIssueもあげないの？
スレの私的利用といい、マジで他人の成果やタイトルに乗っかるだけの寄生虫じゃん
税金も年金もコイツに使うだけ無駄だよ
889: 06/19(木)12:46:46.54 ID:q4l2IwBj(1) AAS
>>885
あんたは早く病院医者板のスレで謝罪しろ
上下前次1-新書関写板覧索設栞歴
ｽﾚ情報 赤ﾚｽ抽出 画像ﾚｽ抽出 歴の未読ｽﾚ
ぬこの手 ぬこTOP 0.399s