Метод, предложенный JD Long, будет извлекать плотность только для точек данных в существующем обучающем наборе.Если вам нужна плотность в других точках (условные или условные переменные), вам нужно использовать функцию predict()
.Следующий код извлекает и строит одномерное распределение плотности, обусловленное годом == 1999, значением, которое не содержится в исходном наборе данных.
Сначала создайте фрейм данных с теми же компонентами, что и набор данных Италии,с регулярным интервалом gdp и с "1999" - упорядоченным множителем.
yr1999<- rep("1999", 100)
gdpVals <-seq(1,35, length.out=100)
nD1999 <- data.frame(year = ordered(yr1999), gdp = gdpVals)
Затем используйте функцию предсказания для извлечения плотностей.
gdpDens1999 <-predict(fhat,newdata = nD1999)
Следующий код строит график плотности.*
plot(gdpVals, gdpDens1999, type='l', col='red', xlab='gdp', ylab = 'p(gdp|yr = 1999)')