# The Old Faithful geyser data
d <- density(faithful$eruptions, bw = "sj")
> head(d$x)
[1] 1.179869 1.188363 1.196857 1.205350 1.213844 1.222338
Я использую функцию density
в {stats}
, и мне интересно, можно ли увидеть плотность при определенных значениях на выходе? Например, в настоящее время у меня есть оценки плотности при значениях извержения [1] 1.179869 1.188363 ...
, но что, если я хочу узнать оценки плотности при значениях извержения 1 2 5 10 ...
? Есть ли способ извлечь эти объекты density
, d
?
interp
для компонентовx
иy
возвращаемого объекта... - person Ben Bolker   schedule 03.08.2017x <- 2; approx(d$x,d$y,xout = x)
? - person AK88   schedule 03.08.2017d$x[3]
для третьего значения. илиd$y[5]
для пятого результата местоположения y, вы можете извлечь любое значение таким образом. - person sconfluentus   schedule 03.08.2017approx
будет в порядке. Даже в документацииstats::density
упоминается, что они будут использоватьapprox
, если вы запросите вычисление плотности более чем в 512 точках. - person AEF   schedule 03.08.2017