A blog for physicists or biologists. Mainly for experimentalists interested in Statistics. R scripts detailed and explained.
Associate Professor at Université de Paris
How to make a very simple linear fit with R (few points, no error bars)
rm(list = ls()) #remove ~ everything in the working environment myX<-1:10 myY<-jitter(1:10) adjust<-lm(myY ~ myX) # Use adjust<-lm(myY ~ 0 + myX) to force the intercept at 0 plot(myX,myY,abline(adjust))