R语言笔记：群体药动学数据探索

msxiaoming

于 2020-05-27 00:12:22 发布

阅读量1.4k

点赞数

分类专栏：群体药代动力学

本文链接：https://blog.csdn.net/m0_47950413/article/details/106368560

版权

本文介绍了使用R语言进行群体药动学数据的探索，包括加载必要的包，设置工作目录，读取和过滤数据，以及对数据的单位和分类因子进行处理，进一步计算均值和方差等统计量。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

用到的包

library(dplyr)
library(ggplot2)
library(gridExtra)

加载资源

source("D:\\工作文档\\R\\col_function.r")

设置工作目录

setwd("C:\\Users\\Administrator.ADMIN-20170618T\\Documents\\demo-data\\PPK\\iv infusion 3cp")

读取数据

#read data----
SDTAB<-read.table("SDTAB1",sep="",skip=1,header = TRUE) 
PKdata<-read.csv("MonkeyPK.csv",header=T,skip=2)

过滤数据

#filter data----
PKdata<-filter(MKPK,MDV==0)
PKdata<-filter(MKPK,C!="C")
INDFIT_PK5<-filter(INDFIT_PK,TIME>829)

单位和分类因子

#add unit----
MKPK$DOSE[MKPK$DOSE==0.5]<-"0.5 mg"
#data group----
MKPK$DOSE_f<-factor(MKPK$DOSE,levels=c("0.5 mg","1.5 mg"))

求均值方差等

#Conc Mean Profile-----
ConcvsTIME<-summarise(group_by(MKPKSAD,TIME,DOSE_f),
                      mean=mean(CONC),
                      sd=sd(CONC),
                      n=sum(!is.na(CONC)),
                      se=sd/sqrt(n))

#---
  ylim(0,500)+
  xlim(0,1200)+
  scale_y_log10()+
  facet_wrap(~DOSE_f,ncol=3)+
  facet_wrap(~ID,ncol=7)+
  guides(color=FALSE)+

  legend.text = el

最低0.47元/天解锁文章