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

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

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

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で変換して表示する。

seminer

今日のゼミは3ページしか進まなかった。 でも、ベイズ推定のところで内容は重要。次回は8月31日14:00から再開。

you need to install the perl-doc package to use this programと言われたので、sudo apt-get install perl-docとすれば、インストールできる。$ perldoc -ml File::Find などとすれば、 /usr/share/perl/5.14/File/Find.pm で、CPAN モジュールの .pm ファ…

射影法

昨日やったゼミのメモで射影法は、 制約条件がでJ(w)を最小にするようなwを求めることに使えるのです。 射影法は、m次元単位球面上に直交射影するというものだが、勘違いをしており、更新前のが更新後のに直交するということと思っていました。しかし、これ…

「はてなサマーインターン2012」に参加させて頂きます!

「はてなサマーインターン2012」にお選びいただいてありがとうございます。「Webアプリケーション開発 基礎コース(2週間)」に参加いたします。昨日連絡を頂いたわけですが、嬉しかったです。本当は、実践コースに参加したかったわけですが、筑波大学は9月か…

IMS-APRM2012

IMS-APRM2012(http://ims-aprm2012.org/)に参加(?)してきた。 スタッフとしてですが。 僕の担当したところは、回帰の話が多かったですかね。他にはカーネル法の話も出てきて、良かったですね。

2次元版、ニュートン法(関数の勾配が0である解を求める)

#loss function1 x=seq(-2,2,length=50) y=seq(-2,2,length=50) lossfunc1=function(x,y){ return(x*y*exp(-x^2-y^2)) } z=matrix(0,50,50) for(i in 1:50){ for(j in 1:50){ z[i,j]=lossfunc1(x[i],y[j]) } } persp(x,y,z,ticktype="detailed",theta = 30, …

1次元版、ニュートン法(関数の勾配が0である解を求める)

func1=function(w){ return(w^2-2) } dffunc1=function(w){ return(2*w) } ddffunc1=function(w){ return(2) } newton102=function(w,times){ cnt=1 for(i in 1:times){ w_before=w w=w-dffunc1(w)/ddffunc1(w) printf("cnt:%d;w=%.8f,error=%.8f\n",cnt,w,a…

2次元版、勾配降下法2(最急降下法)

極値が1個だけだと面白く無いので、次のような関数を損失関数とします。 #loss function1 x=seq(-2,2,length=50) y=seq(-2,2,length=50) lossfunc1=function(x,y){ return(x*y*exp(-x^2-y^2)) } z=matrix(0,50,50) for(i in 1:50){ for(j in 1:50){ z[i,j]=l…

最近、

最近というかパソコンを新しく買ったぐらいから、Google Driveにコードを置くようになった。家や院生室や研究室でも使うので、かなり便利である。それだけ。

書くのが遅くなりましが、PRML 復々習レーンに参加してきた。

渋谷にあるDeNAさんでPRML勉強会があったので、参加してきました。(渋谷には初めて行きました)第一回は時間的余裕がなかったためキャンセルさせて頂きました。僕は3.5(エビデンス近似)と3.5.1(エビデンス関数の評価)の担当になりました。うまくできるかなぁ…

Django

Django×Python (LLフレームワークBOOKS)作者: 露木誠出版社/メーカー: 技術評論社発売日: 2009/02/16メディア: 単行本(ソフトカバー)購入: 10人 クリック: 120回この商品を含むブログ (33件) を見るchapter3のhttps://www.djangoproject.com/download/pyth…

2次元版、勾配降下法1(最急降下法)

#loss function x=seq(-3,3,length=50) y=seq(-3,3,length=50) lossfunc1=function(x,y){ return(x^2+y^2) } z=matrix(0,50,50) for(i in 1:50){ for(j in 1:50){ z[i,j]=pdf(x[i],y[j]) } } persp(x,y,z,theta = 0, phi = 20, expand = 0.5,col = "lightblu…

勾配降下法1(最急降下法)

勾配降下法をやっています。 問題設定:損失関数J(w)を最小にするようなwを求めることが目的です。 今、としておきます。Jのグラフは、のようになります。のとき最小値は0です。勾配はです。グラフは、 です。 func2=function(t){ return(t^4) } dffunc2=fun…

メモ(勾配降下法(最急降下法))

func2=function(t){ return(t^4) } dffunc2=function(t){ return(4*t^3) } grad2=function(x,rho,its){ cnt=1 for(i in 1:its){ xbefore=x x=x-rho*dffunc2(x) if(0==cnt%%1000000){ printf("cnt:%d;x=%f,abs(x-xbefore)=%.13f\n",cnt,x,abs(x-xbefore)) } i…

勾配降下法3(最急降下法)

2つ解がある場合。 x=seq(-10,10,,10000) func4=function(t){ return((t+3)^2*(t-3)^2) } plot(x,func4(x),xlim=c(-10,10),ylim=c(0,100),type="l") x=seq(-10,10,,10000) dffunc4=function(t){ return(2*(t+3)*(t-3)^2+(t+3)^2*(t-3)) } plot(x,dffunc4(x),…