2011/3/11 Albyn Jones <jones@reed.edu>

> but presumably what you really want would be based on a joint confidence
> region for all the proportions.
>


I've had read "On Exact Methods for Testing Equality of Binomial
Proportions" by Akihito Matsuo, but  still, this concept is for me unclear
and I got lost...
We have H0: |pi1-pi2-pi3| = 0.05

n1<-40;n2<-40;n3<-40
s1<-list(1:11);s2<-list(1:17);s3<-list(1:15)
pi1<-max(s1[[1]])/n1;pi2<-max(s2[[1]])/n2;pi3<-max(s3[[1]])/n3
epsilon<-.05
t(c(pi1,pi2,pi3))

T_chi.sq<-sum(sapply(s1,(function(s1){(s1-n1*pi1)^2/n1*pi1*(1-pi1)})))
T_binom<-sum(sapply(s1,function(s1){choose(n1,s1)*(pi1-epsilon)^s1*(pi1+epsilon)^(n1-s1)}))
# Or it's about sum of all s<-list(1:43) and n<-n1+n2+n3 ?

Am I going to right direction?
-- 
Miłego dnia

	[[alternative HTML version deleted]]

