sheephead 200808

旅立つ人

　久しぶりにかまちょさんから電話きたと思ったら、なんとだいちゃんがPDで外国へ行くそうな。マヂっすか！！
　ん～なんかQ大もどんどん知った人がいなくなって寂しい限りだ・・・

　てなわけで、今週の土曜はだいちゃんの壮行会でまたしても福岡に行ってきます。

　先週も逝ったばっかりなんだが。しかも、かぶってるメンツもあるんだが。ま、そんなことはど～でもよいので、ほんとだいちゃんにはこれからも頑張ってもらいたいです。

2008年08月18日 | 日記 | トラックバック：0 | コメント：0

ふつーの人に多次元空間を認識させる方法を教えてください

今のお仕事で、一般の人に自分のやっていること、やりたいことを説明する機会とゆーのが、まあ多々あるのですけど一番苦労するのが、多次元母数な対象をこうやって解析します、みたいな話をするとき。

　ベイズはおろか、重回帰とか言っても会話が成り立たないので、そこは平易な言葉で書くんですけど、そもそも多次元空間を認識できない人間の脳みそにあって、そんな話してもあまり理解してもらえない（本人も、当然理解できてないが）。

　てーわけで、代替として取る戦略としては、そういうややこしいとこをブラックボックスにして、こんなことしたらこんなことができますよーな書き方をしたりとかします。

　ただ、これにもリスクがあって、あまりすっ飛ばしすぎると、なんだこの怪しい話は？詐欺じゃなかろーか？、といぶかしがられたりするので、そこの塩梅が非常に難しい。ま、十分怪しすぎることやってますけどね。

　のようなことを悩みながら、今日もデスクワークに励んでいます。なんか、最近こんなことばっかしやね。

2008年08月15日 | 日記 | トラックバック：0 | コメント：0

最近のCDジャンキー

ここんんところ、CDお買い物レビューをおさぼりしておりましたが、ま、チョコチョコ買っておりまして、久しぶりにレビューです。

まずは、こないだ聴いて感動した大橋トリオの新譜。今回も非常に良いです。今回の方がジャンルの幅が広いかも。カントリーっぽい曲調から、ジャジーな感じの曲までオールラウンド。今回もお薦め度満点の一品です。

THIS IS MUSIC
大橋トリオ

インディーズ・メーカー 2008-07-02
売り上げランキング : 1424

Amazonで詳しく見る by G-Tools

次はジャジスポこと、jazzy sportのコンピシリーズ2作目。1作目も良かったんですけど、今回は看板のgagleなんかに加え、raythoughtなんかも加わり、豪華な感じです。

Pound for Pound 2
GAGLE

PLANETGROOVE 2008-07-23
売り上げランキング : 43880

Amazonで詳しく見る by G-Tools

little tempoの最新盤。今回はアッパーな曲が多い気が。スチールパンが心地よく鳴り響きます。little tempo聴くと、夏だな～って感じが漂うですな。

山と海
	リトル・テンポおすすめ平均これは中毒。 Amazonで詳しく見る by G-Tools

最後は、ino hidefumiの最新盤。今回もフェンダーローズサウンドが渋さを醸し出します。ブランデー片手に飲めば気分は大人。

Living Message
	INO hidefumi おすすめ平均いつも聴いていたい音楽聴かせる音楽キタっ！難しいことはわかりませんが…。素晴らしい音楽。 Amazonで詳しく見る by G-Tools

　今回も、やっぱり大橋トリオにやられました。別にマニアックな曲調でもなくて、ホントに良い曲です。これは、万人におすすめできるミュージック　イン　ダ　ハウス！！
　もっと、がつんと売れても良いと思うんだが・・・

2008年08月14日 | 日記 | トラックバック：0 | コメント：2

marryとかに関する何かのこもごも

　最近、何の前触れもないおめでたいニュースがありまして、まずは大学院の時の同級生k-taro氏が11月に結婚するということと
　個人的に非常にショックだったのは、マイアイドル一十三十一タンまでもが結婚するということorz

　実は、今回結婚する二人は同い年だったりするわけですな。む～何かといろいろとあるお年頃ということなんでしょうな。

　てか、そう思うようにしよう。

　追伸
　最近、自分の腹でもない写真を載せたり、恐ろしげなニシキヘビの写真を載せたりして、なんだこれは的なお叱りうけております。これは、ただ思いつきでGIGAZINE風にflickrから適当な写真を載せてるだけです。

　てことで、私個人とは、全く関係ありませんので悪しからず

2008年08月13日 | 日記 | トラックバック：0 | コメント：0

trapの呪いから抜けられない

どうにか、zipモデルをwinbugsコードで書くことはできたのですが、trapから抜けられそうな気配が全くありません。
　ちなみにcongdon先生バージョンzipモデルはこういう風に書いてみた。

model {
# Likelihood
# Loop for Plot
for (i in 1:n) {

log(output[i]) <- alpha + beta*elv[i] + w[i]
deer[i] ~ dpois(mu.i[i])
mu.i[i] <- m[i,index[i]]
m[i,1] <- 0
m[i,2] <- output[i]
index[i] ~ dcat(theta[i,1:2])
theta[i,1] <- p[i]
theta[i,2] <- 1 - p[i]
logit(p[i]) <- betaz *elv[i]
mu[i] <- 0
}

w[1:n] ~ spatial.exp(mu[],x[],y[],tau,phi,1)

# Parameters and hyper parameters

alpha ~ dnorm(0,0.1)
beta ~ dnorm(0,0.1)
betaz ~ dnorm(0,0.1)
tau ~ dgamma(0.001, 0.1)
phi ~ dunif(0.05, 20)

sigma <- 1/tau

}

kuboさんバージョンzipモデルはこーゆー感じ。

model {
# Likelihood
# Loop for Plot
for (i in 1:n) {
deer[i] ~ dpois(out[i])
out[i] <- (1 - zero[i]) *p.n[i]
log(p.n[i]) <- alpha
zero[i] ~ dbern(p.z[i])
#logit(p.z[i]) <- betaz * elv[i] + w[i]
logit(p.z[i]) <- betaz + w[i]
mu[i] <- 0
}

w[1:n] ~ spatial.exp(mu[],x[],y[],tau,phi,1)

# Parameters and hyper parameters

alpha ~ dnorm(0,0.1)
beta ~ dnorm(0,0.1)
betaz ~ dnorm(0,0.1)
tau ~ dgamma(0.001, 0.1)
phi ~ dunif(0.05, 20)

sigma <- 1/tau

}

　まずは、変数をscale()したら怒られたので、ベイジアンクリギングポアソンモデルでやったときのパラメータ値を初期値にぶっ込んでもみたが駄目。
　次には暗黒面の奥義、事前分布操作にまで手を出したもののこれも駄目。
　次に、変数減らして定数項のみ、つまりモデル自体何の意味もなく、ただMCMC回したいだけ、とゆー本末転倒なコードを流すと２０分ほど回ってまたしてもtrap・・・・
　
　あ゛～～～ぐあ～～～っとＰＣぶっ壊したい衝動に駆られながら、気持ちを落ち着けて冷静に考えてみても次なる手が浮かばない。
　まぢ困った。
　

2008年08月10日 | 日記 | トラックバック：0 | コメント：0

花火見てきた

熊本城で花火大会があったので、逝って来ました。花火見に行くのなんて４年ぶりくらい？
　会場はそんなにごった返しているというわけでもなく、リラックスして見れたです。

　当然、カメラも持って行ったんですが、む～イマイチだ。

　花火写真を撮るときは
　　１．三脚を持って行く手間を惜しむな
　　２．場所命。柔軟に場所を変えられる行動性大事

　ってことか。

http://art9.photozou.jp/bin/photo/11608887/org.bin?size=800

http://art10.photozou.jp/pub/203/61203/photo/11608857_big.jpg

2008年08月08日 | 日記 | トラックバック：0 | コメント：2

Rプログラミングお仕事

以前からお願いされていたRプログラミングをちらっと書いてみる。

　中身は、nls使ったフィッティングで解析自体はややこしくないんですけど、マスターデータから各試験体のデータ部分を個別に抽出しないといけないところで、ちと悩む。

　結局は、各試験体データの間に空白行があったので、いちおそれをキーにして分けるってことでどうにか完成しますた。csvファイルを読み込むときにこうすれば空白行を認識するみたい。

　　read.csv("hogehoge.csv",blank.lines.skip=FALSE)

　read.csv()にそんな引数があったんですな。また一つ賢くなった。ただ、空白セルなのかNAなのかとか諸々の問題で、このあたりはバギーなので、あまりおすすめはできんような気が・・・・

　もちっといぢるか。
　

2008年08月05日 | 日記 | トラックバック：0 | コメント：0

ubuntuでR環境を整える

zipモデルは折れ折れなので、気分かえるためubuntuをいろいろいぢってみました。

　まずは、Rjpwikiを参考にRをapt-get。emacsはすでに入れてるので、ついでにESSも入れてサクッとR実行環境確保。むお～便利すぐる。

　いろいろとパッケージを眺めてみるとRkwardなる便利そうなR実行環境を発見。なんですけど、KDEでしか動かないのか、うまくｲﾝｽｺできんかったです。今、GNOMEだしな。頃合い見てKDE入れた後に試してみよう。

　ま、そんだけ。うん、明日から本気出す。

2008年08月02日 | 日記 | トラックバック：0 | コメント：0

congdon先生バージョンzipモデル

今月の中旬にODさんが来福されるってんで、その時にベイズ談義でも、な話になったのですが、雑務にかまけてシカベイズを放置しっぱなし。

　というわけで、久々に昔自分がやろうとしてたことを復習。
　まずは、ベイジアンクリギングまでは何となくうまくいったんだけど、zero-inflatedなモデルを作らなければいけない、とまではわかった・・・はず。

　と、そこまで来て、brugs使ったMCMCループとか脱線しまくってしまって、収集が付かなくなったんだ。

　うん、まずしないといけないこと整理しとく。

　・zipモデル仕上げる。
　・spatial.pred使ってお絵かきしてみる。
　・線形予測子をいろいろ代えてみる。
　・マップをいろいろと重ねてみる。

　今日は、とりあえずcongdon先生的zipモデルまでは完成させなければ、見通しがたたんな。

　んでまあ、bugsのコードは書けたことは書けたんですが、またしてもtrapにはまる・・・
　ぐぅ、とりあえず今日は頭使いすぎでtrap窓と格闘する気力も残っていないため、ここまでにするか。
　

2008年08月01日 | 日記 | トラックバック：0 | コメント：0

follow myuhe at http://twitter.com

CalendArchive

≪ ≫
-	-	-	-	-	1	2
3	4	5	6	7	8	9
10	11	12	13	14	15	16
17	18	19	20	21	22	23
24	25	26	27	28	29	30
31	-	-	-	-	-	-

最近はこんなの聞いてます

今週のベスト５

Tree-Comment

プロフィール

myuhe

Author:myuhe

物欲の神様
衝動買いの伝道師

mixiは↓
　http://mixi.jp/show_friend.pl?id=1864001
　　
last.fmは
　http://www.last.fm/user/myuhe/
　

最近の晩ご飯