Оценка плотности в R

Я оценил плотность старых верующих, используя:

 library(np)
 library(datasets)
 data(faithful)
 f <- npudens(~eruptions+waiting,data=faithful)

Как найти плотность точки, не соответствующей действительности? Спасибо!


r kernel-density
person Michael Burrows    schedule 17.02.2014    source источник

Ответы (1)


arrow_upward
4
arrow_downward

?npudens говорит, что у npdensity объектов есть predict метод:

predict(f,newdata=data.frame(eruptions=3.5,waiting=70))
## [1] 0.004473444
person Ben Bolker    schedule 17.02.2014
comment
Спасибо. Я не смог найти никакого документа для предсказания().npudens, поэтому я использовал предсказать().lm, в котором не упоминается, что новые данные должны быть именованным аргументом. У меня появляется сильное ощущение, что мне есть чему поучиться в поиске и чтении R doc. - person Michael Burrows; 18.02.2014