ajustement d'une régression par bootstrap et IC des parametres sur logiciel R

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

Code:

---

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:

---

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

---

Mes résultats sont les suivants :

Code:

---

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.

Bonjour,

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


Cordialement,


A.D.

Merci
Cordialement

Citation:

---

Envoyé par manu-111

Merci
Cordialement

---

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

Citation:

---

Envoyé par fafabzh6

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.