Estimates the sample size and the optimum sample size ratio needed for a given width, costs, disease prevalence and significance level under the assumption of binormality.
SS(par1.1, par1.2, par2.1, par2.2=NULL, rho, width,
costs=matrix(c(0, 0, 1, (1-rho)/rho), 2, 2, byrow=TRUE),
R=NULL, var.equal=FALSE, alpha=0.05)an object of class SS which is a list with eight components:
sample size for the diseased group
sample size for the healthy group
sample size ratio between non-diseased and diseased subjects
width of the confidence interval provided by the user
significance level provided by the user
cost matrix provided by the user
R term, the product of the non-disease odds and the cost ratio (for further details about this concept, see References)
disease prevalence provided by the user
healthy population mean.
healthy population standard deviation.
diseased population mean.
diseased population standard deviation. It can be omitted when assuming equal variances (that is, when var.equal=TRUE) and in this situation the common variance is assumed to be equal to par1.2.
disease prevalence.
desired interval width.
cost matrix. Costs should be entered as a 2x2 matrix, where the first row corresponds to the true positive and true negative costs and the second row to the false positive and false negative costs. Default cost values are a combination of costs that yields R=1, which is equivalent to the Youden index method (for details about this concept, see References). It must be set to NULL if the user prefers to set R (see next argument).
if the cost matrix costs is not set, R desired (the algorithm will choose a suitable combination of costs that leads to R). Default, NULL (which leads to R=1 using the default costs).
a logical variable indicating whether to use equal variances. Default, FALSE.
significance level for the confidence interval. Default, 0.05.
Skaltsa K, Jover L, Carrasco JL. (2010). Estimation of the diagnostic threshold accounting for decision costs and sampling uncertainty. Biometrical Journal 52(5):676-697.
par1.1 <- 0
par1.2 <- 1
par2.1 <- 2
par2.2 <- 1
rho <- 0.3
width <- 0.5
SS(par1.1, par1.2, par2.1, par2.2, rho, width, var.equal=TRUE)
Run the code above in your browser using DataLab