Modèle quadratique proc AUTOREG

Bonjour cher réseau,
j'ai une base de données formées de vaches, de leur rang de lactation et de leurs production par rapport au temps(en jour).
vache Rang Milk Jour
1 1 12 1
1 1 14 2
1 1 10 3
1 2 15 1
1 2 20 2
2 6 13 1
2 6 10 2
2 7 12 1
2 7 13 2
Je voudrais appliquer sur mes données un modèle quadratique que j'ai trouvé dans un article scientifique pour avoir un modèle de production de lait cumulatif.
le modèle est le suivant y=bt+at²+e
Y est la production de lait cumulative par rapport aux individus et au rang de lactation.
t est le temps en jours,
a, b sont des constantes et e l'erreur résiduelle.
comment je peux résoudre les paramètres de ce modèle?
et je voudrais calculer la première dérivé:
dy/dt=b+2at
et la dérivé seconde:
d(dy)/d(dt)=2a.
A noter que l'article a mentionné que la procédure utilisée était AUTOREG (sas 1999).
Je suis débutante en sas et j'ai des connaissances superficielles concernant la modélisation, donc je vous de demande de bien vouloir m'aider.
Merci d'avance.
Bien cordialement.