library(gdata)
rawData <- read.xls("平喘藥實(shí)驗(yàn).xlsx",
encoding="GBK", header=F)[c(-1, -2, -3), c(-1)]
dataDF <- data.frame(value=as.numeric(t(rawData)),
time=rep(1:ncol(rawData)-1, times=nrow(rawData)),
group=factor(
rep(c("H", "Iso", "BH", "ACh", "Iso2"),
each=ncol(rawData)), levels=
c("H", "Iso", "BH", "ACh", "Iso2")))
labeli <- function(variable, value){
names_li <- list("H"="組胺", "Iso"="異丙腎上腺素", "BH"="苯海拉明+組胺", "ACh"="乙酰膽堿",
"Iso"="異丙腎上腺素")
return(names_li[value])
}
library(ggplot2)
p <- ggplot(dataDF, aes(x=time, value))+
geom_line(lwd=1.2)+theme_bw() +
facet_wrap(~group, labeller=labeli)+
theme(strip.text.x = element_text(size = 10, face="bold"))+
theme(axis.text.x=element_text(face="bold", size=10),
axis.text.y=element_text(face="bold", size=10),
axis.title.y=element_text(size=10, face="bold"))+
xlab("時(shí)間(min)")+ylab("灌注液流速(ml/min)")
pdf("result.pdf", width=5, height=4)
p
dev.off()
?著作權(quán)歸作者所有,轉(zhuǎn)載或內(nèi)容合作請聯(lián)系作者
【社區(qū)內(nèi)容提示】社區(qū)部分內(nèi)容疑似由AI輔助生成，瀏覽時(shí)請結(jié)合常識與多方信息審慎甄別。
平臺(tái)聲明：文章內(nèi)容（如有圖片或視頻亦包括在內(nèi)）由作者上傳并發(fā)布，文章內(nèi)容僅代表作者本人觀點(diǎn)，簡書系信息發(fā)布平臺(tái)，僅提供信息存儲(chǔ)服務(wù)。
