Trouver les coefficients d'une fonction logistique-débutant

[Anne Didier]

Bonjour,

Je dispose de points expérimentaux dont les X (temps) et Y (masse) sont connus.
J'aimerais faire passer par ces points une fonction logistique du type : Ysim = A / (1+exp(4*Vm*(Pinf-X)/A)) et récupérer les valeurs des coefficients A, Vm et Pinf.
Je voudrais que la logistique démarre forcément à 0.
Pour le moment j'avais travaillé dans Excel, en utilisant le solveur permettant de réduire la somme des carrés des écarts, mais je n'arrive pas à de très bons résultats.

Je ne sais pas vraiment comment m'y prendre pour le faire avec R...

Merci d'avance pour votre aide!

[Stéphane Laurent]

Bonjour,

Avec la fonction nlme() du package nlme. Voir le livre Pinheiro, J.C., and Bates, D.M. (2000) "Mixed-Effects Models in S and S-PLUS".

[Pierre Bady]

bonjour,


fonction nls pour les modèles non-linéaires

un petit document de John Fox sur le sujet:
http://cran.r-project.org/doc/contrib/F ... ession.pdf

la fonction nlme s'utilisent généralement pour les modèles non-linéaires mixtes.


HTH

pierre

[Anne Didier]

Bonjour,

Merci pour vos réponses mais je ne suis pas assez calée en modélisation/stats pour pouvoir comprendre et utiliser seule cette fonction. je vais essayer de me débrouiller autrement.

[Stéphane Laurent]

Pierre a raison, c'est nls().

Anne, tu as tort de chercher à faire autrement, car nls() est la "Rolls Royce" pour faire ce que tu veux. Regarde les exemple de l'aide nls() avec SSlogis() :

Code : Tout sélectionner

> library(nlme)
> ?nls
> ?SSlogis

[jean lobry]

Merci pour vos réponses mais je ne suis pas assez calée en modélisation/stats pour pouvoir comprendre et utiliser seule cette fonction. je vais essayer de me débrouiller autrement.

Bonjour Anne,

pas de panique, tu trouveras bien mieux ici que dans le "solveur" d'excell, avec plein de spécialistes prêts à t'aider sur ce forum. Ce qu'ils ont fait en l'occurrence, mais sans paliers de décompression. Pour une introduction en douceur je te redirige vers un document consistant à estimer les paramètres du modèle logistique à partir des données de son inventeur. Juste pour dire que ce n'est pas si terrible R, et qu'il n'y a que là où l'on puisse s'amuser en "modélisation/stats" pour de vrai, parce que c'est un logiciel libre, et parce que c'est intrinsèquement amusant.

Bon courage pour la suite.

Amicalement,

Jean

[Anne Didier]

Merci pour vos nouvelles réponses,
J'ai pu faire ma logistique avec nls grâce à un peu d'aide supplémentaire!