Description Usage Arguments Details Value Note Author(s) References See Also Examples
View source: R/tTestLnormAltPower.R
Compute the power of a one or twosample ttest, given the sample size, ratio of means, coefficient of variation, and significance level, assuming lognormal data.
1 2 3 
n.or.n1 
numeric vector of sample sizes. When 
n2 
numeric vector of sample sizes for group 2. The default value is the value of

ratio.of.means 
numeric vector specifying the ratio of the first mean to the second mean.
When 
cv 
numeric vector of positive value(s) specifying the coefficient of
variation. When 
alpha 
numeric vector of numbers between 0 and 1 indicating the Type I error level
associated with the hypothesis test. The default value is 
sample.type 
character string indicating whether to compute power based on a onesample or
twosample hypothesis test. When 
alternative 
character string indicating the kind of alternative hypothesis. The possible values
are 
approx 
logical scalar indicating whether to compute the power based on an approximation to
the noncentral tdistribution. The default value is 
If the arguments n.or.n1
, n2
, ratio.of.means
, cv
, and
alpha
are not all the same length, they are replicated to be the same length
as the length of the longest argument.
OneSample Case (sample.type="one.sample"
)
Let \underline{x} = x_1, x_2, …, x_n denote a vector of n
observations from a lognormal distribution with mean
θ and coefficient of variation τ, and consider the null hypothesis:
H_0: θ = θ_0 \;\;\;\;\;\; (1)
The three possible alternative hypotheses are the upper onesided alternative
(alternative="greater"
):
H_a: θ > θ_0 \;\;\;\;\;\; (2)
the lower onesided alternative (alternative="less"
)
H_a: θ < θ_0 \;\;\;\;\;\; (3)
and the twosided alternative (alternative="two.sided"
)
H_a: θ \ne θ_0 \;\;\;\;\;\; (4)
To test the null hypothesis (1) versus any of the three alternatives (2)(4), one might be tempted to use Student's ttest based on the logtransformed observations. Unlike the twosample case with equal coefficients of variation (see below), in the onesample case Student's ttest applied to the logtransformed observations will not test the correct hypothesis, as now explained.
Let
y_i = log(x_i), \;\; i = 1, 2, …, n \;\;\;\;\;\; (5)
Then \underline{y} = y_1, y_2, …, y_n denote n observations from a normal distribution with mean μ and standard deviation σ, where
μ = log(\frac{θ}{√{τ^2 + 1}}) \;\;\;\;\;\; (6)
σ = [log(τ^2 + 1)]^{1/2} \;\;\;\;\;\; (7)
θ = exp[μ + (σ^2/2)] \;\;\;\;\;\; (8)
τ = [exp(σ^2)  1]^{1/2} \;\;\;\;\;\; (9)
(see the help file for LognormalAlt). Hence, by Equations (6) and (8) above, the Student's ttest on the logtransformed data would involve a test of hypothesis on both the parameters θ and τ, not just on θ.
To test the null hypothesis (1) above versus any of the alternatives (2)(4), you
can use the function elnormAlt
to compute a confidence interval for
θ, and use the relationship between confidence intervals and hypothesis
tests. To test the null hypothesis (1) above versus the upper onesided alternative
(2), you can also use
Chen's modified ttest for skewed distributions.
Although you can't use Student's ttest based on the logtransformed observations to test a hypothesis about θ, you can use the tdistribution to estimate the power of a test about θ that is based on confidence intervals or Chen's modified ttest, if you are willing to assume the population coefficient of variation τ stays constant for all possible values of θ you are interested in, and you are willing to postulate possible values for τ.
First, let's rewrite the hypotheses (1)(4) as follows. The null hypothesis (1) is equivalent to:
H_0: \frac{θ}{θ_0} = 1 \;\;\;\;\;\; (10)
The three possible alternative hypotheses are the upper onesided alternative
(alternative="greater"
)
H_a: \frac{θ}{θ_0} > 1 \;\;\;\;\;\; (11)
the lower onesided alternative (alternative="less"
)
H_a: \frac{θ}{θ_0} < 1 \;\;\;\;\;\; (12)
and the twosided alternative (alternative="two.sided"
)
H_a: \frac{θ}{θ_0} \ne 1 \;\;\;\;\;\; (13)
For a constant coefficient of variation τ, the standard deviation of the logtransformed observations σ is also constant (see Equation (7) above). Hence, by Equation (8), the ratio of the true mean to the hypothesized mean can be written as:
R = \frac{θ}{θ_0} = \frac{exp[μ + (σ^2/2)]}{exp[μ_0 + (σ^2/2)]} = \frac{e^μ}{e^μ_0} = e^{μ  μ_0} \;\;\;\;\;\; (14)
which only involves the difference
μ  μ_0 \;\;\;\;\;\; (15)
Thus, for given values of R and τ, the power of the test of the null hypothesis (10) against any of the alternatives (11)(13) can be computed based on the power of a onesample ttest with
\frac{δ}{σ} = \frac{log(R)}{√{log(τ^2 + 1)}} \;\;\;\;\;\; (16)
(see the help file for tTestPower
). Note that for the function
tTestLnormAltPower
, R corresponds to the argument ratio.of.means
,
and τ corresponds to the argument cv
.
TwoSample Case (sample.type="two.sample"
)
Let \underline{x}_1 = x_{11}, x_{12}, …, x_{1n_1} denote a vector of
n_1 observations from a lognormal distribution with mean
θ_1 and coefficient of variaiton τ, and let
\underline{x}_2 = x_{21}, x_{22}, …, x_{2n_2} denote a vector of
n_2 observations from a lognormal distribution with mean θ_2 and
coefficient of variation τ, and consider the null hypothesis:
H_0: θ_1 = θ_2 \;\;\;\;\;\; (17)
The three possible alternative hypotheses are the upper onesided alternative
(alternative="greater"
):
H_a: θ_1 > θ_2 \;\;\;\;\;\; (18)
the lower onesided alternative (alternative="less"
)
H_a: θ_1 < θ_2 \;\;\;\;\;\; (19)
and the twosided alternative (alternative="two.sided"
)
H_a: θ_1 \ne θ_2 \;\;\;\;\;\; (20)
Because we are assuming the coefficient of variation τ is the same for both populations, the test of the null hypothesis (17) versus any of the three alternatives (18)(20) can be based on the Student tstatistic using the logtransformed observations.
To show this, first, let's rewrite the hypotheses (17)(20) as follows. The null hypothesis (17) is equivalent to:
H_0: \frac{θ_1}{θ_2} = 1 \;\;\;\;\;\; (21)
The three possible alternative hypotheses are the upper onesided alternative
(alternative="greater"
)
H_a: \frac{θ_1}{θ_2} > 1 \;\;\;\;\;\; (22)
the lower onesided alternative (alternative="less"
)
H_a: \frac{θ_1}{θ_2} < 1 \;\;\;\;\;\; (23)
and the twosided alternative (alternative="two.sided"
)
H_a: \frac{θ_1}{θ_2} \ne 1 \;\;\;\;\;\; (24)
If coefficient of variation τ is the same for both populations, then the standard deviation of the logtransformed observations σ is also the same for both populations (see Equation (7) above). Hence, by Equation (8), the ratio of the means can be written as:
R = \frac{θ_1}{θ_2} = \frac{exp[μ_1 + (σ^2/2)]}{exp[μ_2 + (σ^2/2)]} = \frac{e^μ_1}{e^μ_2} = e^{μ_1  μ_2} \;\;\;\;\;\; (25)
which only involves the difference
μ_1  μ_2 \;\;\;\;\;\; (26)
Thus, for given values of R and τ, the power of the test of the null hypothesis (21) against any of the alternatives (22)(24) can be computed based on the power of a twosample ttest with
\frac{δ}{σ} = \frac{log(R)}{√{log(τ^2 + 1)}} \;\;\;\;\;\; (27)
(see the help file for tTestPower
). Note that for the function
tTestLnormAltPower
, R corresponds to the argument ratio.of.means
,
and τ corresponds to the argument cv
.
a numeric vector powers.
The normal distribution and
lognormal distribution are probably the two most
frequently used distributions to model environmental data. Often, you need to
determine whether a population mean is significantly different from a specified
standard (e.g., an MCL or ACL, USEPA, 1989b, Section 6), or whether two different
means are significantly different from each other (e.g., USEPA 2009, Chapter 16).
When you have lognormallydistributed data, you have to be careful about making
statements regarding inference for the mean. For the twosample case with
assumed equal coefficients of variation, you can perform the
Student's ttest on the logtransformed observations.
For the onesample case, you can perform a hypothesis test by constructing a
confidence interval for the mean using elnormAlt
, or use
Chen's ttest modified for skewed data.
In the course of designing a sampling program, an environmental scientist may wish
to determine the relationship between sample size, significance level, power, and
scaled difference if one of the objectives of the sampling program is to determine
whether a mean differs from a specified level or two means differ from each other.
The functions tTestLnormAltPower
, tTestLnormAltN
,
tTestLnormAltRatioOfMeans
, and plotTTestLnormAltDesign
can be used to investigate these relationships for the case of
lognormallydistributed observations.
Steven P. Millard ([email protected])
van Belle, G., and D.C. Martin. (1993). Sample Size as a Function of Coefficient of Variation and Ratio of Means. The American Statistician 47(3), 165–167.
Also see the list of references in the help file for tTestPower
.
tTestLnormAltN
, tTestLnormAltRatioOfMeans
,
plotTTestLnormAltDesign
, LognormalAlt,
t.test
, Hypothesis Tests.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103  # Look at how the power of the onesample test increases with increasing
# sample size:
seq(5, 30, by = 5)
#[1] 5 10 15 20 25 30
power < tTestLnormAltPower(n.or.n1 = seq(5, 30, by = 5),
ratio.of.means = 1.5, cv = 1)
round(power, 2)
#[1] 0.14 0.28 0.42 0.54 0.65 0.73
#
# Repeat the last example, but use the approximation to the power instead of the
# exact power. Note how the approximation underestimates the true power for
# the smaller sample sizes:
power < tTestLnormAltPower(n.or.n1 = seq(5, 30, by = 5),
ratio.of.means = 1.5, cv = 1, approx = TRUE)
round(power, 2)
#[1] 0.09 0.25 0.40 0.53 0.64 0.73
#==========
# Look at how the power of the twosample ttest increases with increasing
# ratio of means:
power < tTestLnormAltPower(n.or.n1 = 20, sample.type = "two",
ratio.of.means = c(1.1, 1.5, 2), cv = 1)
round(power, 2)
#[1] 0.06 0.32 0.73
#
# Look at how the power of the twosample ttest increases with increasing
# values of Type I error:
power < tTestLnormAltPower(30, sample.type = "two", ratio.of.means = 1.5,
cv = 1, alpha = c(0.001, 0.01, 0.05, 0.1))
round(power, 2)
#[1] 0.07 0.23 0.46 0.59
#==========
# The guidance document Soil Screening Guidance: Technical Background Document
# (USEPA, 1996c, Part 4) discusses sampling design and sample size calculations
# for studies to determine whether the soil at a potentially contaminated site
# needs to be investigated for possible remedial action. Let 'theta' denote the
# average concentration of the chemical of concern. The guidance document
# establishes the following goals for the decision rule (USEPA, 1996c, p.87):
#
# Pr[Decide Don't Investigate  theta > 2 * SSL] = 0.05
#
# Pr[Decide to Investigate  theta <= (SSL/2)] = 0.2
#
# where SSL denotes the preestablished soil screening level.
#
# These goals translate into a Type I error of 0.2 for the null hypothesis
#
# H0: [theta / (SSL/2)] <= 1
#
# and a power of 95% for the specific alternative hypothesis
#
# Ha: [theta / (SSL/2)] = 4
#
# Assuming a lognormal distribution with a coefficient of variation of 2,
# determine the power associated with various sample sizes for this onesample test.
# Based on these calculations, you need to take at least 6 soil samples to
# satisfy the requirements for the Type I and Type II errors.
power < tTestLnormAltPower(n.or.n1 = 2:8, ratio.of.means = 4, cv = 2,
alpha = 0.2, alternative = "greater")
names(power) < paste("N=", 2:8, sep = "")
round(power, 2)
# N=2 N=3 N=4 N=5 N=6 N=7 N=8
#0.65 0.80 0.88 0.93 0.96 0.97 0.98
#
# Repeat the last example, but use the approximate power calculation instead of
# the exact one. Using the approximate power calculation, you need at least
# 7 soil samples instead of 6 (because the approximation underestimates the power).
power < tTestLnormAltPower(n.or.n1 = 2:8, ratio.of.means = 4, cv = 2,
alpha = 0.2, alternative = "greater", approx = TRUE)
names(power) < paste("N=", 2:8, sep = "")
round(power, 2)
# N=2 N=3 N=4 N=5 N=6 N=7 N=8
#0.55 0.75 0.84 0.90 0.93 0.95 0.97
#==========
# Clean up
#
rm(power)

Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.