INFINITY -数学とかプログラミングとか-

統計とプログラムを使って役に立たせたい

TeX用コマンド入力を支援するための辞書をご利用ください。
sanctuary's blogは,適当なことが書いてあります。

2012-01-01から1年間の記事一覧

最尤法によるICA(メモ)

https://docs.google.com/open?id=0Bzz10UZbU-i5WGctWWtQSERsaTg https://docs.google.com/open?id=0Bzz10UZbU-i5aUNfU0UwWVpXcjQメモ: 勾配法を用いた最尤法による独立成分抽出は,お世辞にも推定できているとはいえない。 最尤法によるFastICAは,まあまあ…

久しぶりに落語を聞きに行きたい

湯屋番柳家三三 立川談志紺屋高尾FullVersion 八代目桂文楽- 明烏明烏が一番好き。

Schedule::Cron

foreachを使ってadd_entry すると指定した時間でちゃんと動かない。(即、jobが実行される) #!/use/bin/perl use strict; use warnings; use Schedule::Cron; sub dispatcher{ print "ID: ", shift, "\n"; print "Args:","@_", "\n"; } sub job{ # 実行した…

デジタル画像処理

デジタル画像処理 (Rで学ぶデータサイエンス 11)作者: 勝木健雄,蓬来祐一郎,金明哲出版社/メーカー: 共立出版発売日: 2011/11/23メディア: 単行本購入: 2人 クリック: 10回この商品を含むブログ (4件) を見る最初、自分でネットで調べてGTK+やImageMagicのイ…

疑問

?

任意の分散の場合は、一様分布が、エントロピー最大の分布である。しかし、分散一定、例えば1のとき、正規分布がエントロピー最大の分布である。一見、図から一様分布のほうが、乱雑さが高く、ある点に密度が集中していず、正規分布のほうが、平均周りで集中…

確率的最急勾配法

R

1つのwのみ。 #確率的最急勾配法 f<-file("rnorm_data.txt","r") data=function(){ a<-readLines(con=f,1) b=unlist(strsplit(a, "\\,")) # 文字 "," で分割 return(as.numeric(b)) } #data()#ファイルの1行目が読み込まれる #data()#ファイルの2行目が読み…

rnorm()によるデータ作成

R

getwd() setwd("作業ディレクトリ") N=1000 out <- file("rnorm_data.txt", "w") # 書き込みモードで開く for(cnt in 1:N){ x=rnorm(5) for(i in 1:length(x)){ if (i < 5) writeLines(paste(x[i]),out,sep=", ") else if(i==5) writeLines(paste(x[i]),out,…

双方な分布のネゲントロピーに関する考察

分布を平均0,分散1にしたときのネゲントロピーを計算し,横軸が0.001刻みでを0から1まで動かし,縦軸がネゲントロピーの近似した図 (ちゃんとやったのに図が間違っているかも?) ネゲントロピーが大きいという事はガウス分布から形状が離れているという事であ…

一般化ガウス族のネゲントロピーに関する考察

#一般化ガウス族(GGF) p=function(xi){ (1/(2*a^(1/a-1)*gamma(1/a)))*exp(-abs(xi)^a/a) } a=2 mean=(a^(1/a))*(gamma(2/a)/gamma(1/a)) var=(a^(2/a))*(gamma(3/a)/gamma(1/a))-((a^(1/a))*(gamma(2/a)/gamma(1/a)))^2 Gauss = function(x){ (1/sqrt(2*pi*…

catalyst

my $row = $c->model('RemaindDB::Contents0')->create({ content => $content }); とかやったときに、データベースにInsertされない。不思議。 原因がよくわからないが、下が気になる。 $ script/remaind_server.pl -r DBIx::Class::Schema::connection(): …

株式会社はてなに、出社しました(*´∀`*)

id:infinity_th4です。念願のはてなインターンシップを、経験してきました。すでに他のインターン生が、ブログでインターンの概要について書いていらっしゃいます。 こちらをご覧頂いたほうが、秀逸な文章が揃っていると思います。 はてなインターン前半まと…

エントロピー近似

が、0になる。 不等式評価とかしなくてよくて、単純にぶち込んで計算する。は、ガウス分布。Hなんちゃらとかは、エルミート多項式。 とする。 少々、面倒な計算を経て, =k_4(\frac{648}{864}k_3^2+\frac{1728}{13824}k_4^2)"/>を得ることができる。

メニューバー

ubuntu で、メニューバーに電源ボタンが表示されなくなったと思ったら、何回か再起動してたら表示された。

glibは2.32.1 pkg-config 0.20はインストール済み。 (http://www.bnote.net/linux/pkg_config.shtml) http://hooktail.org/computer/index.php?GTK%2B-2.0%20%A4%CE%A5%A4%A5%F3%A5%B9%A5%C8%A1%BC%A5%EBを参考にatk-1.0.3をインストールしたいのだが…

R

ubuntu ソフトフェアセンターでr-base-coreをインストール。 r-base-core: Depends: libgfortran3 (>= 4.3) しかし、4.6.3-1ubuntu5 はインストー\ ルされます Depends: libgomp1 (>= 4.2.1) しかし、4.6.3-1ubuntu5 はインストール\ されます Depends: libj…

Emacs

とりあえず、これだけ覚えてれば何とかなる。Cは、Ctlrのこと。C-x h 全て選択Alt-w コピーC-w 切り取りC-y 貼り付けC-s 順方向に検索C-r 逆方向に検索C-x C-s 上書き保存C-x C-c 終了C-x 2 選択されているウィンドウを上下2つに分割する。C-x 3 選択されて…

Gitってやつは

git

sudo apt-get install git-core とし、git からgit-1.7.12をDownloadしてくる $ tar zxvf git-1.7.12.tar.gz $ cd git-1.7.12 $ ./configure $ make prefix=/usr/local all $ sudo make prefix=/usr/local install xxx@xxx:~/Downloads/git-1.7.12$ git conf…

MacBook

Mac

かなり前のMacbookを持っているが、(中古で買った) DVDを入れたら、反応せず、しかも出てこなくなった。http://support.apple.com/kb/HT2286?viewlocale=ja_JP&locale=ja_JPこの辺に書いてあることをひと通りやった。しかし出てこないので、修理に出すとい…

UbuntuでEmacs

UbuntuでEmacsを使うようになったが、Macでは、terminal 上でemacs が起動するのに対し、Ubuntu では、Emacs の別ウィンドウが開きなんだか、不便に感じた。 調べると、 $alias emacs='emacs -nw'とすればよい。

Django×Python

ゼミが終わってから、Django×Python (LLフレームワークBOOKS)作者: 露木誠出版社/メーカー: 技術評論社発売日: 2009/02/16メディア: 単行本(ソフトカバー)購入: 10人 クリック: 120回この商品を含むブログ (33件) を見るDjango×Python (LLフレームワークBO…

GD::Graph的なのをcpanm で入れた時にエラーが起きたら、lib~~~dev的なものを入れれば良い。

Pod::Usageの使い方

DBの取り扱い

mysql -uroot -e "Drop database infinity_th4" mysqladmin -uroot create infinity_th4 mysql -uroot infinity_th4 < db/schema.sql mysql -uroot -e "show tables" infinity_th4 mysql -uroot -e "select * from entry" infinity_th4 mysql -uroot -e "des…

matrix

#データを行列の形にする foreach my $line (@lines){ @tmp=split /,/, $line;#$lineは一行 #warn @tmp; for($j=0;$j

R

data<-read.csv("data.csv",row.names=1) 行の名前をdata.csvの1行目に書いてある時にrow.names=1とする。 setwd("c:/usr") # 作業ディレクトリを変更する getwd() # 現在の作業ディレクトリを確認する

grep

カレントディレクトリ下のファイルで「ようこそようこ」を含むファイル find . -exec grep -H ようこそようこ {} \; カレントディレクトリ下の拡張子plファイルで「ようこそようこ」を含むファイル find . -name "*pl" -exec grep -H ようこそようこ {} \;

エントロピー

R

p=seq(0,1,length=1000) f=function(p){ -p*log(p) } plot(p,f(p),xlab="p",ylab="entropy",xlim=c(0,1),ylim=c(0,0.5),type="l") 確率変数Xが取り得る値が8で確率が、 (1/2,1/4,1/8,1/16,1/64,1/64,1/64,1/64)であるとする。この時、Entropyは、以下のよう…

Hatena Intern 2012

PRML 復々習レーン #4の発表資料

PRML #4 - 3.5,3.5.1 from Hiroshi Tashiro ちなみに、はてなダイアリーでは、iframeが使えないので、http://gimite.net/rails/iframe_gadgetで変換して表示する。