carloslesmes — Sep 10, 2012, 3:30 PM
# uso de algunas funciones del paquete "prob"
# técnicas de conteo
suppressPackageStartupMessages(library(prob)) # cargar el paquete 'prob'
nsamp(n=10, k=5, replace = TRUE, ordered = TRUE) # muestra ordenada con reemplazo
[1] 1e+05
nsamp(n=10, k=5, replace = FALSE, ordered = TRUE) # muestra ordenada sin reemplazo
[1] 30240
nsamp(n=10, k=5, replace = TRUE, ordered = FALSE) # muestra sin orden con reemplazo
[1] 2002
nsamp(n=10, k=5, replace = FALSE, ordered = FALSE) # muestra sin orden sin reemplazo
[1] 252
# aqui están las muestras
urnsamples(1:4, size = 2,ordered=TRUE,replace=TRUE)
X1 X2
1 1 1
2 2 1
3 3 1
4 4 1
5 1 2
6 2 2
7 3 2
8 4 2
9 1 3
10 2 3
11 3 3
12 4 3
13 1 4
14 2 4
15 3 4
16 4 4
urnsamples(1:4, size = 2,ordered=TRUE,replace=FALSE)
Loading required package: combinat
Attaching package: 'combinat'
The following object(s) are masked from 'package:utils':
combn
X1 X2
1 1 2
2 2 1
3 1 3
4 3 1
5 1 4
6 4 1
7 2 3
8 3 2
9 2 4
10 4 2
11 3 4
12 4 3
urnsamples(1:4, size = 2,ordered=FALSE,replace=TRUE)
X1 X2
1 1 1
2 1 2
3 1 3
4 1 4
5 2 2
6 2 3
7 2 4
8 3 3
9 3 4
10 4 4
urnsamples(1:4, size = 2,ordered=FALSE,replace=FALSE)
X1 X2
1 1 2
2 1 3
3 1 4
4 2 3
5 2 4
6 3 4
cards(makespace = TRUE)
rank suit probs
1 2 Club 0.01923
2 3 Club 0.01923
3 4 Club 0.01923
4 5 Club 0.01923
5 6 Club 0.01923
6 7 Club 0.01923
7 8 Club 0.01923
8 9 Club 0.01923
9 10 Club 0.01923
10 J Club 0.01923
11 Q Club 0.01923
12 K Club 0.01923
13 A Club 0.01923
14 2 Diamond 0.01923
15 3 Diamond 0.01923
16 4 Diamond 0.01923
17 5 Diamond 0.01923
18 6 Diamond 0.01923
19 7 Diamond 0.01923
20 8 Diamond 0.01923
21 9 Diamond 0.01923
22 10 Diamond 0.01923
23 J Diamond 0.01923
24 Q Diamond 0.01923
25 K Diamond 0.01923
26 A Diamond 0.01923
27 2 Heart 0.01923
28 3 Heart 0.01923
29 4 Heart 0.01923
30 5 Heart 0.01923
31 6 Heart 0.01923
32 7 Heart 0.01923
33 8 Heart 0.01923
34 9 Heart 0.01923
35 10 Heart 0.01923
36 J Heart 0.01923
37 Q Heart 0.01923
38 K Heart 0.01923
39 A Heart 0.01923
40 2 Spade 0.01923
41 3 Spade 0.01923
42 4 Spade 0.01923
43 5 Spade 0.01923
44 6 Spade 0.01923
45 7 Spade 0.01923
46 8 Spade 0.01923
47 9 Spade 0.01923
48 10 Spade 0.01923
49 J Spade 0.01923
50 Q Spade 0.01923
51 K Spade 0.01923
52 A Spade 0.01923
rolldie(1) # lanzar un dado corriente
X1
1 1
2 2
3 3
4 4
5 5
6 6
rolldie(2,nsides=3) # lanzar dos dados de tres caras
X1 X2
1 1 1
2 2 1
3 3 1
4 1 2
5 2 2
6 3 2
7 1 3
8 2 3
9 3 3
rolldie(3, nsides = 4) # lanzar tres dados de cuatro caras
X1 X2 X3
1 1 1 1
2 2 1 1
3 3 1 1
4 4 1 1
5 1 2 1
6 2 2 1
7 3 2 1
8 4 2 1
9 1 3 1
10 2 3 1
11 3 3 1
12 4 3 1
13 1 4 1
14 2 4 1
15 3 4 1
16 4 4 1
17 1 1 2
18 2 1 2
19 3 1 2
20 4 1 2
21 1 2 2
22 2 2 2
23 3 2 2
24 4 2 2
25 1 3 2
26 2 3 2
27 3 3 2
28 4 3 2
29 1 4 2
30 2 4 2
31 3 4 2
32 4 4 2
33 1 1 3
34 2 1 3
35 3 1 3
36 4 1 3
37 1 2 3
38 2 2 3
39 3 2 3
40 4 2 3
41 1 3 3
42 2 3 3
43 3 3 3
44 4 3 3
45 1 4 3
46 2 4 3
47 3 4 3
48 4 4 3
49 1 1 4
50 2 1 4
51 3 1 4
52 4 1 4
53 1 2 4
54 2 2 4
55 3 2 4
56 4 2 4
57 1 3 4
58 2 3 4
59 3 3 4
60 4 3 4
61 1 4 4
62 2 4 4
63 3 4 4
64 4 4 4
tosscoin(4, makespace = TRUE)
toss1 toss2 toss3 toss4 probs
1 H H H H 0.0625
2 T H H H 0.0625
3 H T H H 0.0625
4 T T H H 0.0625
5 H H T H 0.0625
6 T H T H 0.0625
7 H T T H 0.0625
8 T T T H 0.0625
9 H H H T 0.0625
10 T H H T 0.0625
11 H T H T 0.0625
12 T T H T 0.0625
13 H H T T 0.0625
14 T H T T 0.0625
15 H T T T 0.0625
16 T T T T 0.0625
R <- rolldie(4)
A <- probspace(R)
sim(A,ntrials=3)
Warning: Walker's alias method used: results are different from R < 2.2.0
X1 X2 X3 X4
1 6 5 4 3
2 3 5 5 5
3 4 4 5 6
citation("prob")
To cite package 'prob' in publications use:
G. Jay Kerns (2009). prob: Elementary Probability on Finite
Sample Spaces. R package version 0.9-2.
http://CRAN.R-project.org/package=prob
A BibTeX entry for LaTeX users is
@Manual{,
title = {prob: Elementary Probability on Finite Sample Spaces},
author = {G. Jay Kerns},
year = {2009},
note = {R package version 0.9-2},
url = {http://CRAN.R-project.org/package=prob},
}
ATTENTION: This citation information has been auto-generated from
the package DESCRIPTION file and may need manual editing, see
'help("citation")' .