Je me prépare à exécuter une analyse de dégradation paramétrique par chemin non linéaire.
J'ai un paramètre, n, que je veux adapter globalement, mais un autre paramètre b0 je veux ajuster par ID système
Donc, je veux la deuxième partie de ceci (l'ajustement global à n):

mais le premier paramètre b0, a la forme suivante:

Je me suis moqué de ce que je serais idéalement capable de faire:

où le paramètre b0 est indexé par DUT mais le paramètre n ne l’est pas.
Toute aide serait grandement appréciée.
Bonne idée ... malheureusement, certaines des données en cours d’ajustement contiennent des valeurs négatives.
En fait, j'ai compris que je devais comparer le script des paramètres globaux b0 et n à celui du script avec b0 et n FIT par ID système. Ensuite, j'ai essayé de fusionner les deux scripts ... et cela a fonctionné! (voir ci-dessous)
le rouge est l'endroit où j'ai apporté des modifications ... Tout d'abord, supprimez l'index [DUT] de n, puis réduisez simplement les listes pour qu'elles atteignent les 8 objets traités de b0 et une seule instance de n.
Dégradation(
Y (: shift), Heure (: durée), Label (: DUT), Limite supérieure de spécification (20), Limite de spécification inférieure (-20),
Application (dégradation de mesures répétées), Connecter des marqueurs de données (1), Afficher les lignes ajustées (0), Afficher les limites de spécification (1),
Afficher les courbes médianes (0), Montrer la légende (0), Aucune liste d'onglets (0), Définir la limite de spécification supérieure (20), Définir la limite de spécification inférieure (-20), Définir le temps de censure (. ),
Montrer le tracé résiduel (1), Afficher le tracé de prédiction inverse (1), Intervalle de prédiction inverse (pas d'intervalle), Prédiction Inverse Alpha (0,05),
Spécifications de chemin (
Linéaire simple (
Ajouter une échelle X personnalisée (
{{"Linéaire", Fonction ({x}, x), Fonction ({x}, x)}} ),
Ajouter une échelle Y personnalisée ( {{"Linéaire", Fonction ({x}, x), Fonction ({x}, x)}} ),
Pente (différente), Intercept (différent), Sélectionnez X Scale ("Linear"), Sélectionnez l'échelle Y ("linéaire") ),
Chemin non linéaire ( Ajouter une formule (
Nom de formule ("ShiftModel"), Formule(
Paramètre( {b0 = 1, n = 0,6},
durée ^ n * b0 [DUT] ) ),
Valeurs initiales( [ 0,009, 0,017, 0,026, 0,007, 0,015, 0,025, 0,019, 0,008, 0,600 ] ),
Inférieur( [ -0,016, -0,009, -0,002, -0,017, -0,011, -0,004, -0,008, -0,015, -12,797 ] ),
Plus haut( [ 0,035, 0,043, 0,054, 0,033, 0,041, 0,054, 0,047, 0,032, 13.998 ] ),
Méthode d'adaptation (QuasiNewton BFGS),
Fixé( [ 0, 0, 0, 0, 0, 0, 0, 0, 0 ] ) ),
Sélectionner une formule ("ShiftModel") ) ),
Chemin non linéaire (1), Chemin moyen (1),
SendToReport ( Envoi( {"Superposition", "Spécification du modèle"}, "ShiftModel", TextEditBox, {Taille fixe (1, 102, 16)} ) ) )
Je comprends votre première déclaration et votre image d’un chemin de dégradation non linéaire. Je ne comprends pas votre deuxième exemple, mais ça va.
Que faire si vous utilisez une transformation de journal? Ensuite, votre modèle: Durée ^ n * b0 (DUT) devient linéaire comme suit: log (b0 (DUT)) + n * Log (Durée). Utilisez une interception différente et une option de pente commune pour le chemin linéaire. Puis revenez transformer si vous voulez.