Discussion :

[manu-111]

Bonjour à tous,
Voici les paramètres d'une droite décroissante, chaque xy correspondant a un intervalle:

Code :

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

1
2
x=c(60,59,55,49,43,38,33,28,24,20,17,15,12,10,8,6,5)
y=c(59,55,49,43,38,33,28,24,20,17,15,12,10,8,6,5,4)

Je l'ai ajustée avec une polynomiale d'ordre 3 (pour l'exemple car ce n'est pas le meilleur ajustement...) :

Code :

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

lm(formula = y ~ x + I(x^2) + I(x^3))

Mes résultats sont les suivants :

Code :

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

1
2
3
Coefficients:
(Intercept)            x                       I(x^2)            I(x^3)  
 -1.5558172    1.0738235   -0.0101314    0.0001425

1. Maintenant je voudrais les renvois des valeurs médianes de l'intercept, x, l(x^2) et l(x^3) après boostrap. C'est-à-dire que je veux faire passer 1000 droites (par exemple) entre le premier intervalle 60-59 entre le second 59-55, etc...
2. Je voudrais également le R² ainsi que l'erreur de prédiction pour chacun des coefficients

Merci d'avance pour la routine.

[A. D.]

Bonjour,

Pour répondre au point 2, voir par exemple la fonction summary.


Cordialement,


A.D.

[manu-111]

Merci
Cordialement

[fafabzh6]

Merci
Cordialement

Que veux-tu exactement, le nom d'une fonction R te permettant de faire ça ; ou un pseudo algorithme pour le faire ?

[manu-111]

Que veux-tu exactement, le nom d'une fonction R te permettant de faire ça ; ou un pseudo algorithme pour le faire ?

Merci pour ton intéret Fafa, je voudrais la routine complète à partir de cet exemple pour bien identifier tous les éléments et pouvoir refaire tourner ensuite avec d'autres jeux de données. Je ne suis pas matheux et j'apprends par l'exemple.