Moindres carrés contraints et ajustement exponentielJe cherche à adapter Y ~ 1-exp {- [C1 * (x / Xmax) ^ 1 + ....... CN * (x / Xmax) ^ N]}. La commande peut être n'importe où entre 6 et 12 (valeur maximale de N). Mais avant de générer l'ajustement, je dois estimer le coefficient à l'aide de la LSQ modifiée. J'ai examiné l'ajustement exponentiel (analyser l'ajustement Y -> X) et essayé d'entrer une formule, mais je ne souhaite pas ajuster directement la courbe exponentielle. la nature strictement monotiquement croissante de Y.. Le plus proche que je puisse trouver qui représente mon modèle est https://www.mathworks.com/help/optim/ug/lsqlin.html . On dirait un processus en 2 étapes avec le LSQ st contraintes linéaires, suivies de l'ajustement exponentiel réel. Je ne sais pas s'il existe une approche en 2 étapes ou une approche droite pour l'ajustement pour une limite supérieure donnée sur le C1..CN.
Tous les indicateurs seraient très utiles.
Merci,
Re: Moindres carrés et ajustement exponentielPouvez-vous linéariser la relation avec les journaux: -Log (1 - Y) = c1 * (x / Xmax) + c2 * (x / Xmax) ^ 2 + ... + cN * (x / Xmax) ^ N?
(Cette méthode peut ne pas être valide, en fonction de la contribution des erreurs à l'expression d'origine.)
Re: Moindres carrés et ajustement exponentielJe l'ai essayé, mais les coefficients estimés sont faussés par rapport à certaines références de la littérature. Il me manque peut-être des limites ou des contraintes qui n'ont pas été entrées correctement. Je vais essayer à nouveau.
J'ai également essayé d'ajouter le modèle personnalisé à la bibliothèque, comme indiqué ici. https://www.jmp.com/support/help/14-2/create-a-formula-using-the-model-library.shtml
mais obtenir une erreur I / O. J'ai des droits d'administrateur qui ne savent pas pourquoi cela se produit. Existe-t-il un mode "édition" ou "entrée utilisateur" qui me permettrait d'ajouter mon modèle non linéaire à la bibliothèque?
Re: Moindres carrés et ajustement exponentielLes estimations se trouveront dans l'espace du journal, pas dans les dimensions d'origine. Cette différence expliquerait-elle les résultats inattendus? En outre, ces résultats ne sont pas censés être les meilleures estimations mais servent uniquement de valeurs initiales pour les paramètres avant l'ajustement des données à la plate-forme non linéaire. Vous devez retransformer les valeurs de départ avec la fonction Exp ().
Je ne connais pas d'erreur I / O. Est-il possible que le fichier de script soit verrouillé?
Je vous suggère de créer d'abord votre modèle personnalisé manuellement / interactivement. C'est-à-dire, définissez d'abord le modèle en tant que formule de colonne. Cela permettra de démontrer qu'un tel modèle est réalisable. La bibliothèque de modèles n'est qu'un moyen de s'adapter à un modèle personnalisé. Ce n'est pas obligatoire. De cette façon, vous obtiendrez également l'expression correcte si vous décidez ultérieurement de l'inclure dans la bibliothèque de modèles.
