Deming_Gillard_Fuller_over_groups: Estimate prediction band using Deming regression formulated...

In pernille267/commutability.selectivity: Commutability and non-selectivity evaluation tools

Estimate prediction band using Deming regression formulated by J. Gillard and Fuller over groups

Description

Estimate prediction band using Deming regression formulated by J. Gillard and Fuller over groups

Usage

Deming_Gillard_Fuller_over_groups(
  data,
  groups = "Comparison",
  level = 0.99,
  R = 3,
  Np = 1000,
  evaluated_materials = NULL,
  column_order = c("Comparison", "SampleID", "MP_B", "MP_A", "fit", "lwr", "upr"),
  add_judgement = FALSE,
  judgement_labels = c("no", "yes")
)

Arguments

data	A data frame or data table with format LFDT enclosed with all replicated measurements
groups	The names of the grouping columns of the data
level	A numeric value that captures the overall confidence level of the estimated prediction band
R	An integer signifying the maximum number of replicated measurements performed on each evaluated material
Np	An integer, which captures the number of pointwise prediction intervals making the prediction band across the concentration range. Not relevant if evaluated_materials are !NULL
evaluated_materials	A data frame or data table with format LFDT enclosed with all replicated measurements for evaluated materials such as EQAMs or CRMs. Should be NULL if the PB, that is, pointwise prediction intervals to be estimated
column_order	A vector specifying the order of the columns of the outputted grouped data table. Will be ignored if evaluated_materials is NULL. Note that the length of column_order must be the same as the number of columns of output columns, which is 6 + the number of grouping columns
add_judgement	Should we include a column to the output column stating whether each evaluated material at in each group is inside the estimated prediction interval? Default is FALSE
jugdement_labels	A character vector including the labels for when EQAM is outside estimated PI (FALSE) or inside (TRUE) in that order. For example, judgement_labels = c("Not inside","inside") or judgement_labels = c("No","Yes"). Default is c("no","yes)

Details

This is the prediction interval estimation procedure using parts from Jonathan Gillard's PhD thesis work and Wayne Fuller's book Measurement error models. This is recommended to use over Deming_CLSI_EP14_over_groups if k is most probably larger than 1

Value

A grouped data table enclosed with information regarding the estimated prediction band across the concentration range or for the particular evaluated materials for all groups

Examples

Deming_Gillard_Fuller_over_groups(MS_wise(sampled_cs_measurements), groups = "Comparison", level = 0.95, Np = 1, evaluated_materials = MS_wise(sampled_eqam_measurements))

pernille267/commutability.selectivity documentation built on July 20, 2022, 9:11 a.m.