Comment ajuster un modèle polynomial de degré 2 avec des sous-groupes puis capturer les estimations de paramètresJe voudrais créer le degré polynomial 2 par chacun (SERIAL, P1, P2, P3) et obtenir le coefficient et la constante puis mettre dans chaque colonne en utilisant JSL.
Les données d'exemple ont été jointes dans ce post, mais en fait je vais faire avec les données environ 1 million de lignes.
Re : Comment créer un polynôme de degré 2 et obtenir le coefficient et la constante ?Salut MathStatChem,
Votre solution fonctionne.
Re : Comment créer un polynôme de degré 2 et obtenir le coefficient et la constante ?Salut ian_jmp, Mon objectif est d'obtenir le coefficient et la constante de "y = ax ^ 2 + bx + c" mais le polynôme en forme est " y = a (xd) ^ 2 + bx + c". Cependant, j'ai obtenu la solution de l'utilisateur # MathStatChem par fit special.Re : Comment créer un polynôme de degré 2 et obtenir le coefficient et la constante ?Merci MathStatChem, laissez-moi essayer votre suggestion.Re : Comment créer un polynôme de degré 2 et obtenir le coefficient et la constante ?Merci ian_jmp, laissez-moi essayer.Re : Comment créer un polynôme de degré 2 et obtenir le coefficient et la constante ?Essayez ce script :
_biv=Bivariate(
Y( :y ),
X( :x ),
Fit Special(
Degree( 2 ),
Centered Polynomial( 0 ),
{Line Color( {213, 72, 87} )}
),
By( :SERIAL, :P1, :P2, :P3 )
);
_rbiv=Report(_biv[1]);
_dtParameterEst = _rbiv[Table Box(3)]<< Make Combined Data Table;
_dtParameterEst << Split(
Split By( :Term ),
Split( :Estimate ),
Group( :P3, :P2, :P1, :SERIAL ),
Remaining Columns( Drop All ),
Sort by Column Property
)Re : Comment créer un polynôme de degré 2 et obtenir le coefficient et la constante ?De plus, si vous utilisez la plate-forme Fit Y by X / Bivariate pour générer les ajustements polynomiaux, assurez-vous de faire un "Fit Special" et d'utiliser des polynômes non centrés. Sinon, la quadratique ajustée sera de la forme
f(x) = a + b*x + c*(x-x barre)
où xbar est la moyenne des valeurs x utilisées dans la régression.
Effectuer l'ajustement avec une variable By, puis cliquer avec le bouton droit sur le tableau des estimations des paramètres pour créer un tableau de données combinées est un moyen très efficace de le faire.
Re : Comment créer un polynôme de degré 2 et obtenir le coefficient et la constante ?Il serait possible de développer le script ci-dessous pour obtenir ce dont vous avez besoin :
NamesDefaultToHere(1);
ClearLog();
dt = DataTable("Example data.jmp");
col = Column(dt, "Predicted y By SERIAL By P1 By P2 By P3");
// Get the formula generated by the platform
f = col << getFormula;
// Start to parse this formula expression
for(g = 1, g <= Narg(f) - 1, g = g + 2,
f1 = Arg(f, g);
f2 = Arg(f, g+1);
Print(f1, f2, "\!n");
);
Mais je soupçonne qu'il pourrait être plus simple d'envoyer le message 'MakeCombinedDataTable' à la table appropriée dans la fenêtre de rapport si cela est possible. Regardez dans 'Aide > Guide de script'.
