Displaying 3 results from an estimated 3 matches for "aruike".
2007 Oct 29
2
help please
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.
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.
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-lgam...
2007 Oct 26
5
help
...d 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-lgam...
2007 Oct 25
1
help
...e 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)...