Kosugitti's BLOG

アンドロイドは正規分布の夢を見るか

この画面は、簡易表示です

2015 / 12月

来年の抱負

来年は痩せますよ。宣言通り、人が変わったように、ヘルスサイドに堕ちますよ。

ということで、年末の友達との集いでその話をし、どうすれば痩せられるかについて、いろいろ協議しました。

結果、次の三か条を守る、という約束をしてきました。

  1. ラーメン屋でレンゲを使わない
  2. 夕食時にお酒を飲まない
  3. 夕食後のオヤツは無し

1つ目のやつは、スープを飲むなという意味です。丼鉢つかんで飲むのは論外。麺は食べてもスープを飲まないこと。

2つ目のやつは、食事は食事、お酒はお酒、と分けましょうとのこと。

で、最後のはそのお酒の時に、スナック菓子やピーナッツを食べないということ。

この条件を習慣化したら、4.5キロは痩せられるという。ウソーン、そんなに痩せへんよ、ほなやってみぃ、ということで来年からそうすることにしました!

ちなみに、帰りに食べたラーメン納めのらーめんで、レンゲを使わないトライをしたら、意外といけた。ふむぅ。



入部しました

京都ラーメン部に入部。早速年末の締めのラーメンを食べてきました。クリーミィな方向に言ってるというけど,いやいや,京都ラーメンはいつ食べてもうまし。

 

そのあと,いろいろお買い物,の一環で,ついに万年筆を買ってしまいました。去年の秋口からずっと「いいのが欲しいなー」と思っていて,ゆっくりいろいろ悩んだつもりだけど,結局はテレビで見たものをまんま買いましたな。これで万年筆部にも入部できたかな。

ペリカンのスーベレーンM600です。最初は400かな,と思ったけど書き比べて見てこっちにした。サイズ感で決まりかな。

いいものを買ったなー,とウキウキしています。まるで恋だね。



2015年まとめ

まいど!年末ですね!

25日にお歳暮交換会をして,一足先に仕事納めをしてしまいました。帰省してのんびりしております。

帰省の道中で「今年の俺の10大ニュースってなにかなー」と思って考えていたんですが,なかなか難しい。一つは大人になるにつれて「口外できないこと」みたいな案件がいくつかでてきて,それをこんなところで書けないもんねー,とおもっていたら,去年の俺も書いてなかったw

まあでも,無理してでも出すことに意味があると思ったんです。全部「いろいろあるわな」と言ってしまったら,ますます記憶が薄くなるからね。ということで,頑張って思い出して総括しますよ。

  1. 家のこと
  2. 師匠の退職記念パーティー
  3. イタリア行ってきたよ
  4. 欠番
  5. 丸坊主にしてみました
  6. 欠番
  7. 冷蔵庫が入らないなんて
  8. 新入生合宿研修
  9. ゼミ合宿と学会発表と病院
  10. 40前,さらに仕事が

 

一位の「家のこと」なんですけど,本当はこれだけでランキングをほとんど埋めてしまうんじゃないかしら,ってぐらい。というのも,ほぼ一年がかりでやっていたことだから。土地を探し,うまくいかず,中古マンションを見に行ったりしたものの,結局(まだふみもみず)土地を買うことになり,家を見て回ってローン組んで餅巻いて・・・。一生もののことだけあって,結婚と同じぐらいエネルギーを使う感じですね。例えば7位のことなんかも関係するし,それより何より,本当はいつも一位にしておきたい「家内安全」が今年も叶ったということ。子供たちがすくすくと育ってくれたおかげで,今の狭い賃貸じゃあどうしようもなくなって,うれしい悲鳴というところ。40の前厄ってこともミソで,これで厄落としになる・・・らしいぜ。イベントとしてひとくくりにしてしまっていいのかな。

二位はもうねー。まず,2年ぐらい前から準備しててね,本当にありがたいことにまた気のいい仲間たちとお仕事ができた,イベントができたということが嬉しかったです。イベントの内容的には悲しいことでもあるはずなんだけど,きちんとできた,とか,ちゃんと盛り上がった,というか,いろいろうまくいったことがうれしかったです。よかった。こんなことはもう人生で2度とないでしょう。

第三位は,母とイタリアに行ってきたことです。学会ですけどね。食べ物がなんでも美味しいし,街ゆくイタリア人セニョールはみんなかっこいいし,最後の晩餐も見たりちょっと足を伸ばしたり・・・まあ思いでいっぱいでしたよ。母に仕事している姿も見せられたしねw

第四位,第六位はこれまた守秘義務に関わることでもあるのでー,言えないのでー。

五位に入れたのは珍しく外見に関することです。死ぬまでに一度やってみたかったこと,まわりからやるなやるなと言われていたこと,でも不惑になるまでに惑ってみたかったこと,で,6ミリで丸刈りにするというヤツを体験してみました。頭の触り心地が面白い,周りの反応が面白い,徐々に伸びていく雰囲気が楽しい,いやもういろいろ楽しかったです。個人的にはもう一回ぐらいやりたいんだけど,だめかなー。まあやりたくなったらやります。ただ,顔も含めてあだ名が「丸い」って言われたのはショックだったなw

七位の冷蔵庫云々は,買った冷蔵庫が今の家に運び込まれなかった問題です。「家が狭くて返品して悪かった,だから家を大きくしたよ」と言いながら,秋口にまた買いに行ったんですけどね!このネタができてよかったw

八位,九位は仕事のことですねー。八位の件は,一年生に関することかな。今年から改組されて,一年生チューターにもなったので,いろいろ仕掛けをしてみたんです。個人的にはちゃんとやれたし楽しかったし,ゼミ生が頑張っt稀有れたのも嬉しいですね。九位の件はゼミのことですが,ゼミ合宿と学会発表を一緒にやるってのは効率もいいし面白いな,ってこと。ちょっとしたハプニングもありましたけどね!

 

最後に十位。「30代に来た仕事は断らない」というポリシーを10年間やってきましたが,来年からはちょっとちゃんと考えるぞ!とか思っていたら,駆け込み需要かよ!という感じで大物の仕事が飛び込んできました。偉そうに言いましたが,基本的に仕事をもらえることは嬉しいことです。選ぶつもりも別にありません。誰にも求められなくなったら,この商売も終わりですからね。他人に「あいつ変なことやってんなぁ」「あいつの話を聞いてみたいなあ」と思われるような人間で居られるように,これからも精進していくつもりです。ただ,来年からはそろそろ「仕舞の大仕事」も考えていきたいなってこと。あと20年しか仕事ができにとすると,10年計画で5年でやって5年で仕上げる,ぐらいのつもりがいるからね。ちゃんと研究者人生を全うしてから死にたいと思っています。

でもよくよく考えると,結局この一年は一言で言うと,年齢のことを気にし続けた一年だったと思う。不惑前に,30代最後に,今できること,これからしたいこと,ずーっとそういうことを考えていたような気がする。記録には残りにくいけど,体が徐々に硬くなっていったり動かなくなっていったり,高尿酸値で健康診断で「要検査」が出たり,といろいろ気になることがあったわけで。多分来年来る「40代」ってのをかなり意識しているんだろうな。というのも,俺は昔から兄貴やオヤジのような大人になりたい,と思ってきたからで。

家のローンを組んだ時に,銀行が生命保険を俺にかけるという話で,「ああ世の中そういう仕組みだったのね」ということを思うと同時に,「あれ,これで死ねるな」とも思ったのは事実。もちろん研究者としてはまだ死ねないよ。私人としても,妻や子供を残して死にたいわけではない。ただ,社会的に死ぬことの許可が得られたような,そんな変な感じがした。それは決して不快感ではなく,安心感のような。

まだまだ大人になれないなー,という感じを抱えつつ,今年を収めて,来年からもまだまだやっていきたいと思います。



因子軸の回転の可視化

「因子軸を回転って結局どういうイメージですか」

という質問をされたのですね。
探索的因子分析の後は,解釈をしやすくするために軸を回転しますよ,と説明しているのだけど,そのイメージがわかりにくいということなので,Rでプロットしてみることにしました。

psychパッケージのデータセットbfiを例にあげます。これは五次元のデータなんですけど,仮に2因子だとしましょう。
(後で出てきますが3次元空間の回転行列から角度を求めるのは行列計算的にちょっと面倒なのです。)

因子数を2にし,回転なし,バリマックス回転,プロマックス回転の結果をオブジェクトに代入しておきます。

[crayon-5bcacf4441cc0402412364/]

[crayon-5bcacf4441cd1148452430/]

[crayon-5bcacf4441cd6473950484/]

直交回転の場合

さてここから加工。まず因子負荷行列をmatrix型で抜き出し,データフレーム型に変換。psychパッケージのfaオブジェクトはいろいろ付いているからね。無駄なものをそぎ落とします。

[crayon-5bcacf4441cdb961853083/]

因子名をわかりやすくFx,Fyにしましょう(2次元だから)

[crayon-5bcacf4441ce1703938520/]

回転なしの因子負荷量をつかって因子負荷行列のプロットを描いてみます。

[crayon-5bcacf4441ce6142019603/]

 

unnamed-chunk-4-1

 

回転後の因子負荷行列の散布図も描いてみます。

[crayon-5bcacf4441ceb822955254/]

unnamed-chunk-5-1

うむ,確かにこれでは分かりにくい。

さて,回転は行列をかけることで行っています。座標を変換するんですね。Wikiのページにもあるように,回転行列とはcos,sinをつかった行列です。つ回転行列

この回転行列はrot.matの中に入っています。

[crayon-5bcacf4441cf0745517410/]

1行1列目は$\cos \theta$の値でしたね。この値から逆に角度を求めます。逆cosはacosです。

[crayon-5bcacf4441cf5437967493/]

[crayon-5bcacf4441cfa835606901/]

0.71ラジアンぐらいですね

図にしましょう。まず回転前のプロット。

[crayon-5bcacf4441cff662886442/]

回転軸を引きます。角度を傾きにするため,tanをとっています

[crayon-5bcacf4441d04824419418/]

第2次元は第1次元と直角に交わるので,0.5piを足してtanをとっています。

[crayon-5bcacf4441d08910099689/]

unnamed-chunk-10-1

たしかに,元の青い軸のプロットよりも赤い軸のプロットの方が,点と軸の距離が近いようです。強調した,というのがちょっとは実感できるでしょうか。

斜交回転

promax回転の場合。最初の処理は同じです。

[crayon-5bcacf4441d0e323595690/]

回転行列を抜き出します。

[crayon-5bcacf4441d12795216415/]

角度を抜き出します。

[crayon-5bcacf4441d17387138909/]

さて,斜交回転の場合は軸が相関します。相関行列はPhiに入っています。相関係数は(\cos \theta) の結果なので,逆算(acos)して第一軸と第二軸の角度を求めます。

[crayon-5bcacf4441d1c231777410/]

これでプロット。

[crayon-5bcacf4441d21483873662/]

unnamed-chunk-15-1

こんなもんでどうでしょう。

本当のこというと,適当な2軸を選んで描画し,軸をぐいっと動かして相関係数を見せるようなインタラクティブなR画面を作りたかったけど,それは技量と時間の限界ということで,誰か助けてエロい人。




top