簡單線性迴歸分析

使用包basicTrendline。

安裝包

install.packages("basicTrendline")

查看已經安裝的包

installed.packages("basicTrendline")

調用包

library("basicTrendline")

載入數據

mydata<-read.table("C:/Users/Administrator/Desktop/11.csv",sep=",")

爲了方便調用，使用attach方法

attach(mydata)

線性迴歸

lm.model<-lm(V2~V1+1) #有截距的形式
summary(lm.model) #查看模型

lm.model<-lm(V2~V1-1) #沒有截距的形式，即y=ax+b中，b=0
summary(lm.model) #查看模型

結果

Call:
lm(formula = V2 ~ V1 - 1)

Residuals:
    Min      1Q  Median      3Q     Max 
-2.7126 -0.6409 -0.2093  0.4317  4.4799 

Coefficients:
    Estimate Std. Error t value Pr(>|t|)    
V1 0.1710027  0.0001303    1312   <2e-16 ***
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

Residual standard error: 0.8798 on 1199 degrees of freedom
Multiple R-squared:  0.9993,    Adjusted R-squared:  0.9993 
F-statistic: 1.722e+06 on 1 and 1199 DF,  p-value: < 2.2e-16

使用basicTrendline包的迴歸擬合

• 這裏不能自定義。（或者說還不會）

trendline(x, y, model="line2P", ePos.x = "topleft", summary=TRUE, eDigit=5) #自動添加95%置信區間lines and fill color

trendline(V1, V2, model="line3P", CI.fill = FALSE, CI.color = "black", CI.lty = 2, linecolor = "blue") #只添加95%置信區間的lines，不fill color (set CI.fill = FALSE)

trendline(V1, V2, model="log2P", ePos.x= "top", linecolor = "red", CI.color = NA) #只繪製迴歸曲線，不添加95%置信區間 (set CI.color = NA)

trendline(V1, V2, model="line3P", show.equation = TRUE, show.Rpvalue = FALSE)  #顯示方程，不顯示R值和P值 (set show.Rpvalue = FALSE)

trendline(V1, V2, model="line3P", xname="a", yname=paste(beta^15,b), yhat=FALSE, Rname=1, Pname=0, ePos.x = "bottom") #自定義方程中的參數的名稱:‘xname’, ‘yname’, ‘yhat’, ‘Rname’, ‘Pname’

trendline(V1, V2, model="power2P", ePos.x = "topleft", summary=TRUE, eDigit = 3, eSize = 1.4, text.col = "blue") #改變方程的 小數位，字體顏色，字號大小

trendline(V1, V2, model="power2P",ePos.x = NA) #不顯示方程，只顯示迴歸曲線 (set ePos.x = NA)

#設置繪圖區大小
par(mgp=c(1.5,0.4,0), mar=c(3,3,1,1), tck=-0.01, cex.axis=0.9)
trendline(V1, V2)

• R語言中多使用數據框這種數據結構
• R語言中\代表轉義符，/代表文件夾地址。