從已知回歸方程繪制回歸線

Question

我有以下代碼從所有數據對計算非參數回歸方程（不是 lm）。 在這段代碼中，“mdns”是斜率中值，“mdni”是中值截距，所以我想從散點圖中的數據中繪制點（所有 x，y 對），然后繪制回歸線（1 ) 使用來自原始數據的 lm 和 (2) 使用標准方程 y=mdni+mdns(X) 中計算的非參數回歸值。 我知道如何做所有事情，除了這個在已知方程的圖中畫一條線。

#Calculate pairs of all Slope values
slope<-vector()
 for (i in 1:n){
 for(j in 2:n){ 
   temps<-(data1$y_book[j]-data1$y_book[i])/(data1$x_book[j]-data1$x_book[i])
 if (j>i){  
 slope<-append(slope,temps)}
 }
 }

slope<-slope[!is.na(slope)]
print(slope)
mdns<-median(slope)

#Calculate the Intercept values
inter<-vector()
for (k in 1:n){
    tempi<-(data1$y_book[k] - (mdns*data1$x_book[k]))
    
 inter<-append(inter,tempi)
}
#echo results
print(inter)
mdni<-median(inter)

Answer 1

不知道你在找什么，但有兩種繪制線條的基本方法。

使用 abline()

plot(NULL,ylim=c(-10,10),xlim=c(-10,10)
abline(a=.4,b= 1) # a= slope, b= intercept (see docs for other options)
# Note line extends completely across plot.

使用行（）

plot(NULL,ylim=c(-10,10),xlim=c(-10,10))
x <- seq(from=-9,to=9,length.out=30) 
y <- seq(from=-3,to=3,length.out=30)
lines(cbind(x,y))
# Note line covers only defined extent

請參閱繪圖、線條、abline 文檔。 另見ggplot。