suppressWarnings(library(RODBC))
library(ggplot2)
library(ggpubr)
library(markdown)
library(shiny)
library(shinythemes)
library(tidyverse)
library(magrittr)
library(lubridate)
library(plotly)
library(xts)
library(dygraphs)
library(kableExtra)
library(knitr)
library("readxl")
library(rsconnect)
library(dplyr)
library(summarytools)
library(epiDisplay)
#library(leaflet)
library(haven)
library(epiDisplay)
library("readxl")
library(expss)
library(hrbrthemes)
library(viridis)
library(viridisLite)
library(DescTools)
library(roperators)
library(shinycssloaders)
library(writexl)
library(labelled)
library(tidyverse)
library(haven)
library(readr)
library(sjmisc)
library(WriteXLS)
w <- dataset2006[[6]] %>% attr('labels')
head(w,4)
## Iquique Camiña Colchane Huara
## 1101 1102 1103 1104
#### Aqui se obtienen los nombres de las comunas
a <- names(w)
head(a,4)
## [1] "Iquique" "Camiña" "Colchane" "Huara"
b<-w[[4]]
b
## [1] 1104
# head(a,5)
dataset06 <- read.csv('mydata2006_sub.csv')
a <- dataset06$comuna
b <- dataset06$t4
c <- dataset06$e1
d <- dataset06$sexo
cross_tab= table(a, b, c, d)
cross_tab
## , , c = no, d = hombre
##
## b
## a atacameño aymara diaguita mapuche
## camiña 0 14 0 0
## colchane 0 24 0 0
## huara 0 11 0 0
## iquique 0 0 0 0
## pica 0 0 0 0
## b
## a no pertenece a ningún pueblo indígena quechua rapa nui
## camiña 1 0 0
## colchane 0 1 0
## huara 9 0 0
## iquique 12 0 0
## pica 0 0 0
##
## , , c = sí, d = hombre
##
## b
## a atacameño aymara diaguita mapuche
## camiña 0 175 0 0
## colchane 0 135 0 0
## huara 0 105 0 4
## iquique 1 23 1 8
## pica 0 8 0 1
## b
## a no pertenece a ningún pueblo indígena quechua rapa nui
## camiña 30 0 1
## colchane 5 3 0
## huara 82 13 0
## iquique 392 3 0
## pica 10 0 0
##
## , , c = no, d = mujer
##
## b
## a atacameño aymara diaguita mapuche
## camiña 0 40 0 0
## colchane 0 53 0 0
## huara 0 24 0 0
## iquique 0 0 0 1
## pica 0 1 0 0
## b
## a no pertenece a ningún pueblo indígena quechua rapa nui
## camiña 2 0 0
## colchane 0 1 0
## huara 6 1 0
## iquique 9 0 0
## pica 0 0 0
##
## , , c = sí, d = mujer
##
## b
## a atacameño aymara diaguita mapuche
## camiña 0 143 0 1
## colchane 0 121 0 0
## huara 0 90 0 5
## iquique 3 19 0 8
## pica 0 15 0 0
## b
## a no pertenece a ningún pueblo indígena quechua rapa nui
## camiña 23 0 1
## colchane 8 2 0
## huara 77 14 0
## iquique 391 1 0
## pica 12 0 0
tabla <- as.data.frame(cross_tab)
#head(tabla,10)
head(tabla,10)
## a b c d Freq
## 1 camiña atacameño no hombre 0
## 2 colchane atacameño no hombre 0
## 3 huara atacameño no hombre 0
## 4 iquique atacameño no hombre 0
## 5 pica atacameño no hombre 0
## 6 camiña aymara no hombre 14
## 7 colchane aymara no hombre 24
## 8 huara aymara no hombre 11
## 9 iquique aymara no hombre 0
## 10 pica aymara no hombre 0
ppp <- data_frame()
## Warning: `data_frame()` is deprecated as of tibble 1.1.0.
## Please use `tibble()` instead.
## This warning is displayed once every 8 hours.
## Call `lifecycle::last_warnings()` to see where this warning was generated.
#vamos a recorrer todos los elementos de la tabla de frecuencia, si este valor es existente construimos un arreglo con la fila y lo arreglamos para que sea concistente en la comparacion del if
for( i in 1: nrow(tabla))
{
if(tabla[i,5]!=0)
{
lll<-tabla[i,]
lll
llll<-lll[,1]
llll
sentenceString <- toString(llll)
searchString <- ' '
replacementString <- ''
sentenceString = sub(searchString,replacementString,sentenceString)
sentenceString
for(j in 1: 300){
w <- dataset2006[[6]] %>% attr('labels')
ww<-names(w[j])
vv<-tolower(ww)
if(sentenceString==vv){
lll <- cbind(lll,w[[j]])
lll <- cbind(lll,"2006")
}
}
ppp <-rbind(ppp,lll)
}
}
ppp
## a b c d Freq w[[j]]
## 6 camiña aymara no hombre 14 1102
## 7 colchane aymara no hombre 24 1103
## 8 huara aymara no hombre 11 1104
## 21 camiña no pertenece a ningún pueblo indígena no hombre 1 1102
## 23 huara no pertenece a ningún pueblo indígena no hombre 9 1104
## 24 iquique no pertenece a ningún pueblo indígena no hombre 12 1101
## 27 colchane quechua no hombre 1 1103
## 39 iquique atacameño sí hombre 1 1101
## 41 camiña aymara sí hombre 175 1102
## 42 colchane aymara sí hombre 135 1103
## 43 huara aymara sí hombre 105 1104
## 44 iquique aymara sí hombre 23 1101
## 45 pica aymara sí hombre 8 1105
## 49 iquique diaguita sí hombre 1 1101
## 53 huara mapuche sí hombre 4 1104
## 54 iquique mapuche sí hombre 8 1101
## 55 pica mapuche sí hombre 1 1105
## 56 camiña no pertenece a ningún pueblo indígena sí hombre 30 1102
## 57 colchane no pertenece a ningún pueblo indígena sí hombre 5 1103
## 58 huara no pertenece a ningún pueblo indígena sí hombre 82 1104
## 59 iquique no pertenece a ningún pueblo indígena sí hombre 392 1101
## 60 pica no pertenece a ningún pueblo indígena sí hombre 10 1105
## 62 colchane quechua sí hombre 3 1103
## 63 huara quechua sí hombre 13 1104
## 64 iquique quechua sí hombre 3 1101
## 66 camiña rapa nui sí hombre 1 1102
## 76 camiña aymara no mujer 40 1102
## 77 colchane aymara no mujer 53 1103
## 78 huara aymara no mujer 24 1104
## 80 pica aymara no mujer 1 1105
## 89 iquique mapuche no mujer 1 1101
## 91 camiña no pertenece a ningún pueblo indígena no mujer 2 1102
## 93 huara no pertenece a ningún pueblo indígena no mujer 6 1104
## 94 iquique no pertenece a ningún pueblo indígena no mujer 9 1101
## 97 colchane quechua no mujer 1 1103
## 98 huara quechua no mujer 1 1104
## 109 iquique atacameño sí mujer 3 1101
## 111 camiña aymara sí mujer 143 1102
## 112 colchane aymara sí mujer 121 1103
## 113 huara aymara sí mujer 90 1104
## 114 iquique aymara sí mujer 19 1101
## 115 pica aymara sí mujer 15 1105
## 121 camiña mapuche sí mujer 1 1102
## 123 huara mapuche sí mujer 5 1104
## 124 iquique mapuche sí mujer 8 1101
## 126 camiña no pertenece a ningún pueblo indígena sí mujer 23 1102
## 127 colchane no pertenece a ningún pueblo indígena sí mujer 8 1103
## 128 huara no pertenece a ningún pueblo indígena sí mujer 77 1104
## 129 iquique no pertenece a ningún pueblo indígena sí mujer 391 1101
## 130 pica no pertenece a ningún pueblo indígena sí mujer 12 1105
## 132 colchane quechua sí mujer 2 1103
## 133 huara quechua sí mujer 14 1104
## 134 iquique quechua sí mujer 1 1101
## 136 camiña rapa nui sí mujer 1 1102
## "2006"
## 6 2006
## 7 2006
## 8 2006
## 21 2006
## 23 2006
## 24 2006
## 27 2006
## 39 2006
## 41 2006
## 42 2006
## 43 2006
## 44 2006
## 45 2006
## 49 2006
## 53 2006
## 54 2006
## 55 2006
## 56 2006
## 57 2006
## 58 2006
## 59 2006
## 60 2006
## 62 2006
## 63 2006
## 64 2006
## 66 2006
## 76 2006
## 77 2006
## 78 2006
## 80 2006
## 89 2006
## 91 2006
## 93 2006
## 94 2006
## 97 2006
## 98 2006
## 109 2006
## 111 2006
## 112 2006
## 113 2006
## 114 2006
## 115 2006
## 121 2006
## 123 2006
## 124 2006
## 126 2006
## 127 2006
## 128 2006
## 129 2006
## 130 2006
## 132 2006
## 133 2006
## 134 2006
## 136 2006