Services & Resources / Wolfram Forums
-----
 /
MathGroup Archive
2004
*January
*February
*March
*April
*May
*June
*July
*August
*September
*October
*November
*December
*Archive Index
*Ask about this page
*Print this page
*Give us feedback
*Sign up for the Wolfram Insider

MathGroup Archive 2004

[Date Index] [Thread Index] [Author Index]

Search the Archive

Re: bimodal ditribution form counting signs of Pi digits differences

  • To: mathgroup at smc.vnet.net
  • Subject: [mg51765] Re: bimodal ditribution form counting signs of Pi digits differences
  • From: Roger Bagula <tftn at earthlink.net>
  • Date: Mon, 1 Nov 2004 02:52:52 -0500 (EST)
  • References: <cm22ng$gip$1@smc.vnet.net>
  • Reply-to: tftn at earthlink.net
  • Sender: owner-wri-mathgroup at wolfram.com

Dear Bill Rowe,
I'll try you methods and see what I get.
I appreciate a way to speed up the count
as my method is dreadfully slow and
doesn't always come down.

The null  hypothesis is that the digits of Pi are random.
Texting the same kind of digit array but by a known
pseudorandom method gives a "measure"
of compariason between the two.
That the built in should apprear less random was a unexpected suprise.
I thought I was doing somnething wrong.
It appears I wasn't.
I've been trying to get a bigger sample since without much luck.
Dr. Bob has been very helpful and I'm grateful for both your help and 
that of a friend.
Bill Rowe wrote:

>On 10/30/04 at 3:49 AM, tftn at earthlink.net (Roger Bagula) wrote:
>
>  
>
>>Null hypothesis: the digits of Pi are random. To check make up a
>>random set of base 10 digits using Mathematica using:
>>Random[Integer,{0,9}]
>>I picked a seed off the top of my head. The result is a revolting
>>development. A trimodal distribution of noise pushed way positive.
>>The Pi digits behave more like an "ideal" probablity than the
>>Mathematica random! Anybody see what I did wrong? Dr. Bob you
>>always have an opinion, ha, ha...
>>    
>>
>
>  
>
>>(* random digit array base 10*)
>>SeedRandom[123]
>>a=Table[Random[Integer,{0,9}],{n,1,2000}];
>>b=Table[Sum[Sign[a[[m+1]]-a[[m]]],{m,1,n}],{n,1,Dimensions[a][[1]]-1}];
>>    
>>
>
>Your expression for b is rather convoluted to me but after looking at it, I see you are computing the cumulative sum of the Sign of two pairs of discrete uniform deviates. 
>
>An expression that computes the same result which should run much faster is
>
>Rest@FoldList[Plus, 0, Sign@ListConvolve[{1,-1},a]]
>
>  
>
>>ListPlot[b,PlotJoined->True]
>>(* distribution of the noise that results*)
>>Max[b]
>>Min[b]
>>c=Table[Count[b,m],{m,Min[b]-1,Max[b]+1}]
>>    
>>
>
>a faster expression for c is
>
>Flatten@{0,Length/@Split[Sort@b], 0}
>
>  
>
>>ListPlot[c,PlotJoined->True]
>>    
>>
>
>But other than suggesting ways to speed up your code, I don't know what you are doing "wrong" since I cannot determine what it is you expect nor what you think this code as to do with testing the null hypothesis you posted.
>
>What distribution are you expecting for b?
>
>When you take the sign of pairs of randomly chosen digits from 0 to 9, you clearly map the random numbers Mathematica generates to {-1,0,1}. I would expect the frequency of -1 to be roughly the same as 1 and higher than the frequency of 0. Given this, I would think the distribution of cumulative sums to be bimodal. That is I would expect runs of -1 to occur roughly as often as runs of 1 in b and runs of 0 to be relatively rare. Note, I've not taken the time to rigorously derive the distribution for b. This is simply my rough guess as to what should be expected
>
>Using your seed but increasing the number of samples significantly, i.e.,
>
>SeedRandom[123]; 
>a = Table[Random[Integer,  {0, 9}], {n, 200000}]; 
>ListPlot[b = Rest[FoldList[Plus, 0, 
>              Sign[ListConvolve[{1, -1}, a]]]],
>         PlotJoined -> True, Frame -> True, Axes -> None]; 
>d = ({First[#1], Length[#1]}&)/@Split[Sort[b]]; 
>ListPlot[d, PlotJoined ->True, Frame ->True, Axes -> None]; 
>
>My d here is {x,y} pairs where x is a an element of b and y is the number of times x occurs in b, i.e., ListPlot[d] generates a histogram of b. My d is essentially a rescaled version of your c except I've not force the begining and end values to 0 as you do for your c.
>
>Looking at this last plot, it does appear the distribution of b is bimodal giving me some confidence my rough guess above is correct.
>--
>To reply via email subtract one hundred and four
>
>  
>

-- 
Respectfully, Roger L. Bagula
tftn at earthlink.net, 11759Waterhill Road, Lakeside,Ca 92040-2905,tel: 619-5610814 :
alternative email: rlbtftn at netscape.net
URL :  http://home.earthlink.net/~tftn



  • Prev by Date: Re: Re: New version, new bugs
  • Next by Date: Re: Re: bimodal ditribution form counting signs of Pi digits differences
  • Previous by thread: Re: bimodal ditribution form counting signs of Pi digits differences
  • Next by thread: Re: Re: bimodal ditribution form counting signs of Pi digits differences