2008-09-20から1日間の記事一覧

最急降下法で最適なバンド幅を計算してみた

ナダラヤワトソン推定量は計算できたわけですが、これにはまだバンド幅hを決めるという問題が残っています。バンド幅を最適なものにする*1ための一つの尺度として、最小二乗クロスバリデーション評価関数を最小にするようにhを選択する、という方法がありま…

最小二乗クロスバリデーション評価関数を最小化のためのコード

これをとりあえず最急降下法で出してみることにする。 y.exp <- function(Xj,Yj){ y.exp.intern <- function(Xj,Yj){ return(paste("exp(-1/2 * ((Xi - ",Xj,")/h)^2) * ",Yj,sep="")) } paste(paste(mapply(function(Xj,Yj){y.exp.intern(Xj,Yj)},Xj,Yj),co…

数式微分怖いです

R

ひゃー。