Bonjour, je cherchais à tester la fonction lm de R,
mais je ne retrouve pas ce que je mets en entrée :
Je crée une variable a, issue d'une loi normale de moyenne 5 et de variance 1, que je multiplie par x allant de 0.1 à 10.
Ensuite je lance mon lm à partir des resultats et j'obtiens pour valeur de a : 1
cf le code ci-dessous :
Avez vous des idées ?

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
par(mfrow=c(2,1))
x<-rnorm(100,mean=5,sd=1)*seq(0.1,10,by=0.1)
x
y<-rnorm(100,mean=1,sd=1)*seq(0.1,10,by=0.1)
z<-runif(1,min=4,max=5)*seq(0.1,10,by=0.1)
calcul<-x
plot(calcul~seq(0.1,10,by=0.1),type="l")
lesdata=as.data.frame(x)
a<-lm(calcul~.,data=lesdata)
summary(a)
b=a$Intercept
b
a1=a$coefficients[2]
a1
sequence<-seq(0.1,10,by=0.1)
nouveauCalcul=a1*sequence
plot(nouveauCalcul~sequence)