diff options
Diffstat (limited to 'man')
| -rw-r--r-- | man/mkinresplot.Rd | 16 | ||||
| -rw-r--r-- | man/plot.mkinfit.Rd | 15 |
2 files changed, 22 insertions, 9 deletions
diff --git a/man/mkinresplot.Rd b/man/mkinresplot.Rd index 27e1322f..465b3038 100644 --- a/man/mkinresplot.Rd +++ b/man/mkinresplot.Rd @@ -5,7 +5,8 @@ \title{Function to plot residuals stored in an mkin object} \usage{ mkinresplot(object, obs_vars = names(object$mkinmod$map), xlim = c(0, - 1.1 * max(object$data$time)), xlab = "Time", ylab = "Residual", + 1.1 * max(object$data$time)), standardized = FALSE, xlab = "Time", + ylab = ifelse(standardized, "Standardized residual", "Residual"), maxabs = "auto", legend = TRUE, lpos = "topright", col_obs = "auto", pch_obs = "auto", frame = TRUE, ...) } @@ -18,15 +19,17 @@ the model} \item{xlim}{plot range in x direction.} -\item{xlab}{Label for the x axis. Defaults to "Time [days]".} +\item{standardized}{Should the residuals be standardized by dividing by the +standard deviation given by the error model of the fit?} -\item{ylab}{Label for the y axis. Defaults to "Residual [\% of applied -radioactivity]".} +\item{xlab}{Label for the x axis.} + +\item{ylab}{Label for the y axis.} \item{maxabs}{Maximum absolute value of the residuals. This is used for the scaling of the y axis and defaults to "auto".} -\item{legend}{Should a legend be plotted? Defaults to "TRUE".} +\item{legend}{Should a legend be plotted?} \item{lpos}{Where should the legend be placed? Default is "topright". Will be passed on to \code{\link{legend}}.} @@ -58,7 +61,8 @@ mkinresplot(fit, "m1") } \seealso{ \code{\link{mkinplot}}, for a way to plot the data and the fitted - lines of the mkinfit object. + lines of the mkinfit object, and \code{\link{plot_res}} for a function + combining the plot of the fit and the residual plot. } \author{ Johannes Ranke diff --git a/man/plot.mkinfit.Rd b/man/plot.mkinfit.Rd index 3834eaf5..6824f8a5 100644 --- a/man/plot.mkinfit.Rd +++ b/man/plot.mkinfit.Rd @@ -16,9 +16,13 @@ lpos = "topright", inset = c(0.05, 0.05), show_errmin = FALSE, errmin_digits = 3, frame = TRUE, ...) -plot_sep(fit, show_errmin = TRUE, ...) +plot_sep(fit, show_errmin = TRUE, + show_residuals = ifelse(identical(fit$err_mod, "const"), TRUE, + "standardized"), ...) -plot_res(fit, sep_obs = FALSE, show_errmin = sep_obs, ...) +plot_res(fit, sep_obs = FALSE, show_errmin = sep_obs, + show_residuals = ifelse(identical(fit$err_mod, "const"), TRUE, + "standardized"), ...) plot_err(fit, sep_obs = FALSE, show_errmin = sep_obs, ...) } @@ -54,7 +58,9 @@ corresponding model prediction lines.} \item{show_residuals}{Should residuals be shown? If only one plot of the fits is shown, the residual plot is in the lower third of the plot. Otherwise, i.e. if "sep_obs" is given, the residual plots will be located -to the right of the plots of the fitted curves.} +to the right of the plots of the fitted curves. If this is set to +'standardized', a plot of the residuals divided by the standard deviation + given by the fitted error model will be shown.} \item{show_errplot}{Should squared residuals and the error model be shown? If only one plot of the fits is shown, this plot is in the lower third of @@ -89,6 +95,8 @@ chi2 error percentage.} \item{frame}{Should a frame be drawn around the plots?} \item{\dots}{Further arguments passed to \code{\link{plot}}.} + +\item{standardized}{For} } \value{ The function is called for its side effect. @@ -113,6 +121,7 @@ SFO_SFO <- mkinmod(parent = mkinsub("SFO", "m1", full = "Parent"), fit <- mkinfit(SFO_SFO, FOCUS_2006_D, quiet = TRUE, error_model = "tc") plot(fit) plot_res(fit) +plot_res(fit, standardized = FALSE) plot_err(fit) # Show the observed variables separately, with residuals |