2011-08-19 123 views
-1

我現在對R有一些困難。核心密度估計(R)中的恆定帶寬

我有什麼:

  • 與每個三角形內核
  • 與最大值
  • 與「中間」點向量(每個三角形的挑選應該向量極小值的矢量定)

我需要什麼:

  • 密度估計(三角內核OFC)
+2

授可重複使用的小例子,以及迄今爲止您嘗試的內容。這樣,你有更好的機會收穫(好)答案。 –

+0

一個可重複的例子是這裏解釋:http://stackoverflow.com/questions/5963269/how-to-make-a-great-r-reproducible-example。你需要的是一些幫助文件,比如''density'和''kernel'。 –

回答

0

看看你是否可以建立在不同的帶寬規格,基本密度函數與triagular內核的這個例子:

data("faithful", package = "datasets") 
x <- faithful$waiting 
opar <- par(mfcol= c(2,2)) 
for (bw in c(3,6,12,18)) { 
    hist(x, xlab = "Waiting interval", ylab = "Eruption Freq", 
    probability = TRUE, main = "Density plot: Triangular kernel", 
    border = "gray", sub= bquote("Bandwidth = "*.(bw))) 
    lines(density(x, width = bw, window = "triangular"), lty=3, col="red") 
    } 
par(opar) 

Four bandwidths to triangular kernel density plot