MUESTREO PAPA
Ana María Florián Pulido
8/5/2021
library(clhs)
set.seed(1013687724)
xy_=expand.grid(x=seq(1,60,5),y=seq(1,42,3.5))
sana_=sample(0:1,144,T,prob = c(0.8,0.2))
plot(xy_, pch=15, cex=2, col=sana_+3)
Muestreo aleatorio simple (para estimar una proporción)
\[tamaño~de~muestra~(n)=~ \frac{Np(1-p)}{(N-1)(e/z)^2+p(1-p)}\] \(N=140\) tamaño de la población
\(p_1 = 2/12 ; \ p_2 = 50%\) proporción
\(e<10%\) error de muestreo
$z=1.96 $ valor en la curva normal para un nivel de confianza (\(\alpha=0.95\))
n_aleatorio_=function(N,p,e,z=1.96){
n_=(N*p*(1-p))/(((N-1)*(e/z)^2)+(p*(1-p)))
return(ceiling(n_))
}
n_aleatorio_(N=144,p=0.5,e=0.1)## [1] 58muestra_1_=sample(seq(144),n_aleatorio_(N=144,p=0.5,e=0.1), replace = FALSE)plot(xy_,pch=15,cex=2,col=sana_+3)
points(xy_[muestra_1_,1],xy_[muestra_1_,2],col=3,pch=16)
prev1_=100*table(sana_[muestra_1_])/58;prev1_##
## 0 1
## 75.86207 24.13793prev_real_=100*table(sana_)/144; prev_real_## sana_
## 0 1
## 83.33333 16.66667e_mues_=(abs(prev1_-prev_real_)/prev_real_); e_mues_##
## 0 1
## 0.08965517 0.44827586muestra_2_=sample(seq(144),n_aleatorio_(N=144,p=0.5,e=0.05), replace = FALSE)
plot(xy_,pch=15,cex=2,col=sana_+3)
points(xy_[muestra_2_,1],xy_[muestra_2_,2],col=3,pch=16)
length(muestra_2_)## [1] 105100*table(sana_[muestra_2_])/length(muestra_2_)##
## 0 1
## 83.80952 16.19048muestra_3_ = sample(seq(144),n_aleatorio_(N=144,p=2/12,e=0.1), replace = FALSE)
plot(xy_,pch=15,cex=2,col=sana_+3)
points(xy_[muestra_3_,1],xy_[muestra_3_,2],col=3,pch=16)
length(muestra_3_)## [1] 40100*table(sana_[muestra_3_])/length(muestra_3_)##
## 0 1
## 75 25muestra_4_ = sample(seq(144),n_aleatorio_(N=144,p=2/12,e=0.05), replace = FALSE)
plot(xy_,pch=15,cex=2,col=sana_+3)
points(xy_[muestra_4_,1],xy_[muestra_4_,2],col=3,pch=16)
length(muestra_4_)## [1] 87100*table(sana_[muestra_4_])/length(muestra_4_)##
## 0 1
## 79.31034 20.68966library(clhs)
res_=NULL
for (s in 10:108) {
print(s)
pre_res_ = clhs(xy_,size = s, iter = 1000, progress = FALSE)
res_[[s-9]] = pre_res_
}## [1] 10
## [1] 11
## [1] 12
## [1] 13
## [1] 14
## [1] 15
## [1] 16
## [1] 17
## [1] 18
## [1] 19
## [1] 20
## [1] 21
## [1] 22
## [1] 23
## [1] 24
## [1] 25
## [1] 26
## [1] 27
## [1] 28
## [1] 29
## [1] 30
## [1] 31
## [1] 32
## [1] 33
## [1] 34
## [1] 35
## [1] 36
## [1] 37
## [1] 38
## [1] 39
## [1] 40
## [1] 41
## [1] 42
## [1] 43
## [1] 44
## [1] 45
## [1] 46
## [1] 47
## [1] 48
## [1] 49
## [1] 50
## [1] 51
## [1] 52
## [1] 53
## [1] 54
## [1] 55
## [1] 56
## [1] 57
## [1] 58
## [1] 59
## [1] 60
## [1] 61
## [1] 62
## [1] 63
## [1] 64
## [1] 65
## [1] 66
## [1] 67
## [1] 68
## [1] 69
## [1] 70
## [1] 71
## [1] 72
## [1] 73
## [1] 74
## [1] 75
## [1] 76
## [1] 77
## [1] 78
## [1] 79
## [1] 80
## [1] 81
## [1] 82
## [1] 83
## [1] 84
## [1] 85
## [1] 86
## [1] 87
## [1] 88
## [1] 89
## [1] 90
## [1] 91
## [1] 92
## [1] 93
## [1] 94
## [1] 95
## [1] 96
## [1] 97
## [1] 98
## [1] 99
## [1] 100
## [1] 101
## [1] 102
## [1] 103
## [1] 104
## [1] 105
## [1] 106
## [1] 107
## [1] 108res_## [[1]]
## [1] 69 46 60 4 30 143 99 73 20 125
##
## [[2]]
## [1] 118 79 28 12 143 44 51 122 65 102 95
##
## [[3]]
## [1] 31 8 132 88 111 69 42 133 107 82 26 53
##
## [[4]]
## [1] 61 34 5 105 19 54 74 131 92 111 139 48 64
##
## [[5]]
## [1] 122 91 83 28 140 54 109 17 48 2 81 99 70 117
##
## [[6]]
## [1] 27 37 96 58 134 17 71 7 80 114 93 130 101 136 125
##
## [[7]]
## [1] 128 16 89 36 71 82 21 133 110 2 98 54 144 51 91 38
##
## [[8]]
## [1] 9 133 53 120 15 131 94 32 100 45 62 78 123 115 19 29 140
##
## [[9]]
## [1] 12 27 104 117 38 118 139 64 11 73 92 55 128 53 18 34 50 114
##
## [[10]]
## [1] 1 112 91 140 57 41 132 74 58 97 64 23 52 6 53 27 15 133 85
##
## [[11]]
## [1] 1 24 131 54 29 136 122 109 91 82 108 90 38 63 77 137 93 52 8
## [20] 68
##
## [[12]]
## [1] 14 48 55 142 131 116 25 122 135 100 93 54 74 60 70 10 40 56 72
## [20] 65 64
##
## [[13]]
## [1] 89 140 56 21 35 46 15 7 120 133 118 62 102 80 75 70 28 20 96
## [20] 129 72 103
##
## [[14]]
## [1] 124 61 42 82 22 26 104 59 131 105 47 96 123 6 117 137 39 132 134
## [20] 72 7 75 94
##
## [[15]]
## [1] 121 6 20 110 90 82 39 131 69 50 100 137 103 108 104 84 65 30 94
## [20] 73 140 66 102 59
##
## [[16]]
## [1] 133 12 66 123 95 120 38 35 116 130 23 52 37 142 51 100 76 57 29
## [20] 41 77 54 61 139 128
##
## [[17]]
## [1] 116 22 12 85 83 63 84 38 97 18 127 26 21 123 140 56 96 14 128
## [20] 124 36 49 28 5 130 66
##
## [[18]]
## [1] 37 121 5 66 32 86 136 63 44 15 8 120 51 35 140 69 75 70 115
## [20] 92 40 126 87 98 127 38 84
##
## [[19]]
## [1] 131 121 106 139 18 13 134 22 16 56 109 95 29 71 84 37 15 50 113
## [20] 21 74 98 3 55 100 6 59 66
##
## [[20]]
## [1] 15 143 133 26 51 46 35 126 91 43 36 23 61 8 38 93 120 56 84
## [20] 101 74 14 64 75 67 86 50 68 53
##
## [[21]]
## [1] 134 25 121 3 88 36 13 43 90 82 61 119 102 110 111 140 48 50 5
## [20] 31 104 17 89 23 38 105 62 70 15 69
##
## [[22]]
## [1] 12 123 144 71 126 104 29 94 61 14 15 69 91 97 125 57 66 19 139
## [20] 40 83 117 99 50 89 128 45 18 101 109 7
##
## [[23]]
## [1] 117 14 46 135 28 3 36 110 80 4 29 87 107 116 91 38 67 66 51
## [20] 21 8 72 84 13 6 31 97 18 125 27 101 11
##
## [[24]]
## [1] 19 124 140 27 52 131 45 109 127 94 97 5 125 77 61 86 24 91 132
## [20] 78 55 136 101 41 66 112 69 105 18 51 63 64 68
##
## [[25]]
## [1] 46 84 18 130 81 36 123 73 16 10 138 67 95 110 129 50 68 13 88
## [20] 115 117 17 104 85 106 94 71 96 114 31 11 66 72 39
##
## [[26]]
## [1] 143 31 1 84 87 46 138 111 122 114 33 106 52 101 104 29 45 3 20
## [20] 108 137 34 71 124 88 73 99 69 141 62 55 92 32 40 43
##
## [[27]]
## [1] 21 93 79 111 142 41 119 38 135 59 50 123 58 144 121 112 25 132 105
## [20] 22 33 18 5 9 118 32 34 72 47 31 30 94 62 29 128 137
##
## [[28]]
## [1] 142 37 84 133 29 18 33 45 143 134 99 32 79 21 22 113 118 136 49
## [20] 15 71 14 117 6 119 24 131 23 100 60 27 87 126 38 106 97 5
##
## [[29]]
## [1] 144 26 117 133 16 34 139 106 39 105 95 128 21 8 11 1 125 126 140
## [20] 4 10 109 65 98 124 61 74 18 68 102 91 114 59 71 29 78 55 63
##
## [[30]]
## [1] 109 144 23 135 51 103 1 4 50 27 35 38 143 29 90 68 22 108 45
## [20] 121 7 105 76 2 30 139 11 83 114 106 36 12 61 138 79 130 42 87
## [39] 47
##
## [[31]]
## [1] 8 135 4 96 31 106 25 111 18 128 129 68 82 80 102 118 24 11 48
## [20] 99 2 38 60 74 53 54 21 141 12 26 57 139 130 115 107 89 69 64
## [39] 67 110
##
## [[32]]
## [1] 66 111 43 131 133 63 144 35 122 41 86 71 101 60 59 48 142 120 95
## [20] 33 6 46 76 13 44 138 134 132 53 103 114 55 118 107 85 130 141 8
## [39] 31 116 42
##
## [[33]]
## [1] 115 12 136 1 95 142 112 72 71 38 88 15 17 137 85 40 123 51 68
## [20] 55 5 127 124 118 83 63 125 66 30 100 93 48 75 7 81 42 110 89
## [39] 104 49 39 20
##
## [[34]]
## [1] 35 26 25 98 2 76 32 133 119 65 45 15 137 75 60 48 69 38 9
## [20] 96 125 114 108 71 92 77 115 70 112 123 93 46 66 54 63 127 34 27
## [39] 95 100 1 56 74
##
## [[35]]
## [1] 133 89 52 11 84 2 26 128 135 21 73 141 70 91 88 35 69 140 87
## [20] 127 118 54 29 18 27 7 49 105 126 40 109 80 114 28 50 71 30 115
## [39] 81 5 113 43 124 78
##
## [[36]]
## [1] 128 11 50 123 25 99 85 134 23 33 109 129 91 4 69 28 106 35 89
## [20] 17 6 105 46 84 18 62 39 13 57 96 112 139 74 63 15 14 55 125
## [39] 67 81 29 115 53 68 80
##
## [[37]]
## [1] 49 41 11 77 124 125 32 117 81 84 13 144 106 27 132 92 94 33 122
## [20] 102 44 90 118 74 47 111 34 16 138 39 17 75 8 127 71 45 36 9
## [39] 141 10 88 52 57 55 29 42
##
## [[38]]
## [1] 12 131 62 121 106 1 11 9 114 49 82 6 14 77 17 141 75 83 101
## [20] 81 10 90 29 55 21 42 113 63 115 3 53 7 68 74 136 5 28 34
## [39] 125 140 95 19 66 26 100 59 69
##
## [[39]]
## [1] 131 25 58 23 120 130 19 36 30 101 12 59 87 135 69 37 3 68 49
## [20] 141 74 76 107 22 108 48 91 81 105 124 112 67 78 98 7 39 17 116
## [39] 71 110 106 83 90 8 43 127 118 70
##
## [[40]]
## [1] 54 48 5 75 60 56 120 107 24 70 81 130 76 39 111 8 3 141 31
## [20] 67 34 105 116 85 100 29 137 71 38 32 101 93 64 89 82 6 20 88
## [39] 30 52 123 10 72 99 17 65 83 45 13
##
## [[41]]
## [1] 45 124 32 128 135 59 119 68 74 5 121 17 18 48 84 34 7 15 115
## [20] 81 108 143 116 20 25 14 44 60 139 37 23 113 76 12 93 2 54 140
## [39] 40 91 77 82 100 96 87 22 137 9 38 69
##
## [[42]]
## [1] 26 136 142 71 107 98 80 43 131 58 124 135 13 97 86 99 56 15 39
## [20] 121 25 35 46 29 100 90 60 54 45 112 74 88 5 144 140 111 81 18
## [39] 16 78 44 122 66 106 95 93 55 110 6 113 67
##
## [[43]]
## [1] 11 5 99 133 48 33 143 26 31 49 112 13 34 85 126 80 9 117 32
## [20] 135 105 42 78 118 43 62 139 44 114 142 116 23 141 115 50 52 72 93
## [39] 28 46 53 92 75 58 81 66 106 127 40 102 56 84
##
## [[44]]
## [1] 66 144 23 26 50 34 46 36 121 142 22 130 123 54 86 132 114 89 32
## [20] 10 49 39 76 47 21 14 67 129 138 131 125 128 110 143 42 102 97 79
## [39] 94 92 118 96 4 60 44 72 30 139 109 115 64 77 82
##
## [[45]]
## [1] 40 25 15 1 93 135 144 83 36 48 120 121 11 88 73 13 68 89 43
## [20] 122 50 96 24 116 74 107 57 14 110 29 85 129 32 63 82 53 112 113
## [39] 18 76 77 141 123 91 75 31 102 100 33 56 7 84 114 92
##
## [[46]]
## [1] 39 120 24 107 21 51 58 102 29 40 3 136 1 74 36 106 93 27 88
## [20] 137 122 18 37 109 117 52 128 50 139 35 6 138 65 134 10 4 130 92
## [39] 127 100 104 9 78 90 63 116 34 57 129 79 75 19 81 55 76
##
## [[47]]
## [1] 87 111 13 83 26 23 130 9 121 64 113 60 37 35 135 74 136 110 39
## [20] 138 27 55 119 52 59 54 36 10 107 142 104 102 88 95 22 50 131 106
## [39] 49 62 44 116 71 14 144 8 57 45 76 78 7 32 43 101 61 84
##
## [[48]]
## [1] 107 103 1 84 61 111 7 93 9 51 73 110 44 141 62 139 39 88 13
## [20] 74 142 50 57 49 23 40 33 24 2 105 90 70 6 48 69 80 140 59
## [39] 95 66 137 122 72 65 31 81 56 114 128 99 26 67 71 133 53 124 94
##
## [[49]]
## [1] 144 24 115 119 65 27 39 131 111 58 122 21 109 14 8 110 121 67 96
## [20] 51 43 22 126 71 63 64 50 15 112 37 53 129 31 68 26 69 91 40
## [39] 66 3 41 33 70 92 9 79 23 128 59 100 89 87 38 7 137 13 139
## [58] 47
##
## [[50]]
## [1] 132 48 135 121 1 25 125 28 83 112 21 111 69 4 36 119 114 128 7
## [20] 72 9 139 88 3 24 46 100 92 2 106 136 61 56 85 27 78 123 93
## [39] 65 8 31 51 74 143 80 67 29 115 13 6 63 62 87 81 71 50 73
## [58] 12 16
##
## [[51]]
## [1] 38 37 48 83 132 3 108 128 28 90 109 65 13 118 24 137 53 17 52
## [20] 34 140 33 123 117 60 42 25 85 71 94 77 73 50 11 46 115 23 105
## [39] 21 32 36 76 136 40 119 5 84 107 10 104 78 56 75 79 91 103 70
## [58] 12 129 41
##
## [[52]]
## [1] 15 126 122 47 39 1 134 97 132 89 71 13 70 45 100 143 46 31 124
## [20] 75 66 24 50 29 112 8 91 44 110 49 27 74 56 78 138 142 28 37
## [39] 103 131 20 19 58 121 67 129 73 6 128 62 85 7 53 76 92 5 52
## [58] 18 51 99 48
##
## [[53]]
## [1] 144 32 25 83 132 133 108 114 110 45 40 24 67 87 7 53 86 16 47
## [20] 91 18 103 69 107 35 31 37 135 113 14 43 68 51 109 127 61 119 74
## [39] 112 79 138 54 57 62 27 65 89 71 123 36 30 6 94 10 41 118 72
## [58] 120 115 50 104 56
##
## [[54]]
## [1] 133 132 5 60 128 13 121 144 17 112 75 114 124 141 12 80 38 47 23
## [20] 81 24 63 37 35 131 95 92 125 36 96 29 11 9 16 109 120 78 73
## [39] 118 72 134 111 40 20 53 103 105 79 123 43 108 2 7 25 31 70 119
## [58] 66 94 101 117 107 39
##
## [[55]]
## [1] 133 136 144 9 36 89 51 2 53 21 25 63 44 90 45 7 121 11 73
## [20] 13 29 58 124 72 107 135 65 70 139 78 98 32 54 113 142 17 20 120
## [39] 99 19 50 132 30 86 49 109 115 28 3 138 52 26 82 101 33 129 74
## [58] 92 75 61 87 41 39 69
##
## [[56]]
## [1] 23 128 12 35 127 13 113 118 51 125 137 39 142 68 80 57 124 141 144
## [20] 102 11 37 67 106 38 97 99 28 54 17 5 9 41 96 92 130 90 21
## [39] 43 36 30 7 31 104 101 114 58 62 24 89 71 75 85 8 6 103 26
## [58] 95 140 112 61 65 116 44 94
##
## [[57]]
## [1] 119 39 132 107 12 1 51 60 95 63 43 59 81 83 134 109 24 13 45
## [20] 108 125 25 14 79 139 22 82 11 137 48 130 64 97 121 19 15 10 141
## [39] 44 96 76 120 61 42 116 88 102 126 138 52 65 89 74 47 78 127 85
## [58] 5 7 117 35 115 17 66 31 80
##
## [[58]]
## [1] 22 11 48 142 95 105 65 35 108 121 137 135 102 15 38 111 8 141 136
## [20] 60 110 84 43 99 96 92 124 37 143 62 5 14 27 41 118 34 69 107
## [39] 128 63 46 55 131 7 71 66 72 115 28 127 45 49 31 64 20 33 87
## [58] 1 78 98 103 16 36 93 9 80 94
##
## [[59]]
## [1] 93 55 35 141 96 17 3 133 49 24 106 4 8 60 50 94 74 143 115
## [20] 13 80 142 98 40 132 58 63 15 113 88 20 135 66 37 109 108 62 81
## [39] 54 91 89 99 39 34 114 36 130 138 10 97 41 124 38 27 128 104 125
## [58] 126 48 19 22 7 101 31 77 33 32 140
##
## [[60]]
## [1] 119 1 44 30 142 106 87 22 16 83 14 13 135 105 2 33 127 38 71
## [20] 58 117 144 46 97 85 100 112 8 63 137 59 39 35 109 60 132 73 53
## [39] 64 3 54 111 10 11 5 126 21 102 107 98 122 103 49 123 116 124 86
## [58] 69 138 79 68 75 125 143 114 20 104 18 34
##
## [[61]]
## [1] 54 30 3 136 87 103 75 22 138 81 2 35 45 96 137 53 129 135 84
## [20] 68 110 143 59 73 114 109 29 12 144 13 141 82 80 40 86 117 11 15
## [39] 39 106 23 108 77 92 93 24 118 31 125 128 58 50 88 98 94 55 111
## [58] 34 140 37 89 4 69 72 104 78 9 41 90 65
##
## [[62]]
## [1] 130 34 82 28 111 27 96 109 93 68 6 94 57 60 73 97 112 132 65
## [20] 144 133 80 91 43 69 45 24 89 95 40 56 119 53 131 42 86 62 99
## [39] 110 71 88 3 38 107 36 142 137 143 16 67 10 33 108 100 72 70 83
## [58] 101 11 63 59 128 81 77 15 76 20 79 114 84 19
##
## [[63]]
## [1] 99 123 131 128 28 36 121 144 53 71 13 84 72 133 118 138 2 95 119
## [20] 117 115 97 56 23 54 24 93 114 113 109 120 12 15 31 108 81 89 9
## [39] 62 58 90 30 78 67 49 116 140 126 106 136 40 48 83 85 41 64 14
## [58] 52 94 17 20 80 37 105 66 107 60 98 45 63 32 33
##
## [[64]]
## [1] 64 118 48 107 45 53 128 5 124 108 16 22 142 24 69 40 96 110 139
## [20] 23 104 39 72 123 46 95 141 73 10 67 36 27 103 97 42 140 134 120
## [39] 127 88 34 3 57 91 115 38 8 30 74 78 59 89 55 76 37 135 9
## [58] 94 93 130 131 87 75 50 60 56 83 92 113 7 2 101 15
##
## [[65]]
## [1] 142 70 35 31 1 63 121 46 16 3 18 137 50 82 109 133 8 111 27
## [20] 41 67 54 68 10 139 101 51 65 136 115 22 144 84 24 79 52 74 7
## [39] 120 66 37 85 126 60 32 15 77 40 6 86 129 96 116 19 83 91 36
## [58] 39 125 81 17 75 72 69 5 73 141 53 104 62 103 92 124 80
##
## [[66]]
## [1] 11 118 55 15 43 133 13 36 119 18 1 93 102 134 140 66 76 117 132
## [20] 6 16 127 27 30 24 126 46 90 143 136 109 95 111 108 47 37 10 40
## [39] 75 9 104 77 70 74 8 73 96 105 17 121 4 92 45 79 60 14 68
## [58] 53 82 87 80 128 33 88 31 64 120 32 124 112 49 51 83 135 94
##
## [[67]]
## [1] 124 81 138 71 70 14 45 129 134 87 35 16 5 118 32 121 116 20 115
## [20] 140 100 27 12 91 37 73 19 33 53 25 96 142 21 99 125 143 11 2
## [39] 79 88 123 41 66 60 22 47 43 15 82 92 56 39 90 28 7 31 102
## [58] 44 127 38 120 59 18 36 69 51 139 78 130 84 34 55 68 8 57 72
##
## [[68]]
## [1] 22 83 122 143 25 77 23 11 69 37 108 109 65 54 118 134 113 72 74
## [20] 94 4 93 24 124 32 129 98 141 90 50 20 14 57 96 84 19 99 68
## [39] 102 89 47 13 21 76 95 138 34 73 125 3 46 139 117 88 126 71 49
## [58] 86 87 131 5 15 103 33 130 43 31 16 120 53 70 40 12 66 10 17
## [77] 58
##
## [[69]]
## [1] 131 13 132 12 93 140 107 120 113 29 117 35 110 121 134 6 28 54 83
## [20] 59 129 94 67 86 39 15 76 65 24 115 47 112 81 122 34 5 60 139
## [39] 44 50 88 130 4 58 124 114 57 43 106 87 91 96 105 68 16 48 56
## [58] 20 18 41 101 92 72 85 74 80 79 19 70 26 73 9 46 61 71 141
## [77] 8 143
##
## [[70]]
## [1] 144 4 30 59 109 79 51 11 2 99 134 140 25 106 17 94 131 117 63
## [20] 12 74 40 9 6 65 58 3 101 110 75 123 116 28 49 95 60 136 20
## [39] 47 132 142 32 26 77 113 27 89 120 72 8 87 70 73 81 127 57 5
## [58] 88 52 21 102 97 114 71 91 85 34 78 10 61 128 139 111 64 112 76
## [77] 141 84 33
##
## [[71]]
## [1] 7 22 9 122 103 91 45 50 47 52 141 81 14 27 23 107 72 37 105
## [20] 87 135 98 51 133 4 79 130 116 119 73 132 127 17 59 86 70 136 54
## [39] 64 28 25 112 21 57 19 43 35 142 55 67 71 18 8 69 76 39 61
## [58] 42 41 114 138 78 121 106 102 32 65 49 90 62 109 108 6 85 95 40
## [77] 118 100 92 16
##
## [[72]]
## [1] 101 92 115 136 16 140 121 142 124 9 132 54 128 133 93 5 38 130 67
## [20] 112 66 134 113 110 41 91 39 53 71 46 56 72 137 57 119 59 104 76
## [39] 94 68 30 69 14 64 85 28 34 86 89 126 90 60 37 109 20 62 141
## [58] 31 51 79 80 55 83 49 17 7 138 15 11 139 19 6 3 117 10 50
## [77] 100 13 65 40 21
##
## [[73]]
## [1] 97 105 119 47 69 24 37 61 144 50 63 137 27 15 70 30 66 134 141
## [20] 127 113 36 12 96 90 106 142 54 85 67 98 120 129 3 103 91 26 121
## [39] 133 65 1 136 55 108 107 115 34 117 124 75 32 39 130 100 109 110 5
## [58] 95 58 93 101 99 73 77 44 68 11 135 38 56 94 138 29 80 45 62
## [77] 59 19 126 13 53 92
##
## [[74]]
## [1] 56 70 23 32 13 96 1 118 143 99 112 33 88 47 29 22 15 144 39
## [20] 38 26 134 71 114 28 142 100 58 133 74 18 136 109 130 75 60 120 103
## [39] 137 93 126 64 65 16 67 101 94 17 6 36 119 72 84 35 19 116 128
## [58] 27 63 140 138 79 82 57 25 10 7 43 127 89 9 95 76 21 123 14
## [77] 24 86 125 8 48 30 91
##
## [[75]]
## [1] 20 13 121 12 60 47 66 139 140 134 88 28 23 31 21 79 127 36 26
## [20] 65 16 128 97 87 89 111 25 73 2 38 67 136 44 45 109 92 11 1
## [39] 70 131 27 122 119 110 4 85 8 39 129 144 18 101 104 106 22 55 48
## [58] 43 114 141 53 94 142 51 3 117 17 81 58 80 115 99 100 59 74 69
## [77] 64 118 123 107 68 90 72 52
##
## [[76]]
## [1] 144 23 46 59 32 94 129 24 120 133 121 2 4 88 89 8 56 82 49
## [20] 109 41 139 142 101 97 141 38 127 128 25 74 44 95 140 80 84 29 76
## [39] 6 57 107 90 40 124 123 18 17 131 126 45 102 81 91 36 68 11 110
## [58] 103 79 99 113 20 34 48 19 108 31 10 78 69 130 122 50 1 70 37
## [77] 63 39 54 115 9 119 116 105 104
##
## [[77]]
## [1] 75 63 102 103 22 82 60 54 58 23 131 126 122 132 134 123 1 3 94
## [20] 120 111 20 65 34 9 8 86 47 70 140 66 85 93 27 48 26 24 59
## [39] 62 5 12 52 42 84 13 46 106 41 76 51 143 18 124 44 105 128 40
## [58] 98 91 141 39 81 100 95 136 80 78 114 90 74 32 79 113 144 56 116
## [77] 69 67 7 61 68 17 77 37 88 55
##
## [[78]]
## [1] 9 80 87 116 132 34 142 77 112 136 23 123 56 50 86 22 46 70 72
## [20] 3 137 47 4 138 88 61 5 21 6 143 102 74 16 113 104 117 44 81
## [39] 99 36 20 85 12 107 28 121 18 128 29 39 122 45 140 124 129 126 127
## [58] 115 2 26 35 54 135 7 94 33 1 38 101 84 141 131 83 78 95 31
## [77] 103 40 63 92 51 139 49 15 75 10 73
##
## [[79]]
## [1] 128 16 48 17 120 92 135 144 116 134 73 122 57 84 40 110 58 97 54
## [20] 14 51 129 103 130 31 29 109 133 96 5 25 47 65 56 7 26 38 63
## [39] 71 61 1 127 91 60 59 36 53 131 34 115 123 139 8 19 111 80 82
## [58] 107 88 66 108 49 72 41 4 42 137 35 94 77 15 87 106 24 68 67
## [77] 101 105 100 114 64 12 95 30 86 117 75 55
##
## [[80]]
## [1] 2 16 23 134 112 90 47 110 76 46 109 103 50 132 71 25 125 20 9
## [20] 84 129 141 121 77 35 105 45 128 79 61 10 119 24 85 82 74 140 5
## [39] 51 116 36 32 3 136 21 127 56 120 69 67 115 33 114 30 73 13 72
## [58] 100 95 101 29 4 8 39 133 80 108 49 139 96 83 54 14 102 6 68
## [77] 11 18 113 37 126 142 131 98 78 64 63 42 53
##
## [[81]]
## [1] 103 110 89 25 112 4 113 55 142 67 108 131 135 94 88 133 120 125 75
## [20] 24 82 53 111 117 8 35 39 80 96 62 5 61 93 87 40 37 123 28
## [39] 76 20 101 122 107 52 45 14 100 19 77 119 30 22 114 79 49 42 106
## [58] 51 32 143 115 57 73 38 109 104 66 9 64 31 63 68 59 132 34 48
## [77] 46 137 16 97 71 29 102 126 78 86 98 95 118 105
##
## [[82]]
## [1] 82 95 133 109 48 10 20 13 119 96 98 31 135 107 73 27 81 143 8
## [20] 35 25 120 42 3 118 72 65 134 71 117 121 41 126 84 47 86 18 49
## [39] 9 50 15 127 5 38 70 63 17 124 33 79 103 46 28 56 88 21 30
## [58] 78 142 94 11 57 39 116 141 51 54 43 62 115 66 114 37 83 6 29
## [77] 138 40 53 101 125 58 136 100 130 22 26 12 44 92 61
##
## [[83]]
## [1] 27 92 24 26 115 93 129 72 9 73 109 4 143 119 8 1 126 141 17
## [20] 116 15 83 75 11 106 79 89 136 49 70 55 90 22 124 85 82 35 2
## [39] 6 44 47 60 86 25 88 95 10 103 59 65 120 133 114 105 130 50 14
## [58] 43 97 32 33 20 118 80 98 36 41 117 62 123 3 110 71 142 34 69
## [77] 94 19 78 67 122 81 18 66 77 51 87 84 63 42 29 46
##
## [[84]]
## [1] 24 141 75 90 122 39 26 20 57 9 55 5 139 23 133 2 136 67 81
## [20] 69 94 96 110 36 56 129 35 1 101 105 76 83 28 137 48 54 72 46
## [39] 117 42 143 32 11 135 121 103 113 14 95 74 140 53 62 78 89 63 4
## [58] 88 49 34 17 44 30 19 38 126 45 130 119 93 71 118 86 6 64 29
## [77] 73 37 115 104 58 16 99 112 84 13 18 51 125 80 8 66 27
##
## [[85]]
## [1] 30 47 25 21 46 125 127 108 81 36 135 55 143 142 139 71 97 62 18
## [20] 1 113 131 93 2 112 141 120 90 121 61 110 144 8 98 124 53 69 66
## [39] 65 78 56 111 94 64 92 102 42 59 126 10 87 82 29 88 86 63 73
## [58] 101 34 76 32 68 38 115 134 100 23 74 84 43 31 57 27 106 138 89
## [77] 117 103 75 7 45 19 44 107 95 51 41 14 39 109 40 104 52 136
##
## [[86]]
## [1] 68 117 12 3 109 126 114 49 97 106 15 128 50 19 11 9 84 4 22
## [20] 17 101 6 10 31 140 122 75 82 80 32 51 60 105 36 64 21 35 107
## [39] 38 108 16 28 55 73 121 27 26 133 71 124 134 142 129 103 88 90 65
## [58] 143 96 87 43 8 2 66 33 25 79 69 53 130 135 139 54 52 113 119
## [77] 37 47 72 30 120 40 141 44 86 115 136 74 70 56 18 95 67 42 61
##
## [[87]]
## [1] 60 45 130 135 22 134 30 4 133 68 1 132 23 50 94 107 8 139 136
## [20] 16 40 12 35 5 2 109 39 126 122 56 18 115 114 99 93 78 33 37
## [39] 121 52 83 106 86 142 85 117 74 103 15 128 55 77 111 61 104 120 31
## [58] 25 76 89 101 129 69 95 144 14 67 47 53 140 79 72 59 125 62 48
## [77] 24 96 29 44 137 98 43 116 71 143 51 70 123 49 105 87 65 54 36
## [96] 41
##
## [[88]]
## [1] 57 81 88 135 3 129 61 131 136 96 111 127 36 75 45 84 2 144 133
## [20] 117 124 26 130 27 28 17 105 34 97 121 92 132 68 24 60 7 30 64
## [39] 47 38 67 66 70 100 14 72 118 9 76 103 77 46 59 22 32 123 71
## [58] 142 19 115 139 37 143 54 25 53 86 91 114 35 49 62 51 108 74 69
## [77] 16 55 85 125 23 90 87 5 63 98 48 80 79 20 58 52 31 29 99
## [96] 73 39
##
## [[89]]
## [1] 126 13 132 1 61 112 144 94 8 39 114 96 44 36 119 40 25 81 122
## [20] 53 3 51 87 16 20 117 2 110 38 72 34 82 23 35 116 109 124 58
## [39] 14 92 24 52 98 85 103 63 28 60 21 65 19 59 49 79 123 107 80
## [58] 64 68 102 62 134 115 47 42 130 27 17 88 77 104 7 69 95 90 78
## [77] 125 138 97 127 67 4 93 66 48 75 12 135 29 6 55 111 139 140 137
## [96] 30 26 89
##
## [[90]]
## [1] 54 129 108 49 1 142 119 115 58 100 35 33 72 29 13 136 101 143 107
## [20] 20 16 106 133 78 110 6 135 47 116 111 5 85 43 8 93 38 109 92
## [39] 74 86 112 9 14 62 2 73 63 55 37 114 21 138 120 61 22 71 82
## [58] 52 95 140 10 81 123 4 19 83 44 87 70 68 126 34 76 59 51 139
## [77] 17 125 60 7 69 65 98 11 117 18 134 41 131 91 53 42 122 77 84
## [96] 26 50 79 128
##
## [[91]]
## [1] 115 46 125 51 72 120 13 92 48 38 33 121 3 61 30 25 140 40
## [19] 7 133 63 71 12 16 79 131 14 64 90 89 104 109 129 76 108 130
## [37] 91 67 11 105 126 27 135 24 132 50 26 124 78 43 74 4 100 15
## [55] 32 49 41 73 5 65 144 35 60 9 31 53 98 116 96 88 19 28
## [73] 47 84 17 52 42 23 62 18 117 113 22 114 94 141 56 82 99 138
## [91] 69 21 6 136 8 111 57 66 139 123
##
## [[92]]
## [1] 135 90 73 89 4 2 133 131 102 31 64 45 5 63 87 67 41 11
## [19] 99 36 22 39 51 21 86 97 28 107 75 80 139 132 44 118 100 109
## [37] 112 134 14 55 12 13 127 43 35 138 124 77 130 78 110 58 65 62
## [55] 125 19 95 66 3 113 56 23 128 42 115 60 141 54 48 85 29 7
## [73] 81 1 117 26 27 52 111 126 16 82 15 137 37 88 61 106 129 98
## [91] 119 6 20 84 40 94 96 136 70 30 79
##
## [[93]]
## [1] 10 113 61 26 107 2 130 120 109 131 108 71 37 31 55 125 101 27
## [19] 86 90 112 36 88 35 25 11 133 43 121 48 65 117 3 85 59 82
## [37] 50 97 16 105 81 62 66 136 45 33 5 126 104 18 124 79 115 142
## [55] 13 9 144 84 93 70 135 132 23 14 134 94 40 91 53 38 44 95
## [73] 46 114 24 52 58 138 21 6 51 67 30 140 64 22 128 87 32 15
## [91] 89 12 110 77 68 116 28 100 17 42 141 119
##
## [[94]]
## [1] 2 51 125 11 91 121 103 13 143 79 56 114 41 86 87 46 75 47
## [19] 127 122 12 38 113 32 18 14 50 144 22 7 26 134 136 1 137 40
## [37] 128 57 55 89 61 78 29 142 124 6 116 101 65 31 109 53 115 34
## [55] 49 10 120 67 54 43 59 93 131 62 94 110 68 129 45 15 82 111
## [73] 76 19 3 81 20 42 66 60 17 71 99 98 28 85 5 77 92 132
## [91] 48 117 107 90 21 52 138 30 69 135 73 130 97
##
## [[95]]
## [1] 9 12 122 142 101 22 131 102 66 31 47 106 120 79 143 16 41 126
## [19] 11 73 94 100 136 138 59 1 6 141 35 20 46 117 21 119 40 7
## [37] 5 72 2 90 121 34 15 118 134 132 125 109 53 33 135 89 13 32
## [55] 3 110 37 130 61 137 57 60 27 75 123 49 52 97 68 26 48 88
## [73] 56 65 83 55 28 70 96 128 63 29 54 38 77 71 62 36 78 18
## [91] 129 76 69 64 44 81 51 99 115 17 80 8 10 93
##
## [[96]]
## [1] 131 6 12 13 47 27 26 128 2 59 60 109 68 139 7 64 78 31
## [19] 142 129 90 124 24 132 92 88 133 70 101 16 15 99 72 25 32 103
## [37] 138 143 4 110 106 81 8 20 18 116 21 108 61 125 100 17 10 45
## [55] 135 50 40 85 115 37 137 34 44 118 3 112 38 117 58 94 11 111
## [73] 119 22 123 97 87 39 1 36 19 67 42 69 140 56 93 55 83 46
## [91] 33 120 134 75 130 79 41 57 127 73 86 53 80 54 89
##
## [[97]]
## [1] 121 59 12 25 119 36 126 3 118 51 96 82 24 64 72 133 120 97
## [19] 107 15 131 92 10 23 13 128 29 19 34 140 14 81 11 89 27 94
## [37] 41 122 17 110 108 26 20 7 125 83 142 37 28 39 84 68 61 124
## [55] 112 115 88 95 69 43 91 40 90 75 141 136 62 56 105 2 57 53
## [73] 135 134 76 137 65 101 63 79 104 50 106 70 32 21 71 138 78 9
## [91] 55 80 38 60 67 132 18 127 66 116 139 33 58 109 42 102
##
## [[98]]
## [1] 129 136 2 9 131 12 28 3 62 13 4 105 119 48 47 134 20 6
## [19] 93 5 143 125 111 73 110 44 74 70 11 115 43 102 86 33 25 113
## [37] 38 85 144 94 127 59 133 1 34 26 60 106 39 112 99 132 114 35
## [55] 58 53 80 138 68 56 32 77 107 15 120 98 29 87 21 54 118 18
## [73] 57 141 27 16 135 65 45 79 88 8 19 76 63 24 92 41 51 97
## [91] 90 117 137 14 95 75 122 31 69 67 103 40 22 100 104 36 83
##
## [[99]]
## [1] 34 68 35 2 104 41 99 12 115 60 4 130 140 8 21 108 57 66
## [19] 36 72 111 122 143 24 120 81 137 90 65 70 13 94 110 123 84 118
## [37] 30 46 102 56 59 80 29 124 71 38 77 58 3 127 64 62 25 112
## [55] 44 128 107 93 132 78 51 9 50 16 105 22 31 33 15 17 14 101
## [73] 119 19 95 79 76 75 138 28 73 10 83 91 27 121 54 117 88 106
## [91] 74 133 5 114 126 52 23 100 43 45 92 129 42 7 67 139 47 32prev_= sapply(res_, function(r){100*table(sana_[r])[2]/length(r)})
plot(10:108,prev_,pch=16)
abline(h= 100*table(sana_)[2]/length(sana_), col="red")
plot(xy_,pch=15,cex=2,col=sana_+3)
points(xy_[res_[[41]],], col=0,pch=16)
set.seed(1013687724)
xy_$edad_ = sort(round(runif(144,10,15)))
res2_=clhs(xy_,size = 50, iter = 1000, progress = FALSE)
table(xy_$edad_[res2_])##
## 10 11 12 13 14 15
## 4 9 13 11 11 2plot(xy_[,1:2],pch=15,cex=2,col=sana_+3)
points(xy_[res2_,1:2], col=0,pch=16)
100*table(sana_[res2_])/length(res2_)##
## 0 1
## 80 20