関数内関数とmapplyの付近でなんだかよく分からない挙動になっている。まず、以下のような関数を用意。 hoge <- function(x1,x2){ return(function(y){ cat(x1,x2,fill=TRUE) x1*x2*y}) } 内側のmapplyで関数を生成、外側のmapplyでそれを評価している。 > m…

先週は、1変数でのナダラヤワトソン推定量を計算していました。 オレオレ平滑化をやってみた - Seeking for my unique color. 最急降下法で最適なバンド幅を計算してみた - Seeking for my unique color. 今週は先生に「多変数でやってみな!!」と言われてい…

2変数というとで先週書いた密度トレイスを参考にしました。が、微妙に間違いも発見しました。 カーネル関数を使って密度トレイスを書いてみた→最後のほうに微妙に間違いがあるので注意!! - yasuhisa's blog gauusian関数の仲の1/hにかかる括弧が抜けていた。…

Tsukuba.Rで3次元のplotは練習しておいたから楽勝。とりあえず、MASSパッケージのkde2d関数を自前で実装することに挑戦するのが今の課題。 #2変数のガウス分布のplot double_gaussian <- function(x1,x2,sigma=1,mu1=0,mu2=0){ return(1 / (2*pi*sigma^2) * …