Discussion :

[Invité]

Bonjour,

Je suis actuellement en train d'effectuer une analyse et je voudrais fixer un modèle sous la forme d'une courbe logistique modifié ou d'une courbe de Gompertz. Suite à ce modèle, je voudrais calculer une aire sous la courbe à des pas fixes.

Mon jeu de donnée est de la forme :

variete	parcelle	a	LAI1	LAI2	LAI3	LAI4	LAI5	ST1	ST2	ST3	ST4	ST5
1	113	3.8	2.893193333	2.632766667	2.641633333	2.206786667	0.609266667	212.7	304.2	419.35	502.2	630.6
1	203	4.3	3.13298	3.13298	3.002833333	2.430933333	0.682123333	212.7	304.2	419.35	502.2	630.6

Le tableau est construit de la façon suivante: ST est en abscisse et LAI en ordonnée. J'ai fait un graphique par parcelle (facteur).


L'expression utilisé pour le modèle est

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
 
nls(LAI1+LAI2+LAI3+LAI4+LAI5 ~ a/(1+exp(-K*((ST1+ST2+ST3+ST4+ST5)-m))), start=list(K=0.005,m=600))

Ce que je voudrais faire est d'obtenir pour chaque parcelle la valeur de K et m ainsi que de calculer l'air sous la courbe (avec un pas régulier, ex tout les 50 ST). J'ai pensé à le faire avec une boucle for mais je ne sais pas comment le faire.
Je veux obtenir ce genre de courbe:


Merci pour votre aide d'avance.

Simon