「DV_PRED_IPRED」の編集履歴（バックアップ）一覧はこちら

「DV_PRED_IPRED」(2007/02/03 (土) 10:45:27) の最新版変更点

追加された行は緑色になります。

削除された行は赤色になります。

*DV, PRED, IPRED の重ね描き

 plot.ipred <-
 function(
   tab, ID="ID", DV="DV", IPRE="IPRE",
   layout=c(3, 3),
   xlab="Time (hr)", ylab="Concentration (ng/mL)"
 ) {
   tab.plot <- rbind(
     data.frame(type=rep(1, nrow(tab)), ID=tab[,ID], TIME=tab$TIME, Y=tab$DV),
     data.frame(type=rep(2, nrow(tab)), ID=tab[,ID], TIME=tab$TIME, Y=tab$PRED),
     data.frame(type=rep(3, nrow(tab)), ID=tab[,ID], TIME=tab$TIME, Y=tab[,IPRE])
   )
   
   res <- xyplot(Y ~ TIME | factor(ID), tab.plot,
     type=tab.plot$type,
     layout=layout,
     xlab=xlab, ylab=ylab,
     panel=function(x, y, subscripts, type=type) {
       panel.xyplot(x[type[subscripts]==1], y[type[subscripts]==1], col.symbol=1)
       panel.xyplot(x[type[subscripts]==2], y[type[subscripts]==2], type="l", col.line=1)
       panel.xyplot(x[type[subscripts]==2], y[type[subscripts]==3], type="l", col.line=4)
     },
   )
   res
 }
 
 plot.ipred(tab[tab$MDV==0,], ID="SID", ylab="Plasma concentration (ng/mL)")

*DV, PRED, IPRED の重ね描き

 plot.ipred <-
 function(
   tab, ID="ID", DV="DV", IPRE="IPRE",
   layout=c(3, 3),
   xlab="Time (hr)", ylab="Concentration (ng/mL)",
   ...
 ) {
   tab.plot <- rbind(
     data.frame(type=rep(1, nrow(tab)), ID=tab[,ID], TIME=tab$TIME, Y=tab$DV),
     data.frame(type=rep(2, nrow(tab)), ID=tab[,ID], TIME=tab$TIME, Y=tab$PRED),
     data.frame(type=rep(3, nrow(tab)), ID=tab[,ID], TIME=tab$TIME, Y=tab[,IPRE])
   )
   
   res <- xyplot(Y ~ TIME | factor(ID), tab.plot,
     type=tab.plot$type,
     layout=layout,
     xlab=xlab, ylab=ylab,
     panel=function(x, y, subscripts, type=type) {
       panel.xyplot(x[type[subscripts]==1], y[type[subscripts]==1], col.symbol=1)
       panel.xyplot(x[type[subscripts]==2], y[type[subscripts]==2], type="l", col.line=1)
       panel.xyplot(x[type[subscripts]==2], y[type[subscripts]==3], type="l", col.line=4)
     },
     ...
   )
   res
 }
 
 plot.ipred(tab[tab$MDV==0,], ID="SID", ylab="Plasma concentration (ng/mL)")