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

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

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

リモートデスクトップ

ubuntu からubuntu へどちらも12.04のとき。接続される側でDashホームで、「デスクトップの共有」を起動する。 (ウィンドウがMac的な感じだが) というように設定しておく。vncserverを入れておく 接続される側のアドレスをifconfigで調べておく。 inetアドレ…

Ubuntu でログイン出来なくなったときの解決法

今日というか昨日の夜、院生室で作業をしようと思って、いつも通りUbuntuを起動しようと思ったが、パスワードを入力したら、一瞬(肉眼では確認できないくらい)、文字列 ナンチャラカンチャラ〜 ...... checking battery state ... が表示されて、ログイン画…

独立成分分析について

詳解 独立成分分析―信号解析の新しい世界作者: アーポビバリネン,エルキオヤ,ユハカルーネン,Aapo Hyv¨arinen,Erkki Oja,Juha Karhunen,根本幾,川勝真喜出版社/メーカー: 東京電機大学出版局発売日: 2005/02メディア: 単行本購入: 1人 クリック: 16回この商…

髪切りたい

おはようございます。引越しをして以来、ダンボールに何が入っているのか分からないので、バリカンがどこにあるのかわからない。この前は、冷凍のポテトが食べたくなったのでフライパンを探したのだが、いくら探してもフライパンが見つからなかったので、新…

Reminder Beta版を開始しました。

Reminderは、(http://reminder.asia) ブックマークをしておいても、結局読まないという問題を解決するために作られたサービスです。はてなブックマークで「あとで読む」タグをつけたブックマークが溜まっていませんか? そんな、あなたにReminderはブックマ…

XML::Parserが入らないとき。

XML::Parserが入らないとき。 ログは、Expat.xs:12:19: error: expat.h: No such file or directory。 sudo apt-get install libexpat1-dev で入る。

$useridのブックマークのタグを取得。

#!/usr/bin/perl use strict; use warnings; use HatenaBookmark::Number; #my director use XML::FeedPP; use Data::Dumper; use Array::Uniq; my $userid = "infinity_th4"; my @user = ($userid); sub gettag { my (@user) = @_; my @tags = (); my $alln…

$useridのブックマーク数を取得。

$source = "http://b.hatena.ne.jp/".$userid."/rss"; my $feed = XML::FeedPP->new( $source ); $num = $feed->get('opensearch:totalResults');

それでも僕は、花粉症だと信じていなかった

僕は、1年前、花粉症(杉)だと診断されたが、それを信じていなかった。なぜなら、ほぼ高校3年から365日鼻水が止まらない。1日何度か波のようなものがあり、鼻がつーんをすると、それが合図。鼻炎だと思う。 そして最近になって目がとてもかゆくなった(ここ3日…

screen

screen を使って、sshを切っても処理を実行させたいので、 $ screen $ script/myapp_server.pl -rという感じにする。僕の場合は、script/myapp_server.pl以外に別の.plを走らせる必要があるので、新しくCtrl+Shift+tとして、ssh接続し、 $ screen $ sample.p…

今日という日

強風だったのでTXのダイヤが乱れていた。東京に行くときはやはり、時間に余裕を持っていかなければ、急な出来事に対応できないと思った。千代田線も少し遅れていたようで、不安だった。時間通りにことが進まないと不安になる。なので、余裕のある生活を送り…

主成分分析による手書き数字の復元

手書き数字について、主成分分析を行うことで、特徴抽出、画像の文字の復元することができる。https://www.kaggle.com/c/digit-recognizer/data のtrain.csvデータにheaderにlabel, pixel0 から pixel783まで定義されている。train.csvを用いて、手書き数字…

デジタル画像(Handwritten)

http://bayesianbiologist.com/2012/08/13/the-essence-of-a-handwritten-digit/で使われている https://www.kaggle.com/c/digit-recognizer/data のtrain.csvデータにheaderにlabel, pixel0 から pixel783まで定義されている。 28*28(=784)の手書き画像が42…

Integrate[Exp[-(s^2/2)]*Cosh[s], {s, -Infinity, Infinity}] は、 s=seq(-5,5,,1000) plot(s,exp(-log(2)-2*log(cosh(s))),xlim=c(-5,5),ylim=c(0,0.8),type="l") par(new=T) plot(s,(1/sqrt(2))*exp(-sqrt(2)*abs(s)),xlim=c(-5,5),ylim=c(0,0.8),type="l…

library(RImageBook) #biOpsの場合 f <- readImage(system.file("images/lena2.tif", package="RImageBook")) lenab = readTiff(f) とやったら、 TIFFOpen: 0.886274509803922: Cannot open. 以下にエラー readTiff(f) : Cannot open file. とでて読み込めな…

最尤法による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…