Introducción
En Costa Rica existen siete zonas productoras de caña y trece ingenios azucareros. De acuerdo con el Plan Anual Operativo (PAO) de DIECA (Dirección de Investigación y Extensión de la Caña de Azúcar) en el año 2018, aproximadamente 22% de los ensayos vigentes se encontraban en la Región Brunca. Entre los temas de investigación más importantes que se estudian en esa región están: la evaluación de variedades, la interacción de nutrientes y la dosificación nutricional, los madurantes, los herbicidas y la densidad de siembra.
Las variables objetivo son frecuentemente el rendimiento de campo y el rendimiento agroindustrial. Para este tipo de investigaciones en la región se han utilizado principalmente tres tamaños de parcela: 5 surcos de 10 m con una separación de 1.5 m (75 m2), también se trabaja con 6 surcos de 9 metros con una separación de 1.5 m (81 m2) o bien, 5 surcos de 8 m con una separación de 1.5 m (60 m2), pero nunca se ha validado la eficiencia de esos tamaños de parcela.
Tampoco se ha investigado el tamaño de la parcela experimental con el que se minimiza el error experimental. Por esta razón, el objetivo de este trabajo es encontrar el tamaño óptimo de la parcela experimental con caña de azúcar en la región Brunca de Costa Rica. Sripathi et al. (2017) agrupan, dentro de lo que llaman factores de diseño experimental, al tamaño de la parcela, tamaño del bloque y número de repeticiones y afirman que en los ensayos de rendimiento el error experimental es sensible a esos factores, debido a que la investigación agronómica depende de los datos registrados en los ensayos de campo.
En la definición del tamaño de la parcela experimental influyen muchos factores y los más importantes son: el error experimental y los costos asociados. El tamaño de la parcela debe permitir captar satisfactoriamente toda la heterogeneidad del suelo, considerada la principal fuente de variabilidad de la parcela experimental. Cocco et al. (2009) compararon los resultados que se obtuvieron con distintos tamaños y formas de parcela experimental en el cultivo de fresa sembrada en suelo y bajo el sistema de hidroponía. El cultivo en el suelo presentó mayor variabilidad experimental que el cultivo hidropónico. El desarrollo de los métodos para la estimación del tamaño óptimo de parcela tiene una base importante con el trabajo pionero de Smith (1938), quien encontró una relación asintótica negativa entre la varianza y el tamaño de la parcela. Son muchas las propuestas hechas desde entonces.
Investigaciones actuales fueron en gran parte desarrollados en Brasil y aplicados en diversos cultivos, por ejemplo: crotalaria (Crotalaria juncea L.) (Facco et al., 2017), berenjena (Solanum melongena) (Krysczun et al., 2018), tomate cherry (Solanum lycopersicum L.) (Giacomini y Lúcio, 2018), tomate (Solanum lycopersicum), habichuelas (Phaseolus vulgaris) y calabacín (Cucurbita pepo) (Schwertner et al., 2015), girasol (Helianthus annus L.) (Santos et al., 2015), café (Coffea) (Mendes et al., 2016), taro (Colocasia esculenta) (Silva, 2014) y boniato (Ipomoea batatas) (Rodríguez et al., 2018).
Algunos trabajos recientes también se han realizado en la India, con la mostaza de la India (Brassica juncea L.) (Khan y Tanwar, 2017). En Estados Unidos de América, con sandía (Citrullus lanatus) (Boyhan, 2013). En Costa Rica este tipo de investigación se han realizado para el cultivo de café, a inicios de los sesenta, (Paez, 1962) y para frijol en los años setenta (Mamani, 1971).
Los ensayos más recientes se realizaron en la zona de Bagaces, Guanacaste, para los cultivos de arroz (Oryza sativa) (Vargas y Navarro, 2014) y maíz (Zea mays) (Vargas y Navarro, 2017). Para el cultivo de la caña de azúcar en otros países ya se ha trabajado en el cálculo del tamaño óptimo de la parcela experimental. Por ejemplo: Igue et al. (1991) en Brasil; Bose y Khanna (1939) en la India y Palencia (1965) y Álvarez (1982) en Guatemala.
Materiales y métodos
Modelos de regresión segmentada
Dentro de las propuestas metodológicas para la estimación del tamaño óptimo de la parcela experimental, se encuentran los modelos de regresión segmentada. Paranaiba et al. (2009) proponen ajustar un modelo de regresión lineal con constante (LRP) para la relación entre el coeficiente de variación (CV) del rendimiento y el tamaño de parcela (medido en unidades experimentales básicas (UEB)), tal como se aprecia en la Figura 1, donde para los valores xi ≤ x0 esta relación está descrita por un modelo lineal decreciente hasta cierto punto, después se vuelve constante. El valor CVP que corresponde al punto de inflexión, se obtiene de un proceso iterativo.
Este modelo de regresión lineal segmentada se puede expresar como:
donde: CVx es el coeficiente de variación entre los totales para parcelas con x i unidades básicas; CVP es el coeficiente de variación en el punto donde se unen los dos segmentos; x es el tamaño de parcela medido en unidades básicas; y εi es el error asociado a CVx, supuestamente normal e independientemente distribuido, con media cero y varianza constante.
Dado que los dos segmentos de (1) son iguales en el punto x0 , entonces:
Por lo que el tamaño óptimo de la parcela experimental está dado por:
Otra alternativa dentro de las posibilidades de los modelos segmentados es el método de regresión cuadrática con constante (QRP), descrito por Ferreira (2007) y aplicado por Mendes et al. (2016) para la estimación del tamaño de la parcela experimental. Con respecto al modelo anterior, este método asume una forma polinomial de segundo grado, en lugar de una forma lineal en el primer segmento (Figura 2).
Para valores xi ≤ x0 el modelo es cuadrático y para valores xi > x0 es constante. De igual manera, la intercepción entre los segmentos determina el tamaño óptimo:
4). El punto x0 representa el punto de unión de los dos segmentos y se debe estimar juntamente con los demás parámetros del modelo. Dado que la curva debe ser continua y suavizada, las primeras derivadas con respecto a x en ambos segmentos deben ser la misma para el punto x0 . De acuerdo con Ferreira (2007), esta condición implica que:
Una vez igualada a cero, la ecuación (5) se resuelve para x y después de sustituir x por x0 , se obtiene:
El valor de la constante (CVP) se puede obtener al sustituir este valor en la ecuación (4):
Fuente de información y ensayo de uniformidad
Se utilizaron los datos de un ensayo de uniformidad de caña de azúcar realizado entre 2018 y 2019 en la finca El Porvenir en Pérez Zeledón, Costa Rica, localizada a una altura de 591 msnm. En el año 2017 se registró una precipitación anual de 3 673.8 mm y la temperatura promedio fue de 23.3 °C, con una máxima de 34.5 °C y una mínima de 15.4 °C.
Se plantaron de forma manual 40 surcos de 84 m de largo con una separación entre surco de 1.5 m de la variedad RB 99-381 a una densidad de tres chorros. En total fueron 4 800 m2 de parcela útil, la cual se dividió en 1 600 unidades experimentales básicas (UEB), cada una de 2 m de largo por 1.5 metros de ancho (3 m2).
Las prácticas del cultivo que se realizaron fueron las mismas que se aplican a las plantaciones comerciales de la región. Un mes antes de la siembra se preparó el terreno. La siembra se realizó manualmente colocando tres cañas con un traslape de 15 cm y se utilizan 13 t ha-1 de semilla. Para el control de malezas se empleó una mezcla química de Pendimetalina 50 EC (3 L ha-1) con Terbutilazina 50 SC (2 L ha-1) en preemergencia de la maleza (suelo desnudo). Y en post emergencia temprana, una mezcla química de Hexazinona 75 WG (0.5 kg ha-1) con Diurón 80 WG (2 kg ha-1).
Por ser caña de primer corte se aplicó 140 kg de nitrógeno,140 kg de fósforo,167 kg de potasio, 35 kg de magnesio y 40 kg de azufre por hectárea, distribuidos en tres fertilizaciones. Para la cosecha se consideran características de topografía, tamaño y distribución de las fincas de caña en la región, según lo cual se realiza manual o semi mecanizada. La cosecha se realizó de forma manual los días 6 y 7 de marzo de 2019 a los 10 meses de edad de la plantación. La variable que se midió fue el rendimiento de campo, tomando el peso de cada UEB medido en kilogramos.
Las unidades experimentales básicas se agruparon en unidades secundarias, de distintas formas y tamaños tomando como requisito que esas agrupaciones de parcelas adyacentes emplearan siempre la totalidad del área experimental, tal como lo describen Paranaiba et al. (2009), a cada una de estas unidades secundarias se les calculó la media de la producción, la varianza y el coeficiente de variación. Para la estimación del tamaño óptimo de parcela se emplearon los métodos de regresión lineal con constante (LRP) y regresión cuadrática con constante (QRP). En la estimación de los parámetros y el análisis de la información se empleó el paquete estadístico SAS versión 9.3 y Python versión 3 (2008).
Resultados y discusión
Al momento de la cosecha y pesada del ensayo, de las 1 600 unidades experimentales básicas se obtuvo una producción promedio de 19.9 kg por unidad experimental básica, con una desviación estándar de 3.8 kg. Los valores oscilaron entre 9.5 kg y 38 kg. Tal como se observa en la Figura 3, 41% de las unidades experimentales básicas pesaron entre 18 y 22 kg.
En total se agruparon los datos del ensayo en 20 tamaños que corresponden a 63 formas distintas (Cuadro 1).
Cantidad de parcelas secundarias | Tamaño | Formas | |
UEB | m2 | ||
1 600 | 1 | 3 | (1×1) |
800 | 2 | 6 | (1×2),(2×1) |
400 | 4 | 12 | (1×4),(4×1),(2×2) |
320 | 5 | 15 | (1×5),(5×1) |
200 | 8 | 24 | (1×8),(8×1),(2×4),(4×2) |
160 | 10 | 30 | (1×10),(10×1),(2×5),(5×2) |
100 | 16 | 48 | (4×4),(2×8),(8×2) |
80 | 20 | 60 | (1×20),(20×1),(2×10),(10×2),(4×5),(5×4) |
64 | 25 | 75 | (5×5) |
50 | 32 | 96 | (4×8),(8×4) |
40 | 40 | 120 | (1×40),(40×1),(2×20),(20×2),(4×10),(10×4),(5×8),(8×5) |
32 | 50 | 150 | (5×10),(10×5) |
25 | 64 | 192 | (8×8) |
20 | 80 | 240 | (2×40),(40×2),(4×20),(20×4),(8×10),(10×8) |
16 | 100 | 300 | (5×20),(20×5),(10×10) |
10 | 160 | 480 | (4×40),(40×4),(8×20),(20×8) |
8 | 200 | 600 | (5×40),(40×5),(10×20),(20×10) |
5 | 320 | 960 | (8×40),(40×8) |
4 | 400 | 1 200 | (10×40),(40×10),(20×20) |
2 | 800 | 2 400 | (20×40),(40×20) |
La relación entre el coeficiente de variación de cada de estas formas y el tamaño de la parcela experimental se aprecia en la Figura 4 (para una mejor apreciación se cortó el eje de las abscisas en 400), como se esperaba el coeficiente de variación decrece rápidamente en el segmento de parcelas pequeñas y después, conforme el tamaño aumenta, el coeficiente de variación tiende a disminuir menos que proporcionalmente.
Al aplicar el modelo de regresión lineal segmentado (LRP), se obtuvieron los parámetros del modelo:
El coeficiente de determinación fue de 61.62% y el estadístico F del análisis de varianza para la prueba de significancia de los parámetros estimados también fue significativo (p< 0.01) (Cuadro 2).
Parámetros | LRP(1) | QRP(1) | ||
β 0 | 12.6093 | ** | 14.1172 | ** |
(0.6728) | (0.8055) | |||
β 1 | -0.4338 | ** | -0.7654 | ** |
(0.0537) | (0.1149) | |||
β 2 | 0.0123 | ** | ||
(0.0031) | ||||
n | 63 | 63 | ||
F | 114.3413 | 136.2508 | ||
R2 ajustado | 61.62 | 76.96 | ||
W (2) | 0.9276 | 0.9309 | ||
P | 2.17 | 2.2256 | ||
x0 | 24.0531 | 31.0742 |
*1)= error estándar entre paréntesis; **= significativo al 1%; *= significativo al 5%; (2)= estadísticas de prueba Shapiro-Wilk.
El tamaño óptimo estimado es aproximadamente de 25.05 UEB. Considerado que cada UEB es de 3 m2, el tamaño óptimo equivale a 72.16 m2(Figura 5).
Con base a los datos para ajustar el modelo de regresión cuadrática con constante, se obtuvieron los estimadores:
En este caso, el tamaño óptimo fue de 31.07 UEB, equivalente a parcelas de 93.22 m2 (Figura 6).
Este modelo tiene un coeficiente de determinación de 76.96% y un estadístico F también significativo (p< 0.01) (Cuadro 2). El tamaño de la parcela tiende a ser más alto con el método QRP que con el LRP, debido a la diferencia en el ajuste. Peixoto et al. (2011) en ensayos con maracuyá le atribuyen esta tenencia a la curvatura del modelo. Mendes et al. (2016) encontraron un tamaño de parcela para café de 10.53 UEB más grande con el método QRP respecto al método LRP.
Por otro lado, Silva et al. (2012) utilizando un ensayo de rábano, encontraron una diferencia de 2.57 UEB a favor de método QRP. En ambos casos se aplicó el contraste de normalidad de Shapiro-Wilk (Shapiro y Wilk, 1965) a los errores del modelo (εi), tanto para el modelo LRP como para el QRP, y no se rechaza la hipótesis nula de normalidad con (p= 0.1998) y (p= 0.2181), respectivamente. Mediante los métodos de regresión segmentada se encontraron tamaños de parcela entre 72 y 93 m2.
Álvarez (1982) recomendó que el tamaño de las parcelas de investigación en caña de azúcar oscilara entre 80 y 115 m2 (5 o 6 surcos de ancho de 8 o 9.6 m), con un largo de 10 a 12 m y parcelas de 4 surcos de 12 m de largo (76.8 m2) en los casos en que existiera alguna limitación técnica o práctica (insumos, semilla, terreno, etc.). Ese autor aplicó el análisis de regresión múltiple y el método de curvatura máxima bivariada, sin encontrar grandes diferencias entre los dos métodos y en ambos casos, con coeficientes de determinación superiores 80%.
Aunque los resultados aquí obtenidos se asemejan a los encontrados por Álvarez (1982), la comparación debe tomarse con reservas, debido a que se trata de una investigación que se llevó a cabo en Guatemala y como el autor lo indica en el título de su trabajo, su resultado está sujeto a las condiciones de la finca Bulbuxya, que se ubica en el departamento de Suchitepéquez de la región sur-occidente de la República de Guatemala y la inferencia, más allá de dicha zona, no es válida. Otro aspecto importante es la vigencia de este tipo de investigaciones.
Además, son necesarias las actualizaciones periódicas. Ese trabajo es de inicios de los años ochenta y no hay registros de actualización alguna. En Guatemala, Palencia (1965), empleó el método de Smith (1938) y calculó un tamaño óptimo de 27.24 m2, con lo que redondeó y recomendó una parcela útil de 28 m2 (dos surcos de 8 m de largo con separación de 1.8 a 2 m), mucho más pequeñas que en el caso anterior. La parcela total sería de 4 surcos de 10 m de largo para excluir los dos surcos laterales y un metro en cada extremo de los surcos.
Conclusiones
El tamaño óptimo de la parcela experimental calculado con el método de regresión lineal con constante fue de 72 m2 y con el método de regresión cuadrática con constante de 93 m2. Si se considera la variación que tienen los demás factores que no fueron incluidos en esta estimación, entonces, con el fin de minimizar el error experimental, se recomienda para ensayos de rendimiento que el tamaño de la parcela experimental con caña de azúcar debe estar entre 72 y 93 m2 y ajustar la forma de la parcela según la disponibilidad de espacio.
Se recomienda replicar la metodología aquí desarrollada en las demás regiones cañeras.