Re: Re: RandomReplacement
- To: mathgroup at smc.vnet.net
- Subject: [mg51414] Re: [mg51380] Re: RandomReplacement
- From: sean kim <sean_incali at yahoo.com>
- Date: Sat, 16 Oct 2004 04:20:58 -0400 (EDT)
- Reply-to: sean_incali at yahoo.com
- Sender: owner-wri-mathgroup at wolfram.com
actually, I made a mistake in the last post. i forgot the index for RandomArray. if it's put back in, i get the similar slowing down effect. sorry for any confusions. but i noticed something else peculiar. Please consider the following. without the package loaded, i get In[1]:= <<"Statistics`NormalDistribution`"; In[2]:= Table[RandomArray[NormalDistribution[50,25], 1],{10^4}];//Timing Out[2]= {0.992 Second,Null} In[3]:= Table[RandomArray[NormalDistribution[50,25], 1],{10^5}];//Timing Out[3]= {9.143 Second,Null} In[4]:= Table[RandomArray[NormalDistribution[50,25], 1],{10^6}];//Timing Out[4]= {91.361 Second,Null} In[5]:= Timing[Table[Random[Real, {0, 100}],{10^6}];] Out[5]= {0.671 Second,Null} as Daniel was saying, no secret here. as you increase the number of samples, you increase the time. it take 10x time to pick 10x more numbers. (lol, to think I was wondering why that was...) anyway, obviously using the RandomArray to sample the number from a distribution is less efficient than just using Random[]. but that's not what is kinda peculiar for me. now if we do the same thing with the RandomReplacement loaded... In[1]:= <<"Statistics`NormalDistribution`"; <<RandomReplacement` In[3]:= Table[RandomArray[NormalDistribution[50,25], 1],{10^4}];//Timing Out[3]= {1.221 Second,Null} In[4]:= Table[RandomArray[NormalDistribution[50,25], 1],{10^5}];//Timing Out[4]= {11.867 Second,Null} In[5]:= Table[RandomArray[NormalDistribution[50,25], 1],{10^6}];//Timing Out[5]= {120.953 Second,Null} In[6]:=Timing[Table[Random[Real, {0, 100}],{10^6}];] Out[6]={27.389 Second,Null} so... the differences in timing between loading and not loading the RandomReplacement package is more pronounced in using Random[Real, {}] alone. When I used RandomArray, package slowed things down 1.3 times. (120.953 seconds/91.361 seconds) but using Random[Real], it slows it down by 40 times!( 27.389 seconds /0.671seconds) I wonder if this gives you guys soem ideas to speed things up. hope that helps a bit. sean __________________________________ Do you Yahoo!? Yahoo! Mail - You care about security. So do we. http://promotions.yahoo.com/new_mail