Ajuste de respuestas a niveles de nutrimentos
y análisis de económico

Proyecto Fertilizantes, IBTA-FAO, Cochabamba, Casilla 3299.

 

Introducción

La respuesta de plantas a niveles crecientes de un nutrimento tal como N, P ó K puede describirse con curvas, véase la Gráfica 1).

La finalidad de ajustar las respuestas a niveles de nutrimentos es ante todo la de facilitar un análisis económico y para esta finalidad se destaca la ecuación cuadrática, por la facilidad de hacer dicho análisis.

Planteamiento de la ecuación cuadrática

La ecuación cuadrática es la siguiente:

(1) Y= A + BX + CX2, donde

Y= El rendimiento estimado

A= El rendimiento estimado sin aplicación del nutrimento (también llamado coeficiente de regresión).

X= El nutrimento en estudio

B= Coeficiente de regresión

C= Coeficiente de regresión.

Los coeficientes de regresión A, B Y C pueden resolverse en base a las siguientes tres ecuaciones normales:

(2) An + BΣX + CΣX2 = ΣY

(3) AΣX + BΣX2 + CΣX3 = ΣXY

(4) AΣX2 + BΣX3 + CΣX4 = ΣX2Y

En la ecuación (2), n es igual al número de los diferentes tratamientos. ΣX significa la suma de cada uno de los tratamientos y ΣX2, la suma del cuadrado de cada uno de ellos, etc.

Para poder eliminar el término con el coeficiente A, se múltiplica la ecuación (2) con ΣX y la (3) con n, a saber:

(5) AnΣX + BΣX•ΣX = CΣX2•ΣX = ΣY • ΣX

(6) AnΣX + BnΣX2 + CnΣX3 = nΣXY

Al substraer la ecuación (6) de la (5) se elimina los términos AnΣX y resulta la siguiente ecuación:

(7) B (ΣX * ΣX - nΣX2) + C(ΣX2•ΣX - nΣX3) = ΣYΣX –nΣXY

Multiplicando la ecuación (2) con EX2 y la (4) con n se obtiene:

(8) AnΣX2 + BΣX2 * ΣX + CΣX2 - ΣX2 = ΣYΣX2

(9) AnΣX2 + BnΣX3 + CnΣX4 =nΣX2Y

Se eliminan los términos AnΣX2 substrayendo la ecuación (9) de la (8) y se obtiene:

(10) B(ΣX2. ΣX - nΣX3) + C(ΣX2. ΣX2 - nΣX4) = ΣY•ΣX2 - 2Y

(11) D = ΣX * ΣX - nΣX2 (Ver la ecuación 7)

(12) E = ΣX2 * ΣX - nΣX3 (Verla ecuación 7 y 10)

(13) F = ΣX2 * ΣX2 - nΣX4 (Ver la ecuación 10)

(14) G = ΣY - ΣX - nΣXY (Ver la ecuación 7)

(15) H = ΣY * ΣX2 - nΣX2Y (Ver la ecuación 10)

Ahora las ecuaciones (7) y (10) pueden formularse como sigue:

(16) B • D + C • F = G

(17) B • E + C • F = H

Para poder determinar el coeficiente C se multiplica la ecuación (16) con E y la (17) con D y substrae la (17) de la (16) y se obtiene:

(18) C (E•E - D•F) = G- E - D•H

(19) C = (G * E – D * H) / (E * E – D * F)

Ya se puede convertir la ecuación (16) a la siguiente:

(20) B = (G – C * E) / D

La ecuación (2) puede convertirse a la siguiente:

(21) A = (ΣY – BΣX – CX2) / n

Ejemplo de resolución de la ecuación cuadrática

Datos:

0 1 2 3 (X en códigos)

0

50

100

150

(X en kg N/ha)

2123

2546

2687

2699

(Y en kg/ha)

Para evitar productos y sumas elevadas se hacen primero los cálculos con X en códigos:

X

X2

X3

X4

Y

XY

X2Y

X0

0

0

0

0

2123

0

0

X1

1

1

1

1

2546

2546

2546

X2

2

4

8

16

2687

5374

10748

X3

3

9

27

81

2699

8097

24291

Σ

6

14

36

98

10055

16017

37585

Ahora se pueden calcular las variables D a H [ver las ecuaciones (11) a (15)]:

D = 6x6 - 4x14 = -20

E = 14x6 - 4x36 = -60

F = 14x14 - 4x98 = -196

G = 10055 x 6 - 4 x 16017 = - 3738

H =10055 x 14 - 4x 37585 = - 9570

C =[ (-3738) (-60) - (-20) (-9570) ] [(-60) (-60) - (-20) (-196)

C = 32880 / - 320 = - 102.75

B- [ -3738 - (-102,75) (-60) / - 20 = 495,15

A = 10055 - 495.15x 6 - (-102,75) x 14 / n = 2130.65

Y = 2130.65 + 495.15 X - 102,75 X2 (X en códigos)

Para tener X en kg/ha se divide 495,15 con 50 y - 102.75 con (50X50), o sea:

Y= 2130.65 + 9.903X - 0.0411X2 (X en kgN/ha)

En forma general dicha división se hace con la diferencia entre cada tratamiento y con dicha diferencia al cuadrado, respectivamente. Se observará que resulta más fácil hacer los cálculos manuales si se usa tratamientos equidistantes.

Análisis económico

En la Gráfica 1 se observa la línea q/p - X, donde X es el nivel de nitrógeno, q el costo del nitrógeno en Bs/kg N y p el valor del grano en Bs/kg.

Para encontrar el nivel de X que dé la máxima utilidad se buscará una línea que sea paralela a la q/p X y tangente a la curva BX + CX2, para lo cual se iguala el derivado de la línea a el de la curva:

(22) B + 2CX = q/p

(23) X = (q/p – B) / 2C

(24) X = (B – q/p) / (-2C)

Un agricultor no tendrá necesariamente liquidez para invertir en el insumo para lograr la máxima utilidad y resulta útil el concepto tasa de retorno marginal, R, la que indica la utilidad debida al último boliviano invertido en el insumo.

La línea de costo del insumo se multiplicará con R + 1, véase la Gráfica 1, donde sé indica tanto q/pX como (R + 1) q/p con R=1.

Si en la ecuación a (24) se substituye q/p por (R + 1) q/p se obtiene:

(25) X = (B – (R + 1) q/p) / (-2C)

Habiendo determinado mediante la ejecución de ensayos la respuesta física a niveles de un nutrimento, y mediante cómputos o los cálculos indicados líneas arriba, los coeficientes B y C correspondientes resulta muy fácil para el usuario (el extensionista) actualizar los cálculos económicos.

Primero se definen la tasa de retorno marginal y el costo de fertilizantes. Además se estima el valor de la cosecha. Es conveniente usar precios a nivel de la chacra.

El costo de financiamiento se tratará en una de las siguientes formas

a) Llevar el valor de la cosecha al momento de la siembra, dividiendo el precio, P con (1 + r), donde r es el interés a pagar durante el período vegetativo. Ejemplos: 1) p = 0,40 Bs. el interés anual es de 20 % y el período vegetativo de 5 meses; por consiguiente resulta que:

r = 20/100 X 5/12 = 0,083 y 1 + r = 1,083.

b) Alternativamente se multiplicará el costo de fertilizantes con (1 + r)

Ejemplos de recomendación de niveles de nitrógeno, X:

Supuestos de respuesta física: B = 9,903 y C = -0,0411

a) Supuestos específicos: R = 0, q = 1.60 y p = 0.40

X = (9.903 - 1 x 1.60/0,4) / (-2) (-0.0411) = 71.8

b) R = 0,5,

q =1,60

y p = 0,40

X= (9,903 - 1,5 X 1,60/0,4) / (-2) (-0,0411) = 47,5

c) R = 0.5,

q = 1.60

y p = 0.8

X= (9,903 - 1.5 x 1.60/0,8) / [(-2) (-0,0411)] = 84

Estos tres ejemplos demuestran a la vez la importante influencia que tienen normalmente la tasa de retorno marginal y la relación del costo del insumo y el valor de la cosecha, en las recomendaciones sobre el uso del insumo. Habiendo variaciones entre localidades y entre años tanto en el costo del insumo como en el valor de la cosecha resulta importante que los extensionistas sepan realizar los cálculos según la ecuación (25) y que los investigadores, en lo posible, les entreguen la información sobre la respuesta física en una forma que permita una fácil actualización de los cálculos económicos.

Si, por interés académico se desea conocer el nivel X que dé la máxima producción se pone

R= -1 y el término (R = 1) q/p resulta cero y por consiguiente

X = B / -2C = 120

Formación de promedios

Al tener datos de repeticiones de un ensayo o datos de varios ensayos surge la cuestión si deben usarse estos promedios. El Programa de Fertilizantes ha optado por basar los análisis de regresión en los promedios. Por consiguiente, los cálculos de regresión resultan más simples.

Supongamos que en una zona ecológica determinada se hicieron en diferentes años ensayos con niveles de nitrógeno en el cultivo de trigo, a saber:

PERIODO

Nº DE ENSAYOS

COEFICIENTE DE REGRESIÓN

A

B

C

1989/90

1

2131

9.9030

-0.0411

1990/91

2

1246

11.3401

-0.05530

1991/92

3

1750

8.0556

0.03925

1992/93

4

1675

7.4580

0.02550

Suma

4

6802

36.7567

-0.16115

Promedio

1701

9.1892

-0.04029

La ecuación promedio servirá de base para formular una recomendación general para la zona ecológica en cuestión y se incluye en el promedio la variación debida a diferencias entre años, lo cual se puede justificar, aceptando que no se puede predecir al momento de fertilizar como resultará el clima.

Esquema para resolver una regresión

(Modelo para Técnica de Campo)

Tres niveles equidistantes de una sola variable independiente.

1

2

3

4

5

Y

X

XY

X2

X2Y

0

1

2

0

1

4

ΣY =

ΣXY =

ΣX2Y =

Llene la Columna 1 con los datos de rendimientos correspondientes a los tratamientos respectivos y forme la suma (ΣY). Luego multiplique cada uno de los datos de rendimiento con el X correspondiente de la columna 2, ponga cada producto de la Columna 3 y forme la suma (ΣXY). En la Columna 5, ponga el producto de la Columna 1 por la 4 y forme la suma (ΣX2Y).

Llene el esquema siguiente con los valores indicados:

G = (ΣY) --------------- x3 - (ΣXY) --------------- x3 = ------------

H = (ΣY) --------------- x5 - (ΣX2Y) -------------- x3 = ------------

C = [ (G) ---------- x (-12) - (H) ---------- x (-6) ] / (-12) = = = = =

B = [ (G) ---------- - (C) ------------------- x (-12) ] / (-6) = = = = =

A = [ (ΣY) -------- - (B) -----------x 3 - (C) ------------x 5] / 3 = = = =

Luego convierta B dividiendo con la diferencia entre los niveles, y C, con el cuadrado de la diferencia. (B) ---------- / ---------- = ; (C) ---- / (----) (----) = = = = = = =

Esquema para resolver una regresión

(Modelo para Técnico de campo)

- Cuatro niveles equidistantes de una sola variable independiente

1

2

3

4

5

Y

X

XY

X2

X2Y

0

1

2

3

0

1

4

9

ΣY =

ΣXY =

ΣX2Y =

Llene la Columna 1 con los datos de rendimientos correspondiente a los tratamientos respectivos y forme la suma (ΣY). Luego multiplique cada uno de los datos de rendimiento con el X correspondiente de la Columna 2, ponga cada producto en la Columna 3 y forme la suma (ΣXY). En la columna 5, ponga el producto de la Columna 1 por la 4 y forme la suma (ΣX2Y).

Llene el esquema siguiente con los valores indicados:

G = (ΣY) -------------- x6 - (ΣXY) ---------- x4 =

H = (ΣY) -------------- x14 - (ΣX2Y) ------ x4 =

C = [ (G) ----------(60) - (H) ---------- x (-20) ] / -320 = = = = = =

B = [ (G) --------- - (C) --------------------- x (-60) ] / (-20) = = = = = =

A = [ (ΣY) ------- - (B) -------- x6 - (C) ------------------ x 14 ] / 4 = = = = =

Luego convierta B dividiendo con la diferencia entre los niveles, y C, con el cuadrado de la diferencia. (B) ---------- / ---------- =; (C) ------ / (-----) (-----) = = = = =

Esquema para resolver una regresión

Modelo para técnicos de campo.

Cinco niveles equidistantes de una segunda variable independiente:

1

2

3

4

5

Y

X

XY

X2

X2Y

0

1

2

3

4

0

1

4

9

16

ΣY =

ΣXY =

ΣX2Y =

Llene la columna 1 con los datos de rendimientos correspondientes a los tratamientos respectivos y forme la suma (ΣY). Luego multiplique cada uno de los datos de rendimiento con el X correspondiente de la Columna 2, ponga cada producto en la Columna 3 y forme la suma (ΣXY). En la Columna 5, ponga el producto de la Columna 1 por la 4 y forme la suma (ΣX2Y).

Llene el esquema siguiente con los valores indicados:

G = (ΣY) x 10 - (ΣXY) ---------- x5 = ----------

H = (ΣY) x 30 - (ΣX2Y) ------- x5 = -------

C = [ (G) ---------- (- 200) - (H) ---------- x (- 50) ] / (- 3500) = ======

B =-[ (G) ----------- - (-C) ----------------- x (- 200) ] / (- 50) = =====

A = [ (ΣY) ----------- - (B) -----------x 10 - (C) ------- x 30] / 5 = =======

Luego convierta B dividiendo con la diferencia entre los niveles, Y C, con el cuadrado de la diferencia.

(B) ---------- / ---------- =; (C) ---------- / (---)(---) = =======

Referencias

El programa de Fertilizantes de la FAO dispone de un paquete de programas de análisis de resultados de ensayos de fertilización para usar en computadoras PC, IBM o compatibles con el sistema operativo MS - DOS y con un disco duro.

Para más informaciones sobre este paquete de programas, los mismas que son interactivos, consultar a:

COPE, F., 1985.Statistical design and analysis programme suite for fertilizar field trials. FAO Fertilizer Programe, Roma. s.d.t.

COLWELL, S.A., 1975 Computations for studies of soil fertility and fertilizer requirements. FAO, AGL/M ISC/75/1.

DRAFT EDITION. Roma s.d.t.

Gráfico 1. Respuesta a niveles de nitrógeno y la ecuación cuadrática.