wants <- c("DescTools")
has <- wants %in% rownames(installed.packages())
if(any(!has)) install.packages(wants[!has])
DV <- factor(c("+", "+", "-", "+", "-", "+", "+"), levels=c("+", "-"))
N <- length(DV)
(tab <- table(DV))
DV
+ -
5 2
Exact binomial test
data: tab
number of successes = 5, number of trials = 7, p-value = 0.01288
alternative hypothesis: true probability of success is greater than 0.25
95 percent confidence interval:
0.3412614 1.0000000
sample estimates:
probability of success
0.7142857
Exact binomial test
data: hits and N
number of successes = 10, number of trials = 20, p-value = 0.01704
alternative hypothesis: true probability of success is not equal to 0.25
95 percent confidence interval:
0.2719578 0.7280422
sample estimates:
probability of success
0.5
[1] 0.01703563
est lwr.ci upr.ci
wilson 0.7142857 0.3589345 0.9177811
jeffreys 0.7142857 0.3523383 0.9352717
midp 0.7142857 0.3302059 0.9490209
3-sample test for equality of proportions without continuity
correction
data: hits out of total
X-squared = 50.587, df = 2, p-value = 1.035e-11
alternative hypothesis: two.sided
sample estimates:
prop 1 prop 2 prop 3
0.1462500 0.1220000 0.1796667
R markdown - markdown - R code - all posts