Bonjour,je cherche a faire une régression non linéaire selon une sigmoïde
donc je commence ce script:
Code : Sélectionner tout - Visualiser dans une fenêtre à part
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
 
##################################################
####tentative de fit R stage M2 psychophysique#### 
##################################################
#
# ouvrir les données
data1=read.table('donnearrangec.txt',h=T) 
#
#garder les données actives
attach(data1)
#
#afficher la tableau de données
data1
#
#tentative de fit sur un sujet
x = c(X0[1],X50[1],X70[1],X90[1],X110[1],X130[1],X150[1],X200[1],X250[1],X300[1],X350[1])
y = c(X0[2],X50[2],X70[2],X90[2],X110[2],X130[2],X150[2],X200[2],X250[2],X300[2],X350[2])
plot(y~x)
test =  nls(y ~ a/(b + exp(b-(c*x))))
qui lorsque je le lance me lance : ce message d'erreur
Code : Sélectionner tout - Visualiser dans une fenêtre à part
1
2
3
4
5
6
7
Erreur dans c * x : argument non numérique pour un opérateur binaire
De plus : Message d'avis :
In nls(y ~ a/(b + exp(b - (c * x)))) :
  Pas de valeurs initiales fournies pour certains paramètres.
Initialisation ‘a’, ‘b’ à '1'.
Spécifiez 'start' ou utilisez un model de type 'SelfStart'
J'aimerais savoir comment est la syntaxe,car après en avoir testé plusieurs,je n'arrive toujours à rien...
Comment utiliser un modèle de type selfstart(si il est correct)?
Comment serait la syntaxe correcte (je 'arrangerais pour les paramètres de départ)?
Comment afficher un plot avec le fit associé? et récupérer les valeurs
Merci de l'attention porté à ce post en espérant que je ne vais pas trop exaspérer les connaisseurs