最頻値（モード） mode をRで求める

q<-rnorm(10000)
j<-density(log(q))
qq3[t]<-exp(mean(j$x[which(j$y==max(j$y))]))

plot(density(q))
abline(v=log(qq3[t]),col=3)