高校数学の質問スレ（医者・東大卒専用） Part438 (991ﾚｽ)
上下前次1-新
抽出解除 ﾚｽ栞
ﾘﾛｰﾄﾞ規制です｡10分ほどで解除するので､他のﾌﾞﾗｳｻﾞへ避難してください｡
51: 2024/08/15(木)19:23:19.21 ID:ysyyeRB1(1/2) AAS
j[n_] :=( (* n人でジャンケンして勝者が決まるまでの回数と勝者の数*)
count=0;
Until[Length@Union@a==2,a=RandomChoice[Range[3],n];count++];
b=Sort@Union@a;
If[b=={1,2}, winners=Count[a,2]];
If[b=={2,3}, winners=Count[a,3]];
If[b=={1,3}, winners=Count[a,1]] ;
{winners,count}
)
sim[n_] :=((* 勝者が一人になるまでの回数 *)
For[{winner,counts}=j[n],winner>1k=j[winner];winner=k[[1]];counts=counts+k[[2]]];
counts
)
res9=Table[sim[9],10^5];
Histogram[res9]
Mean[res9]//N
Median[res9]
225: 2024/09/21(土)07:06:35.21 ID:cdm6DP8+(1) AAS
library(binom)
ci=binom.confint(324-300,324)
lu=unlist(ci[11,5:6])
LearnBayes::beta.select(list(p=0.025,x=lu[1]),list(p=0.975,x=lu[2]))
# 信頼区間からβ分布の形状係数を算出し代表値を返す
ci2ab=\(l,u,verbose=FALSE,cl=0.95){ # CI : [l,u], cl : confidence level
if(l==u) return(NA)
options(warn = -1)
HDI=\(InvCDF=qbeta,cred=0.95,...){
opt=optimize(\(p) InvCDF(p+cred,...) - InvCDF(p,...),c(0,1-cred))
lwr=InvCDF(opt$min,...)
upr=lwr+opt$obj
c(lwr,upr)
}
f=\(ab){
LU=HDI(qbeta,cred=cl,shape1=ab[1],shape2=ab[2])
(LU[1]-l)^2 + (LU[2]-u)^2
}
opt=optim(runif(2,1,100),f)
opt=optim(opt$par,f)
par=opt$par
lu=HDI(qbeta,cred=cl,shape1=par[1],shape2=par[2])
if(verbose){
mean=par[1]/sum(par)
median=qbeta(0.50,par[1],par[2])
mode=(par[1]-1)/(sum(par)-2)
cat('α =',round(par[1],2),' β =',round(par[2],2),'\n')
cat('mean =',round(mean,3),' median =',round(median,3))
if(par[1]>1 & par[2]>1) cat(' mode =',round(mode,3))
cat('\nlower =',round(lu[1],3),' upper =',round(lu[2],3),'\n') curve(dbeta(x,par[1],par[2]),type='h',col=2,n=250,bty='l',ann=FALSE,axes=FALSE)
axis(1)
}
options(warn = 0)
invisible(par)
}
ab=ci2ab(lu[1],lu[2])
k=1e5
p=rbeta(k,ab[1],ab[2])
# 検査陽性の事後確率
postp=\(p,s,t) p*s/ (1-t+p*(s+t-1)) # p:事前確率 s:感度 t:特異度
# 検査陰性の事後確率
postn=\(p,s,t) p*(s-1)/(-t+p*(s+t-1)) # p:事前確率 s:感度 t:特異度
# 尿素呼気試験(感度90-100% 特異度80-99%) 便中ピロリ菌抗原 (感度90-98% 特異度87-100%)
abs=ci2ab(0.90,1.00)
abt=ci2ab(0.80,0.99)
s=rbeta(k,abs[1],abs[2])
t=rbeta(k,abt[1],abt[2])
post1=postn(p,s,t)
abs=ci2ab(0.90,0.98)
abt=ci2ab(0.87,1.00)
s=rbeta(k,abs[1],abs[2])
t=rbeta(k,abt[1],abt[2])
post2=postn(post1,s,t)
1/mean(post2)
1/median(post2)
hist(post2,freq=FALSE,breaks='scott',ann=F,axes=F) ; axis(1)
229: 2024/09/28(土)11:06:24.21 ID:1u+u54qN(1) AAS
尿瓶ジジイ都合の悪いレスに噛みついたところで速攻で論派されるのでここでコソコソ書き込むしかできない模様
312: 2024/11/16(土)06:59:43.21 ID:Si+9FQ0m(1) AAS
#
alphabet=c(letters,LETTERS)
a2n=Vectorize(\(a) which(alphabet==a))
w2n=\(w){
y=unlist(strsplit(w,''))
a2n(y)
}
s="supercalifragilisticexpialidocious"
n=w2n(s)
k=1e3
p1=replicate(k,mean(replicate(k,all(diff(sample(n))!=0))))
hist(p1,freq=F,breaks='scott',main='',axes=F,ann=F,col=2) ; axis(1)
mean(p1)
mean(replicate(k^2,all(diff(sample(n))!=0)))
s="Supercalifragilisticexpialidocious"
n=w2n(s)
k=1e3
p1=replicate(k,mean(replicate(k,all(diff(sample(n))!=0))))
hist(p1,freq=F,breaks='scott',main='',axes=F,ann=F,col=3) ; axis(1)
mean(p1)
mean(replicate(k^2,all(diff(sample(n))!=0)))
313(1): 2024/11/17(日)19:57:55.21 ID:WmoqT3ZI(1/2) AAS
数直線上に動点Pがあり、
時刻0には原点に居て、速度が位置の関数としてv(x)= ax+b で与えられているときの
時刻tでのPの位置を求めよ。
画像ﾘﾝｸ
361(1): 2024/12/01(日)16:43:54.21 ID:URtmhVxU(4/7) AAS
(* decimal integer to b-based digits 10進法ｘをｂ進法表記の文字リストに返還*)
i2d[x_,b_:16] :=(
digits=Flatten@{"0","1","2","3","4","5","6","7","8","9",Alphabet[]};
r=List@Mod[x,b];
q=Floor[x/b];
While[q > 0,PrependTo[r,Mod[q,b]];q=Floor[q/b]];
digits[[r+1]]
)
(* b-based digits to decimal integer ｂ進法表記文字列xを10進法の数に返還 *)
d2i[x_,b_:16] :=(
tonum[char_]:=(
digits=Flatten@{"0","1","2","3","4","5","6","7","8","9",Alphabet[]};
Position[digits,char][[1]][[1]]-1
);
ss=StringSplit[x,""];
d10=tonum /@ ss;
le=Length@ss;
Table[(Reverse@d10)[[i]]*b^(i-1),{i,1,le}]//Total
)
(* ハーシャッド数（ハーシャッドすう、英: harshad number）とは、自然数の各位の数字和が元の数の約数に含まれている自然数である。*)
harshadQ[n_,b_] := Divisible[n,Total[d2i[#,b]&/@i2d[n,b]]]
solve[start_:"1",end_:"ff",b_] :=(
li=Select[Range[d2i[start,b],d2i[end,b]],harshadQ[#,b]&];
re=i2d[#,b]& /@ li;
StringJoin /@ re
)
415: 2024/12/11(水)10:04:45.21 ID:qk9aTmTZ(3/6) AAS
"
注射液8本のうち2本が致死的な毒薬である場合
ネズミ1: 注射液1, 2, 3, 4
ネズミ2: 注射液1, 2, 5, 6
ネズミ3: 注射液1, 3, 5, 7
に記載された投与法で実験したら全部のネズミが死亡してしまった。
可能な２つの毒薬の組み合わせをすべて列挙せよ。
この投与法で２個の毒薬が特定できる場合をすべて求めて
そのときの毒薬の組み合わせと合わせて列挙せよ。
"
rm(list=ls())
mat=matrix(
c(1, 2, 3, 4,
1, 2, 5, 6,
1, 3, 5, 7),nrow=3,by=TRUE)
cm=t(combn(8,2)) ; cm
lost=\(y){
tmp1=(apply(mat,1, \(x) y[1] %in% x) | apply(mat,1, \(x) y[2] %in% x))
tmp2=as.integer(tmp1)
tmp2 %*% c(1,2,4)|> as.vector() -> tmp3
tmp3
}
i=apply(cm,1,lost) ; i
cbind(cm,i)
lapply(1:7,\(x) cm[i==x,])
453(1): 2024/12/22(日)12:35:08.21 ID:rWRciXK9(4/8) AAS
>>452
おまけ、R言語でのクズコード
理工系卒ならWolframで書けるはず。
但し、Fラン卒は除く。
r=0.3
source("toolmini.R")
p=0.5+0.5i
q=-0.5+0.5i
r0=-0.5-0.5i
s=0.5-0.5i
fP=\(t) r*exp(-1i*t)+p
fQ=\(t) r*exp(-1i*t)+q
fR=\(t) r*exp(-1i*t)+r0
fS=\(t) r*exp(-1i*t)+s
f=\(tt,verbose=FALSE,...){
P=fP(tt[1])
Q=fQ(tt[2])
R=fR(tt[3])
S=fS(tt[4])
if(verbose){
Plot(-1,1,zero=FALSE)
Polygon(p,q,r0,s,Col=8)
Cir(p,r,col=8)
Cir(q,r,col=8)
Cir(r0,r,col=8)
Cir(s,r,col=8)
Polygon(P,Q,R,S,...)
pta(P);pta(Q);pta(R);pta(S)
}
ABC2S(P,Q,R)+ABC2S(P,S,R)
}
opt=optim(runif(4,-pi,pi),f)
opt=optim(opt$par,f)
tt=opt$par
min=f(tt,TRUE,Col=2)
opt=optim(runif(4,-pi,pi),f,control=list(fnscale=-1))
opt=optim(opt$par,f,control=list(fnscale=-1))
tt=opt$par
max=f(tt,TRUE,Col=4)
c(min,max)
501: 01/01(水)11:12:05.21 ID:vYaE/f7n(1) AAS
>>500
年明けても日本語読めないみたいだな
592(3): 01/21(火)02:36:13.21 ID:eMqBmBFs(2/2) AAS
>>590
ここは医師、東大卒専用に立てられたスレなんだが。
医師＆東大卒でなく医師｜東大卒と解釈して俺は投稿している。
旧二期校の医学部卒なので。
625: 01/24(金)16:26:04.21 ID:WCbmzKUH(1) AAS
>>622
>>624
まんまと乗せられてやがる
885(3): 06/19(木)09:36:40.21 ID:Kb+ol8z3(1) AAS
プログラムが弄れる医者や東大卒なら自力で算出できる問題。
Fランや裏口シリツ医には無理。
直線上にならんだ池が６個ある。１〜６と命名する。
池にはカエルがいる。
観察の結果、カエルは翌日には５０％の確率で隣の池に移る。
隣に２つの池がある場合、どちらの池に移る確率は等しいとする。
（１）池１にカエルが1匹いるとき、このカエルが１００日後はどの池にいる確率が最も高いか。
（２）カエルが池１に1匹、池２に2匹、池3に3匹、池4に4匹、池5に5匹、池6に6匹いるとする。
100日後にもっとも多くのカエルがいる確率が高いのはどの池か。その池のカエルの数の期待値とともに答えよ。
982: 08/08(金)18:36:50.21 ID:sJI18Py3(2/4) AAS
プロトタイプを自分で作ってAIと対話的にデバッグして完成。
丸投げするとこちらがデバッグさせられる羽目になる。
上下前次1-新書関写板覧索設栞歴
ｽﾚ情報 赤ﾚｽ抽出 画像ﾚｽ抽出 歴の未読ｽﾚ
ぬこの手 ぬこTOP 0.063s