たまに必要になるんだけど、毎回書き直しているのでメモっておく。 テストコードの結果はR言語と揃ってるのは確認済。 #include <iostream> #include <cmath> #include <vector> #include <list> #include <numeric> template <template<class T, class Allocator = std::allocator<T> > class Container> double mean(Container<double> & x) { return std::accumula</double></template<class></numeric></list></vector></cmath></iostream>…

この前のエントリーのために作ったグラフが結構気に入ったのでテンプレとして残しておく。例えば対数正規分布するデータを1000個用意して、 x <- exp(rnorm(1000)) これをPLOT。 library(ggplot2) qplot(x, geom="blank") + geom_histogram(aes(y=..density.…