Obtener los datos

library(stringr)
library(stringi)
library(xtable)
library(RCurl)
library(pander)
library(plyr)
library(grid)
library(tidyr)
library(readxl)

### Obtener los datos de rendimiento 
uti <- read.csv('C:/Users/LVARGAS/Google Drive/datosBEM2016/utilidad2016.csv')

#unique(rend$Tipo.de.parcela..testigo.o.innovación.)

### Obtener los datos de caracteristicas de bitacora  y eliminar instituciones diferentes
caracteristicasBitNA <- read_excel('C:/Users/LVARGAS/Google Drive/datosBEM2016/completa2016.xlsx', '01_caracteristicas Bitácora')
valoresNA <- is.na(caracteristicasBitNA[,1])
caraBit <- caracteristicasBitNA[!valoresNA,]
caraBitIns <- caraBit[caraBit$`Nombre de la institución` != 'MASAGRO GUANAJUATO' & caraBit$`Nombre de la institución` != 'PIMAF' & caraBit$`Nombre de la institución` != 'MILPA PENINSULA DE YUCATAN', ] 
 
#names(caraBitIns)
caraBitInsNam <- caraBitIns[, c(1, 17)]

#names(rend)
utilidad <- uti[uti[,1] %in% caraBitIns$`ID de la bitácora (clave primaria)`,]

### Crear nueva variable con los datos de tipo de parcela

#utilidad$nuevoTipoPar <- ifelse(utilidad$nb.Tipo.Parcela == 'Parcela testigo', 'Productor sistema optimizado', ifelse((utilidad$nb.Tipo.Parcela == 'Parcela Área de extensión') & (utilidad$nb.Tipo.Parcela == 'Parcela innovación'), 'Innovacion', 'Innovacion'))

### Unir rendimiento con utilidad para obtener el cultivo

rend <- read.csv('C:/Users/LVARGAS/Google Drive/datosBEM2016/rendimiento2016.csv')

#names(rend)
rendimiento <- rend[rend[,1] %in% caraBitIns$`ID de la bitácora (clave primaria)`,]

rendimiento$nuevoTipoPar <- ifelse(rendimiento$Tipo.de.parcela..testigo.o.innovación. == 'Parcela testigo', 'Productor sistema optimizado', ifelse((rendimiento$Tipo.de.parcela..testigo.o.innovación. == 'Parcela Área de extensión') & (rendimiento$Tipo.de.parcela..testigo.o.innovación. == 'Parcela innovación'), 'Innovacion', 'Innovacion'))

names(rendimiento)[2] <- 'Tipobitacora'

utilidadS <- utilidad[, c(2, 28)]
names(utilidadS)[1] <- 'Tipobitacora'

utilidadRend <- merge(rendimiento, utilidadS, by = 'Tipobitacora', all.x=FALSE)

Procesar y obtener el resultado

vectorCultivo <- sort(unique(utilidadRend$Nombre.del.cultivo.cosechado))
        
##vectorCultivo <- 'Maiz'

for(cultivo in vectorCultivo){
          
        #cultivo <- 'Maiz'
  
        datosCultivo <- utilidadRend[utilidadRend$Nombre.del.cultivo.cosechado == cultivo,]
        
        vectorAo <- sort(unique(datosCultivo$Año))

        #vectorAo <- '2016'
        
        for(ao in vectorAo){
            
                #ao <- '2016'
            
                datosAo <- datosCultivo[datosCultivo$Año == ao,]
                
                VectorCiclo <- sort(unique(datosAo$Ciclo.agronómico))
                
                #VectorCiclo <- 'Primavera-Verano' 
            
                for(ciclo in VectorCiclo){
                  
                        #ciclo <- 'Primavera-Verano'
                  
                        datosCiclo <- datosAo[datosAo$Ciclo.agronómico == ciclo,]
                        
                        vectorTipo <- sort(unique(datosCiclo$tipoProduccion))
                        
                        #vectorTipo <- 'Riego'
                  
                        for(tipo in vectorTipo){
                                
                                #tipo <- 'Riego'
                
                                dataTipo <- datosCiclo[datosCiclo$tipoProduccion == tipo,]
                                names(dataTipo)
                                        
                                
                                dataTipoInn <- dataTipo[dataTipo$nuevoTipoPar == 'Innovacion',]
                                
                                numObInn <- dim(dataTipoInn)[1]

                                #----- Practicas optimizadas ------
                                
                                dataTipoOp <- dataTipo[dataTipo$nuevoTipoPar == 'Productor sistema optimizado',]
                                numObOp <- dim(dataTipoOp)[1]
                                
                                
                                if(numObInn > 0 & numObOp > 0){
                                        
                                        pandoc.header("--------------------------------------------", level = 2)
                                        
                                        etiqueta <- paste(cultivo, ao, ciclo, tipo)
                                        
                                        pandoc.header(etiqueta, level = 3)
                                        
                                                        
                                        #------------ Innovacion ----------
                                        
                                        #dataTipoInn <- dataTipo[dataTipo$nuevoTipoPar == 'Innovacion',]
                                                        
                                        #### Realizar el resumen de rendimiento y conteo de observaciones
                                        
                                        promediosUtiEstadoInn <- with(dataTipoInn, aggregate(RENTABILIDAD....ha., by = list(Estado, Año, Ciclo.agronómico, tipoProduccion), mean))
                                        
                                        names(promediosUtiEstadoInn) <- c('Estado', 'Año', 'Ciclo', 'Tipo producción', 'Utilidad parcela de innovación (Pesos/ha)')
                                        
                                        #Conteo observaciones rendimiento
                                        
                                        numUtiEstadoInn <- with(dataTipoInn, aggregate(RENTABILIDAD....ha., by = list(Estado, Año, Ciclo.agronómico, tipoProduccion), length))
                                        
                                        names(numUtiEstadoInn) <- c('Estado', 'Año', 'Ciclo', 'Tipo producción', 'Número de observaciones parcela de innovación')
                                        
                                        tablaUtiEstadoInn <- join_all(list(promediosUtiEstadoInn, numUtiEstadoInn), by = c ('Estado', 'Año', 'Ciclo', 'Tipo producción'), type = 'full')
                                        
                                        ####-------------------------------
                                        
                                        
                                        #----- Practicas optimizadas ------
                                        
                                        #dataTipoOp <- dataTipo[dataTipo$nuevoTipoPar == 'Productor sistema optimizado',]
                                                        
                                        #### Realizar el resumen de rendimiento y conteo de observaciones
                                        
                                        promediosUtiEstadoOp <- with(dataTipoOp, aggregate(RENTABILIDAD....ha., by = list(Estado, Año, Ciclo.agronómico, tipoProduccion), mean))
                                        
                                        names(promediosUtiEstadoOp) <- c('Estado', 'Año', 'Ciclo', 'Tipo producción', 'Utilidad en parcela del productor optimizando el sistema gracias a las capacitaciones de MasAgro (Pesos/ha)')
                                        
                                        #Conteo observaciones rendimiento
                                        numUtiEstadoOp <- with(dataTipoOp, aggregate(RENTABILIDAD....ha., by = list(Estado, Año, Ciclo.agronómico, tipoProduccion), length))
                                        
                                        names(numUtiEstadoOp) <- c('Estado', 'Año', 'Ciclo', 'Tipo producción', 'Número de observaciones en parcela del productor optimizando el sistema gracias a las capacitaciones de MasAgro')
                                        
                                        tablaUtiEstadoOp <- join_all(list(promediosUtiEstadoOp, numUtiEstadoOp), by = c ('Estado', 'Año', 'Ciclo', 'Tipo producción'), type = 'full')
                                        
                                        tablaUtiEstadoOpS <- tablaUtiEstadoOp[, c(1, 5, 6)]
                                        
                                        ####-------------------------------
                                        
                                        ### Unir las dos tablas  
                                        tablaUtiEstadoInnOp <- merge(tablaUtiEstadoInn,tablaUtiEstadoOpS, by = 'Estado', all.x=TRUE)
                                        
                                        ####-Calcular la tasa de variación-
                                                        
                                        tablaUtiEstadoInnOp$tasaVariacion <- ((tablaUtiEstadoInnOp$`Utilidad parcela de innovación (Pesos/ha)` - tablaUtiEstadoInnOp$`Utilidad en parcela del productor optimizando el sistema gracias a las capacitaciones de MasAgro (Pesos/ha)` ) / tablaUtiEstadoInnOp$`Utilidad en parcela del productor optimizando el sistema gracias a las capacitaciones de MasAgro (Pesos/ha)`) * 100
                                        
                                        names(tablaUtiEstadoInnOp)[9] <- 'Tasa de variación utilidad promedio parcela de innovación vs utilidad en parcela del productor optimizando el sistema gracias a las capacitaciones de MasAgro (%)'
                                        
                                        ####-------------------------------
                                        
                                        carpetaAlmacenamiento <- './utilidad_salida'
                
                                        if(!dir.exists(carpetaAlmacenamiento)){dir.create(carpetaAlmacenamiento)}
                                        
                                        nombreArchivo <- paste(carpetaAlmacenamiento, '/','utilidad_', etiqueta,'.csv')
                                        nombreArchivo <- str_replace_all(nombreArchivo, pattern=" ", repl="")
                                        write.csv(tablaUtiEstadoInnOp, file = nombreArchivo, row.names = FALSE)
                                        
                                        utiTable <- xtable(tablaUtiEstadoInnOp, digits = c(0, 1, 1, 1, 1, 1, 1, 1, 1, 0))
                                                          
                                        print(utiTable, type = 'html', include.rownames = FALSE)
                                        
                                }
                                
                                
                        }       
                        
                }
                 
        }
        
}

——————————————–

Cebada 2016 Otoño-Invierno Riego

Estado Año Ciclo Tipo producción Utilidad parcela de innovación (Pesos/ha) Número de observaciones parcela de innovación Utilidad en parcela del productor optimizando el sistema gracias a las capacitaciones de MasAgro (Pesos/ha) Número de observaciones en parcela del productor optimizando el sistema gracias a las capacitaciones de MasAgro Tasa de variación utilidad promedio parcela de innovación vs utilidad en parcela del productor optimizando el sistema gracias a las capacitaciones de MasAgro (%)
Querétaro 2016 Otoño-Invierno Riego 14335.5 10 13703.2 2 5

——————————————–

Cebada 2016 Primavera-Verano Temporal

Estado Año Ciclo Tipo producción Utilidad parcela de innovación (Pesos/ha) Número de observaciones parcela de innovación Utilidad en parcela del productor optimizando el sistema gracias a las capacitaciones de MasAgro (Pesos/ha) Número de observaciones en parcela del productor optimizando el sistema gracias a las capacitaciones de MasAgro Tasa de variación utilidad promedio parcela de innovación vs utilidad en parcela del productor optimizando el sistema gracias a las capacitaciones de MasAgro (%)
Hidalgo 2016 Primavera-Verano Temporal 6241.2 198 5899.5 13 6
México 2016 Primavera-Verano Temporal 3928.3 19 3059.5 1 28
Oaxaca 2016 Primavera-Verano Temporal 10691.3 3
Tlaxcala 2016 Primavera-Verano Temporal 7146.9 10 5232.4 2 37

——————————————–

Frijol 2016 Otoño-Invierno Riego

Estado Año Ciclo Tipo producción Utilidad parcela de innovación (Pesos/ha) Número de observaciones parcela de innovación Utilidad en parcela del productor optimizando el sistema gracias a las capacitaciones de MasAgro (Pesos/ha) Número de observaciones en parcela del productor optimizando el sistema gracias a las capacitaciones de MasAgro Tasa de variación utilidad promedio parcela de innovación vs utilidad en parcela del productor optimizando el sistema gracias a las capacitaciones de MasAgro (%)
Sinaloa 2016 Otoño-Invierno Riego 23786.4 7 26017.0 2 -9
Sonora 2016 Otoño-Invierno Riego 25620.8 6 19660.2 2 30

——————————————–

Frijol 2016 Primavera-Verano Riego

Estado Año Ciclo Tipo producción Utilidad parcela de innovación (Pesos/ha) Número de observaciones parcela de innovación Utilidad en parcela del productor optimizando el sistema gracias a las capacitaciones de MasAgro (Pesos/ha) Número de observaciones en parcela del productor optimizando el sistema gracias a las capacitaciones de MasAgro Tasa de variación utilidad promedio parcela de innovación vs utilidad en parcela del productor optimizando el sistema gracias a las capacitaciones de MasAgro (%)
Chiapas 2016 Primavera-Verano Riego 3817.5 1 1087.5 1 251
Oaxaca 2016 Primavera-Verano Riego 9270.0 1
Sonora 2016 Primavera-Verano Riego 16813.5 1

——————————————–

Frijol 2016 Primavera-Verano Temporal

Estado Año Ciclo Tipo producción Utilidad parcela de innovación (Pesos/ha) Número de observaciones parcela de innovación Utilidad en parcela del productor optimizando el sistema gracias a las capacitaciones de MasAgro (Pesos/ha) Número de observaciones en parcela del productor optimizando el sistema gracias a las capacitaciones de MasAgro Tasa de variación utilidad promedio parcela de innovación vs utilidad en parcela del productor optimizando el sistema gracias a las capacitaciones de MasAgro (%)
Chiapas 2016 Primavera-Verano Temporal 13086.0 1 10810.0 1 21
Durango 2016 Primavera-Verano Temporal 9868.0 5 4438.0 1 122
Hidalgo 2016 Primavera-Verano Temporal -1329.4 3 -11090.4 1 -88
Oaxaca 2016 Primavera-Verano Temporal 2638.6 38 917.7 8 188
Querétaro 2016 Primavera-Verano Temporal 6991.6 5 3567.2 2 96
Zacatecas 2016 Primavera-Verano Temporal 3691.0 9 3704.1 9 -0

——————————————–

Maiz 2016 Otoño-Invierno Riego

Estado Año Ciclo Tipo producción Utilidad parcela de innovación (Pesos/ha) Número de observaciones parcela de innovación Utilidad en parcela del productor optimizando el sistema gracias a las capacitaciones de MasAgro (Pesos/ha) Número de observaciones en parcela del productor optimizando el sistema gracias a las capacitaciones de MasAgro Tasa de variación utilidad promedio parcela de innovación vs utilidad en parcela del productor optimizando el sistema gracias a las capacitaciones de MasAgro (%)
Sinaloa 2016 Otoño-Invierno Riego 20227.2 77 22233.2 16 -9

——————————————–

Maiz 2016 Otoño-Invierno Temporal

Estado Año Ciclo Tipo producción Utilidad parcela de innovación (Pesos/ha) Número de observaciones parcela de innovación Utilidad en parcela del productor optimizando el sistema gracias a las capacitaciones de MasAgro (Pesos/ha) Número de observaciones en parcela del productor optimizando el sistema gracias a las capacitaciones de MasAgro Tasa de variación utilidad promedio parcela de innovación vs utilidad en parcela del productor optimizando el sistema gracias a las capacitaciones de MasAgro (%)
Campeche 2016 Otoño-Invierno Temporal 2178.6 4 369.7 4 489
Guerrero 2016 Otoño-Invierno Temporal 13272.1 1 7502.0 1 77
Hidalgo 2016 Otoño-Invierno Temporal -1022.7 3
México 2016 Otoño-Invierno Temporal -2210.0 1 -3860.0 1 -43
Oaxaca 2016 Otoño-Invierno Temporal 11179.2 2 7076.5 1 58
Puebla 2016 Otoño-Invierno Temporal -567.4 3 -5110.5 4 -89

——————————————–

Maiz 2016 Primavera-Verano Riego

Estado Año Ciclo Tipo producción Utilidad parcela de innovación (Pesos/ha) Número de observaciones parcela de innovación Utilidad en parcela del productor optimizando el sistema gracias a las capacitaciones de MasAgro (Pesos/ha) Número de observaciones en parcela del productor optimizando el sistema gracias a las capacitaciones de MasAgro Tasa de variación utilidad promedio parcela de innovación vs utilidad en parcela del productor optimizando el sistema gracias a las capacitaciones de MasAgro (%)
Chiapas 2016 Primavera-Verano Riego 8885.5 3 12068.8 2 -26
Chihuahua 2016 Primavera-Verano Riego 17108.0 6 7369.0 1 132
Guanajuato 2016 Primavera-Verano Riego 22834.1 1
Hidalgo 2016 Primavera-Verano Riego 26356.6 3 22794.6 3 16
Jalisco 2016 Primavera-Verano Riego 27078.0 2 24597.0 1 10
México 2016 Primavera-Verano Riego 18481.9 8 11416.4 5 62
Michoacán de Ocampo 2016 Primavera-Verano Riego 23216.1 31 22994.6 9 1
Oaxaca 2016 Primavera-Verano Riego 15447.9 21 10277.5 7 50
Puebla 2016 Primavera-Verano Riego 6490.4 6 4360.0 2 49
Querétaro 2016 Primavera-Verano Riego 18850.1 20 11577.1 7 63
Sinaloa 2016 Primavera-Verano Riego 28341.0 1
Sonora 2016 Primavera-Verano Riego 9038.0 1

——————————————–

Maiz 2016 Primavera-Verano Temporal

Estado Año Ciclo Tipo producción Utilidad parcela de innovación (Pesos/ha) Número de observaciones parcela de innovación Utilidad en parcela del productor optimizando el sistema gracias a las capacitaciones de MasAgro (Pesos/ha) Número de observaciones en parcela del productor optimizando el sistema gracias a las capacitaciones de MasAgro Tasa de variación utilidad promedio parcela de innovación vs utilidad en parcela del productor optimizando el sistema gracias a las capacitaciones de MasAgro (%)
Campeche 2016 Primavera-Verano Temporal 7528.7 49 4865.0 17 55
Chiapas 2016 Primavera-Verano Temporal 3748.0 201 7542.5 34 -50
Chihuahua 2016 Primavera-Verano Temporal 6090.0 2 6090.0 2 0
Durango 2016 Primavera-Verano Temporal 4625.9 45 2015.4 17 130
Guanajuato 2016 Primavera-Verano Temporal 12530.1 2
Guerrero 2016 Primavera-Verano Temporal 6287.0 1243 5546.4 127 13
Hidalgo 2016 Primavera-Verano Temporal 8850.6 56 7657.6 11 16
Jalisco 2016 Primavera-Verano Temporal 18432.5 47 15665.9 22 18
México 2016 Primavera-Verano Temporal 10363.4 21 14223.3 7 -27
Michoacán de Ocampo 2016 Primavera-Verano Temporal 13391.3 25 11663.9 7 15
Morelos 2016 Primavera-Verano Temporal 11871.6 386 9615.8 19 23
Oaxaca 2016 Primavera-Verano Temporal 4312.7 206 2994.0 37 44
Puebla 2016 Primavera-Verano Temporal 6166.9 20 344.9 5 1688
Querétaro 2016 Primavera-Verano Temporal 7521.8 47 4699.4 14 60
San Luis Potosí 2016 Primavera-Verano Temporal 2083.0 9 1936.0 3 8
Tabasco 2016 Primavera-Verano Temporal 6308.8 4 4724.3 3 34
Tlaxcala 2016 Primavera-Verano Temporal 11318.0 24 6666.6 8 70
Veracruz de Ignacio de la Llave 2016 Primavera-Verano Temporal 3583.8 30 1417.9 16 153
Yucatán 2016 Primavera-Verano Temporal 3830.1 13 2878.9 3 33
Zacatecas 2016 Primavera-Verano Temporal 3639.2 7 1873.3 8 94

——————————————–

Sorgo 2016 Otoño-Invierno Temporal

Estado Año Ciclo Tipo producción Utilidad parcela de innovación (Pesos/ha) Número de observaciones parcela de innovación Utilidad en parcela del productor optimizando el sistema gracias a las capacitaciones de MasAgro (Pesos/ha) Número de observaciones en parcela del productor optimizando el sistema gracias a las capacitaciones de MasAgro Tasa de variación utilidad promedio parcela de innovación vs utilidad en parcela del productor optimizando el sistema gracias a las capacitaciones de MasAgro (%)
Chiapas 2016 Otoño-Invierno Temporal -653.5 2 592.5 1 -210

——————————————–

Sorgo 2016 Primavera-Verano Riego

Estado Año Ciclo Tipo producción Utilidad parcela de innovación (Pesos/ha) Número de observaciones parcela de innovación Utilidad en parcela del productor optimizando el sistema gracias a las capacitaciones de MasAgro (Pesos/ha) Número de observaciones en parcela del productor optimizando el sistema gracias a las capacitaciones de MasAgro Tasa de variación utilidad promedio parcela de innovación vs utilidad en parcela del productor optimizando el sistema gracias a las capacitaciones de MasAgro (%)
Michoacán de Ocampo 2016 Primavera-Verano Riego 22034.5 1
Sinaloa 2016 Primavera-Verano Riego 17455.0 1 16360.0 1 7
Sonora 2016 Primavera-Verano Riego 2196.7 1 796.8 1 176
Zacatecas 2016 Primavera-Verano Riego 18838.2 2

——————————————–

Sorgo 2016 Primavera-Verano Temporal

Estado Año Ciclo Tipo producción Utilidad parcela de innovación (Pesos/ha) Número de observaciones parcela de innovación Utilidad en parcela del productor optimizando el sistema gracias a las capacitaciones de MasAgro (Pesos/ha) Número de observaciones en parcela del productor optimizando el sistema gracias a las capacitaciones de MasAgro Tasa de variación utilidad promedio parcela de innovación vs utilidad en parcela del productor optimizando el sistema gracias a las capacitaciones de MasAgro (%)
Chiapas 2016 Primavera-Verano Temporal 4582.0 2
Durango 2016 Primavera-Verano Temporal 5060.0 7 3085.0 1 64
Michoacán de Ocampo 2016 Primavera-Verano Temporal 12965.0 1
Puebla 2016 Primavera-Verano Temporal 17942.3 2 4375.9 3 310
San Luis Potosí 2016 Primavera-Verano Temporal 3057.9 4 -1345.0 1 -327

——————————————–

Trigo 2016 Otoño-Invierno Riego

Estado Año Ciclo Tipo producción Utilidad parcela de innovación (Pesos/ha) Número de observaciones parcela de innovación Utilidad en parcela del productor optimizando el sistema gracias a las capacitaciones de MasAgro (Pesos/ha) Número de observaciones en parcela del productor optimizando el sistema gracias a las capacitaciones de MasAgro Tasa de variación utilidad promedio parcela de innovación vs utilidad en parcela del productor optimizando el sistema gracias a las capacitaciones de MasAgro (%)
Baja California 2016 Otoño-Invierno Riego 11903.2 6 10329.8 6 15
Sinaloa 2016 Otoño-Invierno Riego 12663.9 5 8468.8 2 50
Sonora 2016 Otoño-Invierno Riego 11167.6 40 9258.1 20 21

——————————————–

Trigo 2016 Otoño-Invierno Temporal

Estado Año Ciclo Tipo producción Utilidad parcela de innovación (Pesos/ha) Número de observaciones parcela de innovación Utilidad en parcela del productor optimizando el sistema gracias a las capacitaciones de MasAgro (Pesos/ha) Número de observaciones en parcela del productor optimizando el sistema gracias a las capacitaciones de MasAgro Tasa de variación utilidad promedio parcela de innovación vs utilidad en parcela del productor optimizando el sistema gracias a las capacitaciones de MasAgro (%)
Michoacán de Ocampo 2016 Otoño-Invierno Temporal 6998.3 8 10452.5 2 -33

——————————————–

Trigo 2016 Primavera-Verano Temporal

Estado Año Ciclo Tipo producción Utilidad parcela de innovación (Pesos/ha) Número de observaciones parcela de innovación Utilidad en parcela del productor optimizando el sistema gracias a las capacitaciones de MasAgro (Pesos/ha) Número de observaciones en parcela del productor optimizando el sistema gracias a las capacitaciones de MasAgro Tasa de variación utilidad promedio parcela de innovación vs utilidad en parcela del productor optimizando el sistema gracias a las capacitaciones de MasAgro (%)
Hidalgo 2016 Primavera-Verano Temporal 2703.7 1
México 2016 Primavera-Verano Temporal 6256.4 2 6995.0 1 -11
Oaxaca 2016 Primavera-Verano Temporal 10236.2 13 -3715.0 1 -376
Tlaxcala 2016 Primavera-Verano Temporal 6699.4 2 9109.6 1 -26
Zacatecas 2016 Primavera-Verano Temporal 5190.0 1 3630.0 1 43