臨床統計もおもしろいですよ、その3

1 ：卵の名無しさん：2020/03/05(木) 20:17:05.03 ID:naSB8128.net

　
　内科認定医受験の最低限の知識、
　製薬会社の示してくる臨床データ、
　論文の考察、
　論文を書くときの正当性、
　というのが、臨床統計の今までの目的の大きい部分でしたが、
　
　AI=機械学習の基本も、結局は統計学と確率に支配されます。
　そういう雑多な話をするスレです。
　
※前スレ
臨床統計もおもしろいですよ、その１
https://egg.2ch.net/test/read.cgi/hosp/1493809494/
臨床統計もおもしろいですよ、その2
https://egg.5ch.net/test/read.cgi/hosp/1540905566/

239 ：卵の名無しさん：2020/12/29(火) 00:31:19.09 ID:RkgAFAj8.net

エセ医者がようやるわ。

240 ：卵の名無しさん：2020/12/29(火) 07:35:14.53 ID:v/EM6r6h.net

>>239
では、同業者の意見を聞いてみましょう。

https://egg.5ch.net/test/read.cgi/hosp/1592662437/73

73 卵の名無しさん sage 2020/06/23(火) 13:24:47.79 ID:riQXI/fH
宮廷卒だけど、一括りに医師免許と言ってるが、私大卒など医者とは思ってへんよ
私大入学というインチキを経由したイシャモドキが、あんま調子のんなや

241 ：卵の名無しさん：2020/12/29(火) 09:55:53.10 ID:RkgAFAj8.net

自演w

242 ：卵の名無しさん：2020/12/29(火) 09:56:56.13 ID:RkgAFAj8.net

面白いねたしかに
AIみたいに定型文かコピペしか送ってこないんだもん

243 ：卵の名無しさん：2020/12/29(火) 10:03:17.22 ID:RkgAFAj8.net

じゃあそんな君にこれをプレゼント！

うりゅうひろゆきの本スレにようこそ！
このウンコスレは基本的にさらしアゲて使用します！
ココは医師真似事務員ジイさんが、医師妬みの挙句、テメェで勝手に自称で医科歯科大を卒業した医師って設定で
妄想願望日記を発表するスレです
お医者さんのマトモなレスなど全く不要！
ゴミ箱まがいの公衆便所掃き溜め痰壷下品スレとして使用しましょうね
そんなわけで、いまだ懲りない恥さらし（自称）医科歯科事務員ジジイ登場〜〜

│　　　　　　彡川川川三三三ミ〜 　プウゥ〜ン
｜　　　　　　川｜川／　　＼｜〜　ポワ〜ン 　　　＿＿＿＿＿＿＿＿
｜　　　　　‖｜‖　◎---◎｜〜 　　　　　　　／
｜　　　　　川川‖　　　　3　　ヽ〜　　　　　　＜ 僕はうりゅうひろゆき
｜　　　　　川川　　　∴）д（∴）〜 　　　　　　＼＿＿＿＿＿＿＿＿
｜　　　　　川川　　　　　　〜　/〜　カタカタカタ
｜　　　　　川川‖　　　　〜　/‖　＿＿＿＿＿
｜　　　　川川川川＿＿＿／‖　 |　　|￣￣＼　＼
｜　　　　　　／　　　　　　　＼＿＿|　　|　　　　|￣￣|
｜　　　　　/　　＼医科歯科 　　　 |　　|　　　　|＿＿|
｜　　　　　|　＼　　　　　　|つ　　　|＿_|＿＿／　／
｜　　　　 /　医師うらやましいぜ・・　|￣￣￣￣|　　〔￣￣〕
　　　　　　　よし、医科歯科の医者と称して
　　　　　　　私立のお医者様を攻撃だ！なんたって、医科歯科の看板借りてんだからな
　　　　　　　え？証明？？？自称だから医科歯科の学生証、卒業証書あるわけねえよ
　　　　　　　同期の医師あげてみろ、ったて、そもそも医師じゃねえし居ねえよバカ
　　　　　　　でも病院でもらった、医科歯科の封筒はもってま〜す！これを証拠にしてやるぜ
　　　　　　　専門医？欲しいけど・・いいや！そんなの必要ねえ！ 俺は事務員だし

244 ：卵の名無しさん：2021/01/28(木) 15:59:37.22 ID:RbXTChKG.net

某院でn人の患者感染。死亡率pとすると
予想される死亡者数とその95%信頼区間はいくらになるか？
"
p=0.20
n=100
(CI=HDInterval::hdi(qbinom,size=n,prob=p))
y=dbinom(0:n,n,p)
plot(0:n,y,type='h',bty='l')
sum(dbinom(CI[1]:CI[2],n,p))

245 ：卵の名無しさん：2021/02/01(月) 18:02:03.90 ID:e24XNIzt.net

「Googleドライブ」フォルダをDドライブにする

http://win10faq.blog.jp/archives/2208311.html

246 ：卵の名無しさん：2021/02/11(木) 08:40:29.92 ID:Sfff+JMX.net

ウリュウのジジイ、5chでも社会でもまるで相手にされず！

247 ：卵の名無しさん：2021/05/03(月) 20:55:00.69 ID:ukTrFNeY.net

1 meter * 1 meter * 1 milli meter (water) = 1kg

1 mmHg柱/m2 = 13.6kg重/m2
13.6*9.8 N/m2 = 133.3 Pa =1.33hPa
1 hPa = 1/1.33 mmHg柱 = 0.752 mmHg柱

1 mmH20 = 1 kg重/m2
1*9.8 N/m2 = 9.8 Pa = 0.98 hPa
1 hPa = 1/0.98 = 1.02 mmH2O

248 ：卵の名無しさん：2021/05/17(月) 17:24:47.71 ID:vLTAoLk6.net

https://news.mynavi.jp/article/20160317-a313/

about:config

privatebrowsing.autostartと入力して設定対象を絞り込み

249 ：卵の名無しさん：2021/06/16(水) 09:36:16.08 ID:Iaydx6qX.net

ほんと庭なしバカ安達はくそ犬はうるせーわ道路でばかでかい声で会話するわやりたい放題
家族みんなデブだし

250 ：卵の名無しさん：2021/06/25(金) 06:15:51.74 ID:9PeLNEbV.net

【HKEY_CURRENT_USER > Control Panel > Personazation > Desktop Slideshow】

251 ：卵の名無しさん：2021/06/25(金) 17:47:33.56 ID:OGI6KuDR.net

尿瓶ジジイ＝トケジ＝ウリュウ＝ニセ内視鏡医　笑

252 ：卵の名無しさん：2021/07/03(土) 07:35:22.04 ID:C912OJj6.net

tryCatch has a slightly complex syntax structure. However, once we understand the 4 parts which constitute a complete tryCatch call as shown below, it becomes easy to remember:

expr: [Required] R code(s) to be evaluated

error : [Optional] What should run if an error occured while evaluating the codes in expr

warning : [Optional] What should run if a warning occured while evaluating the codes in expr

finally : [Optional] What should run just before quitting the tryCatch call, irrespective of if expr ran successfully, with an error, or with a warning

tryCatch(
expr = {
# Your code...
# goes here...
# ...
},
error = function(e){
# (Optional)
# Do this if an error is caught...
},
warning = function(w){
# (Optional)
# Do this if an warning is caught...
},
finally = {
# (Optional)
# Do this at the end before quitting the tryCatch structure...
}
)

253 ：卵の名無しさん：2021/07/03(土) 07:36:30.21 ID:C912OJj6.net

log_calculator <- function(x){
tryCatch(
expr = {
message(log(x))
message("Successfully executed the log(x) call.")
},
error = function(e){
message('Caught an error!')
print(e)
},
warning = function(w){
message('Caught an warning!')
print(w)
},
finally = {
message('All done, quitting.')
}
)
}

254 ：卵の名無しさん：2021/07/04(日) 06:05:02.42 ID:MvYVZ7Z5.net

Windows Registry Editor Version 5.00

[HKEY_CURRENT_USER\Control Panel\Desktop\WindowMetrics]

255 ：卵の名無しさん：2021/07/06(火) 19:26:40.83 ID:7oMfliWM.net

古い低スペックPCに、Insider Programの正規版でWindows 11を入れてみました。

↓ Windows Insider Program公式サイト
https://insider.windows.com/ja-jp/

＜レジストリの変更＞
HKEY_LOCAL_MACHINE\SOFTWARE\Microsoft\WindowsSelfHost\UI\Selection
UIBranch ← Dev
UIContentType ←Mainline
UIRing ← External
HKEY_LOCAL_MACHINE\SOFTWARE\Microsoft\WindowsSelfHost\Applicability
BranchName ← Dev
ContentType ←Mainline
Ring ← External
HKEY_LOCAL_MACHINE\SYSTEM\Setup
新しいキー LabConfig を追加
LabConfigに以下のDWORD(32ビット値)を追加
BypassTPMCheck ← 1
BypassSecureBootCheck ← 1

Windows Updateでエラー0cx190012eが発生した場合、
C:\Users\Default\AppData\Local\Microsoft\Windows\WSUS
からsetupconfig.iniを削除します。

AppraiserRes.dllのダウンロードサイト
https://github.com/CodeProf14/Fix-TPM
C:\$WINDOWS.~BT\SourcesのDLLと置き換える

256 ：卵の名無しさん：2021/08/08(日) 22:23:50.74 ID:DBdOU9fb.net

'%IN%' <- function(v,M) any(rowSums(M == v[col(M)]) == ncol(M))
c(4,5,6) %IN% matrix(c(1,2,3,
4,5,6,
7,8,9),
nrow=3,ncol=3,byrow=TRUE)

257 ：卵の名無しさん：2021/08/26(木) 18:32:55.73 ID:qkTk8DPf.net

>>1
【全スレに注意連絡】

医療空間で連携して悪質な行為をする創.価.学.会.員.に注意してください。

【創.価.学.会.員.が連携して特定の患者をいじめたり、嫌がらせ、ガスライティング、不審な診察や治療、病院に受診できなくさせたり、たらい回し】

そのようなことが絶対にないように、みんなで不審な学.会.員の医療従事者の行動を見守りましょう。

そしておかしなことがあれば、そのまま放置せず、各所に報告や通報、ネットなどにも声をあげていきましょう。

隠蔽されやすい連携した嫌がらせは、そのように表沙汰にすることで弱体化していきます。
みんなで声をあげて悪質な連携を阻止していきましょう。

関連スレッド↓
医療法人アイクレセント　いなだ小児科・アレルギー科 【佐賀】２
https://itest.5ch.net/egg/test/read.cgi/hosp/1624113085
集団ストーカー自称おおせき
https://lavender.5ch.net/test/read.cgi/minor/1622387124/l50
佐賀大学医学部附属病院
https://egg.5ch.net/test/read.cgi/hosp/1622027754/l50

情報交換などはこのスレでもできます↓
◆◆◆◆◆◆◆ 医療空間での創.価.学.会　◆◆◆◆◆◆◆
https://egg.5ch.net/test/read.cgi/hosp/1628326388/l50

情報のリークはこのスレにも出来ます↓
誰にも言えない事を吐くスレ
https://egg.5ch.net/test/read.cgi/hosp/1627796592/l50 hdkdl

258 ：卵の名無しさん：2021/08/29(日) 23:14:57.19 ID:RCoITSyi.net

[R] retrieving object names passed indirectly to a function

https://stat.ethz.ch/pipermail/r-help/2010-September/252840.html

foo <- function(p1, p2) {
p1name <- deparse(substitute(p1))
p2name <- deparse(substitute(p2))

# The real function runs a simulation, but
# here we just confirm the object names
# and values passed as args
cat("passed", p1name, "with value", p1, "\n")
cat("passed", p2name, "with value", p2, "\n")
}

259 ：卵の名無しさん：2021/10/30(土) 16:47:02.51 ID:k7cvh2/f.net

"
レミフェンタニルの添付文書には，BMIが25を超える肥満患者には，標準体重を用いた投与量の調整を行うことが推奨されている．
高齢者:血圧低下などの副作用があらわれやすいため，開始用量の半減など，
患者の全身状態を観察しながら投与量を決定
"
# 2mg/20mL lwr=0.25μg/kg/min upr=0.05/kg/min ideal BMI=22
Ultiva<-function(cm,kg,age,lwr=0.25,upr=0.50,conc=2000/20,bmi=22){
BMI=kg/(cm/100)^2
LBM=ifelse(BMI>25,bmi*(cm/100)^2,kg) # lean body mass
L=lwr*LBM*60/conc
U=upr*LBM*60/conc
Lower=ifelse(age<75,L,L/2) # mL/h
Upper=ifelse(age<75,U,U/2)
round(c(BMI=BMI,LBM=LBM,'L(mL/h)'=Lower,'U(mL/h)'=Upper),2)
}
Ultiva(167.5,77.5,33)
Ultiva(137.2,50.8,87)

"

260 ：卵の名無しさん：2021/10/30(土) 16:47:07.75 ID:k7cvh2/f.net

持続投与量の目安は、セボフルランでは3-4 μg/kg/min、プロポフォールでは7.5 μg/kg/minと言われる。
添付文書:持続注入により投与する場合は、7μg/kg/分の投与速度で持続注入を開始する。
"

Eslax<-function(kg,sevo=FALSE,lwr=0.6,upr=0.9,conc=50e3/5){ # 50mg/5mL
precurarization=kg*0.03*1e3/conc
l=lwr*1e3*kg/conc # mL bolus
u=upr*1e3*kg/conc
L=ifelse(sevo,3*kg*60/conc,7*kg*60/conc) # on label
U=ifelse(sevo,4*kg*60/conc,7.5*kg*60/conc) # off label
cat('precararization(mL) =',precurarization,'\n')
cat('bolus(mL) =',l,'-',u,'\n')
cat('continuous(mL/h) =',L,'-',U,'\n')
}
Eslax(77.5,sevo=FALSE)
Eslax(50.8,sevo=TRUE)

261 ：卵の名無しさん：2021/10/31(日) 21:39:49.32 ID:asgKs/eN.net

# Practically loading kg(mL/h) for 10 min and maintenance kg/10(ml/h)
Precedex<-function(kg,conc=200/50,lwr=0.2,upr=0.7,load=6){
I=kg*load/conc # 初期(10分)負荷(添付文書)
L=kg*lwr/conc　# 維持下限(添付文書)
U=kg*upr/conc # 維持上限(添付文書)
i=kg*4/conc # 推奨初期負荷
m=kg*0.4/conc # 推奨維持
d=kg*0.1/conc # 譫妄予防
cat('sedation(mL/h) =',round(L,2),'-',round(U,2),'\n')
cat('delirium prevention(mL/h) =',round(d,2),'\n\n')
invisible(list(初期負荷=I,維持下限=L,維持上限=U,推奨初期負荷=i,推奨維持=m,譫妄予防=d))
}
Precedex(50)

> Precedex(50)
sedation(mL/h) = 2.5 - 8.75
delirium prevention(mL/h) = 1.25

262 ：卵の名無しさん：2021/10/31(日) 21:40:12.06 ID:asgKs/eN.net

# pdfからcdfの逆関数を作ってhdiを表示させて逆関数を返す
# 両端での演算を回避 ∫[0,1]は∫[1/nxxx,1-1/nxx]で計算
pdf2hdi <- function(pdf,xMIN=0,xMAX=1,cred=0.95,Print=TRUE,nxx=1001){
xx=seq(xMIN,xMAX,length=nxx)
xx=xx[-nxx]
xx=xx[-1]
xmin=xx[1]
xmax=xx[nxx-2]
AUC=integrate(pdf,xmin,xmax)$value
PDF=function(x)pdf(x)/AUC
cdf <- function(x) integrate(PDF,xmin,x)$value
ICDF <- function(x) uniroot(function(y) cdf(y)-x,c(xmin,xmax))$root
hdi=HDInterval::hdi(ICDF,credMass=cred)
print(c(hdi[1],hdi[2]),digits=5)
if(Print){
par(mfrow=c(3,1))
plot(xx,sapply(xx,PDF),main='pdf',type='h',xlab='x',ylab='Density',col='lightgreen',bty='l')
legend('top',bty='n',legend=paste('HDI:',round(hdi,3)))
plot(xx,sapply(xx,cdf),main='cdf',type='h',xlab='x',ylab='Probability',col='lightblue',bty='l')
pp=seq(0,1,length=nxx)
pp=pp[-nxx]
pp=pp[-1]
plot(pp,sapply(pp,ICDF),type='l',xlab='p',ylab='x',main='ICDF',bty='l')
par(mfrow=c(1,1))
}
invisible(ICDF)
}

263 ：卵の名無しさん：2021/10/31(日) 21:41:36.08 ID:asgKs/eN.net

# ベータ分布の差の分布

beta.diff<-function(r1,r2,n1,n2,a=1,b=1){
f <- function(x,y) dbeta(x+y,a+r1,b+n1-r1)*dbeta(y,a+r2,b+n2-r2)
vf=Vectorize(f,vectorize.args = 'y')
pdf <- function(x) integrate(function(y) vf(x,y),-1,1)$value
pdf=Vectorize(pdf)
invisible(pdf)
}
pdf=beta.diff(r1,r2,n1,n2)

264 ：卵の名無しさん：2021/11/03(水) 13:49:42.44 ID:WyXoRruB.net

MAC2Sevo=function(MAC,age,N2O=0,MAC40sevo=1.8){
GOS=function(
age=85, # 年齢
FE_sevo=0.5, # 呼気セボフルラン濃度(%)
FE_N20=66.6){ # 呼気笑気濃度(%)

MAC40_sevo=MAC40sevo # 40歳セボフルランMAC(1.71,1.8,2.1と諸説あり)
MAC40_N2O=104 # 笑気MAC
MACage_sevo=MAC40_sevo*10^(-0.00269*(age-40))
MACage_N2O=MAC40_N2O*10^(-0.00269*(age-40))
MACage_total = FE_sevo/MACage_sevo + FE_N20/MACage_N2O
return(MACage_total) # 混合ガスのMAC
}
f=function(x) GOS(age,x,N2O) - MAC
f=Vectorize(f)
c('FEsevo(%)'=round(uniroot(f,c(0,5))$root,2))
}
MAC2Sevo(1.5,57)

265 ：卵の名無しさん：2021/11/04(木) 06:15:33.44 ID:6I+CHuG0.net

# 年齢を設定してMACを笑気とセボフルレン濃度でMACをグラフ化
Age=87
GOS=function(
age=85, # 年齢
FE_sevo=0.5, # 呼気セボフルラン濃度(%)
FE_N20=66.6, # 呼気笑気濃度(%)
MAC40sevo=1.8){ # 40歳セボフルランMAC(1.71,1.8,2.1と諸説あり)
MAC40_sevo=MAC40sevo
MAC40_N2O=104 # 笑気MAC
MACage_sevo=MAC40_sevo*10^(-0.00269*(age-40))
MACage_N2O=MAC40_N2O*10^(-0.00269*(age-40))
MACage_total = FE_sevo/MACage_sevo + FE_N20/MACage_N2O
return(MACage_total) # 混合ガスのMAC
}
GOS=Vectorize(GOS)
sev=seq(0,3,by=0.1)
N2O=seq(0,70,by=1)
MAC=outer(sev,N2O,function(x,y)GOS(Age,x,y))
contour(sev,N2O,MAC,nlevels=50,
bty='l',xlab='Sevoflurane',ylab='nitrous oxide',main='MAC for age 87')
contour(sev,N2O,MAC,level=0.75,lwd=2,add=T)

266 ：卵の名無しさん：2022/04/06(水) 10:23:07.52 ID:HkRaCX9n.net

f=\(n){
dat=c(0,rep(100,n-1))
m=mean(dat)
s=sd(dat)
50+(0-m)/s
}

nn=1000:5000
plot(nn,sapply(nn,f))

uniroot(f,c(1000,10000))
f(2501)
f(2502)

267 ：卵の名無しさん：2022/04/07(木) 06:36:23.98 ID:Ou/qCTfb.net

calc = \(L,m,n){
f1 = \(x){
y = as.character(x)
strsplit(y,'') |> unlist() |> as.numeric()
}

f2 = \(x){
y = f1(x)
x%%n==0 & all(y==sort(y)) & sum(y)==m
}
x=((10^L-1)/9):(10^L-1)
x[sapply(x,f2)]
}

calc(3,16,2)
calc(5,24,3)

268 ：卵の名無しさん：2022/04/09(土) 06:30:01.59 ID:mZjVxAsM.net

FF 83 F8 03 72

FF 83 F8 03 EB

269 ：卵の名無しさん：2022/04/09(土) 16:45:12.95 ID:wXi20OuD.net

pH=7.32
pCO2=23
HCO3=16

pH_l=7.36
pH_u=7.44
pCO2_l=36
pCO2_u=44
HCO3_l=21
HCO3_u=27

if(pH<pH_l & pCO2<pCO2_l){
primary="metabolic acidosis"
expected_pCO2=40 + 1.2*(HCO3-24)
if(pCO2<expected_pCO2) additional= "secondary respiratory alkalosis")
if(pCO2>expected_pCO2) additional= "secondary respiratory acidosis")
}

if(pH>pH_u & pCO2>pCO2_u){
primary="metabolic alkalosis"
expected_pCO2=40 + 0.7*(HCO3-24)
if(pCO2<expected_pCO2) additional= "secondary respiratory alkalosis")
if(pCO2>expected_pCO2) additional= "secondary respiratory acidosis")
}

270 ：卵の名無しさん：2022/04/09(土) 16:45:22.36 ID:wXi20OuD.net

if(pH<pH_l & pCO2>pCO2_u){
primary="respiratory aciosis"
status=NULL
if(HCO3_l<HCO3 & HCO3<HCO3_u) status="acute" else status="chronic"
if(status=="acute") expected_HCO3=24+0.1*(pCO2-40)
if(status=="chronic") expected_HCO3=24+0.4*(pCO2-40)
if(HCO3<expectedHCO3) additional="secondary metabolic acidosis")
if(HCO3>expectedHCO3) additional="secondary metabolic alkalosis")
}

if(pH>pH_u & pCO2<pCO2_l){
primary="respiratory alkalosis"
status=NULL
if(HCO3_l<HCO3 & HCO3<HCO3_u) status="acute" else status="chronic"
if(status=="acute") expected_HCO3=24+0.2*(40-pCO2)
if(status=="chronic") expected_HCO3=24+0.4*(40-pCO2)
if(HCO3<expectedHCO3) additional="secondary metabolic acidosis")
if(HCO3>expectedHCO3) additional="secondary metabolic alkalosis")
cat(paste(status,primary),'\n')
cat(additional)

271 ：卵の名無しさん：2022/04/09(土) 16:45:58.04 ID:wXi20OuD.net

if(pH<pH_l & pCO2>pCO2_u){
primary="respiratory aciosis"
status=NULL
if(HCO3_l<HCO3 & HCO3<HCO3_u) status="acute" else status="chronic"
if(status=="acute") expected_HCO3=24+0.1*(pCO2-40)
if(status=="chronic") expected_HCO3=24+0.4*(pCO2-40)
if(HCO3<expectedHCO3) additional="secondary metabolic acidosis")
if(HCO3>expectedHCO3) additional="secondary metabolic alkalosis")
}

if(pH>pH_u & pCO2<pCO2_l){
primary="respiratory alkalosis"
status=NULL
if(HCO3_l<HCO3 & HCO3<HCO3_u) status="acute" else status="chronic"
if(status=="acute") expected_HCO3=24+0.2*(40-pCO2)
if(status=="chronic") expected_HCO3=24+0.4*(40-pCO2)
if(HCO3<expectedHCO3) additional="secondary metabolic acidosis")
if(HCO3>expectedHCO3) additional="secondary metabolic alkalosis")
cat(paste(status,primary),'\n')
cat(additional)

272 ：卵の名無しさん：2022/04/09(土) 17:37:32.78 ID:wXi20OuD.net

Astrap <- \(pH,pCO2,HCO3){
pH_l=7.36
pH_u=7.44
pCO2_l=36
pCO2_u=44
HCO3_l=22
HCO3_u=26
if(pH_l<pH & pH<pH_u){
if(pCO2>pCO2_u){
primary="mixed respiratory acidosis"
additional="and metabolic alkalosis"
}
if(pCO2<pCO2_l){
primary="mixed respiratory alkalosis"
additional="and metabolic acidosis"
}
}

if(pCO2_l<pCO2 & pCO2<pCO2_u){
if(pH>pH_u){
primary="mixed metaboic alkalosis"
additional="and respiratory acidosis"
}
if(pH<pH_l){
primary="mixed metabolic acidosis"
additional="and respiratory alkalosis"
}
}

273 ：卵の名無しさん：2022/04/09(土) 17:38:10.02 ID:wXi20OuD.net

if(pH<pH_l & pCO2<pCO2_l){
primary="primary metabolic acidosis"
expected_pCO2=40 + 1.2*(HCO3-24)
if(pCO2<expected_pCO2) additional= "secondary respiratory alkalosis"
if(pCO2>expected_pCO2) additional= "secondary respiratory acidosis"
}
if(pH>pH_u & pCO2>pCO2_u){
primary="primary metabolic alkalosis"
expected_pCO2=40 + 0.7*(HCO3-24)
if(pCO2<expected_pCO2) additional= "secondary respiratory alkalosis"
if(pCO2>expected_pCO2) additional= "secondary respiratory acidosis"
}
status=NULL
if(pH<pH_l & pCO2>pCO2_u){
primary="primary respiratory aciosis"
if(HCO3_l<HCO3 & HCO3<HCO3_u) status="acute" else status="chronic"
if(status=="acute") expected_HCO3=24+0.1*(pCO2-40)
if(status=="chronic") expected_HCO3=24+0.4*(pCO2-40)
if(HCO3<expected_HCO3) additional="secondary metabolic acidosis"
if(HCO3>expected_HCO3) additional="secondary metabolic alkalosis"
}
if(pH>pH_u & pCO2<pCO2_l){
primary="primary respiratory alkalosis"
if(HCO3_l<HCO3 & HCO3<HCO3_u) status="acute" else status="chronic"
if(status=="acute") expected_HCO3=24+0.2*(40-pCO2)
if(status=="chronic") expected_HCO3=24+0.4*(40-pCO2)
if(HCO3<expected_HCO3) additional="secondary metabolic acidosis"
if(HCO3>expected_HCO3) additional="secondary metabolic alkalosis"
}
cat(status,primary,'\n')
cat(additional)
}

274 ：卵の名無しさん：2022/04/09(土) 17:38:59.63 ID:wXi20OuD.net

> Astrap(pH=7.23,pCO2=23,HCO3=16)
primary metabolic acidosis
secondary respiratory alkalosis
> Astrap(pH=7.54,pCO2=23,HCO3=17)
chronic primary respiratory alkalosis
secondary metabolic acidosis
> Astrap(pH=7.40, pCO2=23, HCO3=20)
mixed respiratory alkalosis
and metabolic acidosis
> Astrap(pH=7.40, pCO2=23, HCO3=35)
mixed respiratory alkalosis
and metabolic acidosis
> Astrap(pH=7.40, pCO2=23, HCO3=17)
mixed respiratory alkalosis
and metabolic acidosis

275 ：卵の名無しさん：2022/04/26(火) 16:51:59.52 ID:aUc6q9el.net

>>1
【全スレ注意連絡】※告発内容はｺﾋﾟｰ自由ｶｸｻﾝ希望


【裏で集団ストーカー犯罪を指示してる医師、稲田成安に気をつけて】


http://www.inada-kids.com/wordpress/wp-content/uploads/greeting.jpg

佐賀県のいなだ小児科・アレルギー科 の小児科医、稲田成安（いなだ　しげやす）に気をつけてください

自分に批判的な医療従事者や患者、その他の人たちを、探偵やソウカガッカイの集団ストーカーを使って調べ上げて、嫌がらせをやらせて仕返し、弾圧してきます
あらゆる嫌がらせで自殺や自滅に追い込む工作を仕掛けてきます

（※ ５ちゃんねるでは「ソウカガッカイ」という言葉を漢字四字で投稿するとアクセス禁止になる事があるのでカタカナ表記にしてあります）

もし稲田による集団ストーカー犯罪にあったら、証拠を取り、周囲や各機関に声を上げていきましょう。

そうする事で集団ストーカー組織は弱体化していきます

詳細は、↓のスレ
【ソウカ在日】稲田成安（いなだ　しげやす）
https://egg.5ch.net/test/read.cgi/hosp/1644311065/

いなだ小児科・アレルギー科　【佐賀】
https://egg.5ch.net/test/read.cgi/hosp/1622365146/-100

黒猫世（くろねこよ）
https://rio2016.5ch.net/test/read.cgi/police/1648800949/l50 gyhjoooo


「　　　いなだ小児科・アレルギー科　　　集団ストーカー犯罪　　」で検索

276 ：卵の名無しさん：2022/06/12(日) 19:47:04.91 ID:48n0Max4.net

rm(list=ls())
library(RcppAlgos)
pm=permuteGeneral(v=1:3,freqs=c(3,3,4))
n=nrow(pm)

f=\(x){
t1=(1:10)[x==1]
t2=(1:10)[x==2]
t3=(1:10)[x==3]
if(t1[1]>t2[1]){
temp=t1
t1=t2
t2=temp
}
list(t1,t2,t3)
c(paste(LETTERS[t1],collapse=''),paste(LETTERS[t2],collapse=''),
paste(LETTERS[t3],collapse=''))
}
ans=unique(t(apply(pm,1,f)))
print(head(ans,15),quote=FALSE)
print(tail(ans,15),quote=FALSE)

277 ：卵の名無しさん：2022/06/12(日) 20:19:02.97 ID:48n0Max4.net

"
10人の子供(A,B,C,....,J)を３チームに割り当てる。
各チームは１人以上の子供が割り当てられる。
何通りの分け方があるか？

例.
(A) (J) (B,C,D,E,F,G,H,I)
(A,C) (H,I,J) (B,D,E,F,G)

"
rm(list=ls())
gr=expand.grid(1:8,1:8)
gr1=gr[rowSums(gr)<10,]
f=\(x) 10-sum(x)
gr2=cbind(gr1,apply(gr1,1,f))
gr3=unique(t(apply(gr2,1,sort)))

gr3
ans=apply(gr3,1,\(x) choose(10,x[1])*choose(10-x[1],x[2])/ length(unique(x)) )
sum(ans)
ans

278 ：卵の名無しさん：2022/06/12(日) 20:25:08.35 ID:48n0Max4.net

"10人の子供(A,B,C,....,J)を３人、３人、４人のチームに割り当てる。
何通りの分け方があるか、列挙せよ"

rm(list=ls())
library(RcppAlgos)
pm=permuteGeneral(v=1:3,freqs=c(3,3,4))

f=\(x){
t1=(1:10)[x==1]
t2=(1:10)[x==2]
t3=(1:10)[x==3]
if(t1[1]>t2[1]){
temp=t1
t1=t2
t2=temp
}
list(t1,t2,t3)
c(paste(LETTERS[t1],collapse=''),paste(LETTERS[t2],collapse=''),
paste(LETTERS[t3],collapse=''))
}
ans=unique(t(apply(pm,1,f)))
print(head(ans,20),quote=FALSE)
print(tail(ans,20),quote=FALSE)

279 ：卵の名無しさん：2022/06/12(日) 20:25:12.57 ID:48n0Max4.net

"
10人の子供(A,B,C,....,J)を３チームに割り当てる。
各チームは１人以上の子供が割り当てられる。
何通りの分け方があるか？

例.
(A) (J) (B,C,D,E,F,G,H,I)
(A,C) (H,I,J) (B,D,E,F,G)

"
rm(list=ls())
gr=expand.grid(1:8,1:8)
gr1=gr[rowSums(gr)<10,]
f=\(x) 10-sum(x)
gr2=cbind(gr1,apply(gr1,1,f))
gr3=unique(t(apply(gr2,1,sort)))

ans=apply(gr3,1,\(x) choose(10,x[1])*choose(10-x[1],x[2])/ length(unique(x)) )
sum(ans)
res=cbind(gr3,ans)
colnames(res)=NULL
rownames(res)=NULL
res

280 ：卵の名無しさん：2022/08/15(月) 21:02:59.07 ID:WCPJCNZ4.net

多摩南部地域病院に入院した事あるけど看護師がクソ過ぎ
患者のことはどうでもよく自分達の仕事さえまわせればそれでいいって感じ
底意地の悪いクソ女看護師の清○桃○は要注意
そのうえ入院環境もクソ過ぎ
狭いトイレ、備品等もボロいものばかりで配慮が足りなさすぎる
病院自体のシステムも不親切だし二度と行きたくない
こんな病院は無いほうがいい
開業医もこんな病院紹介しないでほしい

281 ：卵の名無しさん：2022/10/07(金) 07:26:02.86 ID:i93G7UfB.net

ganrikinto<-function(A=NULL,I=NULL,N=NULL,x=NULL){ # A:Principal　I:Annual Interest N:months of repayment x:monthly repayment
if(is.null(I)){
I=uniroot(\(m) x-A*m/(1-(1+m)^(-N)),c(1e-12,1-1e-12))$root*12
return(c('Annual Interest'=I))
}
m=I/12
if(is.null(x)){
x=A*m/(1-(1+m)^(-N))
return(c('monthly repayment'=round(x,2)))
}
if(is.null(N)){
N=ifelse(x>A*m,log(x/(x-A*m))/log(m+1),Inf)
return(c('repayment period in months'=round(N,2)))
}
if(is.null(A)){
A=x/m*(1-(1+m)^(-N))
return(c('Principal'=round(A)))
}

}

282 ：卵の名無しさん：2024/01/04(木) 05:47:47.73 ID:s7ncVxDB.net

f=\(n){
n1=n%%10
i=1
flg <- (n1*n1^i)%%10==n1
if(flg){ p=i
}else{
while(!flg & i<2024){
i=i+1
flg <- (n1*n1^i)%%10==n1
}
p=i
}
j=ifelse(2024%%p,2024%%p,p)
(n1^p)%%10
}
sapply(1:2024,f) |> table()

283 ：卵の名無しさん：2024/01/04(木) 11:13:27.39 ID:mx4lFG6x.net

2024＾2024
π^2024
の先頭の数字を算出

臨床医の嗜み：R言語での一般解を算出するコード
calc=function(m,n) floor(10^(n*log10(m) -floor(n*log10(m))))
calc(2024,2024)
calc(pi,2024)

284 ：卵の名無しさん：2024/01/05(金) 07:22:34.13 ID:lKTSD30o.net

"
地球の公転周期は秒単位だと365日5時間48分46秒であるという
5時間48分46秒＝20926秒
閏年は400年に97回
1日＝86400秒
"
p=20926/86400 ; p
97/400 - p
8/33 - p
fn=\(x){
xp=x*p
fl=floor(xp)
ce=ceiling(xp)
dfl=abs(fl/x-p)
dce=abs(ce/x-p)
c(min(dfl,dce),ifelse(dfl<dce,fl,ce),x)
}
dat=t(sapply(1:40,fn))
rbind(dat[which.min(dat[,1]),])

dat=t(sapply(1:10000,fn))
rbind(dat[which.min(dat[,1]),])

n=seq(4,400,4)
dat=t(sapply(n,fn))
rbind(dat[which.min(dat[,1]),])

285 ：卵の名無しさん：2024/01/05(金) 19:36:46.56 ID:Nm2/P+r6.net

"
短径20cm　長径24cmの楕円のピザを楕円弧の長さが等しいように３つの楕円扇形に分割する。
楕円扇形の面積が最大になるようにするにはどのように切ればよいか？
そのときの面積はいくつか？（答は整数値でよい）。

あらゆるリソースを用いてよい。


"
a=12
b=10
n=3
S=\(t1,t2) a*b/2*(t2-t1)
L=\(t1,t2) integrate(\(x) sqrt((-a*sin(x))^2+(b*cos(x))^2),t1,t2,abs.tol = 1e-24)$value
Lt=L(-pi,pi) ; Lt # 全周長
t12t2=\(t1) uniroot(\(x) L(t1,x)-Lt/n,c(t1,t1+pi),tol=1e-24)$root
t12S=\(t1) S(t1,t12t2(t1))
t12S=Vectorize(t12S)
curve(t12S,0,pi)
opt=optimize(t12S,c(-pi/2,pi/2),maximum = TRUE,tol=1e-24) ; opt
source('toolmini.R')
Plot(-a,a,zero=FALSE)
Ellip(0i,a,b,lty=3)
t1=opt$maximum ; t2=t12t2(t1)
tt=seq(t1,t2,length=100)
lines(a*cos(tt),b*sin(tt),lwd=2)
T1=a*cos(t1)+1i*b*sin(t1)
T2=a*cos(t2)+1i*b*sin(t2)
seg(0i,T1,lty=2,lwd=2)
seg(0i,T2,lty=2,lwd=2)
Angle(T1,0i,T2)

(pi*a*b - opt$obj)/2
S(t2,pi)

286 ：卵の名無しさん：2024/01/06(土) 16:06:08.99 ID:YqxPaf+q.net

rm(list=ls())
par(bty='l')

calc=function(m=2024,n=2024) floor(10^(n*log10(m) -floor(n*log10(m))))
calc(2024,2024)
calc(pi,2024)
calc(exp(1),2024)

"
nを1以上2024以下の整数とする。
n^2024の先頭に現れる数字はなにが一番多いか？
頻度順に並べよ。
"
calc.top =\(n) sapply(1:n,calc) |> table() |> sort() |> rev()
re=calc.top(2024)
re
p=re/sum(re)
calc.top(1000)
sapply(1:1000,\(x) as.numeric(names(calc.top(x)[1])))

k=1:9
plot(k,p,pch=16)
lines(k,log10(1+1/k),type='h',col=2) # ベンフォードの法則

287 ：卵の名無しさん：2024/01/12(金) 16:58:50.71 ID:NVcGRVeI.net

医者板では全く相手にされないから数学板で高校生相手にイキれると思ってる尿瓶チンパンジジイ滑稽だね

288 ：卵の名無しさん：2024/04/08(月) 05:59:31.45 ID:0nrZk2Or.net

これが現実

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


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

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