Discussion :

[remissssss]

Bonjour a tous,
Je suis nouveau dans le forum et j'espère que je serai suffisament clair pour exposer mon problème !

Voila le topo :

Mon jeu de données est tres simple, je n'ai que 2 vecteurs à gerer, l'un correspond a des temperatures (x) , l'autre (y) à une variable dépendante de la température :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
y<- c(0.06659643, 0.42096700, 0.80442651, 0.61664031, 0.06694109, 0.09041552 ,0.08000000)
x<- c(5,10,15,20,25,30,35)

L'idée est de savoir comment evolue la variable y en fonction de x.
En plottant y~x, j'obtiens une courbe en cloche, pas symétrique, qui à mon avis ressemble beaucoup à une fonction de densité de la loi Gamma (dgamma). Je souhaite donc fitter mes données avec cette fonction là.
Pour trouver les coeff de ma regression, j'ai voulu utiliser la fonction nls(), en lui precisant un modele gamma, je pensais à un truc comme ca :

nls(y~dgamma(x,alpha,beta)), avec alpha et beta les paramètres du modèle.

Mais, je ne sais pas comment fixer alpha et beta a priori.
Avez-vous des solutions?

Merci.

[AMDBG]

Bonjour,

Tu peux utiliser le package fitdistrplus
http://cran.r-project.org/web/packag...lus/index.html

Ensuite la fonction fitdist fait tout le travail

fitdist( y , "gamma")

[remissssss]

Merci,
j'y avais pensé, mais ça n'a pas l'air de donner les bons résultats.

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
> model<-fitdist( y , "gamma")
> summary(model)
Fitting of the distribution ' gamma ' by maximum likelihood 
Parameters : 
      estimate Std. Error
shape 1.100299  0.5221874
rate  3.589157  2.1386908
Loglikelihood:  1.295946   AIC:  1.408108   BIC:  1.299928 
Correlation matrix:
          shape      rate
shape 1.0000000 0.7964532
rate  0.7964532 1.0000000

D'apres ça, ma variable y suit une loi dgamma de paramètre shape = 1.1 et rate = 3.5 .
Donc quand je fais dgamma(15,shape=1.1,rate=3.5), je devrais retrouver ma valeur de y pour x=15, à savoir 0.80442651, mais ce n'est pas le cas :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
 
> dgamma(15,shape=1.1,rate=3.5)
[1] 8.655317e-23

[Sengar]

Bonjour,

La somme des y est supérieure à 1 donc déjà c'est mal parti pour dire que c'est une densité (que ce soit gamma ou une autre).

D'autre part, je n'ai jamais vu un modèle du type y=densitéLoi(x, paramètres).

On fait plutôt des modèles du genre:
y~Loi(paramètres=function(x))+residus

Tu devrais changer de méthode.



Edit: en fait en relisant tes messages j'ai l'impression que tu confonds densité et valeurs prises par une loi.

[remissssss]

salut,
en effet, je ne veux pas que mes valeurs representent une densité, je veux simplement trouver une formule mathematique en y=f(x), qui corresponde au mieux à mes données.
A l'allure de mon plot(x,y), ca m'a fait pensé à la fonction densité de la loi gamma, c'est pour ca.
Quelle est alors la méthode pour trouver la fonction qui s'ajuste au mieux à ma courbe de données?

[Sengar]

Vu ton problème et la forme de ta courbe (peu de points, et correspond pas franchement à une fonction connue) j'aurais dis un spline ou une approximation polynomiale.

Ca tombe bien ca en parle juste le topic d'en dessous x)

http://www.developpez.net/forums/d13...-polynomiales/