library(gdata)
rawData <- read.xls("平喘药实验.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("时间(min)")+ylab("灌注液流速(ml/min)")
pdf("result.pdf", width=5, height=4)
p
dev.off()
