diff options
| author | jranke <jranke@edb9625f-4e0d-4859-8d74-9fd3b1da38cb> | 2010-03-30 19:49:44 +0000 |
|---|---|---|
| committer | jranke <jranke@edb9625f-4e0d-4859-8d74-9fd3b1da38cb> | 2010-03-30 19:49:44 +0000 |
| commit | 12934520d7ec3218ce1505787b6066334a24a562 (patch) | |
| tree | e3d3439a63260dfe8fc01037ef297e09e876cc67 /R/kinresplot.R | |
Initial import of the kinfit package developed from 2008-07 to 2010-03 at
Harlan Laboratories Ltd (former RCC Ltd).
Supports fitting of parent data with the usual FOCUS kinetic models.
git-svn-id: svn+ssh://svn.r-forge.r-project.org/svnroot/kinfit/pkg/kinfit@2 edb9625f-4e0d-4859-8d74-9fd3b1da38cb
Diffstat (limited to 'R/kinresplot.R')
| -rw-r--r-- | R/kinresplot.R | 14 |
1 files changed, 14 insertions, 0 deletions
diff --git a/R/kinresplot.R b/R/kinresplot.R new file mode 100644 index 0000000..be0a85d --- /dev/null +++ b/R/kinresplot.R @@ -0,0 +1,14 @@ +kinresplot <- function(kinobject, kinmodel,
+ xlab = "Time [days]", ylab = "Residual [% of applied radioactivity]",
+ maxabs = "auto")
+{
+ m <- kinobject$fits[[kinmodel]]
+ t <- m$model$t
+ residuals <- residuals(m)
+ if (maxabs == "auto") maxabs = max(abs(residuals))
+ plot(t, residuals,
+ xlab = xlab,
+ ylab = ylab,
+ ylim = c( -1.2 * maxabs, 1.2 * maxabs))
+ title(paste("Residuals of", kinmodel, "fit"), font.main = 1)
+}
|