Proyecto Estadísticaaa
Doménica Obando
2025-06-04
setwd("/cloud/project/PROYECTO/DATOS")
datos<-read.csv("petroo.csv",
header= TRUE, sep= ",", dec=".")
setwd("/cloud/project/PROYECTO/DATOS/")
datos2<-read.csv("DISCRETAS.csv",
header= TRUE, sep= ",", dec=".")ESTADÍSTICA DESCRIPTIVA
VARIABLES NOMINALES
Tipo de pozo
Tabla de distribución de frecuencias
Tipo<-datos$WELL_TYPE
Tipo<-na.omit(Tipo)
TDFTipo<-table(Tipo)
TablaTipo<-as.data.frame(TDFTipo)
hi_Tipo<-TablaTipo$Freq/sum(TablaTipo$Freq)
hi_Tipo<-hi_Tipo*100
TablaTipo<-data.frame(TablaTipo,hi_Tipo)Crear agrupación de Tipo de pozo
for (i in 1:27) {
tipo <- as.character(TablaTipo$Tipo[i])
if (tipo %in% c("Oil Well", "Natural Gas Well", "Oil and Gas Well", "Private Gas Well", "Historical Oil Well")) {
TablaTipo$GrupoTipo[i] <- "Productivo (Oil/Gas)"
} else if (tipo %in% c("Brine Well", "Brine Well Oil Well", "Brine Well Natural Gas Well", "Solution Mining Well", "Source Well")) {
TablaTipo$GrupoTipo[i] <- "Productivo (Brine/Solution Mining)"
} else if (tipo %in% c("Injection Well", "Historical Injection Well", "Natural Gas Storage Well", "Cavern Storage Well", "Compressed Air Energy Storage Well", "Disposal Well")) {
TablaTipo$GrupoTipo[i] <- "Inyección/Almacenamiento"
} else if (tipo %in% c("Dry Hole", "Gas Show", "Oil Show", "Oil and Gas Show", "Stratigraphic Test")) {
TablaTipo$GrupoTipo[i] <- "Exploratorio/Evaluación"
} else {
TablaTipo$GrupoTipo[i] <- "Observación/Otro"
}
}
GrupoTipo<-table(TablaTipo$GrupoTipo)
colnames(TablaTipo)<-c("Tipo de pozo","ni","hi","Grupo")
TablaTipo## Tipo de pozo ni hi
## 1 Brine Well 34 0.127503188
## 2 Brine Well Natural Gas Well 4 0.015000375
## 3 Brine Well Oil Well 5 0.018750469
## 4 Cavern Storage Well 124 0.465011625
## 5 Compressed Air Energy Storage Well 2 0.007500188
## 6 Disposal Well 129 0.483762094
## 7 Dry Hole 4546 17.047926198
## 8 Gas Show 2728 10.230255756
## 9 Gas Well Oil Show 56 0.210005250
## 10 Historical Injection Well 49 0.183754594
## 11 Historical Oil Well 1210 4.537613440
## 12 Injection Well 235 0.881272032
## 13 Licensed 2 0.007500188
## 14 Location 936 3.510087752
## 15 Natural Gas Storage Well 333 1.248781220
## 16 Natural Gas Well 9355 35.082127053
## 17 Observation Well 176 0.660016500
## 18 Oil and Gas Show 795 2.981324533
## 19 Oil and Gas Well 253 0.948773719
## 20 Oil Show 977 3.663841596
## 21 Oil Well 2102 7.882697067
## 22 Oil Well Gas Show 53 0.198754969
## 23 Other 19 0.071251781
## 24 Private Gas Well 1057 3.963849096
## 25 Solution Mining Well 204 0.765019125
## 26 Source Well 10 0.037500938
## 27 Stratigraphic Test 1272 4.770119253
## Grupo
## 1 Productivo (Brine/Solution Mining)
## 2 Observación/Otro
## 3 Productivo (Brine/Solution Mining)
## 4 Inyección/Almacenamiento
## 5 Inyección/Almacenamiento
## 6 Inyección/Almacenamiento
## 7 Exploratorio/Evaluación
## 8 Exploratorio/Evaluación
## 9 Observación/Otro
## 10 Inyección/Almacenamiento
## 11 Productivo (Oil/Gas)
## 12 Inyección/Almacenamiento
## 13 Observación/Otro
## 14 Observación/Otro
## 15 Observación/Otro
## 16 Productivo (Oil/Gas)
## 17 Observación/Otro
## 18 Exploratorio/Evaluación
## 19 Productivo (Oil/Gas)
## 20 Exploratorio/Evaluación
## 21 Productivo (Oil/Gas)
## 22 Observación/Otro
## 23 Observación/Otro
## 24 Observación/Otro
## 25 Productivo (Brine/Solution Mining)
## 26 Productivo (Brine/Solution Mining)
## 27 Exploratorio/EvaluaciónGDF versión local grupos de tipo de pozo
nombres_Tipo <- names(TDFTipo)
nuevos_nombres_Tipo<- gsub("/", "/\n", nombres_Tipo)
nombres_Grupo_Tipo <- names(GrupoTipo)
nuevos_nombres_Grupo_Tipo<- gsub("/", "/\n", nombres_Grupo_Tipo)
barplot(GrupoTipo,
main="Gráfica No.1:
Distrubución de grupos de tipo de pozo",
xlab="Tipo de pozo",
ylab="Cantidad",
col="darkblue",
names.arg = nuevos_nombres_Grupo_Tipo,
cex.main = 1.5,
cex.lab = 1,
cex.axis = 0.7,
cex.names = 0.7)
GDF versión global de agrupación de Tipo de pozo
barplot(GrupoTipo,
main="Gráfica No.2:
Distrubución de grupos de tipo de pozo",
xlab="Tipo de pozo",
ylab="Cantidad",
col="darkblue",
names.arg = nuevos_nombres_Grupo_Tipo,
cex.main = 1.5,
cex.lab = 1,
cex.axis = 0.7,
cex.names = 0.7,
ylim=c(0,sum(GrupoTipo)))
GDF porcentaje de agrupación de Tipo de pozo
hi_grupo_Tipo<-GrupoTipo/sum(GrupoTipo)
hi_grupo_Tipo<-hi_grupo_Tipo*100
barplot(hi_grupo_Tipo,
main="Gráfica No.3:
Distrubución de grupos de tipo de pozo",
xlab="Tipo de pozo",
ylab="Porcentaje",
col="darkblue",
names.arg = nuevos_nombres_Grupo_Tipo,
cex.main = 1.5,
cex.lab = 1,
cex.axis = 0.7,
cex.names = 0.7,
ylim=c(0,100))
Diagrama círcular de Tipo de pozo
circular_tipo<-etiqueta_Tipo<-paste(round(hi_grupo_Tipo,2),"%")
pie(hi_grupo_Tipo, main="Gráfica No.4:
Distribución de grupos de Tipos de pozo",
radius = 0.9,
col=colores<-c(rev(heat.colors(5))),
labels = etiqueta_Tipo)
legend("bottomleft", legend= names(hi_grupo_Tipo),
title="Leyenda",
cex = 0.5,
fill=colores<-c(rev(heat.colors(5))))
Variable Condado
Tabla de distribución de frecuencias
Condado<-datos$COUNTY
Condado<-na.omit(Condado)
TDFcondado<-table(Condado)
TablaCondado<-as.data.frame(TDFcondado)
hiCondado<-TablaCondado$Freq/sum(TablaCondado$Freq)
hiCondado<-hiCondado*100
TablaCondado<-data.frame(TablaCondado,hiCondado)Agrupación de Condados
for (i in 1:44) {
condado<-TablaCondado$`Condado`[i]
if (condado %in% c("Algoma", "Cochrane", "Kenora", "Kenora/Patricia Portion", "Manitoulin", "Timiskaming")) {
TablaCondado$GrupoCondado[i] <- "Norte de Ontario"
} else if (condado %in% c("Dundas", "Glengarry","Hastings", "Ottawa/Carleton", "Prescott/Russell", "Stormont/Dundas/Glengarry", "Prince Edward")) {
TablaCondado$GrupoCondado[i] <- "Este de Ontario"
} else if (condado %in% c("Dufferin", "Lanark","Northumberland", "Peterborough", "Simcoe","Victoria")) {
TablaCondado$GrupoCondado[i] <- "Centro de Ontario"
} else {
TablaCondado$GrupoCondado[i] <- "Sur de Ontario"
}
}
GrupoCondado<-table(TablaCondado$GrupoCondado)
colnames(TablaCondado)<-c("Condado","ni","hi","Grupo")
TablaCondado## Condado ni hi Grupo
## 1 2 0.007500188 Sur de Ontario
## 2 Algoma 3 0.011250281 Norte de Ontario
## 3 Brant 461 1.728793220 Sur de Ontario
## 4 Bruce 123 0.461261532 Sur de Ontario
## 5 Cochrane 194 0.727518188 Norte de Ontario
## 6 Dufferin 14 0.052501313 Centro de Ontario
## 7 Dundas 4 0.015000375 Este de Ontario
## 8 Durham 34 0.127503188 Sur de Ontario
## 9 Elgin 1905 7.143928598 Sur de Ontario
## 10 Essex 1115 4.181354534 Sur de Ontario
## 11 Glengarry 15 0.056251406 Este de Ontario
## 12 Grey 82 0.307507688 Sur de Ontario
## 13 Haldimand 6014 22.553063827 Sur de Ontario
## 14 Halton 66 0.247506188 Sur de Ontario
## 15 Hastings 21 0.078751969 Este de Ontario
## 16 Huron 256 0.960024001 Sur de Ontario
## 17 Kenora 19 0.071251781 Norte de Ontario
## 18 Kenora/Patricia Portion 1 0.003750094 Norte de Ontario
## 19 Kent 3117 11.689042226 Sur de Ontario
## 20 Lambton 5712 21.420535513 Sur de Ontario
## 21 Lanark 2 0.007500188 Centro de Ontario
## 22 Lennox/Addington 5 0.018750469 Sur de Ontario
## 23 Lincoln 519 1.946298657 Sur de Ontario
## 24 Manitoulin 152 0.570014250 Norte de Ontario
## 25 Middlesex 1148 4.305107628 Sur de Ontario
## 26 Norfolk 2693 10.099002475 Sur de Ontario
## 27 Northumberland 10 0.037500938 Centro de Ontario
## 28 Ottawa/Carleton 30 0.112502813 Este de Ontario
## 29 Oxford 486 1.822545564 Sur de Ontario
## 30 Peel 28 0.105002625 Sur de Ontario
## 31 Perth 34 0.127503188 Sur de Ontario
## 32 Peterborough 3 0.011250281 Centro de Ontario
## 33 Prescott/Russell 38 0.142503563 Este de Ontario
## 34 Prince Edward 28 0.105002625 Este de Ontario
## 35 Simcoe 53 0.198754969 Centro de Ontario
## 36 Stormont/Dundas/Glengarry 16 0.060001500 Este de Ontario
## 37 Timiskaming 1 0.003750094 Norte de Ontario
## 38 Toronto 24 0.090002250 Sur de Ontario
## 39 Victoria 20 0.075001875 Centro de Ontario
## 40 Waterloo 3 0.011250281 Sur de Ontario
## 41 Welland 1892 7.095177379 Sur de Ontario
## 42 Wellington 51 0.191254781 Sur de Ontario
## 43 Wentworth 253 0.948773719 Sur de Ontario
## 44 York 19 0.071251781 Sur de OntarioGDF de agrupación por Condado versión local
barplot(GrupoCondado,
main = "Gráfica No.1:
Distribución por grupos de Condados de Ontario",
xlab = "Condados",
ylab = "Cantidad",
cex.main = 1.5,
cex.lab = 1,
cex.axis = 0.7,
cex.names = 0.7,
col = "darkgreen")
GDF de agrupación por Condado versión global
barplot(GrupoCondado,
main="Gráfica No.2:
Distrubución por grupos de Condados de Ontario",
xlab="Condados",
ylab="Cantidad",
col="darkgreen",
cex.names = 0.7,
ylim=c(0,sum(GrupoCondado)))
GDF porcentaje de agrupación de condado
hi_grupoCondado<-GrupoCondado/sum(GrupoCondado)
hi_grupoCondado<-round((hi_grupoCondado*100),2)
barplot(hi_grupoCondado,
main="Gráfica No.3:
Distrubución de grupos de Condados",
xlab="Condados",
ylab="Cantidad",
col="darkgreen",
cex.axis = 0.7,
cex.names = 0.7,
ylim=c(0,100))
Diagrama círcular de agrupación de condados
etiquetaCondado<-paste(hi_grupoCondado,"%")
pie(hi_grupoCondado, main="Gráfica No.4:
Distribución de grupos de Condados de Ontario",
radius = 1,
col=colores<-c(rev(heat.colors(4))),
labels = etiquetaCondado)
legend("bottomright", legend= names(hi_grupoCondado),
title="Leyenda",
cex = 0.75,
fill=colores<-c(rev(heat.colors(4))))
Variable Municipio
municipio<-datos$TOWNSHIP
municipio<-na.omit(municipio)
TDFmunicipio<-table(municipio)
TablaMunicipio<-as.data.frame(TDFmunicipio)
hiMunicipio<-TablaMunicipio$Freq/sum(TablaMunicipio$Freq)
hiMunicipio<-round(hiMunicipio*100,2)
TablaMunicipio<-data.frame(TablaMunicipio,hiMunicipio)Agrupación de Municipio
for (i in 1:293) {
frecuencia3 <- TablaMunicipio$Freq[i]
if (frecuencia3<100) {
TablaMunicipio$Grupo_Actividad_Municipio[i] <- "Baja actividad"
} else if (frecuencia3 >= 100 & frecuencia3 < 700) {
TablaMunicipio$Grupo_Actividad_Municipio[i] <- "Actividad media"
} else {
TablaMunicipio$Grupo_Actividad_Municipio[i] <- "Alta actividad"
}
}
Grupo_Actividad_Municipio<-table(TablaMunicipio$Grupo_Actividad_Municipio)
colnames(TablaMunicipio)<-c("Municipio","ni","hi","Grupo")
TablaMunicipio## Municipio ni hi Grupo
## 1 3 0.01 Baja actividad
## 2 Adelaide 86 0.32 Baja actividad
## 3 Adjala 4 0.02 Baja actividad
## 4 Albemarle 13 0.05 Baja actividad
## 5 Albion 1 0.00 Baja actividad
## 6 Aldborough 732 2.75 Alta actividad
## 7 Alfred 2 0.01 Baja actividad
## 8 Allan 1 0.00 Baja actividad
## 9 Amabel 38 0.14 Baja actividad
## 10 Amaranth 1 0.00 Baja actividad
## 11 Ameliasburgh 1 0.00 Baja actividad
## 12 Ancaster 8 0.03 Baja actividad
## 13 Anderdon 99 0.37 Baja actividad
## 14 Artemesia 6 0.02 Baja actividad
## 15 Arthur 14 0.05 Baja actividad
## 16 Ashfield 14 0.05 Baja actividad
## 17 Assiginack 30 0.11 Baja actividad
## 18 Athol 4 0.02 Baja actividad
## 19 Barton 6 0.02 Baja actividad
## 20 Bayham 296 1.11 Actividad media
## 21 Bentinck 1 0.00 Baja actividad
## 22 Bertie 544 2.04 Actividad media
## 23 Beverly 15 0.06 Baja actividad
## 24 Bexley 1 0.00 Baja actividad
## 25 Biddulph 7 0.03 Baja actividad
## 26 Bidwell 17 0.06 Baja actividad
## 27 Binbrook 161 0.60 Actividad media
## 28 Blandford 73 0.27 Baja actividad
## 29 Blanshard 5 0.02 Baja actividad
## 30 Blenheim 132 0.50 Actividad media
## 31 Bosanquet 58 0.22 Baja actividad
## 32 Brant 4 0.02 Baja actividad
## 33 Brantford 49 0.18 Baja actividad
## 34 Brock 2 0.01 Baja actividad
## 35 Brooke 261 0.98 Actividad media
## 36 Bruce 12 0.05 Baja actividad
## 37 Burford 51 0.19 Baja actividad
## 38 Caistor 318 1.19 Actividad media
## 39 Caledon 7 0.03 Baja actividad
## 40 Caledonia 2 0.01 Baja actividad
## 41 Camden 249 0.93 Actividad media
## 42 Campbell 1 0.00 Baja actividad
## 43 Canborough 392 1.47 Actividad media
## 44 Caradoc 41 0.15 Baja actividad
## 45 Carden 2 0.01 Baja actividad
## 46 Carnarvon 4 0.02 Baja actividad
## 47 Charlottenburgh 9 0.03 Baja actividad
## 48 Charlotteville 458 1.72 Actividad media
## 49 Chatham 230 0.86 Actividad media
## 50 Chinguacousy 7 0.03 Baja actividad
## 51 City of Oshawa 1 0.00 Baja actividad
## 52 City of Owen Sound 1 0.00 Baja actividad
## 53 City of Toronto 24 0.09 Baja actividad
## 54 City of Windsor 5 0.02 Baja actividad
## 55 Clarence 4 0.02 Baja actividad
## 56 Clarke 2 0.01 Baja actividad
## 57 Clinton 5 0.02 Baja actividad
## 58 Cockburn Island 4 0.02 Baja actividad
## 59 Colborne 7 0.03 Baja actividad
## 60 Colchester North 9 0.03 Baja actividad
## 61 Colchester South 72 0.27 Baja actividad
## 62 Collingwood 24 0.09 Baja actividad
## 63 Cornwall 3 0.01 Baja actividad
## 64 Cramahe 4 0.02 Baja actividad
## 65 Crowland 161 0.60 Actividad media
## 66 Culross 4 0.02 Baja actividad
## 67 Cumberland 4 0.02 Baja actividad
## 68 Darlington 16 0.06 Baja actividad
## 69 Dawn 620 2.33 Actividad media
## 70 Delaware 21 0.08 Baja actividad
## 71 Dereham 145 0.54 Actividad media
## 72 District of Cochrane 194 0.73 Actividad media
## 73 District of Kenora 19 0.07 Baja actividad
## 74 Dover 202 0.76 Actividad media
## 75 Dunn 374 1.40 Actividad media
## 76 Dunwich 201 0.75 Actividad media
## 77 East Flamborough 2 0.01 Baja actividad
## 78 East Gwillimbury 2 0.01 Baja actividad
## 79 East Hawkesbury 1 0.00 Baja actividad
## 80 East Nissouri 7 0.03 Baja actividad
## 81 East Oxford 4 0.02 Baja actividad
## 82 East Wawanosh 4 0.02 Baja actividad
## 83 East Williams 8 0.03 Baja actividad
## 84 East Zorra 29 0.11 Baja actividad
## 85 Eastnor 24 0.09 Baja actividad
## 86 Egremont 4 0.02 Baja actividad
## 87 Ekfrid 91 0.34 Baja actividad
## 88 Eldon 3 0.01 Baja actividad
## 89 Elma 6 0.02 Baja actividad
## 90 Emily 3 0.01 Baja actividad
## 91 Enniskillen 2154 8.08 Alta actividad
## 92 Eramosa 3 0.01 Baja actividad
## 93 Esquesing 37 0.14 Baja actividad
## 94 Essa 3 0.01 Baja actividad
## 95 Euphemia 178 0.67 Actividad media
## 96 Euphrasia 1 0.00 Baja actividad
## 97 Fenelon 1 0.00 Baja actividad
## 98 Flos 1 0.00 Baja actividad
## 99 Fullarton 1 0.00 Baja actividad
## 100 Gainsborough 167 0.63 Actividad media
## 101 Glanford 55 0.21 Baja actividad
## 102 Glenelg 1 0.00 Baja actividad
## 103 Gloucester 13 0.05 Baja actividad
## 104 Goderich 59 0.22 Baja actividad
## 105 Gordon 12 0.05 Baja actividad
## 106 Gosfield North 30 0.11 Baja actividad
## 107 Gosfield South 172 0.65 Actividad media
## 108 Grantham 3 0.01 Baja actividad
## 109 Greenock 2 0.01 Baja actividad
## 110 Grey 3 0.01 Baja actividad
## 111 Grimsby 14 0.05 Baja actividad
## 112 Guelph 9 0.03 Baja actividad
## 113 Hallowell 19 0.07 Baja actividad
## 114 Harris 1 0.00 Baja actividad
## 115 Harvey 1 0.00 Baja actividad
## 116 Harwich 120 0.45 Actividad media
## 117 Hay 23 0.09 Baja actividad
## 118 Hibbert 1 0.00 Baja actividad
## 119 Hilton 1 0.00 Baja actividad
## 120 Hope 5 0.02 Baja actividad
## 121 Houghton 150 0.56 Actividad media
## 122 Howard 95 0.36 Baja actividad
## 123 Howland 6 0.02 Baja actividad
## 124 Hullett 18 0.07 Baja actividad
## 125 Humberstone 448 1.68 Actividad media
## 126 Huron 5 0.02 Baja actividad
## 127 Indian Reserve no. 26 43 0.16 Baja actividad
## 128 Indian Reserve no. 44 1 0.00 Baja actividad
## 129 Jocelyn 1 0.00 Baja actividad
## 130 Kaladar/Anglesea/Effingham 1 0.00 Baja actividad
## 131 Kenyon 2 0.01 Baja actividad
## 132 Keppel 21 0.08 Baja actividad
## 133 Killarney 1 0.00 Baja actividad
## 134 Kincardine 6 0.02 Baja actividad
## 135 King 1 0.00 Baja actividad
## 136 Kinloss 2 0.01 Baja actividad
## 137 Lake Erie 2107 7.90 Alta actividad
## 138 Lake Ontario 2 0.01 Baja actividad
## 139 Lake St. Clair 11 0.04 Baja actividad
## 140 Lancaster 9 0.03 Baja actividad
## 141 Laxton/Digby/Longford 1 0.00 Baja actividad
## 142 Lindsay 4 0.02 Baja actividad
## 143 Lobo 11 0.04 Baja actividad
## 144 Lochiel 2 0.01 Baja actividad
## 145 Logan 7 0.03 Baja actividad
## 146 London 11 0.04 Baja actividad
## 147 Louth 3 0.01 Baja actividad
## 148 Maidstone 25 0.09 Baja actividad
## 149 Malahide 195 0.73 Actividad media
## 150 Malden 77 0.29 Baja actividad
## 151 Mara 4 0.02 Baja actividad
## 152 Mariposa 2 0.01 Baja actividad
## 153 Markham 6 0.02 Baja actividad
## 154 Maryborough 1 0.00 Baja actividad
## 155 Matilda 5 0.02 Baja actividad
## 156 McGillivray 61 0.23 Baja actividad
## 157 McKillop 8 0.03 Baja actividad
## 158 Medonte 1 0.00 Baja actividad
## 159 Melancthon 6 0.02 Baja actividad
## 160 Mersea 357 1.34 Actividad media
## 161 Metcalfe 136 0.51 Actividad media
## 162 Middleton 261 0.98 Actividad media
## 163 Minto 5 0.02 Baja actividad
## 164 Mono 5 0.02 Baja actividad
## 165 Moore 1032 3.87 Alta actividad
## 166 Morris 4 0.02 Baja actividad
## 167 Mosa 604 2.27 Actividad media
## 168 Moulton 503 1.89 Actividad media
## 169 Mulmur 2 0.01 Baja actividad
## 170 Murray 11 0.04 Baja actividad
## 171 Nassagaweya 8 0.03 Baja actividad
## 172 Nelson 4 0.02 Baja actividad
## 173 Nepean 7 0.03 Baja actividad
## 174 Niagara 9 0.03 Baja actividad
## 175 Nichol 1 0.00 Baja actividad
## 176 North Cayuga 656 2.46 Actividad media
## 177 North Dorchester 21 0.08 Baja actividad
## 178 North Dumfries 1 0.00 Baja actividad
## 179 North Easthope 3 0.01 Baja actividad
## 180 North Elmsley 2 0.01 Baja actividad
## 181 North Gwillimbury 1 0.00 Baja actividad
## 182 North Norwich 6 0.02 Baja actividad
## 183 North Oxford 32 0.12 Baja actividad
## 184 North Plantagenet 3 0.01 Baja actividad
## 185 North Walsingham 364 1.37 Actividad media
## 186 Nottawasaga 26 0.10 Baja actividad
## 187 Oakland 4 0.02 Baja actividad
## 188 Oneida 422 1.58 Actividad media
## 189 Onondaga 183 0.69 Actividad media
## 190 Ops 1 0.00 Baja actividad
## 191 Orford 141 0.53 Actividad media
## 192 Orillia 2 0.01 Baja actividad
## 193 Osgoode 5 0.02 Baja actividad
## 194 Osprey 4 0.02 Baja actividad
## 195 Otonabee 2 0.01 Baja actividad
## 196 Peel 2 0.01 Baja actividad
## 197 Pelee 23 0.09 Baja actividad
## 198 Pelham 21 0.08 Baja actividad
## 199 Pickering 7 0.03 Baja actividad
## 200 Pilkington 1 0.00 Baja actividad
## 201 Plympton 163 0.61 Actividad media
## 202 Proton 1 0.00 Baja actividad
## 203 Puslinch 9 0.03 Baja actividad
## 204 Rainham 626 2.35 Actividad media
## 205 Raleigh 429 1.61 Actividad media
## 206 Richmond 4 0.02 Baja actividad
## 207 Robinson 8 0.03 Baja actividad
## 208 Rochester 62 0.23 Baja actividad
## 209 Romney 421 1.58 Actividad media
## 210 Roxborough 1 0.00 Baja actividad
## 211 Russell 24 0.09 Baja actividad
## 212 Saltfleet 2 0.01 Baja actividad
## 213 Sandfield 3 0.01 Baja actividad
## 214 Sandwich East 1 0.00 Baja actividad
## 215 Sandwich South 2 0.01 Baja actividad
## 216 Sandwich West 68 0.26 Baja actividad
## 217 Sarawak 7 0.03 Baja actividad
## 218 Sarnia 426 1.60 Actividad media
## 219 Saugeen 3 0.01 Baja actividad
## 220 Scott 1 0.00 Baja actividad
## 221 Seneca 660 2.48 Actividad media
## 222 Seymour 1 0.00 Baja actividad
## 223 Sheguiandah 16 0.06 Baja actividad
## 224 Sherbrooke 71 0.27 Baja actividad
## 225 Sidney 1 0.00 Baja actividad
## 226 Sombra 542 2.03 Actividad media
## 227 Somerville 2 0.01 Baja actividad
## 228 South Cayuga 391 1.47 Actividad media
## 229 South Dorchester 9 0.03 Baja actividad
## 230 South Dumfries 2 0.01 Baja actividad
## 231 South Easthope 9 0.03 Baja actividad
## 232 South Fredericksburgh 1 0.00 Baja actividad
## 233 South Marysburgh 4 0.02 Baja actividad
## 234 South Norwich 32 0.12 Baja actividad
## 235 South Walsingham 151 0.57 Actividad media
## 236 Southwold 28 0.11 Baja actividad
## 237 St. Edmunds 2 0.01 Baja actividad
## 238 St. Joseph 1 0.00 Baja actividad
## 239 St. Vincent 10 0.04 Baja actividad
## 240 Stamford 16 0.06 Baja actividad
## 241 Stanley 26 0.10 Baja actividad
## 242 Stephen 61 0.23 Baja actividad
## 243 Sullivan 2 0.01 Baja actividad
## 244 Sydenham 3 0.01 Baja actividad
## 245 Tay 1 0.00 Baja actividad
## 246 Tecumseth 4 0.02 Baja actividad
## 247 Tehkummah 2 0.01 Baja actividad
## 248 Thorah 1 0.00 Baja actividad
## 249 Thorold 67 0.25 Baja actividad
## 250 Thurlow 5 0.02 Baja actividad
## 251 Tilbury East 407 1.53 Actividad media
## 252 Tilbury North 19 0.07 Baja actividad
## 253 Tilbury West 38 0.14 Baja actividad
## 254 Tiny 2 0.01 Baja actividad
## 255 Torbolton 1 0.00 Baja actividad
## 256 Toronto 13 0.05 Baja actividad
## 257 Tosorontio 1 0.00 Baja actividad
## 258 Town of Fort Erie 6 0.02 Baja actividad
## 259 Town of Whitchurch-Stouffville 2 0.01 Baja actividad
## 260 Townsend 196 0.74 Actividad media
## 261 Trafalgar 17 0.06 Baja actividad
## 262 Tuckersmith 11 0.04 Baja actividad
## 263 Turnberry 1 0.00 Baja actividad
## 264 Tuscarora 172 0.65 Actividad media
## 265 Tyendinaga 4 0.02 Baja actividad
## 266 Usborne 5 0.02 Baja actividad
## 267 Vaughan 7 0.03 Baja actividad
## 268 Verulam 4 0.02 Baja actividad
## 269 Wainfleet 334 1.25 Actividad media
## 270 Wallace 2 0.01 Baja actividad
## 271 Walpole 1291 4.84 Alta actividad
## 272 Warwick 276 1.04 Actividad media
## 273 Waterloo 1 0.00 Baja actividad
## 274 West Flamborough 4 0.02 Baja actividad
## 275 West Garafraxa 3 0.01 Baja actividad
## 276 West Gwillimbury 4 0.02 Baja actividad
## 277 West Hawkesbury 2 0.01 Baja actividad
## 278 West Luther 3 0.01 Baja actividad
## 279 West Nissouri 4 0.02 Baja actividad
## 280 West Oxford 17 0.06 Baja actividad
## 281 West Wawanosh 12 0.05 Baja actividad
## 282 West Williams 30 0.11 Baja actividad
## 283 West Zorra 9 0.03 Baja actividad
## 284 Westminster 17 0.06 Baja actividad
## 285 Whitby 2 0.01 Baja actividad
## 286 Williamsburg 2 0.01 Baja actividad
## 287 Willoughby 212 0.80 Actividad media
## 288 Wilmot 1 0.00 Baja actividad
## 289 Winchester 2 0.01 Baja actividad
## 290 Windham 199 0.75 Actividad media
## 291 Woodhouse 398 1.49 Actividad media
## 292 Yarmouth 33 0.12 Baja actividad
## 293 Zone 402 1.51 Actividad mediaGDF de agrupación por Municipio versión local
barplot(Grupo_Actividad_Municipio,
main = "Gráfica No.1:
Distribución por grupos de municipios según su actividad",
xlab = "Municipios",
ylab = "Cantidad",
cex.names = 0.7,
col = "red")
GDF de agrupación por Municipio versión global
barplot(Grupo_Actividad_Municipio,
main="Gráfica No.2:
Distrubución por grupos de municipios de Ontario",
xlab="Municipios",
ylab="Cantidad",
col="red",
cex.names = 0.7,
ylim=c(0,sum(Grupo_Actividad_Municipio)))
GDF porcentaje de agrupación de Municipio
hi_grupoMunicipio<-Grupo_Actividad_Municipio/sum(Grupo_Actividad_Municipio)
hi_grupoMunicipio<-round((hi_grupoMunicipio*100),2)
barplot(hi_grupoMunicipio,
main="Gráfica No.3:
Distrubución de grupos de municipios de Ontario",
xlab="Municipios",
ylab="Cantidad",
col="red",
cex.axis = 0.7,
cex.names = 0.7,
ylim=c(0,100))
Diagrama círcular de agrupación de municipios
etiquetaMunicipio<-paste(hi_grupoMunicipio,"%")
pie(hi_grupoMunicipio, main="Gráfica No.20:
Distribución de grupos de municipios de Ontario",
radius = 1,
col=colores<-c(rev(heat.colors(3))),
labels = etiquetaMunicipio)
legend("bottomright", legend= names(hi_grupoMunicipio),
title="Leyenda",
cex = 0.75,
fill=colores<-c(rev(heat.colors(3))))
Variable Concesión
Tabla de distribución de frecuencias
concesion<-datos$CONCESSION
concesion<-na.omit(concesion)
TDFconcesion<-table(concesion)
Tabla_Concesion<-as.data.frame(TDFconcesion)
hi_Concesion<-Tabla_Concesion$Freq/sum(Tabla_Concesion$Freq)
hi_Concesion<-hi_Concesion*100
Tabla_Concesion<-data.frame(Tabla_Concesion,hi_Concesion)Agrupación de Concesiones
for (i in 1:390) {
frecuencia <- Tabla_Concesion$Freq[i]
if (frecuencia < 100) {
Tabla_Concesion$Grupo_Concesion[i] <- "Baja frecuencia"
} else if (frecuencia >= 100 & frecuencia < 700) {
Tabla_Concesion$Grupo_Concesion[i] <- "Frecuencia media"
} else {
Tabla_Concesion$Grupo_Concesion[i] <- "Alta frecuencia"
}
}
Grupo_Concesion<-table(Tabla_Concesion$Grupo_Concesion)
colnames(Tabla_Concesion)<-c("Concesión","ni","hi","Grupo")
Tabla_Concesion## Concesión ni hi Grupo
## 1 2629 9.858996475 Alta frecuencia
## 2 34 0.127503188 Baja frecuencia
## 3 1C 6 0.022500563 Baja frecuencia
## 4 A 176 0.660016500 Frecuencia media
## 5 AB 6 0.022500563 Baja frecuencia
## 6 ABF 20 0.075001875 Baja frecuencia
## 7 AC 2 0.007500188 Baja frecuencia
## 8 ACC 12 0.045001125 Baja frecuencia
## 9 AED 5 0.018750469 Baja frecuencia
## 10 ARB 11 0.041251031 Baja frecuencia
## 11 ASTR 1 0.003750094 Baja frecuencia
## 12 B 60 0.225005625 Baja frecuencia
## 13 BA 13 0.048751219 Baja frecuencia
## 14 BAC 1 0.003750094 Baja frecuencia
## 15 BB 35 0.131253281 Baja frecuencia
## 16 BBF 12 0.045001125 Baja frecuencia
## 17 BF 12 0.045001125 Baja frecuencia
## 18 BFC 232 0.870021751 Frecuencia media
## 19 BFCNR 8 0.030000750 Baja frecuencia
## 20 BFCWR 1 0.003750094 Baja frecuencia
## 21 BFGR 3 0.011250281 Baja frecuencia
## 22 BFLE 28 0.105002625 Baja frecuencia
## 23 BFTR 6 0.022500563 Baja frecuencia
## 24 BI 1 0.003750094 Baja frecuencia
## 25 BK 1 0.003750094 Baja frecuencia
## 26 BLC 2 0.007500188 Baja frecuencia
## 27 BRT 2 0.007500188 Baja frecuencia
## 28 BUB 1 0.003750094 Baja frecuencia
## 29 C 28 0.105002625 Baja frecuencia
## 30 CBF 1 0.003750094 Baja frecuencia
## 31 CC 6 0.022500563 Baja frecuencia
## 32 COB 6 0.022500563 Baja frecuencia
## 33 CT 41 0.153753844 Baja frecuencia
## 34 D 13 0.048751219 Baja frecuencia
## 35 DB 3 0.011250281 Baja frecuencia
## 36 DEB 5 0.018750469 Baja frecuencia
## 37 DT 28 0.105002625 Baja frecuencia
## 38 E 1 0.003750094 Baja frecuencia
## 39 EBDR 6 0.022500563 Baja frecuencia
## 40 EFC 23 0.086252156 Baja frecuencia
## 41 ELR 3 0.011250281 Baja frecuencia
## 42 ENR 25 0.093752344 Baja frecuencia
## 43 EPR 1 0.003750094 Baja frecuencia
## 44 ETR 1 0.003750094 Baja frecuencia
## 45 F 4 0.015000375 Baja frecuencia
## 46 FC 89 0.333758344 Baja frecuencia
## 47 FCC 23 0.086252156 Baja frecuencia
## 48 FCED 94 0.352508813 Baja frecuencia
## 49 FCEMM 73 0.273756844 Baja frecuencia
## 50 FCNPA 1 0.003750094 Baja frecuencia
## 51 FCR 6 0.022500563 Baja frecuencia
## 52 FCSPA 1 0.003750094 Baja frecuencia
## 53 FCT 5 0.018750469 Baja frecuencia
## 54 FCWD 29 0.108752719 Baja frecuencia
## 55 FCWMM 27 0.101252531 Baja frecuencia
## 56 FGREM 2 0.007500188 Baja frecuencia
## 57 FGRWM 3 0.011250281 Baja frecuencia
## 58 FI 5 0.018750469 Baja frecuencia
## 59 FLH 8 0.030000750 Baja frecuencia
## 60 FRNLR 3 0.011250281 Baja frecuencia
## 61 FRSLR 2 0.007500188 Baja frecuencia
## 62 FRSUR 2 0.007500188 Baja frecuencia
## 63 FT 2 0.007500188 Baja frecuencia
## 64 G 12 0.045001125 Baja frecuencia
## 65 GA 6 0.022500563 Baja frecuencia
## 66 GB 11 0.041251031 Baja frecuencia
## 67 GO 90 0.337508438 Baja frecuencia
## 68 H 1 0.003750094 Baja frecuencia
## 69 HRC 2 0.007500188 Baja frecuencia
## 70 HUT 30 0.112502813 Baja frecuencia
## 71 I 2142 8.032700818 Alta frecuencia
## 72 IBF 1 0.003750094 Baja frecuencia
## 73 ICC 46 0.172504313 Baja frecuencia
## 74 ICEHS 1 0.003750094 Baja frecuencia
## 75 ICLE 48 0.180004500 Baja frecuencia
## 76 ICNC 54 0.202505063 Baja frecuencia
## 77 IE 13 0.048751219 Baja frecuencia
## 78 IEBR 15 0.056251406 Baja frecuencia
## 79 IECR 3 0.011250281 Baja frecuencia
## 80 IED 3 0.011250281 Baja frecuencia
## 81 IEPLR 58 0.217505438 Baja frecuencia
## 82 IGC 11 0.041251031 Baja frecuencia
## 83 IGCH 5 0.018750469 Baja frecuencia
## 84 IHRS 2 0.007500188 Baja frecuencia
## 85 II 2395 8.981474537 Alta frecuencia
## 86 IIC 13 0.048751219 Baja frecuencia
## 87 IICC 52 0.195004875 Baja frecuencia
## 88 IICLE 12 0.045001125 Baja frecuencia
## 89 IIE 7 0.026250656 Baja frecuencia
## 90 IIEBR 12 0.045001125 Baja frecuencia
## 91 IIECR 2 0.007500188 Baja frecuencia
## 92 IIED 9 0.033750844 Baja frecuencia
## 93 IIEPLR 23 0.086252156 Baja frecuencia
## 94 IIFTRT 3 0.011250281 Baja frecuencia
## 95 IIGC 3 0.011250281 Baja frecuencia
## 96 IIGCH 9 0.033750844 Baja frecuencia
## 97 III 1372 5.145128628 Alta frecuencia
## 98 IIIC 12 0.045001125 Baja frecuencia
## 99 IIICC 24 0.090002250 Baja frecuencia
## 100 IIICFE 6 0.022500563 Baja frecuencia
## 101 IIICLE 49 0.183754594 Baja frecuencia
## 102 IIIDC 2 0.007500188 Baja frecuencia
## 103 IIIE 26 0.097502438 Baja frecuencia
## 104 IIIEBR 11 0.041251031 Baja frecuencia
## 105 IIIED 8 0.030000750 Baja frecuencia
## 106 IIIEPLR 15 0.056251406 Baja frecuencia
## 107 IIIFTRT 5 0.018750469 Baja frecuencia
## 108 IIIG 1 0.003750094 Baja frecuencia
## 109 IIIGC 2 0.007500188 Baja frecuencia
## 110 IIIGCH 9 0.033750844 Baja frecuencia
## 111 IIIHRS 2 0.007500188 Baja frecuencia
## 112 IIILES 38 0.142503563 Baja frecuencia
## 113 IIILRS 2 0.007500188 Baja frecuencia
## 114 IIIMT 3 0.011250281 Baja frecuencia
## 115 IIIN 3 0.011250281 Baja frecuencia
## 116 IIINER 7 0.026250656 Baja frecuencia
## 117 IIINR 39 0.146253656 Baja frecuencia
## 118 IIINTR 5 0.018750469 Baja frecuencia
## 119 IIIREMPR 1 0.003750094 Baja frecuencia
## 120 IIIS 59 0.221255531 Baja frecuencia
## 121 IIISER 30 0.112502813 Baja frecuencia
## 122 IIISHR 1 0.003750094 Baja frecuencia
## 123 IIISR 43 0.161254031 Baja frecuencia
## 124 IIISSCR 35 0.131253281 Baja frecuencia
## 125 IIISTR 106 0.397509938 Frecuencia media
## 126 IIITRS 1 0.003750094 Baja frecuencia
## 127 IIIW 32 0.120003000 Baja frecuencia
## 128 IIIWBS 2 0.007500188 Baja frecuencia
## 129 IIIWCR 9 0.033750844 Baja frecuencia
## 130 IIIWD 1 0.003750094 Baja frecuencia
## 131 IIIWFC 2 0.007500188 Baja frecuencia
## 132 IIIWPR 13 0.048751219 Baja frecuencia
## 133 IILES 39 0.146253656 Baja frecuencia
## 134 IILRS 1 0.003750094 Baja frecuencia
## 135 IIMT 1 0.003750094 Baja frecuencia
## 136 IIN 1 0.003750094 Baja frecuencia
## 137 IINDS 1 0.003750094 Baja frecuencia
## 138 IINER 15 0.056251406 Baja frecuencia
## 139 IINGR 1 0.003750094 Baja frecuencia
## 140 IINR 32 0.120003000 Baja frecuencia
## 141 IINTR 30 0.112502813 Baja frecuencia
## 142 IINWSCR 3 0.011250281 Baja frecuencia
## 143 IIOR 1 0.003750094 Baja frecuencia
## 144 IIRGR 45 0.168754219 Baja frecuencia
## 145 IIRR 1 0.003750094 Baja frecuencia
## 146 IIRSHR 2 0.007500188 Baja frecuencia
## 147 IIRWMPR 1 0.003750094 Baja frecuencia
## 148 IIS 6 0.022500563 Baja frecuencia
## 149 IISDR 2 0.007500188 Baja frecuencia
## 150 IISER 22 0.082502063 Baja frecuencia
## 151 IISR 64 0.240006000 Baja frecuencia
## 152 IISSCR 54 0.202505063 Baja frecuencia
## 153 IISTR 243 0.911272782 Frecuencia media
## 154 IIW 18 0.067501688 Baja frecuencia
## 155 IIWBS 5 0.018750469 Baja frecuencia
## 156 IIWCR 4 0.015000375 Baja frecuencia
## 157 IIWD 3 0.011250281 Baja frecuencia
## 158 IIWFC 19 0.071251781 Baja frecuencia
## 159 ILES 26 0.097502438 Baja frecuencia
## 160 ILRS 1 0.003750094 Baja frecuencia
## 161 IN 7 0.026250656 Baja frecuencia
## 162 INCP 12 0.045001125 Baja frecuencia
## 163 INDR 1 0.003750094 Baja frecuencia
## 164 INDS 1 0.003750094 Baja frecuencia
## 165 INER 17 0.063751594 Baja frecuencia
## 166 INR 48 0.180004500 Baja frecuencia
## 167 INTR 203 0.761269032 Frecuencia media
## 168 INWSCR 14 0.052501313 Baja frecuencia
## 169 IOTRT 9 0.033750844 Baja frecuencia
## 170 IR 11 0.041251031 Baja frecuencia
## 171 IRGR 16 0.060001500 Baja frecuencia
## 172 IRSER 1 0.003750094 Baja frecuencia
## 173 ISDS 1 0.003750094 Baja frecuencia
## 174 ISER 17 0.063751594 Baja frecuencia
## 175 ISR 74 0.277506938 Baja frecuencia
## 176 ISSCR 77 0.288757219 Baja frecuencia
## 177 ISTR 211 0.791269782 Frecuencia media
## 178 IV 1272 4.770119253 Alta frecuencia
## 179 IVC 18 0.067501688 Baja frecuencia
## 180 IVCC 6 0.022500563 Baja frecuencia
## 181 IVE 45 0.168754219 Baja frecuencia
## 182 IVEBR 4 0.015000375 Baja frecuencia
## 183 IVED 3 0.011250281 Baja frecuencia
## 184 IVFTRT 1 0.003750094 Baja frecuencia
## 185 IVG 1 0.003750094 Baja frecuencia
## 186 IVGC 3 0.011250281 Baja frecuencia
## 187 IVGCH 15 0.056251406 Baja frecuencia
## 188 IVN 2 0.007500188 Baja frecuencia
## 189 IVNER 13 0.048751219 Baja frecuencia
## 190 IVNR 4 0.015000375 Baja frecuencia
## 191 IVNTR 1 0.003750094 Baja frecuencia
## 192 IVS 111 0.416260407 Frecuencia media
## 193 IVSER 47 0.176254406 Baja frecuencia
## 194 IVSFCWD 5 0.018750469 Baja frecuencia
## 195 IVSR 77 0.288757219 Baja frecuencia
## 196 IVSSCR 47 0.176254406 Baja frecuencia
## 197 IVSTR 11 0.041251031 Baja frecuencia
## 198 IVW 7 0.026250656 Baja frecuencia
## 199 IVWBS 3 0.011250281 Baja frecuencia
## 200 IVWCR 9 0.033750844 Baja frecuencia
## 201 IVWD 2 0.007500188 Baja frecuencia
## 202 IVWGL 1 0.003750094 Baja frecuencia
## 203 IW 13 0.048751219 Baja frecuencia
## 204 IWBR 1 0.003750094 Baja frecuencia
## 205 IWBS 1 0.003750094 Baja frecuencia
## 206 IWCR 10 0.037500938 Baja frecuencia
## 207 IWD 3 0.011250281 Baja frecuencia
## 208 IWFC 20 0.075001875 Baja frecuencia
## 209 IWPLR 8 0.030000750 Baja frecuencia
## 210 IX 836 3.135078377 Alta frecuencia
## 211 IXE 3 0.011250281 Baja frecuencia
## 212 IXED 1 0.003750094 Baja frecuencia
## 213 IXGC 10 0.037500938 Baja frecuencia
## 214 IXLES 5 0.018750469 Baja frecuencia
## 215 IXNR 26 0.097502438 Baja frecuencia
## 216 IXWD 5 0.018750469 Baja frecuencia
## 217 JTESGR 30 0.112502813 Baja frecuencia
## 218 JTWSGR 55 0.206255156 Baja frecuencia
## 219 KI 7 0.026250656 Baja frecuencia
## 220 KT 1 0.003750094 Baja frecuencia
## 221 LP 1 0.003750094 Baja frecuencia
## 222 LR 9 0.033750844 Baja frecuencia
## 223 LRCN 4 0.015000375 Baja frecuencia
## 224 LRCS 2 0.007500188 Baja frecuencia
## 225 LRE 6 0.022500563 Baja frecuencia
## 226 LROE 2 0.007500188 Baja frecuencia
## 227 LROW 2 0.007500188 Baja frecuencia
## 228 LRW 4 0.015000375 Baja frecuencia
## 229 LT 1 0.003750094 Baja frecuencia
## 230 MC 8 0.030000750 Baja frecuencia
## 231 MT 1 0.003750094 Baja frecuencia
## 232 NB 1 0.003750094 Baja frecuencia
## 233 NBC 23 0.086252156 Baja frecuencia
## 234 NBR 4 0.015000375 Baja frecuencia
## 235 NFR 51 0.191254781 Baja frecuencia
## 236 NG 2 0.007500188 Baja frecuencia
## 237 NLR 22 0.082502063 Baja frecuencia
## 238 NMR 44 0.165004125 Baja frecuencia
## 239 NSTREB 5 0.018750469 Baja frecuencia
## 240 NSTRNB 1 0.003750094 Baja frecuencia
## 241 NT 102 0.382509563 Frecuencia media
## 242 NTR 80 0.300007500 Baja frecuencia
## 243 OT 8 0.030000750 Baja frecuencia
## 244 R1GR 33 0.123753094 Baja frecuencia
## 245 R2 10 0.037500938 Baja frecuencia
## 246 R3 1 0.003750094 Baja frecuencia
## 247 R4 11 0.041251031 Baja frecuencia
## 248 R5 5 0.018750469 Baja frecuencia
## 249 RE 2 0.007500188 Baja frecuencia
## 250 RECTR 6 0.022500563 Baja frecuencia
## 251 REPR 26 0.097502438 Baja frecuencia
## 252 RF 40 0.150003750 Baja frecuencia
## 253 RI 1 0.003750094 Baja frecuencia
## 254 RII 3 0.011250281 Baja frecuencia
## 255 RIII 2 0.007500188 Baja frecuencia
## 256 RIIIS 11 0.041251031 Baja frecuencia
## 257 RIIN 35 0.131253281 Baja frecuencia
## 258 RIIS 29 0.108752719 Baja frecuencia
## 259 RIN 38 0.142503563 Baja frecuencia
## 260 RIS 31 0.116252906 Baja frecuencia
## 261 RIV 4 0.015000375 Baja frecuencia
## 262 RIVS 13 0.048751219 Baja frecuencia
## 263 RP 1 0.003750094 Baja frecuencia
## 264 RP338 9 0.033750844 Baja frecuencia
## 265 RP35 7 0.026250656 Baja frecuencia
## 266 RP384 6 0.022500563 Baja frecuencia
## 267 RP459 1 0.003750094 Baja frecuencia
## 268 RR 297 1.113777844 Frecuencia media
## 269 RV 1 0.003750094 Baja frecuencia
## 270 RVS 6 0.022500563 Baja frecuencia
## 271 RW 4 0.015000375 Baja frecuencia
## 272 RWCTR 6 0.022500563 Baja frecuencia
## 273 RWPR 12 0.045001125 Baja frecuencia
## 274 SBC 13 0.048751219 Baja frecuencia
## 275 SECP 3 0.011250281 Baja frecuencia
## 276 SFR 23 0.086252156 Baja frecuencia
## 277 SG 2 0.007500188 Baja frecuencia
## 278 SMR 48 0.180004500 Baja frecuencia
## 279 SR 1 0.003750094 Baja frecuencia
## 280 SRSUR 2 0.007500188 Baja frecuencia
## 281 SSTREB 2 0.007500188 Baja frecuencia
## 282 SSTRNB 2 0.007500188 Baja frecuencia
## 283 STR 80 0.300007500 Baja frecuencia
## 284 TR 9 0.033750844 Baja frecuencia
## 285 TRNLR 2 0.007500188 Baja frecuencia
## 286 TRS 9 0.033750844 Baja frecuencia
## 287 TRW 393 1.473786845 Frecuencia media
## 288 V 1474 5.527638191 Alta frecuencia
## 289 VE 16 0.060001500 Baja frecuencia
## 290 VEBR 1 0.003750094 Baja frecuencia
## 291 VECR 1 0.003750094 Baja frecuencia
## 292 VED 2 0.007500188 Baja frecuencia
## 293 VFTRT 6 0.022500563 Baja frecuencia
## 294 VGC 7 0.026250656 Baja frecuencia
## 295 VI 1335 5.006375159 Alta frecuencia
## 296 VIE 2 0.007500188 Baja frecuencia
## 297 VIED 2 0.007500188 Baja frecuencia
## 298 VIGC 19 0.071251781 Baja frecuencia
## 299 VII 1159 4.346358659 Alta frecuencia
## 300 VIIE 7 0.026250656 Baja frecuencia
## 301 VIIEBER 1 0.003750094 Baja frecuencia
## 302 VIIED 1 0.003750094 Baja frecuencia
## 303 VIIG 2 0.007500188 Baja frecuencia
## 304 VIIGC 14 0.052501313 Baja frecuencia
## 305 VIII 851 3.191329783 Alta frecuencia
## 306 VIII WBR 1 0.003750094 Baja frecuencia
## 307 VIIIE 12 0.045001125 Baja frecuencia
## 308 VIIIGC 11 0.041251031 Baja frecuencia
## 309 VIIILES 4 0.015000375 Baja frecuencia
## 310 VIIINR 15 0.056251406 Baja frecuencia
## 311 VIIIW 1 0.003750094 Baja frecuencia
## 312 VIIIWCR 1 0.003750094 Baja frecuencia
## 313 VIIIWD 11 0.041251031 Baja frecuencia
## 314 VIILES 3 0.011250281 Baja frecuencia
## 315 VIINER 2 0.007500188 Baja frecuencia
## 316 VIINR 16 0.060001500 Baja frecuencia
## 317 VIIW 1 0.003750094 Baja frecuencia
## 318 VILES 6 0.022500563 Baja frecuencia
## 319 VIN 1 0.003750094 Baja frecuencia
## 320 VINER 1 0.003750094 Baja frecuencia
## 321 VINR 18 0.067501688 Baja frecuencia
## 322 VIRR 2 0.007500188 Baja frecuencia
## 323 VISER 59 0.221255531 Baja frecuencia
## 324 VIW 3 0.011250281 Baja frecuencia
## 325 VIWCR 3 0.011250281 Baja frecuencia
## 326 VN 1 0.003750094 Baja frecuencia
## 327 VNCA 15 0.056251406 Baja frecuencia
## 328 VNER 14 0.052501313 Baja frecuencia
## 329 VNR 20 0.075001875 Baja frecuencia
## 330 VNTR 5 0.018750469 Baja frecuencia
## 331 VOR 1 0.003750094 Baja frecuencia
## 332 VS 115 0.431260782 Frecuencia media
## 333 VSA 1 0.003750094 Baja frecuencia
## 334 VSER 115 0.431260782 Frecuencia media
## 335 VSR 40 0.150003750 Baja frecuencia
## 336 VW 5 0.018750469 Baja frecuencia
## 337 VWBS 5 0.018750469 Baja frecuencia
## 338 VWD 1 0.003750094 Baja frecuencia
## 339 WBDR 1 0.003750094 Baja frecuencia
## 340 WBT 1 0.003750094 Baja frecuencia
## 341 WNR 22 0.082502063 Baja frecuencia
## 342 WPR 3 0.011250281 Baja frecuencia
## 343 WTR 2 0.007500188 Baja frecuencia
## 344 X 996 3.735093377 Alta frecuencia
## 345 XE 5 0.018750469 Baja frecuencia
## 346 XGC 15 0.056251406 Baja frecuencia
## 347 XI 932 3.495087377 Alta frecuencia
## 348 XIE 5 0.018750469 Baja frecuencia
## 349 XIGC 8 0.030000750 Baja frecuencia
## 350 XII 931 3.491337283 Alta frecuencia
## 351 XIIE 3 0.011250281 Baja frecuencia
## 352 XIIG 1 0.003750094 Baja frecuencia
## 353 XIIGC 8 0.030000750 Baja frecuencia
## 354 XIII 546 2.047551189 Frecuencia media
## 355 XIIIE 6 0.022500563 Baja frecuencia
## 356 XIIIGC 39 0.146253656 Baja frecuencia
## 357 XIIILES 1 0.003750094 Baja frecuencia
## 358 XIIINR 39 0.146253656 Baja frecuencia
## 359 XIIIWD 3 0.011250281 Baja frecuencia
## 360 XIILES 1 0.003750094 Baja frecuencia
## 361 XIINR 50 0.187504688 Baja frecuencia
## 362 XIIW 1 0.003750094 Baja frecuencia
## 363 XIIWD 1 0.003750094 Baja frecuencia
## 364 XILES 2 0.007500188 Baja frecuencia
## 365 XINR 39 0.146253656 Baja frecuencia
## 366 XIV 458 1.717542939 Frecuencia media
## 367 XIVE 3 0.011250281 Baja frecuencia
## 368 XIVED 1 0.003750094 Baja frecuencia
## 369 XIVGC 14 0.052501313 Baja frecuencia
## 370 XIVNR 33 0.123753094 Baja frecuencia
## 371 XIWD 2 0.007500188 Baja frecuencia
## 372 XIX 11 0.041251031 Baja frecuencia
## 373 XNR 35 0.131253281 Baja frecuencia
## 374 XV 147 0.551263782 Frecuencia media
## 375 XVE 3 0.011250281 Baja frecuencia
## 376 XVI 35 0.131253281 Baja frecuencia
## 377 XVIE 1 0.003750094 Baja frecuencia
## 378 XVII 10 0.037500938 Baja frecuencia
## 379 XVIII 10 0.037500938 Baja frecuencia
## 380 XVINR 19 0.071251781 Baja frecuencia
## 381 XVNR 42 0.157503938 Baja frecuencia
## 382 XWD 3 0.011250281 Baja frecuencia
## 383 XX 4 0.015000375 Baja frecuencia
## 384 XXI 9 0.033750844 Baja frecuencia
## 385 XXII 12 0.045001125 Baja frecuencia
## 386 XXIV 1 0.003750094 Baja frecuencia
## 387 XXV 2 0.007500188 Baja frecuencia
## 388 XXVI 2 0.007500188 Baja frecuencia
## 389 XXVIII 1 0.003750094 Baja frecuencia
## 390 YT 51 0.191254781 Baja frecuenciaGDF de agrupación por Concesión versión local
barplot(Grupo_Concesion,
main = "Gráfica No.21:
Distribución por grupos de concesiones de Ontario según su actividad",
xlab = "Concesión",
ylab = "Cantidad",
cex.names = 0.7,
col = "orange")
GDF de agrupación por Concesión versión global
barplot(Grupo_Concesion,
main="Gráfica No.22:
Distrubución por grupos de concesiones de Ontario según su actividad",
xlab="Concesión",
ylab="Cantidad",
col="orange",
cex.names = 0.7,
ylim=c(0,sum(Grupo_Concesion)))
GDF porcentaje de agrupación de concesión
hi_grupoConcesion<-Grupo_Concesion/sum(Grupo_Concesion)
hi_grupoConcesion<-round((hi_grupoConcesion*100),2)
barplot(hi_grupoConcesion,
main="Gráfica No.23:
Distrubución de grupos de concesiones de Ontario según su actividad",
xlab="Municipios",
ylab="Cantidad",
col="orange",
cex.axis = 0.7,
cex.names = 0.7,
ylim=c(0,100))
Diagrama círcular de agrupación de condados
etiquetaConcesion<-paste(hi_grupoConcesion,"%")
pie(hi_grupoConcesion, main="Gráfica No.24:
Distribución de grupos de concesiones de Ontario según su actividad",
radius = 1,
col=colores<-c(rev(heat.colors(3))),
labels = etiquetaMunicipio)
legend("bottomright", legend= names(hi_grupoMunicipio),
title="Leyenda",
cex = 0.75,
fill=colores<-c(rev(heat.colors(3))))
Variable Mes de inicio de perforación
Tabla de distribución de frecuencias
inicio_perforacion_mes<-datos2$SPUD_DATE_MES
inicio_perforacion_mes<-na.omit(inicio_perforacion_mes)
TDFinicio_perforacion_mes<-table(inicio_perforacion_mes)
Tabla_inicio_perforacion_mes<-as.data.frame(TDFinicio_perforacion_mes)
hi_inicio_perforacion_mes<- (Tabla_inicio_perforacion_mes$Freq/sum(Tabla_inicio_perforacion_mes$Freq)*100)
sum(hi_inicio_perforacion_mes)## [1] 100Tabla_inicio_perforacion_mes<-data.frame(Tabla_inicio_perforacion_mes, round(hi_inicio_perforacion_mes,2))
Niasc_inicio_perforacion_mes<-cumsum(Tabla_inicio_perforacion_mes$Freq)
Hiasc_inicio_perforacion_mes<-cumsum(hi_inicio_perforacion_mes)
Nidsc_inicio_perforacion_mes<-rev(cumsum(rev(Tabla_inicio_perforacion_mes$Freq)))
Hidsc_inicio_perforacion_mes<-rev(cumsum(rev(hi_inicio_perforacion_mes)))
Tabla_inicio_perforacion_mes_Final<-data.frame(Tabla_inicio_perforacion_mes$inicio_perforacion_mes,
Tabla_inicio_perforacion_mes$Freq,
round(hi_inicio_perforacion_mes,2),
Niasc_inicio_perforacion_mes,
round(Hiasc_inicio_perforacion_mes,2),
Nidsc_inicio_perforacion_mes,
round(Hidsc_inicio_perforacion_mes,2))
colnames(Tabla_inicio_perforacion_mes_Final)<-c("Mes de inicio de perforación","ni","hi","Ni asc","Hi asc","Ni dsc","Hi dsc")
Tabla_inicio_perforacion_mes_Final## Mes de inicio de perforación ni hi Ni asc Hi asc Ni dsc Hi dsc
## 1 1 7411 27.80 7411 27.80 26663 100.00
## 2 2 1288 4.83 8699 32.63 19252 72.20
## 3 3 1365 5.12 10064 37.75 17964 67.37
## 4 4 1199 4.50 11263 42.24 16599 62.25
## 5 5 2048 7.68 13311 49.92 15400 57.76
## 6 6 2250 8.44 15561 58.36 13352 50.08
## 7 7 2252 8.45 17813 66.81 11102 41.64
## 8 8 2166 8.12 19979 74.93 8850 33.19
## 9 9 2022 7.58 22001 82.52 6684 25.07
## 10 10 1780 6.68 23781 89.19 4662 17.48
## 11 11 1552 5.82 25333 95.01 2882 10.81
## 12 12 1330 4.99 26663 100.00 1330 4.99GDF
mes_perforacion<-as.numeric(Tabla_inicio_perforacion_mes_Final$ni)
barplot(mes_perforacion,
main="Gráfica No.1:
Distribución de meses de inicio de perforación",
xlab="Meses del año",
ylab="Cantidad",
col="blue",
names.arg=Tabla_inicio_perforacion_mes$inicio_perforacion_mes,
ylim=c(0,length(inicio_perforacion_mes)))
boxplot(inicio_perforacion_mes,
horizontal=T,
col="darkblue",
main="Gráfica No.2:
Distribución de meses de inicio de perforación",
xlab = "Meses del año")
# OJIVAS ASCENDENTES
x_ni_asc<- as.numeric(Tabla_inicio_perforacion_mes_Final$`Mes de inicio de perforación`)
y_ni_asc<- Tabla_inicio_perforacion_mes_Final$`Ni asc`
plot(x_ni_asc,
y_ni_asc,
type="p",
main="Gráfica No.30:
Ojivas asc de meses de inicio de perforación",
ylab="Cantidad",
xlab = "Meses del año",
col="blue")
#OJIVAS DESCENDENTES
x_ni_dsc<- as.numeric(Tabla_inicio_perforacion_mes_Final$`Mes de inicio de perforación`)
y_ni_dsc<- Tabla_inicio_perforacion_mes_Final$`Ni dsc`
plot(x_ni_dsc,y_ni_dsc,
type="p",
main="Gráfica No.31:
Ojivas dsc de meses de inicio de perforación",
ylab="Cantidad",
xlab = "Meses del año",
col="red")
#OJIVAS
plot(x_ni_asc,y_ni_asc,
type="p",
main="Gráfica No.32:
Ojivas de meses de inicio de perforación",
ylab="Cantidad",
xlab = "Meses del año",
col="blue")
lines(x_ni_dsc,Tabla_inicio_perforacion_mes_Final$`Ni dsc`, col="red", type="p")
Variable Mes profundidad alcanzada
profundidad_total_mes<-datos2$TOTAL_DEPTH_REACHED_DATE_MES
profundidad_total_mes<-na.omit(profundidad_total_mes)
TDFprofundidad_total_mes<-table(profundidad_total_mes)
Tabla_profundidad_total_mes<-as.data.frame(TDFprofundidad_total_mes)
hi_profundidad_total_mes<- (Tabla_profundidad_total_mes$Freq/sum(Tabla_profundidad_total_mes$Freq)*100)
Tabla_profundidad_total_mes<-data.frame(Tabla_profundidad_total_mes, round(hi_profundidad_total_mes,2))
Niasc_profundidad_total_mes<-cumsum(Tabla_profundidad_total_mes$Freq)
Hiasc_profundidad_total_mes<-cumsum(hi_profundidad_total_mes)
Nidsc_profundidad_total_mes<-rev(cumsum(rev(Tabla_profundidad_total_mes$Freq)))
Hidsc_profundidad_total_mes<-rev(cumsum(rev(hi_profundidad_total_mes)))
Tabla_profundidad_total_mes_Final<-data.frame(Tabla_profundidad_total_mes$profundidad_total_mes,
Tabla_profundidad_total_mes$Freq,
round(hi_profundidad_total_mes,2),
Niasc_profundidad_total_mes,
round(Hiasc_profundidad_total_mes,2),
Nidsc_profundidad_total_mes,
round(Hidsc_profundidad_total_mes,2))
colnames(Tabla_profundidad_total_mes_Final)<-c("Mes de profundidad alcanzada","ni","hi","Ni asc","Hi asc","Ni dsc","Hi dsc")
Tabla_profundidad_total_mes_Final## Mes de profundidad alcanzada ni hi Ni asc Hi asc Ni dsc Hi dsc
## 1 1 7411 27.80 7411 27.80 26663 100.00
## 2 2 1288 4.83 8699 32.63 19252 72.20
## 3 3 1365 5.12 10064 37.75 17964 67.37
## 4 4 1199 4.50 11263 42.24 16599 62.25
## 5 5 2048 7.68 13311 49.92 15400 57.76
## 6 6 2250 8.44 15561 58.36 13352 50.08
## 7 7 2252 8.45 17813 66.81 11102 41.64
## 8 8 2166 8.12 19979 74.93 8850 33.19
## 9 9 2022 7.58 22001 82.52 6684 25.07
## 10 10 1780 6.68 23781 89.19 4662 17.48
## 11 11 1552 5.82 25333 95.01 2882 10.81
## 12 12 1330 4.99 26663 100.00 1330 4.99#GDF
mes_profundidad<-as.numeric(Tabla_profundidad_total_mes_Final$ni)
barplot(mes_profundidad,
main="Gráfico No.1:
Distribución de meses de de profundidad total alcanzada",
xlab="Meses del año",
ylab="Cantidad",
col="red",
names.arg=Tabla_profundidad_total_mes$profundidad_total_mes,
ylim=c(0,length(profundidad_total_mes)))
boxplot(profundidad_total_mes,
horizontal=T,
col="red",
main="Gráfico No.2:
Distribución de meses de profundidad total alcanzada",
xlab = "Meses del año")
# OJIVAS ASCENDENTES
x_ni_asc_t<- as.numeric(Tabla_profundidad_total_mes_Final$`Mes de profundidad alcanzada`)
y_ni_asc_t<- Tabla_profundidad_total_mes_Final$`Ni asc`
plot(x_ni_asc_t,
y_ni_asc_t,
type="p",
main="Gráfica No.3 Ojivas asc de meses de profundidad total alcanzada",
ylab="Cantidad", xlab = "Meses del año", col="blue")
#OJIVAS DESCENDENTES
x_ni_dsc_t<- as.numeric(Tabla_profundidad_total_mes_Final$`Mes de profundidad alcanzada`)
y_ni_dsc_t<- Tabla_inicio_perforacion_mes_Final$`Ni dsc`
plot(x_ni_dsc_t,
y_ni_dsc_t,
type="p",
main="Gráfica No.4:
Ojivas dsc de meses de profundidad total alcanzada",
ylab="Cantidad",
xlab = "Meses del año",
col="red")
#OJIVAS
plot(x_ni_asc_t,y_ni_asc_t,
type="p",
main="Gráfica No.5:
Ojivas de meses de profundidad total alcanzada",
ylab="Cantidad",
xlab = "Meses del año",
col="blue")
lines(x_ni_dsc_t,
Tabla_profundidad_total_mes_Final$`Ni dsc`,
col="red",
type="p")
Variables Continuas
# VARIABLE X
X<-datos$X
X<-na.omit(X)
# GDF
min(X)## [1] -88.79041max(X)## [1] -74.36656R_X<-max(X)-min(X)
k_X<-1+(3.3*log10(length(X)))
k_X<-floor(k_X)
A_X<-R_X/k_X
HistoSturges_X<-hist(X, main="Gráfica No.1:
Distribución de coordenadas X",
breaks = seq(min(X),max(X),A_X),
xlab="Coordenadas X",
ylab="Cantidad",
col="red")
HistoSturges_X$counts## [1] 2 3 6 5 12 241 7756 4871 4377 8048 1180 24 36 58 47HistogramaX<-hist(X, main="Gráfica No. :",
xlab="Coordenadas X",
ylab="Cantidad",
col="red")
limites_X<-HistogramaX$breaks
liminf_X<-limites_X[1:15]
limsup_X<-limites_X[2:16]
MC_X<-(liminf_X+limsup_X)/2
# TABLA COORDENADAS X
ni_X<-c()
for(i in 1:k_X) {
if(1==15)
ni_X[i]<- length(subset(X, X>=liminf_X[i] & X<=limsup_X[i]))
else
ni_X[i]<- length(subset(X, X>=liminf_X[i]& X<limsup_X[i]))
}
hi_X<-(ni_X/length(X))*100
Niasc_X<-cumsum(ni_X)
Hiasc_X<-cumsum(hi_X)
Nidsc_X<-rev(cumsum(rev(ni_X)))
Hidsc_X<-rev(cumsum(rev(hi_X)))
Tabla_X_Final<-data.frame(liminf_X,
limsup_X,
MC_X,
ni_X,
round(hi_X,2),
Niasc_X,
round(Hiasc_X,2),
Nidsc_X,
round(Hidsc_X,2))
colnames(Tabla_X_Final)<-c("Xi-1","Xi+1","Xi (MC)","ni","hi","Ni ↓","Hi ↓","Ni ↑","Hi ↑")
Tabla_X_Final## Xi-1 Xi+1 Xi (MC) ni hi Ni ↓ Hi ↓ Ni ↑ Hi ↑
## 1 -89 -88 -88.5 2 0.01 2 0.01 26666 100.00
## 2 -88 -87 -87.5 3 0.01 5 0.02 26664 99.99
## 3 -87 -86 -86.5 6 0.02 11 0.04 26661 99.98
## 4 -86 -85 -85.5 5 0.02 16 0.06 26655 99.96
## 5 -85 -84 -84.5 12 0.05 28 0.11 26650 99.94
## 6 -84 -83 -83.5 261 0.98 289 1.08 26638 99.89
## 7 -83 -82 -82.5 8224 30.84 8513 31.92 26377 98.92
## 8 -82 -81 -81.5 4530 16.99 13043 48.91 18153 68.08
## 9 -81 -80 -80.5 5646 21.17 18689 70.09 13623 51.09
## 10 -80 -79 -79.5 7630 28.61 26319 98.70 7977 29.91
## 11 -79 -78 -78.5 182 0.68 26501 99.38 347 1.30
## 12 -78 -77 -77.5 57 0.21 26558 99.59 165 0.62
## 13 -77 -76 -76.5 6 0.02 26564 99.62 108 0.41
## 14 -76 -75 -75.5 67 0.25 26631 99.87 102 0.38
## 15 -75 -74 -74.5 35 0.13 26666 100.00 35 0.13# OJIVAS CONTINUAS COORDENADAS X
plot(limsup_X, Niasc_X, col="blue", type="b",
main="Gráfica No.2:
Distribución de coordenadas X",
ylab = "Cantidad",
xlab = "Coordenadas X",)
lines(liminf_X, Nidsc_X, col="red", type="b")
# OJIVAS CONTINUAS PORCENTAJE COORDENADAS X
plot(limsup_X, Hiasc_X, col="blue", type="b",
main="Gráfica No.3:
Distribución de coordenadas X",
ylab = "Porcentaje",
xlab = "Coordenadas X",)
lines(liminf_X, Hidsc_X, col="red", type="b")
# GDF CAJA Y BIGOTES
boxplot(X,
horizontal=T,
col="darkgreen",
xlab= "Coordenadas X",
main="Gráfica No. 4:
Distribución de coordenadas X")
# VARIABLE Y
Y<-datos$Y
Y<-na.omit(Y)
# GDF
min(Y)## [1] 41.68497max(Y)## [1] 58.50058R_Y<-max(Y)-min(Y)
k_Y<-1+(3.3*log10(length(Y)))
k_Y<-floor(k_Y)
A_Y<-R_Y/k_Y
HistoSturges_Y<-hist(Y,
main="Gráfica No.1:
Distribución de coordenadas Y",
breaks = seq(min(Y),max(Y),A_Y),
xlab="Coordenadas Y",
ylab="Cantidad",
col="red")
HistogramaY<-hist(Y,
main="Gráfica No.:",
xlab="Coordenadas Y",
ylab="Cantidad",
col="red")
HistoSturges_Y$counts## [1] 11942 13867 387 255 1 1 3 118 66 16 0 7
## [13] 0 2 1limites_Y<-HistogramaY$breaks
liminf_Y<-limites_Y[1:15]
limsup_Y<-limites_Y[2:16]
MC_Y<-(liminf_Y+limsup_Y)/2
# TABLA COORDENADAS Y
ni_Y<-c()
for(i in 1:k_Y) {
if(1==k_Y)
ni_Y[i]<- length(subset(Y, Y>=liminf_Y[i] & Y<=limsup_Y[i]))
else
ni_Y[i]<- length(subset(Y, Y>=liminf_Y[i]& Y<limsup_Y[i]))
}
hi_Y<-(ni_Y/length(Y))*100
Niasc_Y<-cumsum(ni_Y)
Hiasc_Y<-cumsum(hi_Y)
Nidsc_Y<-rev(cumsum(rev(ni_Y)))
Hidsc_Y<-rev(cumsum(rev(hi_Y)))
Tabla_Y_Final<-data.frame(liminf_Y,
limsup_Y,
MC_Y,
ni_Y,
round(hi_Y,2),
Niasc_Y,
round(Hiasc_Y,2),
Nidsc_Y,
round(Hidsc_Y,2))
colnames(Tabla_Y_Final)<-c("Xi-1","Xi+1","Xi (MC)","ni","hi","Ni ↓","Hi ↓","Ni ↑","Hi ↑")
Tabla_Y_Final## Xi-1 Xi+1 Xi (MC) ni hi Ni ↓ Hi ↓ Ni ↑ Hi ↑
## 1 41 42 41.5 128 0.48 128 0.48 26663 99.99
## 2 42 43 42.5 22335 83.76 22463 84.24 26535 99.51
## 3 43 44 43.5 3379 12.67 25842 96.91 4200 15.75
## 4 44 45 44.5 352 1.32 26194 98.23 821 3.08
## 5 45 46 45.5 255 0.96 26449 99.19 469 1.76
## 6 46 47 46.5 3 0.01 26452 99.20 214 0.80
## 7 47 48 47.5 1 0.00 26453 99.20 211 0.79
## 8 48 49 48.5 0 0.00 26453 99.20 210 0.79
## 9 49 50 49.5 3 0.01 26456 99.21 210 0.79
## 10 50 51 50.5 154 0.58 26610 99.79 207 0.78
## 11 51 52 51.5 35 0.13 26645 99.92 53 0.20
## 12 52 53 52.5 11 0.04 26656 99.96 18 0.07
## 13 53 54 53.5 0 0.00 26656 99.96 7 0.03
## 14 54 55 54.5 7 0.03 26663 99.99 7 0.03
## 15 55 56 55.5 0 0.00 26663 99.99 0 0.00# OJIVAS CONTINUAS
plot(limsup_Y, Niasc_Y, col="blue", type="b",
main="Gráfica No.2 : Distribución de coordenadas Y",
ylab = "Cantidad",
xlab = "Coordenadas Y",)
lines(liminf_Y, Nidsc_Y, col="red", type="b")
# OJIVAS CONTINUAS PORCENTAJE
plot(limsup_Y, Hiasc_Y, col="blue", type="b",
main="Gráfica No.3 :Distribución de coordenadas Y",
ylab = "Porcentaje",
xlab = "Coordenadas Y",)
lines(liminf_Y, Hidsc_Y, col="red", type="b")
# GDF CAJA Y BIGOTES
boxplot(Y,horizontal=T,
col="darkgreen",
xlab= "Coordenadas X",
main="Gráfica No.4:
Distribución de coordenadas Y")
# VARIABLE LATITUD EN SUPERFICIE
latitud_superficie<-datos$SURFACE_LATITUDE_83
latitud_superficie<-na.omit(latitud_superficie)
# GDF
min(latitud_superficie)## [1] 41.68496max(latitud_superficie)## [1] 58.50057R_latitud_superficie<-max(latitud_superficie)-min(latitud_superficie)
k_latitud_superficie<-1+(3.3*log10(length(latitud_superficie)))
k_latitud_superficie<-floor(k_latitud_superficie)
A_latitud_superficie<-R_latitud_superficie/k_latitud_superficie
HistoSturges_latitud_superficie<-hist(latitud_superficie,
main="Gráfica No.1:Distribución de Latitud en superficie",
breaks = seq(min(latitud_superficie),max(latitud_superficie),A_latitud_superficie),
xlab="Latitud en superficie", ylab="Cantidad", col="red")
Histograma_latitud_superficie<-hist(latitud_superficie, main="Gráfica No. :",
xlab="Latitud superficie", ylab="Cantidad", col="red")
limites_latitud_superficie<-Histograma_latitud_superficie$breaks
liminf_latitud_superficie<-limites_latitud_superficie[1:15]
limsup_latitud_superficie<-limites_latitud_superficie[2:16]
MC_latitud_superficie<-(liminf_latitud_superficie+limsup_latitud_superficie)/2
# TABLA LONGITUD EN SUPERFICIE
ni_latitud_superficie<-c()
for(i in 1:k_latitud_superficie) {
if(1==k_latitud_superficie)
ni_latitud_superficie[i]<- length(subset(latitud_superficie, latitud_superficie>=liminf_latitud_superficie[i] & latitud_superficie<=limsup_latitud_superficie[i]))
else
ni_latitud_superficie[i]<- length(subset(latitud_superficie, latitud_superficie>=liminf_latitud_superficie[i]& latitud_superficie<limsup_latitud_superficie[i]))
}
hi_latitud_superficie<-(ni_latitud_superficie/length(latitud_superficie))*100
Niasc_latitud_superficie<-cumsum(ni_latitud_superficie)
Hiasc_latitud_superficie<-cumsum(hi_latitud_superficie)
Nidsc_latitud_superficie<-rev(cumsum(rev(ni_latitud_superficie)))
Hidsc_latitud_superficie<-rev(cumsum(rev(hi_latitud_superficie)))
Tabla_latitud_superficie_Final<-data.frame(liminf_latitud_superficie,
limsup_latitud_superficie,
MC_latitud_superficie,
ni_latitud_superficie,
round(hi_latitud_superficie,2),
Niasc_latitud_superficie,
round(Hiasc_latitud_superficie,2),
Nidsc_latitud_superficie,
round(Hidsc_latitud_superficie,2))
colnames(Tabla_latitud_superficie_Final)<-c("Xi-1","Xi+1","Xi (MC)","ni","hi","Ni ↓","Hi ↓","Ni ↑","Hi ↑")
Tabla_latitud_superficie_Final## Xi-1 Xi+1 Xi (MC) ni hi Ni ↓ Hi ↓ Ni ↑ Hi ↑
## 1 41 42 41.5 128 0.48 128 0.48 26663 99.99
## 2 42 43 42.5 22335 83.76 22463 84.24 26535 99.51
## 3 43 44 43.5 3379 12.67 25842 96.91 4200 15.75
## 4 44 45 44.5 352 1.32 26194 98.23 821 3.08
## 5 45 46 45.5 255 0.96 26449 99.19 469 1.76
## 6 46 47 46.5 3 0.01 26452 99.20 214 0.80
## 7 47 48 47.5 1 0.00 26453 99.20 211 0.79
## 8 48 49 48.5 0 0.00 26453 99.20 210 0.79
## 9 49 50 49.5 3 0.01 26456 99.21 210 0.79
## 10 50 51 50.5 154 0.58 26610 99.79 207 0.78
## 11 51 52 51.5 35 0.13 26645 99.92 53 0.20
## 12 52 53 52.5 11 0.04 26656 99.96 18 0.07
## 13 53 54 53.5 0 0.00 26656 99.96 7 0.03
## 14 54 55 54.5 7 0.03 26663 99.99 7 0.03
## 15 55 56 55.5 0 0.00 26663 99.99 0 0.00# OJIVAS CONTINUAS LATITUD EN SUPERFICIE
plot(limsup_latitud_superficie, Niasc_latitud_superficie, col="blue", type="b",
main="Gráfica No.54 : Distribución de Latitud en superficie",
ylab = "Cantidad",
xlab = "Latitud en superficie",)
lines(liminf_latitud_superficie, Nidsc_latitud_superficie, col="red", type="b")
# OJIVAS CONTINUAS PORCENTAJE LATITUD EN SUPERFICIE
plot(limsup_latitud_superficie, Hiasc_latitud_superficie, col="blue", type="b",
main="Gráfica No.55 : Distribución de Latitud en superficie",
ylab = "Porcentaje",
xlab = "Latitud en superficie",)
lines(liminf_latitud_superficie, Hidsc_latitud_superficie, col="red", type="b")
boxplot(latitud_superficie,horizontal=T,
col="darkgreen",
xlab= "Latitud en superficie",
main="Gráfica No.3:
Distribución de Latitud en superficie")