Regression un peu bizarre
Mon problème est le suivant:
J'ai des observations xi dans l'intervalle [0,1] et yi réels, i = 1,...,n.
Je veux faire une fonction sous R qui ajuste le modèle suivant (en utilisant le maximum de vraisemblance et supposant que l'erreur est normale) aux données (xi,yi):
Code:
Y <- exp(a + b*sin(2*pi*(X)))
Sachant que cette fonction doit retourner les paramètres a, b et sigma, les valeurs ajustées et les résidus, est-ce que quelqu'un a une petite idée de ce à quoi ça pourrait rassembler ?
Merci par avance,
