Articulo de referencia

Función logística generalizada

A=M=0, K=C=1, B=3, ν=0,5, Q=0,5 Efecto de la variación del parámetro A. Todos los demás parámetros son 1. Efecto de la variación del parámetro B. A = 0, todos los demás parámetr...

A=M=0, K=C=1, B=3, ν=0,5, Q=0,5
Efecto de la variación del parámetro A. Todos los demás parámetros son 1.
Efecto de la variación del parámetro B. A = 0, todos los demás parámetros son 1.
Efecto de la variación del parámetro C. A = 0, todos los demás parámetros son 1.
Efecto de la variación del parámetro K. A = 0, todos los demás parámetros son 1.
Efecto de la variación del parámetro Q. A = 0, todos los demás parámetros son 1.
Efecto de la variación del parámetroν{\displaystyle \nu }. A = 0, todos los demás parámetros son 1.

La función o curva logística generalizada es una extensión de las funciones logística o sigmoide . Desarrollada originalmente para modelar el crecimiento, permite obtener curvas en forma de S más flexibles. A veces se la denomina curva de Richards, en honor a F. J. Richards , quien propuso la forma general para esta familia de modelos en 1959.  

Definición

La curva de Richards tiene la siguiente forma:

Y(t)=A+KA(do+QmiBt)1/ν{\displaystyle Y(t)=A+{K-A \over (C+Qe^{-Bt})^{1/\nu }}}

dóndeY{\displaystyle Y}= peso, altura, tamaño, etc., yt{\displaystyle t}= tiempo. Tiene seis parámetros:

  • A{\displaystyle A}: la asíntota horizontal izquierda;
  • K{\displaystyle K}: la asíntota horizontal derecha cuandodo=1{\displaystyle C=1}. SiA=0{\displaystyle A=0}ydo=1{\displaystyle C=1}entoncesK{\displaystyle K}se denomina capacidad de carga ;
  • B{\displaystyle B}: la tasa de crecimiento;
  • ν>0{\displaystyle \nu >0} : afecta cerca de qué punto se produce el crecimiento máximo asintótico.
  • Q{\displaystyle Q}: está relacionado con el valorY(0){\displaystyle Y(0)}
  • do{\displaystyle C}: normalmente toma un valor de 1. De lo contrario, la asíntota superior esA+KAdo1/ν{\displaystyle A+{K-A \over C^{\,1/\nu }}}

La ecuación también se puede escribir:

Y(t)=A+KA(do+miB(tMETRO))1/ν{\displaystyle Y(t)=A+{K-A \over (C+e^{-B(t-M)})^{1/\nu }}}

dóndeMETRO{\displaystyle M}puede considerarse como un momento de inicio, en el queY(METRO)=A+KA(do+1)1/ν{\displaystyle Y(M)=A+{K-A \over (C+1)^{1/\nu }}}. Incluyendo ambosQ{\displaystyle Q}yMETRO{\displaystyle M}puede ser conveniente:

Y(t)=A+KA(do+QmiB(tMETRO))1/ν{\displaystyle Y(t)=A+{K-A \over (C+Qe^{-B(t-M)})^{1/\nu }}}

Esta representación simplifica la configuración tanto de un tiempo de inicio como del valor deY{\displaystyle Y}En ese tiempo.

La función logística , con tasa de crecimiento máxima en el tiempoMETRO{\displaystyle M}, es el caso dondeQ=ν=1{\displaystyle Q=\nu =1}.

Ecuación diferencial logística generalizada

Un caso particular de la función logística generalizada es:

Y(t)=K(1+Qmiαν(tt0))1/ν{\displaystyle Y(t)={K \over (1+Qe^{-\alpha \nu (t-t_{0})})^{1/\nu }}}

que es la solución de la ecuación diferencial de Richards (EDR):

Y(t)=α(1(YK)ν)Y{\displaystyle Y^{\prime }(t)=\alpha \left(1-\left({\frac {Y}{K}}\right)^{\nu }\right)Y}

con condición inicial

Y(t0)=Y0{\displaystyle Y(t_{0})=Y_{0}}

dónde

Q=1+(KY0)ν{\displaystyle Q=-1+\left({\frac {K}{Y_{0}}}\right)^{\nu }}

siempre queν>0{\displaystyle \nu >0}yα>0{\displaystyle \alpha >0}

La ecuación diferencial logística clásica es un caso particular de la ecuación anterior, conν=1{\displaystyle \nu =1}, mientras que la curva de Gompertz puede recuperarse en el límiteν0+{\displaystyle \nu \rightarrow 0^{+}}siempre que:

α=O(1ν){\displaystyle \alpha =O\left({\frac {1}{\nu }}\right)}

De hecho, para los pequeñosν{\displaystyle \nu }es

Y(t)=Yr1exp(νln(YK))νrYln(YK){\displaystyle Y^{\prime }(t)=Yr{\frac {1-\exp \left(\nu \ln \left({\frac {Y}{K}}\right)\right)}{\nu }}\approx rY\ln \left({\frac {Y}{K}}\right)}

El modelo RDE representa numerosos fenómenos de crecimiento que surgen en campos como la oncología y la epidemiología.

Gradiente de la función logística generalizada

Al estimar parámetros a partir de datos, a menudo es necesario calcular las derivadas parciales de la función logística con respecto a los parámetros en un punto de datos determinado.t{\displaystyle t}(véase [ 1 ] ). Para el caso en quedo=1{\displaystyle C=1},

YA=1(1+QmiB(tMETRO))1/νYK=(1+QmiB(tMETRO))1/νYB=(KA)(tMETRO)QmiB(tMETRO)ν(1+QmiB(tMETRO))1ν+1Yν=(KA)ln(1+QmiB(tMETRO))ν2(1+QmiB(tMETRO))1νYQ=(KA)miB(tMETRO)ν(1+QmiB(tMETRO))1ν+1YMETRO=(KA)QBmiB(tMETRO)ν(1+QmiB(tMETRO))1ν+1{\displaystyle {\begin{aligned}\\{\frac {\partial Y}{\partial A}}&=1-(1+Qe^{-B(t-M)})^{-1/\nu }\\\\{\frac {\partial Y}{\partial K}}&=(1+Qe^{-B(t-M)})^{-1/\nu }\\\\{\frac {\partial Y}{\partial B}}&={\frac {(K-A)(t-M)Qe^{-B(t-M)}}{\nu (1+Qe^{-B(t-M)})^{{\frac {1}{\nu }}+1}}}\\\\{\frac {\partial Y}{\partial \nu }}&={\frac {(K-A)\ln(1+Qe^{-B(t-M)})}{\nu ^{2}(1+Qe^{-B(t-M)})^{\frac {1}{\nu }}}}\\\\{\frac {\partial Y}{\partial Q}}&=-{\frac {(K-A)e^{-B(t-M)}}{\nu (1+Qe^{-B(t-M)})^{{\frac {1}{\nu }}+1}}}\\\\{\frac {\partial Y}{\partial M}}&=-{\frac {(K-A)QBe^{-B(t-M)}}{\nu (1+Qe^{-B(t-M)})^{{\frac {1}{\nu }}+1}}}\\\end{aligned}}}

Casos especiales

Las siguientes funciones son casos específicos de las curvas de Richards:

Notas a pie de página

  1. Fekedulegn, Desta; Mairitin P. Mac Siurtain; Jim J. Colbert (1999). "Estimación de parámetros de modelos de crecimiento no lineales en silvicultura" (PDF) . Silva Fennica . 33 (4): 327–336 . doi : 10.14214/sf.653 . Archivado del original (PDF) el 29 de septiembre de 2011. Recuperado el 31 de mayo de 2011 .

Referencias

  • Richards, FJ (1959). "Una función de crecimiento flexible para uso empírico". Journal of Experimental Botany . 10 (2): 290– 300. doi : 10.1093/jxb/10.2.290 .
  • Pella, JS; Tomlinson, PK (1969). "Un modelo generalizado de producción de stock". Bull. Inter-Am. Trop. Tuna Comm . 13 : 421– 496.
  • Lei, YC; Zhang, SY (2004). "Características y derivadas parciales del modelo de crecimiento de Bertalanffy-Richards en silvicultura" . Análisis no lineal: modelado y control . 9 (1): 65–73 . doi : 10.15388/NA.2004.9.1.15171 .