いまだに概収束と確率収束の違いがよく分かってなくって、ぐぐっていると「概収束と確率収束の違いがよく分からないんだよねー」という自分のBlogが検索結果に出てきてイライラするのでまとめてみることにしました。あってるか知らん。以前のちょっとまとめ…

K近傍法による密度推定はいまいちな気がしたので、K近傍法によるクラス分類を行なうことにしました。多クラスでもいいんだけど、まあ2クラスから始めます。クラスに帰属する事後確率が(P123の(2.256式))であることを利用しています。自分で適当にデータ生成…

PRMLの2章の最後の付近にK近傍法というのが出てきました。 www.yasuhisay.info で、僕が卒研でたぶんやるカーネル使ったやつに比べてやったら考え方がシンプルなので、とりあえず実装してみることにしました。密度トレイスということなので、1次元の球を考え…

実は、RでK近傍法を用いたクラス分類をするためのコードを書いてみた - yasuhisa's blogのやつは小さな正方形に区間をくぎって、そこでぬりつぶす、という作業をやっています。Rでぬりつぶし、といったらpolygon関数なんですが、あんまり使い方覚えてないの…

場所から内容からなんかすごいぞ、これは。 Tsukuba.R#3 - Tsukuba.R - Seesaa Wiki（ウィキ） 広報用のテンプレをid:wakutekaが作ってくれたので、貼りますー。 Tsukuba.R #3 〜つくば周辺(半径2万キロくらい)をターゲットとしたRの勉強会〜◆日 時： 2008年…

Perlでデバッグをやろうとしていたら、cpanでつっこんだモジュールを探せないと怒られた。terminalとかでやるとちゃんと動くのに、Emacs上でデバッグすると言われる。悩んでたら、@elimさんに教えてもらった。 https://twitter.com/elim/statuses/955405133 …

yasnippetってまず何って話なんですけど、なんかテンプレートを挿入してくれるelispらしい。ifとか書いてtabとかを押すと if () { } という風に展開してくれます。分かんねーよっていう人は下のやつ見るといいよ。 http://youtube.com/watch?v=vOj7btx3ATg …

Emacsで編集してて、文法エラーみたいな単純なやつは実行とかコンパイルする前に指摘して欲しい。そんな欲求を満たしてくれるのがflymakeだそうです。エラーがあると下みたいに指摘してくれる。 Eclipseみたいなやつですね。で、flymakeはPerlとかRubyとかに…

@wakuteka先生の一言により、Tsukuba.R#3が開催されることになりました。Tsukuba.Rはつくば周辺(半径2万キロくらい)をターゲットとしたRの勉強会みたいなのをやっているものなんですが、Tsukuba.R#2でも宣言した通り「つくば→東京→つくば→東京→…」という感じ…

今週で二章が終わりました。確率分布の章はわりとスムーズに進んだような気がしなくもない。上巻が5章くらいあるので、そういう意味では40%を読んだことになります。わっしょい。混合ガウス分布正規分布の和の分布と、混合正規分布が昔ごっちゃになっていた…

3週間くらい空いてしまって久し振りな感じです。続けますよー。継続重要。あとでログを上げる。 加算集合 加算集合の性質 連続の濃度、非可算集合 羃集合の濃度

Limechatとかで検索ができなさげなので、もうbotを作ることにした。channelにfiltterをかけたchannelみたいなのってすぐに作れないもんなのかなー?tig.rbでtwitterの発言とかはrieceのlogにたまってきている。だから、定期的にそのlogを特定userとかでgrepみ…

結構色々あって面白い。自分が昔書いたやつが必要なくなる感じだ。 www.yasuhisay.infoインスールなど install.packages("sets") library(sets) #library(help=sets) 集合の基本演算Rのbaseパッケージ(かな)にも集合演算はあるんだけど、listを使わないとい…

メモメモ。てか、RubyKaigiとかと比べても発表数めっちゃ多いんだが、これ全部発表あったのかな? useR! 2008: Abstracts and Slides 適当に見ていったやつ Good Relations with R 集合とか扱えてるっぽい。自前でやろうとしていたw。ファジィーなのとかもあ…

以前カーネル単回帰のところで最小二乗クロスバリデーション評価関数 を最小にするhを決めるというのをやりました。「カーネル重回帰でも最小二乗クロスバリデーション評価関数を最小にするのがある意味(hogehogeに確率収束する、など)においてよい」という…

関数内関数と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) * …

8月から始めたskypeによるパターン認識と機械学習 上 - ベイズ理論による統計的予測の読書会も、こつこつとやってきて今回が第十回目です。予習とか全くなしで、2時間くらいやって10Pとかいう進行度合いですが、10回もやると100Pくらいに到達しました。この…

たてこみまくっている。。。 卒論 カーネル重回帰のところで最適なhのベクトルを決めるところの問題は、とりあえず本に載っていなかったので今週はやらなくてもよさげになった。まあ、いつかはやらないといけないんですが。。。で、たぶん最小二乗CV評価関数…

今日は「Hedging Strateries Using Furures」ということで先物を使ったヘッジ戦略というところらしいですです。やっべ、来週俺だわwww。 longとshortでのhedge 価格を狙いさだめることができる hedgeのいい点 会社がリスクを最小化したい 色んな変動のリスク…

「それLimechatでできるよ」なんだけど、Rieceでやってみた。まず、~/.riece/initファイルを編集。 (setq riece-server "ircnet" riece-server-alist '(("ircnet" :host "irc.tokyo.wide.ad.jp") ("freenode" :host "chat.freenode.net" :coding utf-8) ("tw…

環境はleopardです。ハイライトとかじゃなくって、検索キーワードとかを投げたら、そこのwindowで一致するやつだけ表示されるとか。logのファイルのありかとか、それはテキストファイルってことは分かっているので、 perl -ne 'print $_ if$_ =~ /T_Hash/;' …

まあ、誰か書いてるだろうけど。 require 'rubygems' require 'mechanize' url = "http://see.stanford.edu/see/materials/icspmcs106a/handouts.aspx" agent = WWW::Mechanize.new source = Hpricot(agent.get_file(url)) (source/"a[@href^='http://course…

今日はchaper 2の「mechanics of futures markets」っていうところ。 広い資産がやり取り可能 取り引き所取り引き settled daily 値洗い 重要らしい 日々決算をする形になるらしい margins 証拠金 現金か有価証券 プレイヤーがブローカーに払うやつ brokerっ…

「@reboot」を付けると、立ち上げたときだけ実行というのができるらしい。 /Users/yasuhisa% crontab -l HOME=/Users/yasuhisa @reboot /opt/local/bin/ruby $HOME/net-irc/examples/gmail.rb --debug @reboot /opt/local/bin/ruby $HOME/net-irc/examples/t…

この前書いてたやつに関して。applyでデータの列ごとに処理していて、列の名前も含めたいよねという状況を考えている。 applyを使ってると、data.frameのカラム名を触れない - yasuhisa's blog なんか素晴らしくよく分からないことやっているような気がする…

なんだかよく分かんなくなってきはじめたところ。極限操作に微妙に自信が持ててないところとか、線形代数忘れまくってるところとかやばいですね。Zoutendijkの条件「ズーテンダイク」と読むらしいです。ズーテンダイク!!(何。 リプシッツ連続 リプシッツ定数…

昔Rubyではてダラのファイル検索みたいなのを作ったんですが、なんかもっと高機能にしたいと思った。正規表現とか使いたいじゃないですか!!そんなわけでPerlとかでやろうかなあと考えてたんですが、「データ構造がこういう感じになってたら、mapplyをかませ…