2009-02-18から1日間の記事一覧

台形公式で積分

R

数値計算を勉強したことがないのがバレバレですねwww。 my.integrate <- function(f,a,b,n){ ak <- seq(from=a,to=b,length.out=n) sum <- 0 for(i in 2:n){ sum <- sum + (ak[i] - ak[i-1]) * (f(ak[i-1]) + f(ak[i])) / 2 } return(sum) } よしよし。 my.i…

Confidence interval for kernel density estimator

R

お、なんか大分間が空いた気がする。修論スタート、ということでノンパラな推定量に関しての信頼区間の付近を準備しています。とりあえず今は一番簡単そうなカーネル密度推定の1変数のバージョンの信頼区間。週末とかに「分からねー」を連発しながらやってた…