実用性があるのか分からない機能を紹介？ progs.sas dm 'log; clear; output; clear'; proc datasets kill; run; option linesize = 130 pagesize = 500; %let execpath = " "; %let Path = " "; %macro setexecpath; %let execpath = %sysfunc(getoption(sy…

リサンプリングとかブートストラップで統計量の分布を作って検定とか推定をするのに使えます。 /* ダミーデータの作成 */ data mult; call streaminit(070605); do y = 0, 1; do i = 1 to 100; x = rand('Normal') + y; output; end; end; run; /* multtest …

StatLib---Applied Statistics algorithms Applied Statistics(JRSS C)に載ったアルゴリズムのソースがみられる。 ほとんどFortran？

2〜20次元関数のAdaptive Quadrature(適応的求積法)ができるパッケージ。 library(adapt) # Example of p - dimensional normal distribution mu <- rbind(c(0, 0, 0)) sigma <- rbind(c(1, 0, 0), c(0, 1, 0), c(0, 0, 1)) pdfnm <- function(z) { 1 / ( (2…