R为 n 维观测值向量（比如：R = 31, 25, 52, 36, .....）

X Y 均为已知的 n 维向量

T = at2 +bt+c （t=1,2，.....，n； a、b、c为待求的系数，均为标量）

Rg= f(X,Y,T) + E (E也是n维向量，是残差项,服从均值为0，标准差为 d 的正态分布；Rg是R的估计值)

已知 Rg 服从均值为 f(X,Y,T) ，标准差为 d 的正态分布。

问题是：要用极大似然法求 a b c d 的值。请问，R 语言里如何实现？