23/03/13 16:59:51.10 G3xv/LvTH.net
# プロタノールL注0.2mgを全量20mLに希釈して体重50kgで0.02～0.20μg/kg/minで投与するときの速度とグラフを出力
gamma2hourlyVol=function(mg,mL,kg,gamma){ # mcg/kg/min -> mL/h
# mg/1000mcg * kg/body * 60min/h * mL/mg * gammma
coef=60/1000*mL/mg*kg
velocity=coef*gamma
mat=cbind(gamma,velocity)
colnames(mat)=c('mcg/kg/min','mL/h')
L=ifelse(length(gamma)==1,0,min(gamma))
U=ifelse(length(gamma)==1,5*gamma,max(gamma))
main=paste('concenteration =',mg,'mg /',mL,'mL\nweight =',kg,'kg')
curve(coef*x,L,U,bty='l',xlab='mcg/kg/min',ylab='mL/h',main=main)
n=21
X=seq(L,U,length=n)
Y=coef*X
for(i in 1:n){
for(j in 1:n){
segments(X[i],Y[1],X[i],Y[n],col=c(8,2)[i%%2+1],lty=3)
segments(X[1],Y[j],X[n],Y[j],col=c(8,2)[j%%2+1],lty=3)
}
}
points(mat,pch=16)
mat
}
# isoproterenol 0.2mg/20mL at 0.02-0.20mcg/kg/min
gamma2hourlyVol(mg=0.2,mL=20,kg=50,gamma=c(0.02,0.20))