Programming the statistical procedures from SAS

fix for indeterminate 0^0

Reply
Occasional Contributor
Posts: 6

fix for indeterminate 0^0

Hello, the "unprotect" fix that seems to be recommended for eliminating the indeterminate 0^0 message doesn't seem to work in my case.  Advice please?

In[308]:= Unprotect[Power];

Power[0,0]=1

Out[309]= 1


(*My purpose is to calculate probability of order statistics ranges.

y5 and y1 are the max and min of a set of 5 order stats based on uniform distribution. In general, the min can be 0 and the max 1, I believe. *)

(*the pdf for the range, y1 to y5.  Note if y1=0, there is a term that is 0^0  *)

In[288]:= f[y1_, y5_] := (5!/(0!* 3!* 0!)) *    y1 ^0 * (y5 - y1) ^ 3 * (1 - y5)^0 * 1 * 1

In[289]:= f[y1, y5]  (*just to get display of the pdf, simplified *)

Out[289]= 20 (y5-y1)^3

In[310]:= Table[f[y1, y5], {y1, 0., 0.9, 0.1}]

During evaluation of In[310]:= Power::indet: Indeterminate expression 0.^0 encountered. >>

Out[310]= {Indeterminate,20. (y5-0.1)^3,20. (y5-0.2)^3,20. (y5-0.3)^3,20. (y5-0.4)^3,20. (y5-0.5)^3,20. (y5-0.6)^3,20. (y5-0.7)^3,20. (y5-0.8)^3,20. (y5-0.9)^3}

Ask a Question
Discussion stats
  • 0 replies
  • 123 views
  • 0 likes
  • 1 in conversation