Displaying 2 results from an estimated 2 matches for "u_0".
Did you mean:
u10
2009 Mar 21
5
Fisher test problem
Hi, I noted a discrepancy between R and openepi when I ran a fisher test
with the same matrix. In R:
> a=matrix(c(1,2,6,17), nrow=2)
> a
[,1] [,2]
[1,] 1 6
[2,] 2 17
> fisher.test(a, conf.int=T)
Fisher's Exact Test for Count Data
data: a
p-value = 1
alternative hypothesis: true odds ratio is not equal to 1
95 percent confidence interval:
0.02061498
2004 Jul 10
1
Exact Maximum Likelihood Package
...*(1-s)^2 + 6*(1-p)*t^2*(1-t)^2;
f3 = 4*p*s^3*(1-s) + 4*(1-p)*t^3*(1-t);
f4 = p*s^4 + (1-p)*t^4;
The polynomial f_i represents the probability of seeing i successes.
Suppose we repeat this experiment 1000 times, and u_i is the number of
times we saw i successes.
The likelihood of this event is f_0^u_0*f_1^u_1*f_2^u_2*f_3^u_3*f_4^u_4,
and we seek to find those parameter values for s,t,p which maximize the
likelihood.
My Singular package has as input the 5 polynomials and a data vector u.
For the particular example of u = (3,5,7,11,13). The output are the
following
four roots (p,s,t) and the corr...