probability - Simulate experiment with R Programming -

-

i new in r programming. need solve 1 problem in r. need simulate following experiment in r. poker hand consists of 5 cards dealt conventional pack of 52 cards, order of cards not being important. find probability given hand has @ least 1 king , @ least 1 queen.

i know how find atleast 1 king not @ least 1 king , @ least 1 queen. atleast 1 king code : 

deck<- rep(1:13 , each=4) #here j=11 , q=12, k=13 nhands <- 1000 xk<- c(rep(-1, nhands)) for( in 1:nhands){ hand <- sample( deck , 5 , replace= false) numberofk<-0 for( j in 1:5){ # count kings if( hand[j] == 13){ numberofk <- numberofk +1 } } #print(numberofk) xk[i] <-numberofk #print(hand) } table(xk) /nhands

can please me in coding required 1.. thanks

the probability of hand of 5 cards contain @ least 1 king or 1 queen can written following, sample(deck, 5) gives hand of 5 cards while any(c(12, 13) ...) checks whether king or queen within hand , sum counts how many times such case happens within 1000 simulation:

set.seed(10) sum(sapply(1:100000, function(i) { any(c(12, 13) %in% sample(deck, 5))}))/100000 # [1] 0.58365

theoretically, probability of such case be:

(choose(52, 5) - choose(44, 5))/choose(52, 5) # [1] 0.5821375

which pretty close.

and on other hand, if indeed and means @ least king , queen, simulation gives:

set.seed(10) sum(sapply(1:100000, function(i) { all(c(12, 13) %in% sample(deck, 5))}))/100000 # [1] 0.09932

and theoretically:

(choose(52, 5) - choose(44, 5) - (2*(choose(48, 5) - choose(44, 5))))/choose(52, 5) # [1] 0.1001785

and number matches closely.


Comments

Post a Comment

Popular posts from this blog

matlab - error with cyclic autocorrelation function -

-
Image
i'm trying code cyclic autocorrelation function in matlab follows: t=0:(n-1); t=t*te; i_alpha=0; tau=-n2*te:te:n2*te; alpha=-1/2:1/n:1/2; ryy_cl=zeros(length(alpha),length(tau)); alpha=alpha/te ind_tau=0; i_alpha=i_alpha+1; k=tau ind_tau=ind_tau+1; if k>0 % ryy_cl(i_alpha,ind_tau)=1/length(sig_bin_syl_mod(1:end- k))*sum((sig_bin_syl_mod(1:end-k)).*sig_bin_syl_mod(k+1:end)).*exp(1i*2*pi*alpha*t(1+k:length(sig_bin_syl_mod(1:end)))); ryy_cl(i_alpha,ind_tau)=(1/n)*sum((sig_bin_syl_mod(1:end-k)).*sig_bin_syl_mod(k+1:end)).*exp(-1i*2*pi*alpha*t(1:end-k)); else ryy_cl(i_alpha,ind_tau)=(1/n)*sum((sig_bin_syl_mod(1-k:end)).*sig_bin_syl_mod(1:end+k)).*exp(-1i*2*pi*alpha*t(1-k:end)); end end end i'm getting errors, , doesn't show expecting. according below formulae, how can fix it? the code provided incomplete, difficult confidently reproduce error. in particular, haven't de...
Read more

django - (fields.E300) Field defines a relation with model 'AbstractEmailUser' which is either not installed, or is abstract -

-
i trying create customuser django project email username , add radio button truck , company. in registration process email-id registered per truck or company. mentioned radio button 'tag' , add manytomany field tag in emailuser model. when makemigrations, raising error : (fields.e300) field defines relation model 'abstractemailuser' either not installed, or abstract. i quite new django , not sure whether have created correct code wants. please me solving this. here code, models.py: import django django.contrib.auth.models import ( abstractbaseuser, baseusermanager, permissionsmixin) django.core.mail import send_mail django.db import models django.utils import timezone django.utils.translation import ugettext_lazy _ class emailusermanager(baseusermanager): """custom manager emailuser.""" def _create_user(self, email, password, is_staff, is_superuser, **extra_fields): """create , save emailus...
Read more

c# - What is a good .Net RefEdit control to use with ExcelDna? -

-
i using exceldna , c# build suite of tools. part of forms i'd want user select cells/ranges workbook. i've got application.inputbox route working...but doesn't cool... i haven't read many things excel's refedit control, , see many have in past written workarounds, in (shudder) vba or vb.net. however, none of these seem recent... best approach/control refedit functionality on form run via exceldna (if applicable)? the refedit sample project of github explores different ways of dealing threading when showing form excel add-in, main issue making refedit control .net add-in. check page links various other implementations: https://github.com/excel-dna/exceldna/wiki/links-about-refedit
Read more