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)
Partager