Group Abstract Group Abstract

Message Boards Message Boards

0
|
10.5K Views
|
11 Replies
|
1 Total Like
View groups...
Share
Share this post:

Histogram problem

Mikel Antoñana
Posted 11 years ago

I want to plot a Histogram, 

 Mu = Mean[Hamerdif];
    Sigma = StandardDeviation[Hamerdif];
    ir1 = Histogram[Hamerdif, {\[Sigma]/2^6}, "PDF"];

But sometimes I have the next error. I do not understand why.

Power::infy: Infinite expression 1/0 encountered. >>
Infinity::indet: Indeterminate expression 0. ComplexInfinity encountered. >>
POSTED BY: Mikel Antoñana
11 Replies
Sort By:
Mikel Antoñana
Mikel Antoñana
Posted 11 years ago

Finally I have noticed that something is wrong with Mathematica 10.2 version. For example with Mathematica 9.0 the execution is right.
POSTED BY: Mikel Antoñana

Histogram does a lot, but I've found it can be a bit wonky with some data sets. As noted scaling may be an issue, but I have no idea what the underlying algorithm scheme is. Just for fun try to Histogram the following set of numbers: {1, 6^-8, -3^-8, -6^-8, -8^-8, 3*^-8}. On my machine (Mac, 10.2.0.0) it hangs until I abort it. This is a pared down version of a case I actually ran into - confused me for a while.
POSTED BY: Michael Burns

My Mathematica does not hang (Mavericks with 10.2.0) but the computation takes ~18 seconds roughly...
POSTED BY: Sander Huisman
Mikel Antoñana
Mikel Antoñana
Posted 11 years ago

I think a normalization of the numbers could be a solution. Thanks for your help.
POSTED BY: Mikel Antoñana

Perhaps your numbers are too small. If I rescale them I get no errors:

Hamerdif = 10^17*{-7.75205766163486/10^17, -5.629097590588154/10^17, 
         -7.30880056733321/10^18, 8.376594318634389/10^17, 
    7.959982479649063/10^18, 
         -5.850425470324006/10^17, -7.563430648620345/10^17, 
    8.99610086221389/10^17, 
         8.047050888271066/10^17, 4.95873023031791/10^17}; 
Sigma = StandardDeviation[Hamerdif];
ir1 = Histogram[Hamerdif, {Sigma/2^6}, "PDF"]
POSTED BY: Gianluca Gorni

Strange. On my system (Mma 10.2 and Mac) I get a good histogram and no errors.

Attachments:
POSTED BY: Gianluca Gorni
Mikel Antoñana
Mikel Antoñana
Posted 11 years ago

I send you my values again. I think it is better to probe with the same "InputForm" values, Hamerdif = {-7.75205766163486^-17, -5.629097590588154^-17, -7.30880056733321^-18, 8.376594318634389^-17, 7.959982479649063^-18, -5.850425470324006^-17, -7.563430648620345^-17, 8.99610086221389^-17, 8.047050888271066^-17, 4.95873023031791^-17}; Sigma = StandardDeviation[Hamerdif]; ir1 = Histogram[Hamerdif, {Sigma/2^6}, "PDF"]
POSTED BY: Mikel Antoñana
David Keith
David Keith
Posted 11 years ago
POSTED BY: David Keith

What is Hamerdif? If it is a list of numbers I get no errors:

Hamerdif = {1, 2, 3, 4, 5};
Sigma = StandardDeviation[Hamerdif];
ir1 = Histogram[Hamerdif, {Sigma/2^6}, "PDF"]
POSTED BY: Gianluca Gorni
Mikel Antoñana
Mikel Antoñana
Posted 11 years ago

I get error for the next example: Hamerdif = {-7.7520610^-17, -5.629110^-17, -7.308810^-18, 8.3765910^-17, 7.9599810^-18, -5.8504310^-17, 7.5634310^-17, 8.996110^-17, 8.0470510^-17, 4.9587310^-17}; Sigma = StandardDeviation[Hamerdif]; ir1 = Histogram[Hamerdif, {Sigma/2^6}, "PDF"]

POSTED BY: Mikel Antoñana
David Keith
David Keith
Posted 11 years ago
POSTED BY: David Keith
Reply to this discussion
Community posts can be styled and formatted using the Markdown syntax.
Tag limit exceeded
Note: Only the first five people you tag will receive an email notification; the other tagged names will appear as links to their profiles.
Reply Preview
Attachments
Remove
or Discard