pairwiseCImethodsProp: Confidence intervals for two sample comparisons of binomial...
In pairwiseCI: Confidence Intervals for Two Sample Comparisons
Description
Usage
Arguments
Details
Value
References
See Also
Examples
Description
For the comparison of two independent samples of binomial observations, confidence intervals for
the difference (RD), ratio (RR) and odds ratio (OR) of proportions are implemented.
Usage
1
2
3
4
5
6
7
8
Arguments
x
observations of the first sample: either a vector with number of successes and failures, or a data.frame with two columns (the successes and failures))
y
observations of the second sample: either a vector with number of successes and failures, or a data.frame with two columns (the successes and failures))
alternative
character string, either "two.sided", "less" or "greater"
conf.level
the comparisonwise confidence level of the intervals, where 0.95 is the default
CImethod
a single character string, see below for details
...
further arguments to be passed to the individual methods, see details
Details
Generally, the input are two vectors x and y giving the number of successes and failures in the two samples,
or, alternatively, two data.frames x and y each containing one column for the successes
and one column for the failures, and the rows containing repeated observations from the same treatment.
Please note, that the confidence intervals available in this function assume counts of successes and failures from a binomial distribution and thus do NOT APPROPRIATELY account for extra-binomial variability between repeated observations for the same treatment!
When there are repeated observations (input as a data.frame with several rows), intervals are calculated based on the sums over the rows of success and failure!
The following methods are available for the risk difference:
-
Prop.diff: asymptotic continuity corrected confidence interval for the difference of proportions by calling prop.test in package stats, see ?prop.test for details,
-
Prop.diff with CImethod="AC": asymptotic Wald-type interval after adding one pseudo-observation to each cell (Agresti and Caffo, 2000).
-
Prop.diff with CImethod="NHS": asymptotic Newcombes Hybrid Score Interval (Newcombe, 1998).
For the risk ratio:
-
Prop.ratio with CImethod="Score": asymptotic Score interval for the ratio of proportions (by iteratively inverting a Chi-Squared test) according to Koopman (1984), following the description by Gart and Nam (1988). This method does NOT involve the skewness correction or extensions to stratification described in Gart and Nam (1988).
-
Prop.ratio with CImethod="MNScore": asymptotic Score interval for the ratio of proportions (by iteratively inverting a Chi-Squared test) according to Miettinen and Nurminen (1985), following the description by Gart and Nam (1988).
-
Prop.ratio with CImethod="MOVER": asymptotic method of variance estimate recovery for ratios (Donner and Zou, 2012; Fagerland and Newcombe, 2013), relying on the Wilson Score interval to obtain the confidence limits for the single proportions.
-
Prop.ratio with CImethod="GNC": crude normal approximation on the log-scale after adding 0.5 to the observed number of successes and the samples sizes.
For the odds ratio:
-
Prop.or with CImethod="Woolf" asymptotic adjusted Woolf confidence interval for the odds ratio of proportions (normal approximation after adding 0.5 to each cell count), as, e.g., described in Lawson (2005).
-
Prop.or with CImethod="Exact" calculates the exact confidence interval for the odds ratio of proportions corresponding to Fishers exact test, by calling to fisher.test in stats. For details, see ?fisher.test.
Value
A list containing:
conf.int
a vector containing the lower and upper confidence limit, and the methods name as an attribute
estimate
a single named value
quantile
the quantile used for constructing the interval
conf.level
the confidence level
References
-
Agresti A and Caffo B (2000): Simple and effective confidence intervals for proportions and differences of proportions result from adding two successes and two failures. American Statistician 54 (4), 280-288.
-
Donner A and Zou GY (2012): Closed-form confidence intervals for functions of the normal mean and standard deviation.Statistical Methods in Medical Research 21(4):347-359.
-
Fagerland MW, Newcombe RG (2013): Confidence intervals for odds ratio and relative risk based on the inverse hyperbolic sine transformation. Statistics in Medicine, DOI: 10.1002/sim.5714
-
Gart JJ and Nam J (1988): Approximate interval estimation of the ratio of binomial parameters: A review and corrections for skewness. Biometrics 44, 323-338.
-
Koopman, PAR (1984): Confidence Intervals for the Ratio of Two Binomial Proportions. Biometrics 40(2), 513-517.
-
Lawson R (2005): Small sample confidence intervals for the odds ratio. Communication in Statistics Simulation and Computation, 33, 1095-1113.
-
Miettinen O and Nurminen M (1985): Comparative Analysis of Two Rates. Statistics in Medicine 4, 213-226.
-
Newcombe RG (1998): Interval Estimation for the Difference Between Independent Proportions: Comparison of Eleven Methods. Statistics in Medicine 17, 873-890.
-
Venables WN and Ripley BD (2002). Modern Applied Statistics with S. Fourth Edition. Springer New York.
See Also
An alternative implementation of the Score interval for the risk ratio in package propCIs, function riskscoreci.
Examples
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
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
# The rooting data.
data(rooting)
# the first comparison should be the same as:
Age5_PosB_IBA0 <- subset(rooting,
Age=="5" & Position=="B" & IBA=="0")[,c("root", "noroot")]
Age5_PosB_IBA0.5 <- subset(rooting,
Age=="5" & Position=="B" & IBA=="0.5")[,c("root", "noroot")]
Age5_PosB_IBA0
Age5_PosB_IBA0.5
Prop.diff(x=Age5_PosB_IBA0, y=Age5_PosB_IBA0.5)
Prop.ratio(x=Age5_PosB_IBA0, y=Age5_PosB_IBA0.5)
Prop.or(x=Age5_PosB_IBA0, y=Age5_PosB_IBA0.5)
# is the same as input two vectors x,y each containing
# the count of successes and the count of failures
colSums(Age5_PosB_IBA0)
colSums(Age5_PosB_IBA0.5)
Prop.diff(x=c(16,32),y=c(29,19))
Prop.ratio(x=c(16,32),y=c(29,19))
Prop.or(x=c(16,32),y=c(29,19))
# # #
# Comparison with original publications:
# I. Risk difference:
# Continuity corrected interval:
# 1.Comparison with results presented in Newcombe (1998),
# Table II, page 877, 10. Score, CC
# column 1 (a): 56/70-48/80: [0.0441; 0.3559]
Prop.diff(x=c(56,70-56),y=c(48,80-48), alternative="two.sided",
conf.level=0.95, CImethod="CC")
# Risk difference, NHS
# Newcombes Hybrid Score interval:
# 1.Comparison with results presented in Newcombe (1998),
# Table II, page 877, 10. Score, noCC
# column 1 (a): 56/70-48/80: [0.0524; 0.3339]
Prop.diff(x=c(56,70-56),y=c(48,80-48), alternative="two.sided",
conf.level=0.95, CImethod="NHS")
Prop.diff(x=c(56,70-56),y=c(48,80-48), alternative="greater",
conf.level=0.975, CImethod="NHS")
Prop.diff(x=c(56,70-56),y=c(48,80-48), alternative="less",
conf.level=0.975, CImethod="NHS")
# 2.Comparison with results presented in Newcombe (1998),
# Table II, page 877, 10. Score, noCC
# column 2 (b): 9/10-3/10: [0.1705; 0.8090]
Prop.diff(x=c(9,1),y=c(3,7), alternative="two.sided",
conf.level=0.95, CImethod="NHS")
# 3.Comparison with results presented in Newcombe (1998),
# Table II, page 877, 10. Score, noCC
# column 2 (h): 10/10-0/10: [0.6075; 1.000]
Prop.diff(x=c(10,0),y=c(0,10), alternative="two.sided",
conf.level=0.95, CImethod="NHS")
# II. Risk ratio,
# Score interval according to Koopman(1984), Gart and Nam (1988)
# 1.Comparison with results presented in Gart and Nam (1998),
# Section 5 (page 327), Example 1
# x1/n1=8/15 x0/n0=4/15:
# Log: [0.768, 4.65]
# Score: [0.815; 5.34]
# Log (GNC)
Prop.ratio(x=c(8,7),y=c(4,11), alternative="two.sided",
conf.level=0.95, CImethod="GNC")
# Score (Score)
Prop.ratio(x=c(8,7),y=c(4,11), alternative="two.sided",
conf.level=0.95, CImethod="Score")
Prop.ratio(x=c(8,7),y=c(4,11), alternative="less",
conf.level=0.975, CImethod="Score")
Prop.ratio(x=c(8,7),y=c(4,11), alternative="greater",
conf.level=0.975, CImethod="Score")
# 2.Comparison with results presented in Gart and Nam (1998),
# Section 5 (page 328), Example 2
# x1/n1=6/10 x0/n0=6/20:
# Crude Log: [0.883, 4.32]
# Score: [0.844; 4.59]
Prop.ratio(x=c(6,4),y=c(6,14), alternative="two.sided",
conf.level=0.95, CImethod="GNC")
Prop.ratio(x=c(6,4),y=c(6,14), alternative="two.sided",
conf.level=0.95, CImethod="Score")
# Koopman (1984), page 517
# x1=36, n1=40, x0=16, n0=80: [2.94, 7.15]
Prop.ratio(x=c(36, 4), y=c(16, 64), CImethod="Score")$conf.int
# Miettinen, Nurminen (1985) p. 217 (Example 6):
# x1=10, n1=10, x0=20, n0=20: [0.72, 1.20]
Prop.ratio(x=c(10, 0), y=c(20, 0), CImethod="MNScore")$conf.int
# MOVER-R Wilson in Newcombe and Fagerland, 2013, Table VIII:
# x1=24, n1=73,x0=53, n0=65: [0.282, 0.563]
Prop.ratio(x=c(24, 49), y=c(53, 12), CImethod="MOVER")$conf.int
# x1=29, n1=55, x0=11, n0=11: [0.398, 0.761]
Prop.ratio(x=c(29, 26), y=c(11,0), CImethod="MOVER")$conf.int
# x1=7, n1=18, x0=1, n0=18: [(1.27, 40.8)]
Prop.ratio(x=c(7, 11), y=c(1, 17), CImethod="MOVER")$conf.int
pairwiseCI documentation built on May 1, 2019, 6:51 p.m.
R Package Documentation
Browse R Packages
We want your feedback!
Note that we can't provide technical support on individual packages. You should contact the package authors for that.
Description Usage Arguments Details Value References See Also Examples
Description
For the comparison of two independent samples of binomial observations, confidence intervals for the difference (RD), ratio (RR) and odds ratio (OR) of proportions are implemented.
Usage
1 2 3 4 5 6 7 8 |
Arguments
x |
observations of the first sample: either a vector with number of successes and failures, or a data.frame with two columns (the successes and failures)) |
y |
observations of the second sample: either a vector with number of successes and failures, or a data.frame with two columns (the successes and failures)) |
alternative |
character string, either "two.sided", "less" or "greater" |
conf.level |
the comparisonwise confidence level of the intervals, where 0.95 is the default |
CImethod |
a single character string, see below for details |
... |
further arguments to be passed to the individual methods, see details |
Details
Generally, the input are two vectors x and y giving the number of successes and failures in the two samples, or, alternatively, two data.frames x and y each containing one column for the successes and one column for the failures, and the rows containing repeated observations from the same treatment.
Please note, that the confidence intervals available in this function assume counts of successes and failures from a binomial distribution and thus do NOT APPROPRIATELY account for extra-binomial variability between repeated observations for the same treatment! When there are repeated observations (input as a data.frame with several rows), intervals are calculated based on the sums over the rows of success and failure!
The following methods are available for the risk difference:
-
Prop.diff: asymptotic continuity corrected confidence interval for the difference of proportions by calling prop.test in package stats, see ?prop.test for details,
-
Prop.diff with CImethod="AC": asymptotic Wald-type interval after adding one pseudo-observation to each cell (Agresti and Caffo, 2000).
-
Prop.diff with CImethod="NHS": asymptotic Newcombes Hybrid Score Interval (Newcombe, 1998).
For the risk ratio:
-
Prop.ratio with CImethod="Score": asymptotic Score interval for the ratio of proportions (by iteratively inverting a Chi-Squared test) according to Koopman (1984), following the description by Gart and Nam (1988). This method does NOT involve the skewness correction or extensions to stratification described in Gart and Nam (1988).
-
Prop.ratio with CImethod="MNScore": asymptotic Score interval for the ratio of proportions (by iteratively inverting a Chi-Squared test) according to Miettinen and Nurminen (1985), following the description by Gart and Nam (1988).
-
Prop.ratio with CImethod="MOVER": asymptotic method of variance estimate recovery for ratios (Donner and Zou, 2012; Fagerland and Newcombe, 2013), relying on the Wilson Score interval to obtain the confidence limits for the single proportions.
-
Prop.ratio with CImethod="GNC": crude normal approximation on the log-scale after adding 0.5 to the observed number of successes and the samples sizes.
For the odds ratio:
-
Prop.or with CImethod="Woolf" asymptotic adjusted Woolf confidence interval for the odds ratio of proportions (normal approximation after adding 0.5 to each cell count), as, e.g., described in Lawson (2005).
-
Prop.or with CImethod="Exact" calculates the exact confidence interval for the odds ratio of proportions corresponding to Fishers exact test, by calling to fisher.test in stats. For details, see ?fisher.test.
Value
A list containing:
conf.int |
a vector containing the lower and upper confidence limit, and the methods name as an attribute |
estimate |
a single named value |
quantile |
the quantile used for constructing the interval |
conf.level |
the confidence level |
References
-
Agresti A and Caffo B (2000): Simple and effective confidence intervals for proportions and differences of proportions result from adding two successes and two failures. American Statistician 54 (4), 280-288.
-
Donner A and Zou GY (2012): Closed-form confidence intervals for functions of the normal mean and standard deviation.Statistical Methods in Medical Research 21(4):347-359.
-
Fagerland MW, Newcombe RG (2013): Confidence intervals for odds ratio and relative risk based on the inverse hyperbolic sine transformation. Statistics in Medicine, DOI: 10.1002/sim.5714
-
Gart JJ and Nam J (1988): Approximate interval estimation of the ratio of binomial parameters: A review and corrections for skewness. Biometrics 44, 323-338.
-
Koopman, PAR (1984): Confidence Intervals for the Ratio of Two Binomial Proportions. Biometrics 40(2), 513-517.
-
Lawson R (2005): Small sample confidence intervals for the odds ratio. Communication in Statistics Simulation and Computation, 33, 1095-1113.
-
Miettinen O and Nurminen M (1985): Comparative Analysis of Two Rates. Statistics in Medicine 4, 213-226.
-
Newcombe RG (1998): Interval Estimation for the Difference Between Independent Proportions: Comparison of Eleven Methods. Statistics in Medicine 17, 873-890.
-
Venables WN and Ripley BD (2002). Modern Applied Statistics with S. Fourth Edition. Springer New York.
See Also
An alternative implementation of the Score interval for the risk ratio in package propCIs, function riskscoreci.
Examples
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 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136 137 138 139 140 141 142 143 144 | # The rooting data.
data(rooting)
# the first comparison should be the same as:
Age5_PosB_IBA0 <- subset(rooting,
Age=="5" & Position=="B" & IBA=="0")[,c("root", "noroot")]
Age5_PosB_IBA0.5 <- subset(rooting,
Age=="5" & Position=="B" & IBA=="0.5")[,c("root", "noroot")]
Age5_PosB_IBA0
Age5_PosB_IBA0.5
Prop.diff(x=Age5_PosB_IBA0, y=Age5_PosB_IBA0.5)
Prop.ratio(x=Age5_PosB_IBA0, y=Age5_PosB_IBA0.5)
Prop.or(x=Age5_PosB_IBA0, y=Age5_PosB_IBA0.5)
# is the same as input two vectors x,y each containing
# the count of successes and the count of failures
colSums(Age5_PosB_IBA0)
colSums(Age5_PosB_IBA0.5)
Prop.diff(x=c(16,32),y=c(29,19))
Prop.ratio(x=c(16,32),y=c(29,19))
Prop.or(x=c(16,32),y=c(29,19))
# # #
# Comparison with original publications:
# I. Risk difference:
# Continuity corrected interval:
# 1.Comparison with results presented in Newcombe (1998),
# Table II, page 877, 10. Score, CC
# column 1 (a): 56/70-48/80: [0.0441; 0.3559]
Prop.diff(x=c(56,70-56),y=c(48,80-48), alternative="two.sided",
conf.level=0.95, CImethod="CC")
# Risk difference, NHS
# Newcombes Hybrid Score interval:
# 1.Comparison with results presented in Newcombe (1998),
# Table II, page 877, 10. Score, noCC
# column 1 (a): 56/70-48/80: [0.0524; 0.3339]
Prop.diff(x=c(56,70-56),y=c(48,80-48), alternative="two.sided",
conf.level=0.95, CImethod="NHS")
Prop.diff(x=c(56,70-56),y=c(48,80-48), alternative="greater",
conf.level=0.975, CImethod="NHS")
Prop.diff(x=c(56,70-56),y=c(48,80-48), alternative="less",
conf.level=0.975, CImethod="NHS")
# 2.Comparison with results presented in Newcombe (1998),
# Table II, page 877, 10. Score, noCC
# column 2 (b): 9/10-3/10: [0.1705; 0.8090]
Prop.diff(x=c(9,1),y=c(3,7), alternative="two.sided",
conf.level=0.95, CImethod="NHS")
# 3.Comparison with results presented in Newcombe (1998),
# Table II, page 877, 10. Score, noCC
# column 2 (h): 10/10-0/10: [0.6075; 1.000]
Prop.diff(x=c(10,0),y=c(0,10), alternative="two.sided",
conf.level=0.95, CImethod="NHS")
# II. Risk ratio,
# Score interval according to Koopman(1984), Gart and Nam (1988)
# 1.Comparison with results presented in Gart and Nam (1998),
# Section 5 (page 327), Example 1
# x1/n1=8/15 x0/n0=4/15:
# Log: [0.768, 4.65]
# Score: [0.815; 5.34]
# Log (GNC)
Prop.ratio(x=c(8,7),y=c(4,11), alternative="two.sided",
conf.level=0.95, CImethod="GNC")
# Score (Score)
Prop.ratio(x=c(8,7),y=c(4,11), alternative="two.sided",
conf.level=0.95, CImethod="Score")
Prop.ratio(x=c(8,7),y=c(4,11), alternative="less",
conf.level=0.975, CImethod="Score")
Prop.ratio(x=c(8,7),y=c(4,11), alternative="greater",
conf.level=0.975, CImethod="Score")
# 2.Comparison with results presented in Gart and Nam (1998),
# Section 5 (page 328), Example 2
# x1/n1=6/10 x0/n0=6/20:
# Crude Log: [0.883, 4.32]
# Score: [0.844; 4.59]
Prop.ratio(x=c(6,4),y=c(6,14), alternative="two.sided",
conf.level=0.95, CImethod="GNC")
Prop.ratio(x=c(6,4),y=c(6,14), alternative="two.sided",
conf.level=0.95, CImethod="Score")
# Koopman (1984), page 517
# x1=36, n1=40, x0=16, n0=80: [2.94, 7.15]
Prop.ratio(x=c(36, 4), y=c(16, 64), CImethod="Score")$conf.int
# Miettinen, Nurminen (1985) p. 217 (Example 6):
# x1=10, n1=10, x0=20, n0=20: [0.72, 1.20]
Prop.ratio(x=c(10, 0), y=c(20, 0), CImethod="MNScore")$conf.int
# MOVER-R Wilson in Newcombe and Fagerland, 2013, Table VIII:
# x1=24, n1=73,x0=53, n0=65: [0.282, 0.563]
Prop.ratio(x=c(24, 49), y=c(53, 12), CImethod="MOVER")$conf.int
# x1=29, n1=55, x0=11, n0=11: [0.398, 0.761]
Prop.ratio(x=c(29, 26), y=c(11,0), CImethod="MOVER")$conf.int
# x1=7, n1=18, x0=1, n0=18: [(1.27, 40.8)]
Prop.ratio(x=c(7, 11), y=c(1, 17), CImethod="MOVER")$conf.int
|
R Package Documentation
Browse R Packages
We want your feedback!
Note that we can't provide technical support on individual packages. You should contact the package authors for that.
Embedding an R snippet on your website
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.