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.