MBESS (version 4.8.0)

mr.smd: Minimum risk point estimation of the population standardized mean difference

Description

A function for the sequential estimation of the standardized mean difference with minimum risk. The function implements the ideas of Chattopadhyay and Kelley (submitted, Psychological Methods), which considers study cost and accuracy of the estimated standardized mean difference simultaniously. This is important to specify that mr.smd.R was developed under the assumption of normally distributed data with equal sample size and equal cost of sampling per observation for each group.

Usage

mr.smd(A, structural.cost, epsilon, d, n, sampling.cost, pilot = FALSE, m0 = 4, 
gamma = 0.49)

Arguments

A

is the price one is willing to pay in order to have a maximum allowable difference of \(epsilon^2\) between the estimate of the standardized mean difference and its corresponding parameter.

structural.cost

epsilon

The maximum desired difference between the estimated standardized mean difference and the population value)

d

the current estimate of the standardized mean difference

n

current sample size per group (thus total sample size is \(2n\)); requires equal sample size per group.

sampling.cost

The sampling cost to collect an additional observation. For example, if each survey costs 10 dollars to distribute and score, sampling.cost would be 10 dollars per additional observation.

pilot

TRUE or FALSE based on whether the users is using the function to plan a pilot sample size (TRUE) or if it is being used to assess if the optimization criterion has been satisfied (FALSE)

m0

the minimum bound on the initial pilot sample size

gamma

A correction factor in which we suggest .49; see the two Chattopadhyay & Kelley articles for more details (ignorable for most users).

Value

Risk

Per group sample size (this simply repeats what was supplied to the function)

n1

Sample size for group 1 (echos the input value)

n1

Sample size for group 2 (echos the input value)

d

Observed value of the standardized mean difference (i.e., d; echos the input value)

Is.Satisfied?

A TRUE or FALSE statement of that evaluates a stopping rule using the risk function to determine if the optimation criterion has been satisfied (based on the goals of the researcher and current information available)

Details

The standardized mean difference is a widely used measure effect size. In this article, we developed a general theory for estimating the population standardized mean difference by minimizing both the mean square error of the estimator and the total sampling cost. This function implements our ideas discussed in Chattopadhyay and Kelley (submitted). See also Kelley and Rausch (2006) for additional information on the standardized mean difference.

References

Chattopadhyay, B., & Kelley, K. (submitted, minor revision requested). Estimating the standardized mean difference with minimum risk: Maximizing accuracy and minimizing cost with sequential estimation. Psychological Methods, X, X--X.

Chattopadhyay, B., & Kelley, K. (in press). Estimation of the Coefficient of Variation with Minimum Risk: A Sequential Method for Minimizing Sampling Error and Study Cost. Multivariate Behavioral Research, X, X--X.

Kelley, K., & Rausch, J. R. (2006). Sample size planning for the standardized mean difference: Accuracy in Parameter Estimation via narrow confidence intervals. Psychological Methods, 11, 363--385.

See Also

ci.smd, mr.cv

Examples

Run this code
# NOT RUN {
# To obtain pilot sample size in a situation in which A=10000. Note that 'A' is 
# 'structural.cost' divided by the square of 'epsilon'.

# From Chattopadhyay and Kelley (submitted, minor revision requested)
mr.smd(pilot=TRUE, A=10000, sampling.cost=2.4, gamma=.49)

High.SLS <- c(11, 7, 22, 13, 6, 9, 11, 16, 12, 17, 14, 8, 16)
Low.SLS  <- c(3, 6, 10, 8, 14, 5, 12, 10, 6, 8, 13, 5, 9)

mr.smd(d=1.021484, n=13, A=10000, sampling.cost=2.40, gamma=.49)

# Or, using the smd() function:
mr.smd(d=smd(Group.1=High.SLS, Group.2=Low.SLS), n=13, A=10000, sampling.cost=2.40, gamma=.49)

# Here, for this situation, the stopping rule is satisfied:
mr.smd(d=1.00, n=75, A=10000, sampling.cost=2.40, gamma=.49)

# }

Run the code above in your browser using DataCamp Workspace