สร้างกราฟเส้นตรงด้วย R (How to make linear graph in R)

    กราฟเส้นตรงมีประโยชน์หลายอย่าง เช่น การสร้างกราฟมาตรฐาน (calibration curve) ในการวิเคราะห์ด้วยเครื่องมือ สมมติว่าเราวัดค่าการดูดกลืนแสงสารละลายกลูโคสความเข้มข้นต่างๆ ด้วย spectrophotometer ที่ความยาวคลื่นแสงค่าหนึ่ง โดยวัด 2 ซ้ำ ได้ค่าดังนี้

บน R console

> file <-file.choose()
> data <- read.csv(file)
> data
     conc   Abs
1   16.64 0.163
2   33.28 0.328
3   66.56 0.677
4   99.84 1.004
5  133.13 1.350
6   16.64 0.167
7   33.28 0.351
8   66.56 0.675
9   99.84 0.998
10 133.13 1.347

> plot(data, xlab="Concentration [mg/L]",ylab="Absorbance",las=1)

   สร้างสมการเส้นตรงด้วย linear regression

> linear <- lm(Abs~conc,data=data)
> linear

Call:
lm(formula = Abs ~ conc, data = data)

Coefficients:
(Intercept)         conc 
 -0.0002884    0.0101057  

ลากเส้นตรงจากสมการ

> abline(linear$coef)

   แสดงสมการบนกราฟ

> txt <-paste("Y= ",format(linear$coef[1],digits=3)," + ",format(linear$coef[2],digits=3),"x","\n", "R-squared = ",format(summary(linear)$r.squared,digits=4))
> legend (55,0.4,txt,bty="n")