Python中LOWESS的置信區間

Question

我如何計算Python中LOWESS迴歸的置信區間？我想將這些作爲陰影區域添加到使用下面的代碼創建的LOESS plot（其他軟件包比statsmodels也好）。

import numpy as np 
import pylab as plt 
import statsmodels.api as sm 

x = np.linspace(0,2*np.pi,100) 
y = np.sin(x) + np.random.random(100) * 0.2 
lowess = sm.nonparametric.lowess(y, x, frac=0.1) 

plt.plot(x, y, '+') 
plt.plot(lowess[:, 0], lowess[:, 1]) 
plt.show() 

我已經添加與來自webblog Serious Stats（它是使用ggplot創建在R）以下的置信區間的示例曲線圖。

Answer 1

黃土不具有標準錯誤的明確概念。在這種情況下，這並不意味着什麼。既然這樣，你用堅強的方法堅持下去。

引導您的數據。您將要將LOESS曲線擬合到自舉數據。看到這個頁面的中間部分，你可以找到你所做的一些漂亮的照片。 http://statweb.stanford.edu/~susan/courses/s208/node20.html

一旦你有你大量不同的黃土曲線，你可以找到頂部和底部的X個百分點。