책에 있는 two_comp.csv로 one-compartment model, two-compart model 둘을 fit한후 그림을 그려봤습니다.
> two.comp=read.csv("two_comp.csv")
> two=nls(conc~SSbiexp(time,A1,lrc1,A2,lrc2),data=two.comp)
> one=nls(conc~C0*exp(-K*time),start=list(C0=0.8,K=0.42),data=two.comp)
> time=rep(two.comp$time,2)
> conc=c(fitted(one),fitted(two))
> model=rep(c('one','two'),c(11,11))
# 그래프
> library(lattice)
> xyplot(conc~time,groups=model,type="b",col=c("red","blue"))
다음검색