Re: How to simulate random samples from crooked coin toss?

*To*: mathgroup at smc.vnet.net*Subject*: [mg58517] Re: How to simulate random samples from crooked coin toss?*From*: "Jens-Peer Kuska" <kuska at informatik.uni-leipzig.de>*Date*: Tue, 5 Jul 2005 04:52:25 -0400 (EDT)*Organization*: Uni Leipzig*References*: <dad7vt$41$1@smc.vnet.net>*Sender*: owner-wri-mathgroup at wolfram.com

Hi, say a 1 mean a head and a 0 mean a tail than HeadOrTail[]:=If[Random[]<0.6,1,0] will generate one realisation of your random variable and with sample = Table[HeadOrTail[], {1000}]; you get a sample of the random variable, here is the number of "heads" in the sample Count[sample, 1] and here the number of tails Count[sample, 0] Regards Jens "DAN BIKLE" <dbikle at gmail.com> schrieb im Newsbeitrag news:dad7vt$41$1 at smc.vnet.net... > Math People, > > I'm new to Mathematica so please bear with me. > > I have a simple problem I want to solve and I'm > looking > for some ideas to get started on it. > > I want to simulate 100 tosses of a crooked coin > where > probability of Heads is 0.6 and Tails is 0.4. > > After I figure the coin out, I want to move on > to crooked > dice and crooked card decks. > > So, really what I'm looking for has more to do > with developing > a method to simulate the collection of random > samples. > > I looked around in Help and found a promising > sentence: > > "Random[distribution] gives a random number with > the specified > statistical distribution." > > Then I found a page here which has some clues: > http://mathworld.wolfram.com/DiscreteDistribution.html > > Unfortunately, I'm weak with math notation so > the above > page makes little sense to me. > > So, let me put it to you as simply as I can: > > I want to define a function that could be called > like this: > > numberOfHeads = coinTosser[infoAboutCoin, > numberOfTosses, head] > numberOfTails = coinTosser[infoAboutCoin, > numberOfTosses, tail] > > Once I figure out how to write the above > function, > I could template off of it to build something > like this: > > numberOfOnes = diceRoller[infoAboutThisDie, > numberOfRolls, aOne] > > Do any of you have thoughts on how to approach > the > writing of the above functions? > > -Dan > -- > dbikle at gmail.com >