Análisis de Datos del Estudio de Semilla a través de PGIS V03 Participativo 2017-2022

Author

Omar Benites Alfaro

1 Descripción del estudio

Se realizó un estudio de encuestas de GIS participativo (PGIS) que se llevo a cabo en la región Huanuco en 3 distritos en los años distintos 2017 y 2022 . El estudio cuenta con 441 observaciones y con 45 variables. A continuación un resumen.

Además, para el presente estudio se tuvieron las siguientes consideraciones:

  • Se consideraron los datos relacionados a tubérculos y raices, cereales hortalizas y leguminosas. En un segundo caso tambien se considero dividir los datos en cultivos para huerto y chacra.

  • Se consideraron los datos relacionados a las familias que se repiten en ambos estudios para poder realizar una comparación más precisa y no alterar demasiado las muestra. Las demás que no se repetin podrían considerarse para otro estudio a nivel provincial.

  • Se consideró los datos solo de rehuso de semilla, es decir las familias que adquirieron semillas en sus parcelas y las rehusaron o mantuvieron para una siguiente evento.

    Debajo una tabla resumen con las diferentes variables consideradas en el estudio, conjuntamente con sus estadisticos, tanto variables cuantitativas (numéricas) y cualitativas (caractér)

Var.

Tipo

Media

SD

Categorias

zona

Character

Amarilis:78, Huánuco-Quisqui:129, Molino-Umari:234

comunidad

Character

auchi:6, auragshay:46, aurimachay:50, challhuayog:27, colpa alta:15, cucho:28, huallacayan:6, huancapallac:8, la linda:3, malconga:22, marabamba:7, marcajananda:13, matibamba:5, mitocucho:10, mitotambo:16, ninacocha:3, nueva florida:2, palmapampa:17, picahuay:39, pichao:2, pucuchinche:8, punchao:8, rodeo de margos:29, san jose de paucar:13, san pablo de lanjas:15, san pedro de cani:18, santo domingo de nauyan:4, shismay:9, yaca:12

cultivo

Character

aji:1, alfalfa:1, arracacha:10, arveja:9, arverja:8, avena:1, brocoli:1, camote:3, cebada:10, cebolla:2, frijol:91, haba:42, maiz:174, mashua:3, oca:6, olluco:3, papa:40, quinua:2, tara:1, tarwi:21, tomate:2, trigo:6, yacon:1, zapallo:3

lugar_siembra

Character

parcela:441

duracion

Character

Anual:421, No_aplica:1, Perenne:19

categoria_varietal

Character

-:17, Mejorada:97, Mejorada y nativa:1, Nativa:162

f_siembra

Character

01 de Febrero:1, 15 de Enero:1, 20 de Diciembre:1, 20 de Enero:1, 20 de Octubre:3, 28 de Junio:1, Abril:4, Agosto:2, Diciembre:21, Diciembre 2015:1, Diciembre 2016:2, Enero:15, Febrero:14, Hace 1 año:1, Julio:11, Junio:1, Marzo:10, Mayo:3, Noviembre:43, Noviembre 2016:12, Octubre:31, Octubre 2016:6, Setiembre:9, Setiembre 2016:1, todo el año:1, Todo el año:2

f_cosecha

Character

01 Junio:1, 05 de Marzo:1, 1° semana de Mayo:1, 10 de Junio:1, 15 Julio:1, 20 de Julio:1, 25 de Marzo:1, 6 Mayo:1, 8 Abril:1, Abril:12, Agosto:38, Diciembre:1, Diciembre 2016:1, Enero:4, Fines de Abril:1, Julio:10, Julio 2017:53, Junio:31, Junio - Julio:1, Junio 2017:1, Marzo:3, Mayo:10, Mayo - Junio:1, mensual:1, Octubre:7, Quincena de Junio:1, Setiembre:10, Todo el año:3

f_siembra_normal

Character

-:41, Abril:10, agosto:1, Agosto:2, Agosto - Diciembre:1, Agosto - Octubre:1, Diciembre:50, Diciembre - Enero:1, Diciembre - Marzo:1, Enero:11, Enero - Marzo:1, Febrero:10, Julio:9, Junio:5, Marzo:12, Mayo:9, Noviembre:108, Noviembre - Diciembre:1, Octubre:68, Octubre - Diciembre:5, Octubre - Noviembre:2, Setiembre:29, Todo el año:54

f_cosecha_normal

Character

-:42, Abril:9, Abril - Mayo:1, Agosto:79, Diciembre:4, Enero:1, Enero - Marzo:1, Febrero:2, Febrero - Marzo:1, Julio:152, Julio - Setiembre:2, Junio:36, Junio - Agosto:1, Junio - Julio:6, Marzo:3, mayo:1, Mayo:10, Mayo - Junio:2, Noviembre:3, Octubre:13, Septiembre:5, Setiembre:5, Todo el año:53

aplico_quimico

Character

no:329, si:112

fuente_semilla

Character

comprado:69, propio:372

transaccion

Character

consumo:301, consumo y venta:131, venta:9

destino_produccion

Character

semilla:441

cant_semilla_sembro

Numeric

29.2321759

263.563646

nombre_zona_sector

Character

Auragshay:12, Aurimachay:32, Av._Peru:2, Ayrampu:1, Baja Villca:1, Bella Vista:4, Campanilla_Ragra:6, Cani:10, Chacualoma:2, Challhuayog:22, Colpa alta:2, Colpa_alta:3, Colquigocha:1, Cucho:9, Cucho barrio alto:5, Cucho_Borrio_Alto:9, Curaupampa:11, Dominio Llanupunaga-Lugmapata:1, Dos Claveles:19, Dos_esquinas:9, Era Pata:2, Gallay paccha:3, Gichga pata:2, Gloriapata:4, Huaga_Molle_Pampa:10, Huagamulli:4, Huallacayan:6, Huancanio:3, Huarico:1, Huaylaspampa:1, Huillca Baja:8, Huillca_baja:4, Iscupata:5, La linda:2, Lanjas:8, Limon pampa:1, Lugmapata:5, Lurin:2, Malconga:12, Marabamba:1, Maray Pata:3, Maray_pata:2, Maraypampa:1, Marcajananda:13, Matibamba:5, Mision:4, Mitocucho:6, Mitotambo:6, Monte azul:1, Ninacocha:2, Niñacocha:1, Nueva_Florida:2, Obaspampa:1, obaspampa alta:1, other_sector:33, Pachamuña:2, Palmapampa:4, Pampa:2, Paucar:5, Paucar_Grande:4, Picahuay:13, Pichao:2, Pichipampa:3, Pichipampa Chico:1, Pucochinche:4, Puitush:5, Pustuypunta:3, Quitasol:4, Racco:1, Rashu:3, Rodeo de Margos:6, Rodeo de Margus:2, Rodeo_de_Margos:5, Runas loma:3, San Pablo de Lanjas:1, Santo Domingo de Nauyan:1, Santo_Domingo_de_Nauyan:3, Shismay:5, Shogosh:3, Shoguscucho:2, tambo alto:1, Tranca:11, umpa:1, Varios alta cucho:2, Vista Alegre:4, Yaca:8, Yawargocha:1

medio_transporte_parcela_casa

Character

a pie:434, auto:4, moto:3

zona_ecologica

Character

altura con mañay:39, altura sin mañay:90, bajera:80, media:232

calendario

Character

Campana_Chica_tuberculos:14, Campana_grande_tuberculos:7, Campaña chica (tuberculos):16, Campaña grande (tuberculos):19, No corresponde:13, No_aplica:18, Tardio (leguminosas y granos):97, Tardio_leguminosas_y_granos:104, Temprano (leguminosas y granos):53, Temprano_leguminosas_y_granos:98

sistema_de_siembra

Character

Barbecho (yunta):134, Barbecho_yunta:130, Labranza minima, Chacmeo:51, Labranza minima, Chacmeo y Barbecho (yunta):2, Labranza_minima_Chacmeo:100, No_aplica:7, Otro:2, Otros:8, Tractor:7

propiedad

Character

Al_partir:6, Alquilado:25, Familiar:410

parte_sector_comunal

Character

No:441

irrigacion

Character

Regado (bomba y/o asperción):61, Regado (canal, graavedad):145, Secano (sin irrigacion):235

pendiente

Character

Fuertemente inclinado:118, Ligeramente inclinado:288, Plano:33

mano_de_obra

Character

ayni minka:6, contratada:62, familiar:359, familiar y ayni o minka:7, familiar y contratada:3, familiar, concentrada y ayni o minka:3

rango_100

Character

1900-2000:5, 1900 - 2000:1, 2000-2100:14, 2000 - 2100:12, 2100 - 2200:3, 2200-2300:1, 2200 - 2300:2, 2400-2500:7, 2400 - 2500:6, 2500-2600:8, 2500 - 2600:14, 2600-2700:28, 2600 - 2700:38, 2700-2800:29, 2700 - 2800:41, 2800-2900:37, 2800 - 2900:39, 2900-3000:18, 2900 - 3000:34, 3000-3100:14, 3000 - 3100:26, 3100-3200:2, 3100 - 3200:7, 3200-3300:5, 3200 - 3300:4, 3300-3400:4, 3300 - 3400:6, 3400-3500:3, 3400 - 3500:3, 3500-3600:3, 3600 - 3700:4, 3700-3800:3, 3700 - 3800:3, 3800-3900:14, 3900-4000:3

rango_200

Character

1800-2000:5, 1800 - 2000:1, 2000-2200:14, 2000 - 2200:15, 2200-2400:1, 2200 - 2400:2, 2400-2600:15, 2400 - 2600:20, 2600-2800:57, 2600 - 2800:79, 2800-3000:55, 2800 - 3000:73, 3000-3200:16, 3000 - 3200:33, 3200-3400:9, 3200 - 3400:10, 3400-3600:6, 3400 - 3600:3, 3600-3800:3, 3600 - 3800:7, 3800-4000:17

anual

Numeric

2,019.7551020

2.489775

clasificacion

Character

cereal:193, hortalizas:9, leguminosas:173, tuberculos_y_raices:66

units_cant_semilla_sembro

Character

:1, . Arrobas:1, . kg:1, . Kg:19, . sacos:1, / kg:5, / Kg:1, / lata:2, / saco:1, Arrobas:1, costal:1, gr:6, kg:65, Kg:284, lata:3, latas:3, mallas:1, mantadas:1, paquetes:1, plantas:5, saco:3, Saco_o_costal:8, sacos:16, unidades:1, Unidades:10

sca_cant_semilla_sembro

Numeric

-0.1095618

0.300183

2 Análisis descriptivo

2.1 Número de observaciones por distrito y comunidad

  • Tabla del número de observaciones por distrito

Distrito

Año

N

% Reg.

% Reg.-Año

Amarilis

2,017

36

46.2

8.2

Amarilis

2,022

42

53.8

9.5

Huánuco-Quisqui

2,017

63

48.8

14.3

Huánuco-Quisqui

2,022

66

51.2

15.0

Molino-Umari

2,017

99

42.3

22.4

Molino-Umari

2,022

135

57.7

30.6
  • Distribución de las observaciones por distrito

2.2 Área de estudio

La área de estudio esta comprendida en la región Huánuco ubicado en el centro norte del Perú. Se consideraron 3 distritos Quisqui, Amarilis y Molino-Umari. Los datos georeferencias, latitud y longitud, fueron tomados en dos años distintos: 2017 y 2022.

El número de observaciones en el 2017 esta distribuido en Amarilis (36), Huánuco-Quisqui (63) y Molino-Umari (99).

Luego, en el 2022, las observaciones estan distribuidas en Amarilis (42), Huánuco-Quisqui (66) y Molino-Umari (135).

2.3 Comparación según fuente de semilla en el 2017-2022

  • Tabla de variación en el uso de semilla para cultivos entre 2017-2022

Fuente de semilla

Año

N

%

Año

N

%

Propio

2017

154

77.8

2022

218

89.7

Comprado

2017

44

22.2

2022

25

10.3

2.3.1 Gráfico de variación porcentual según fuente de semilla en el periodo 2017-2022.

2.4 Variacion porcentual en rehuso de semilla entre 2017 y 2022

Tabla de variación porcentual en rehuso de semilla entre 2017 y 2022

anual

destino_produccion

n

pct_1

2,017

semilla

198

44.9

2,022

semilla

243

55.1

2.4.1 Gráfico de variación en el uso de semilla según destino entre 2017-2022 según clasificación

2.5 Variacion porcentual en el rehuso de semilla según la región

zona

destino_produccion

n

pct_1

Amarilis

semilla

78

17.7

Huánuco-Quisqui

semilla

129

29.3

Molino-Umari

semilla

234

53.1

2.6 Variacion porcentual en el rehuso de semilla considerando la región y año

Zona

Año

N

%

Año

N

%

Molino-Umari

2,017

99

42.3

2,022

135

57.7

Huánuco-Quisqui

2,017

63

48.8

2,022

66

51.2

Amarilis

2,017

36

46.2

2,022

42

53.8

2.7 Variación porcentual según la fuente de semilla y la clasificación

2.7.1 Tabla de la variación según la fuente de semilla y la clasificación entre 2017-2022

Fuente de semilla

Clasificación

Año

N

%

Año

N

%

comprado

cereal

2017

16

8.1

2022

7

2.9

comprado

hortalizas

2017

6

3.0

2022

comprado

leguminosas

2017

10

5.1

2022

9

3.7

comprado

tuberculos_y_raices

2017

12

6.1

2022

9

3.7

propio

cereal

2017

66

33.3

2022

104

42.8

propio

hortalizas

2017

2

1.0

2022

1

0.4

propio

leguminosas

2017

63

31.8

2022

91

37.4

propio

tuberculos_y_raices

2017

23

11.6

2022

22

9.1

2.7.2 Gráfico de variación en el uso de semilla para cultivos entre 2017-2022 según clasificación

2.8 Distribución de fuente de semilla segun rehuso y altitud entre el 2017 y 2022

Año

Fuente de semilla

Media

D.Est.

2,017

comprado

2,761.114

470.7357

2,017

propio

2,898.935

455.5229

2,022

comprado

2,938.758

237.3024

2,022

propio

2,799.399

326.1671

Diagrama

Diagrama

2.9 Distribución según fuente de semilla dividido año y por zona

2.9.1 Tabla de fuentes de semilla agrupado por año

Zona

Fuente de semilla

Año

N

%

Año

N

%

Amarilis

propio

2017

32

88.9

2022

40

95.2

Amarilis

comprado

2017

4

11.1

2022

2

4.8

Huánuco-Quisqui

propio

2017

43

68.3

2022

51

77.3

Huánuco-Quisqui

comprado

2017

20

31.7

2022

15

22.7

Molino-Umari

propio

2017

79

79.8

2022

127

94.1

Molino-Umari

comprado

2017

20

20.2

2022

8

5.9

2.10 Distribución de semillas segun cultivo de chacra y huerto

id

Año

Fuente de semilla

Sub clasificacion

N

%

Año

N

%

concat

1

2,017

comprado

chacra

46

79.3

2,022

25

100.0

chacra:comprado

2

2,017

comprado

huerto

12

20.7

huerto:comprado

3

2,017

propio

chacra

154

83.2

2,022

219

85.5

chacra:propio

4

2,017

propio

huerto

31

16.8

2,022

37

14.5

huerto:propio

Gráfico con la comparación

3 Análisis estadístico inferencial

3.1 Análisis de proporciones considerando rehuso de semilla

  • Prueba Chi-cuadrado según fuente de semilla y año

    Pearson's Chi-squared test

data:  table_clasif
X-squared = 11.773, df = 1, p-value = 0.0006008

Para un \(p<0.05\), se afirma que existe diferencias significativas entre las proprociones de las categorias de fuente de semilla entre el 2017 y 2022.

  • Prueba de Chi-cuadrado por clasificacion de cultivo y fuente de semilla

    Pearson's Chi-squared test

data:  table_clasif
X-squared = 35.713, df = 3, p-value = 8.609e-08

Para un \(p<0.05\), se afirma que existe diferencias significativas entre las proprociones de las clasificaciones de semilla y las categorias de fuente de semilla.

  • Prueba de Chi-cuadrado por clasficacion de cultivo y año

    Pearson's Chi-squared test

data:  table_anual_clasif
X-squared = 9.7681, df = 3, p-value = 0.02064

Se obtiene tambien que existen diferencias en el uso semilla según su clasificación tomando en cuenta que fueron medidos en 2017 (antes de pandemia) y 2021 (despúes de pandemia) con \(p<0.05\)

  • Gráfico de visualización de cambios en la proporción segun fuente de semilla, clasificación y año

3.1.1 Modelo lineal generalizado (GLM) de Poisson considerando la interacción año, clasificación y fuente de semilla bajo un contexto de rehuso de semilla.

Se contruyó el modelo generalizado de Poisson con covariables fuente de semilla, año y clasificación y como variable de respuesta la frecuencia o conteo de los grupos. Esto con el fin de observar las interacciones y que tan significativas fueron.

Realizamos un anova para ver efectos de interaccion y efectos principales son significativos al 95% de confianza:

ID Termino Gl Devianza Gl. Resid Dev.Resid p-valor
1 NULL NA NA 15 542.23224 NA
2 clasificacion 3 259.195819 12 283.03642 0.0000000
3 anual 1 4.599839 11 278.43659 0.0319750
4 fuente_semilla 1 228.781958 10 49.65463 0.0000000
5 clasificacion:anual 3 10.445614 7 39.20901 0.0151344
6 clasificacion:fuente_semilla 3 27.788721 4 11.42029 0.0000040
7 anual:fuente_semilla 1 7.342892 3 4.07740 0.0067328
8 clasificacion:anual:fuente_semilla 3 4.077400 0 0.00000 0.2532268

  • Los efectos principales de clasificación, año y fuente de semilla es significativo con un \(p<0.05\)

  • El efecto de interacción del año y fuente de semilla es significativo con un \(p<0.05\).

  • El efecto de interacción l de la clasificación y fuente de semilla es significativo con un \(p<0.05\).

  • El efecto de interacción de clasificacion y año es tambien significativo con \(p<0.05\).

  • El efecto de interacción triple clasificacion, año y fuente de semilla no es significativo \(p>0.05\).

3.1.2 Modelo lineal mixto (LMM) de Poisson considerando la interacción año, clasificación y fuente de semilla bajo un contexto de rehuso de semilla.

Usando la tabla anterior podemos

Se planteó 2 modelos, uno con solamente efectos principales aditivos y otro con efectos de interacción

Realizamos un análisis de varianza (ANOVA) para comparar que modelo resultó significativamente diferente. Con lo que se obtuvo que el modelo con interacción es significativo (\(p<0.05\)) y con un menor \(AIC=136.7\) con respecto al modelo sin interacción

npar AIC BIC logLik deviance Chisq Df Pr(>Chisq)
mod_rehuso_clas_fuente2 6 141.2009 145.8365 -64.60046 129.2009 NA NA NA
mod_rehuso_clas_fuente1 9 136.7553 143.7086 -59.37765 118.7553 10.44561 3 0.0151344

Se obtuvo que el menor AIC=136.76 para el modelo que incluye la interacción, por lo tanto es el modelo que escogemos.

Gráficamos los efectos mixtos relacionados a fuente de semilla, el cual nos dio

                                       2.5 %      97.5 %
.sig01                            0.39199087  3.53707367
(Intercept)                       1.40389249  5.48828030
clasificacionhortalizas          -3.81031583 -2.45762352
clasificacionleguminosas         -0.31516715  0.09560212
clasificaciontuberculos_y_raices -1.35940002 -0.79946221
anual2022                         0.01761303  0.39317451

3.2 Modelos mixto considerando interacción entre fuente de semilla y año, con intercepto en la clasificación.

npar AIC BIC logLik deviance Chisq Df Pr(>Chisq)
mod_rehuso_anual_clas2 4 148.1957 151.2861 -70.09787 140.1957 NA NA NA
mod_rehuso_anual_clas1 5 138.4298 142.2927 -64.21488 128.4298 11.76598 1 0.0006032
contrast estimate SE df z.ratio p.value
comprado anual2017 - propio anual2017 -1.2527634 0.1705698 Inf -7.344578 0.0000000
comprado anual2017 - comprado anual2022 0.5653127 0.2499104 Inf 2.262061 0.1421618
comprado anual2017 - propio anual2022 -1.6003068 0.1649120 Inf -9.704003 0.0000000
propio anual2017 - comprado anual2022 1.8180761 0.2151554 Inf 8.450060 0.0000000
propio anual2017 - propio anual2022 -0.3475434 0.1050362 Inf -3.308797 0.0056219
comprado anual2022 - propio anual2022 -2.1656195 0.2106983 Inf -10.278298 0.0000000