# [Rd] help

From: ezeh osita <ezehosita_at_yahoo.com>
Date: Thu, 25 Oct 2007 03:22:29 -0700 (PDT)

hello,

please can anyone help me out. Am a new user of R program. Am having problem with this code below, not getting the expected results.

1. Each m, the cumulative sum should be 1.000 but the 2nd and 3rd m returned 2.000 and 3.000 instead of 1.000.
2. to get the LCL(m) and UCL(m) for each m base on these instructions if out.cum > 0.025 then LCL(m)= y-1 if out.cum >0.975 then ucl(m)= y-1 how do I code these instructions into this code.

thanks
Aruike

pp=function(x,n,M){z=1.0;a=2.3071430;b=7.266064;H=3

out.h=c()

out.y=c()

out.m=c()

out.prob=c()

for(h in 1:H){

for(m in 1:M){

for(y in
0:m){

g=lgamma(m+z)+lgamma(n[h]+a+b)+lgamma(x[h]+y+a)+lgamma(n[h]+m+b-x[h]-y)

g=g-lgamma(y+z)-lgamma(m-y+z)-lgamma(x[h]+a)-lgamma(n[h]+b-x[h])- lgamma(n[h]+m+a+b)

out.h=c(out.h,h)

out.y=c(out.y,y)

out.m=c(out.m,m)

out.prob=c(out.prob, exp(g))

out.cum=c(cumsum(out.prob))

Result=data.frame(out.h,out.y,out.m,out.prob,out.cum)

}}}}
Kings=pp(x=c(19,20,30), n=c(52,60,80),3)  Kings

out.h out.y out.m
out.prob out.cum

1 1 0
1 0.65395431 0.6539543

2 1 1 1
0.34604569 1.0000000

3 1 0
2 0.43127277 1.4312728

4 1 1
2 0.44536308 1.8766358

5 1 2
2 0.12336415 2.0000000

6 1 0
3 0.28672775 2.2867277

7 1 1
3 0.43363507 2.7203628

8 1 2 3 0.23440955 2.9547724

9 1 3
3 0.04522764 3.0000000

Send instant messages to your online friends http://uk.messenger.yahoo.com

R-devel_at_r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-devel Received on Thu 25 Oct 2007 - 10:25:28 GMT

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.2.0, at Tue 30 Oct 2007 - 03:46:12 GMT.

Mailing list information is available at https://stat.ethz.ch/mailman/listinfo/r-devel. Please read the posting guide before posting to the list.