octave-maintainers
[Top][All Lists]
Advanced

[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

betai -> betainc


From: Paul Kienzle
Subject: betai -> betainc
Date: Fri, 9 Aug 2002 19:21:46 -0400

Paul Kienzle <address@hidden>

        * statistics/distributions/beta_cdf.m: Replace betai with betainc.
        * statistics/distributions/binomial_cdf.m: Likewise.
        * statistics/distributions/f_cdf.m: Likewise.
        * statistics/distributions/t_cdf.m: Likewise.


Index: statistics/distributions/beta_cdf.m
===================================================================
RCS file: /cvs/octave/scripts/statistics/distributions/beta_cdf.m,v
retrieving revision 1.4
diff -c -r1.4 beta_cdf.m
*** statistics/distributions/beta_cdf.m 2002/05/01 06:48:35     1.4
--- statistics/distributions/beta_cdf.m 2002/08/09 19:31:01
***************
*** 57,63 ****
  
    k = find ((x > 0) & (x < 1) & (a > 0) & (b > 0));
    if (any (k))
!     cdf (k) = betai (a(k), b(k), x(k));
    endif
  
    cdf = reshape (cdf, r, c);
--- 57,63 ----
  
    k = find ((x > 0) & (x < 1) & (a > 0) & (b > 0));
    if (any (k))
!     cdf (k) = betainc (x(k), a(k), b(k));
    endif
  
    cdf = reshape (cdf, r, c);
Index: statistics/distributions/binomial_cdf.m
===================================================================
RCS file: /cvs/octave/scripts/statistics/distributions/binomial_cdf.m,v
retrieving revision 1.4
diff -c -r1.4 binomial_cdf.m
*** statistics/distributions/binomial_cdf.m     2002/05/01 06:48:35     1.4
--- statistics/distributions/binomial_cdf.m     2002/08/09 19:31:01
***************
*** 60,66 ****
            & (p >= 0) & (p <= 1));
    if (any (k))
      tmp = floor (x(k));
!     cdf(k) = 1 - betai (tmp + 1, n(k) - tmp, p(k));
    endif
  
    cdf = reshape (cdf, r, c);
--- 60,66 ----
            & (p >= 0) & (p <= 1));
    if (any (k))
      tmp = floor (x(k));
!     cdf(k) = 1 - betainc (p(k), tmp + 1, n(k) - tmp);
    endif
  
    cdf = reshape (cdf, r, c);
Index: statistics/distributions/f_cdf.m
===================================================================
RCS file: /cvs/octave/scripts/statistics/distributions/f_cdf.m,v
retrieving revision 1.5
diff -c -r1.5 f_cdf.m
*** statistics/distributions/f_cdf.m    2002/05/01 06:48:35     1.5
--- statistics/distributions/f_cdf.m    2002/08/09 19:31:01
***************
*** 57,63 ****
  
    k = find ((x > 0) & (x < Inf) & (m > 0) & (n > 0));
    if (any (k))
!     cdf(k) = 1 - betai (n(k) / 2, m(k) / 2, 1 ./ (1 + m(k) .* x(k) ./ n(k)));
    endif
  
    cdf = reshape (cdf, r, c);
--- 57,63 ----
  
    k = find ((x > 0) & (x < Inf) & (m > 0) & (n > 0));
    if (any (k))
!     cdf(k) = 1 - betainc (1 ./ (1 + m(k) .* x(k) ./ n(k)), n(k) / 2, m(k) / 
2);
    endif
  
    cdf = reshape (cdf, r, c);
Index: statistics/distributions/t_cdf.m
===================================================================
RCS file: /cvs/octave/scripts/statistics/distributions/t_cdf.m,v
retrieving revision 1.5
diff -c -r1.5 t_cdf.m
*** statistics/distributions/t_cdf.m    2002/05/01 06:48:36     1.5
--- statistics/distributions/t_cdf.m    2002/08/09 19:31:01
***************
*** 56,62 ****
  
    k = find ((x > -Inf) & (x < Inf) & (n > 0));
    if (any (k))
!     cdf(k) = betai (n(k) / 2, 1 / 2, 1 ./ (1 + x(k) .^ 2 ./ n(k))) / 2;
      ind = find (x(k) > 0);
      if (any (ind))
        cdf(k(ind)) = 1 - cdf(k(ind));
--- 56,62 ----
  
    k = find ((x > -Inf) & (x < Inf) & (n > 0));
    if (any (k))
!     cdf(k) = betainc (1 ./ (1 + x(k) .^ 2 ./ n(k)), n(k) / 2, 1 / 2) / 2;
      ind = find (x(k) > 0);
      if (any (ind))
        cdf(k(ind)) = 1 - cdf(k(ind));



reply via email to

[Prev in Thread] Current Thread [Next in Thread]