R/precision_plot.R In MethodCompare: Bias and Precision Plots

Documented in precision_plot

```#' Precision plot used to compare the standard deviation of a new measurement
#' method with that of a reference standard with possibly heteroscedastic errors
#'
#' This plot allows the visual comparison of the precision (i.e. standard deviation)
#' of the new measurement method with that of a reference standard by creating
#' a scatter plot of the estimated standard deviation against the best linear
#' prediction (BLUP) of the latent variable x.
#'
#'
#' @param object an object retunred by a call to \link{measure_compare}
#' @author  Mingkai Peng
#' @importFrom stats aggregate confint fitted na.exclude na.omit quantile residuals coef
#' @export
#' @examples
#' data(data1)
#' ### analysis
#' measure_model <- measure_compare(data1)
#' ### Precision plot
#' precision_plot(measure_model)
#'

precision_plot <- function(object){
data_old <- object\$Ref
data_new <- object\$New
min_y <- min(min(data_old\$sig_resid_y2,na.rm=TRUE),min(data_new\$sig_resid_y1_corrected,na.rm=TRUE))
max_y <- max(max(data_old\$sig_resid_y2,na.rm=TRUE),max(data_new\$sig_resid_y1_corrected,na.rm=TRUE))
min_y <- floor(min_y)
range = max_y - min_y
max_y <- ceiling(max_y+range*0.2)
par(mar=c(3.5,3.5,2,2)+0.1)
plot(data_old\$y2_hat,data_old\$sig_resid_y2,xlab="",
ylab="",axes = F,
cex=0.8,ylim=c(min_y,max_y))
title(main="Precision plot",cex.main=0.9)
axis(2,col="black",las=1)
mtext("Standard deviation of measurement errors",side = 2,line=2,cex=0.9)
box(col="black")
axis(1)
mtext("BLUP of x",side=1,col="black",line=2,cex=0.9)
points(data_new\$y2_hat,data_new\$sig_resid_y1_corrected,
cex=0.8,pch=19,col="blue")
legend("topleft",legend=c("Reference method (y2)","New method(corrected y1)"),
pch=c(1,19),col=c("black","blue"),xpd=T,horiz = F,
box.lwd=0,cex=0.7)
}
```

Try the MethodCompare package in your browser

Any scripts or data that you put into this service are public.

MethodCompare documentation built on Jan. 5, 2022, 9:06 a.m.