Displaying 1 result from an estimated 1 matches for "meantheta1".
2005 Apr 12
1
Cumulative Points and Confidence Interval Manipulation in barplot2
...ng the test data I am working on
(xzero):
xzero <- table(factor(WWNEW[HUNTTYPE=="DOVEONLY"], levels=0:12))
> xzero
0 1 2 3 4 5 6 7 8 9 10 11 12
179 20 9 2 2 0 1 0 0 0 0 0 0
> n <- sum(xzero)
> k <- sum(table(xzero))
> meantheta1 <-((2*xzero + 1)/(2*n + k))
> vartheta1
<-((2*(((2*n)+k)-((2*xzero)+1)))*((2*xzero)+1))/((((2*n)+k)^2)*(((2*n)+k)+2))
> stderr <- sqrt(vartheta1)
> cl.l <- meantheta1-(stderr*2)#Fake CI: Test
> cl.u <- meantheta1+(stderr*2)#Fake CI: Test
> barplot2(meantheta1, xlab=&q...