Посмотрите, можете ли вы опираться на этот пример изменяющихся спецификаций полосы пропускания для функции базовой плотности с триагулярным ядром:
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)