廿TT

譬如水怙牛過窓櫺 頭角四蹄都過了 因甚麼尾巴過不得

カーネル密度推定で推定された密度からのリサンプリング

R です。乱数です。

何の意味があるのかわかりませんけれども。

一次元です。

set.seed(1)
dens1 <-density(faithful$eruptions)
r1 <-sample(dens1$x,1000,prob=dens1$y,replace = TRUE)
hist(r1,freq = FALSE)
lines(dens1,lwd=2)

f:id:abrahamcow:20170909004252p:plain

二次元です。

library(MASS)
set.seed(2)
dens2 <-kde2d(faithful$eruptions,faithful$waiting)
densdf <-con2tr(dens2)
r2 <-densdf[sample.int(nrow(densdf),1000,replace = TRUE,densdf$z),1:2]
contour(dens2)
points(r2,col=rgb(0,0,0,0.3),pch=16)

f:id:abrahamcow:20170909004312p:plain