読者です 読者をやめる 読者になる 読者になる

300億円欲しい

メジャーリーグのデータ解析します

俳句で形態素解析をしてみたい

R

序論 友人から17,000の俳句を渡されました. せっかくなので形態素解析してみました. 参考文献 データマイニング入門解析の結果を利用して, 雰囲気だけそれっぽい俳句が作れそうな気がします 面白かった俳句 パンツ脱ぐ 遠くの少年 泳ぐのか 分析が明確です泳…

メジャーリーグのデータ解析をしたい (安打製造機かどうかなんて分からん)

前々回に試した手法が全然ダメそうだ, という内容です. 序論 こんな記事を書きました メジャーリーグのデータ解析をしたい (イチローは本当に安打製造機なのか) - 300億円欲しい内容をまとめると, サロゲートデータ法を使って, イチローの打席結果に系列相関…

効率的にデータフレームの処理がしたい ( doBy vs plyr vs dplyr)

R

メジャーリーグのデータを弄って遊んでいます. 100Mbオーダーのデータフレームの処理をしています. 序論 今まではplyrのddplyで計算させていましたが, 遅いです. doByパッケージで大体同じことができて, しかも速く計算出来ました. どう使い分ければいいんで…

イチローは本当に安打製造機なのか!?

この記事はR Advent Calendar 2013 : ATND , 1月3日の配当記事です. よろしくお願いします. 序論 イチロー. 凄いですよね. MLBでの成績をのせます. YEAR GAME ATBAT HIT AVG 1 2001 157 692 242 0.3497110 2 2002 157 647 208 0.3214838 3 2003 159 679 212 …

効率的にデータフレームの処理がしたい(doByパッケージとは)

序論 Rでデータフレームを効率的に弄りたいです. plyrパッケージはとても優秀です. 使いましょう. doByパッケージもとても優秀です. 使いましょう. irisデータの処理 データフレームを処理する道具として, plyrパッケージを紹介します. みんな大好きirisデー…

メジャーリーグのデータ解析をしたい (pitch f/xって何)

pitch f/x とは また野球の話です. こんな記事がありました. 【MLB】 PITCHf/xとかいう球速測定装置wwwwww - キャッチャーライナー http://blog.livedoor.jp/catcherliner/archives/7557031.html複数の場所からボールを撮影することで投球の軌道を追い…

メジャーリーグのデータ解析をしたい( 犠牲フライくらい打てよ )

序論 Rを使ってデータ解析します. またメジャーリーグの話です.0アウトか1アウトでランナーが3塁にいると, 「おっしゃ1点はとったな」 と思いますよね.外野フライでもいいし, 深めの内野ゴロでもいいです. 高く弾んでもいいです. 何か少し手間のかかるところ…

メジャーリーグのデータ解析をしたい (ダルビッシュがすごい)

序論 Rを使ってメジャーリーグのデータ解析がしたいです.全ての試合結果データが公開されています. retrosheetという名前で, 無料で使えます. http://www.retrosheet.org/game.htm せっかくなので使います. 試合結果データを整形してcsvファイルにしたあとで…

メジャーリーグのデータ解析をしたい (K/BBの計算したら上原が凄い)

序論 Rを使ってメジャーリーグのデータ処理をします. 今回は前処理に役立つplyrパッケージを使いながら, レッドソックスの上原浩治選手の成績を確認します. 今回使うのは, plyrパッケージです. データのまとめにはとても便利です. ddplyの使い方がわかると思…

メジャーリーグのデータ解析をしたい ( 21世紀限定HRランキング)

序論 ある程度まとまったデータの解析をしたいです. でも, 興味のないデータを漫然と弄っても楽しくないです. 今回は, みんなが大好きなメジャーリーグのデータで遊びます. 野球に詳しくなって, さらにRの関数の使い方も練習できます. さっそく, 21世紀のデ…

beamerで補足スライドにジャンプしたい

発表内容を補足するようなスライドを作ったとします. 本編のスライドから補足スライドまでリンクを貼って, ジャンプしたいですね. PowerPointなら簡単にできますが, 貧乏なので持ってません. 無料で使えるbeamerでも, 簡単にジャンプできます. \hyperlink コ…

beamerでスライドを再掲したい

少し前に出したスライドを再掲したい時ってありますよね. たとえば, 目次. 発表の進行状況と流れを意識させるために, 何回も再掲したいです. beamerを使っているときは, アホなのでframeの内容をコピペしてました. こんな感じです. \documentclass[10pt, dvi…

beamerでonslideコマンドを使いたい

序論 TeXでスライドを作れるbeamer. 某M社のPowerPointほど自由な操作はできません.その分, 簡潔なスライドは作りやすいと思います. beamer力を上げて表現力を高めましょう.今回は, onslideコマンドについてです. onslideコマンド こういうスライドを簡単に…

AtCoder016B 音楽ゲーム

ARC016です. http://arc016.contest.atcoder.jp/tasks/arc016_2 プロコンに興味はないです. プログラミングの勉強です. 問題文 楽器を演奏できない高橋君は、音楽ゲームが大好きです。 このゲームでは 9 個のボタンが存在し、曲に合わせてボタンをタイミング…

何を食べたらいいのかを教えてくれるbotを作りたい

序論 お腹が空いて, 何食べようかな...と考えても, 何も思い浮かばないときってありますよね.友達や家族と一緒にいるなら, 「何食べたい?」と尋ねればいいと思います.. しかし, 僕は基本的に1人でいます. どうすればいいのでしょう.こんな状況を打開したいの…

Vimでノーマルモードに入るときに英数入力へ切り替えたい

vim

Vimあるあるとして, ノーマルモードでの全角入力があげられると思います.こんな感じ. ( `Д´)フォオオオオオオオオオ! ノ#`Д´)ノ⌒┻━┻ Insertモードで全角入力のままEscapeするとこうなります. この現象の対策としては, 1. 日本語を使わない 英語だけで生きる 英語だけ使…

Rで研究者の業績を調べたい

R

Rは日々パッケージが更新されています. レポジトリには数千のパッケージがあるのでフォローしきれません. 面白そうなものは使っていきたいです.今回は "scholar"パッケージを使います.参考文献 http://www.r-bloggers.com/new-r-package-scholar/ Package "s…

Vimrcを晒す

vim

keymapとvimrcを整理しました これで500行くらい. まだまだ初心者です. " まずはNeoBundleの設定 if has('vim_starting') set runtimepath+=~/.vim/bundle/neobundle.vim/ endif call neobundle#rc(expand('~/.vim/bundle/')) " Let NeoBundle manage NeoBun…

Rで作った地図に統計情報を載せたい

R

アメリカの犯罪者状況を可視化したい Rのデフォルトで使えるデータフレームで, アメリカの州別の暴力犯罪発生率が入っているものがあります. >data(USArrests) >head(USArrests) Murder Assault UrbanPop Rape Alabama 13.2 236 58 21.2 Alaska 10.0 263 48 …

Rで地図を描きたい

R

地図上にデータをプロットしたいときってありますよね.そのための第一段階として, 地図を描きます.ggplot2パッケージを利用します. library(ggplot2) world <- map_data("world") g <- ggplot() + geom_polygon(data = world, aes(x = long, y=lat, group = …

Rで作ったsnowクラスタの操作

このブログが凄いhttp://d.hatena.ne.jp/MikuHatsune/並列計算のコードを拝借します. すみません.clusterCallを使いたい. 高水準関数を使わずに並列計算します.結論から言うと, 面倒なのでもうやりたくないです.snowクラスタを作って操作するのは難しいです.…

beamerのテンプレ

特に意味は無いです \documentclass[10pt, dvipdfm]{beamer} % \AtBeginDvi{\special{pdf:tounicode 90ms-RKSJ-UCS2}} \usetheme{Frankfurt} \usepackage{color} \usepackage{amsmath} \usefonttheme{professionalfonts} \renewcommand{\kanjifamilydefault}…

Rで計算時間の比較をしたい(その1)

R

Rは工夫してコードを書かないと遅いことがあります. どの手法が遅いのか. 計算時間を比較して, 最速技法を身につけたいです. 計算時間を知りたければ, time <- system.time( 何かの関数 ) とすればOKです. 計算時間をグラフに表して比較したいです.同じこと…

Rで計算の進行状況を確認したい

R

Rでfor文を回すと遅いです. これを, 高速化とは別の方法で解決したいです. Rでプログレスバーの表示 待ち時間が暇だからよくないのです.コンソールにプログレスバーを表示しましょう. 進捗状況が分かると, 精神が安定しますし.計算の進行に応じて, こんな感…

Rのforeach関数を使って並列計算をしたい (その1)

R

foreach関数とは Rでfor文を回すことはご法度らしいです. 計算が遅いので. Rの高速化のために様々な手段が考えられます. その1つが並列計算です.Rのforeach関数は, Rで並列計算をする際にしばしば登場します. この関数で何ができるのかよく分かりません.なの…

TeX文書で, 本文は2段, 概要は1段で表示させたい

目的 理系なのでTeXで文書を作ります. アブストラクトは真ん中, 本文は2段にしたいです. その方法を調べると2通り出てきました. 実際に文書を作って比べてみました. 結果 \twocolumn[ hoge ]とする方法と, \usepackage{multicolumn}する方法があります. TeX…

Vimrcを晒す

vim

Vimrcを晒します 好きなプラグイン NeoBundle NeoBundleだけインストールしておけば, 後はvimrcにちょちょっと書くだけです. 全てはNeoBundleから始まります. http://kazyury.hatenadiary.jp/entry/2013/04/30/235944 quickrun コード書いて即実行するのに超…

三重対角行列の固有値と固有ベクトルを計算したい

3重対角行列の固有値と固有ベクトルの計算をします. ググってみても, 詳しい計算がまとまった記事が見つかりません. ここでは詳しく書きます. 問題 $n$次の3重対角行列$A$を考えます. \[ A = \left( \begin{array}{cccccc} b & c & 0 & \ldots & 0 & 0 \\ a …

Gershgorinの定理で遊びたい

行列の固有値に関する定理, Gerschgorinの定理があります. スペル分かりません.Gerschgorin? Gershgorin? 2通り観測されます. 正しいほうが分かり次第, 修正します.内容を整理して, 図示して動画も作りたいと思います. Gershgorinの定理 弱い形から. Gershgo…

Rのggplot2でグラフを作りたい

R

Rおじさんだよ 描画が綺麗にできるRパッケージ, ggplot2を使いたいです. 紹介のために, ローレンツ方程式の数値解を描画します. ローレンツ方程式とは. http://en.wikipedia.org/wiki/Lorenz_system ローレンツ方程式はカオス的な振る舞いをする有名な問題で…

quickrun.vimとR

R vim

Macの端末でvimを使っています. .rファイルをquickrunで実行したかったのですが動きません.辛かったので, quickrun.vimを覗いて ’r'のところを \ 'r': { \ 'exec': 'sh -c "R -q --no-save --slave < %s"', \ }, としてしまえば動くようにはなりましたが, こ…

勝利確率で見る巨人戦 (vs. 日本ハム 6月6日)

元ネタ http://blog.livedoor.jp/guusoku/archives/6361018.html野球について, "勝利確率 (WPA)" という指標を紹介します. 勝利確率とは 選手のプレイを評価したいです. ヒットの価値. ファインプレーの価値. どの程度勝利を近づけたのかを知りたいです. 既…

6月6日 巨人 vs. 日本ハム 後半戦

続きです. 6回表 投手交代で笠原が登板. 先頭はアブレイユ. ホームランじゃなければOKです. おそらく首脳陣からも同様の指示が出ていたはず. 選手 打席結果 アウト 走者 スコア 勝利確率 アブレイユ ホームラン 0 なし 2-1 中田 三振 2 なし 2-1 小谷野 フォ…

6月6日 巨人 vs. 日本ハム 前半戦

野球を観に行きました. ちゃんと友だちと行きました. 思ったことを書きます. 試合開始前 オーダー 日本ハム 打順 選手 1 陽岱鋼 2 西川 3 アブレイユ 4 中田 5 大谷 6 小谷野 7 大引 8 鶴岡 9 木佐貫 日本ハム打線. 3番にHRキングのアブレイユ. この人の前に…

飽きた

非負定値対称行列を非負の実数のように扱いたいのです.いずれ書きたいこと. $x > 1$なら$x > \frac{1}{x}$ですよね. 同様に, $A > I$なら \[ A > A^{-1} \]です. $y = \frac{1}{x}$は$x>0$なら下に凸です. 下に凸とは. $0 \[ f(ax + (1-a)y) \] なることです…

半正定値対称行列で不等式とは

半正定値対称行列を非負の実数のように扱いたい...という話その1です.まずは大小関係を定義します. $\geq$という記号を行列に用いて, $A\geq O$を$A$が半正定値対称行列であることを表すことにします.対称行列$A$, $B$を持って来て, $A-B \geq O$とすると何…

半正定値対称行列の和と固有値を計算したい

Matrix Analysisを読んでいます. 間違っていたら直します. 主張 対称行列$A$, $B$の固有値を小さい順に並べた時, 先頭から$i$番目の固有値を$\lambda_i (A)$, $\lambda_i (B)$などと表すことにする.このとき, n次の半正定値対称行列$A$, $B$について, \[ \la…

Cramerの公式を証明したい

Matrix Analysis読んでます. Cramerの公式 行列$A$の$i$列目を列ベクトル$x$に置き換える操作を \[ A\leftarrow ^i x \] と表すことにする.このとき, 正方行列$A$とベクトル$b$に関する連立一次方程式$Ax = b$の解$x$について, \[ x_i = \frac{1}{\det A} \d…

正定値対称行列のHadamard積は正定値対称行列

またHadamard積の話です. 経済周りの線形代数の本で出て来ました 主張 正定値対称行列のHadamard積は正定値対称行列. つまり, Hadamard積$\odot$を, \[ (A \odot B)_{ij} = a_{ij} b_{ij} \] と表すことにして, 正定値対称行列$A$, $B$について, そのHadamar…

ブロック行列の公式を使いたい

覚えたばかりの公式を使いたいので, 使えそうな問題を探しました. Hadamardの不等式 $n$次の正方行列$A$について, 列ベクトルを考えて$A = [a_1 \ a_2 \ a_3 \cdots \ a_n]$とすると, \[ {|}\det A {|} \leq \|a_1\| \|a_2\| \cdots \|a_n\| \] が成り立つ. …

線形数理要論 第4回

問題 $m \times n$行列$A$と$n\times p$行列Bに関して, \[ \mathrm{rank} AB \geq \mathrm{rank} A + \mathrm{rank}B - n \] であることの証明 解答 ブロック行列を使った証明の練習です. 次のようなブロック行列を考えます. $I$は$n$次の単位行列です. \[ \…

線形数理要論 第2回

授業の最後に出てきた式を補足します. 間違っているところは直します. 問題 可逆な正方行列$A$とする。$A$の逆行列の$(i,j)$成分について, \begin{align} (A^{-1})_{ij} = \frac{\Delta_{ji}}{\det A} \end{align} が成り立つことの証明. 解答 次のようなブ…

qtmplsel.vimを導入

Vimプラグインの紹介をします。 序文 TeXで何かを書くときにはテンプレートを利用すると思います。 TeXのプリアンプルは長いですからね。 一例。 \documentclass[dvipdfmx]{beamer} \AtBeginDvi{\special{pdf:tounicode 90ms-RKSJ-UCS2}} \usetheme{Cambridg…

線形数理要論予備試験 第1問

問題 $m \times n $行列$A$と$n \times m $行列$B$に関して, \begin{align} \det \left( I_m + AB \right) = \det \left( I_n + BA\right) \end{align} であることを示す. 解答 簡単な場合から考えます 正方行列の場合 $m = n$のとき, つまり$A$, $B$が正方…

Vimrcを晒す

Vimrc Vimfiler, Vimshell, Unite は入れただけで殆ど使っていません。 これから勉強します。ごめんなさい。 Bundleのままです。NeoBundleに移行していません。ごめんなさい。 " 色々 set set number set clipboard=unnnamed set clipboard=autoselect set s…

ARC005 A問題 大好き高橋くん

Haskellは光です。 問題 高橋君は自分が大好きなので、自分の評判を調べるために自分の名前が入っているツイートを調べたいと考えました。しかし彼のフォロー数は多いので、タイムラインを 1 つずつ確認して自分に関係あるつぶやきを探すのは面倒です。 そこ…

ARC003 問題B 逆さま辞書

やっぱHaskellって凄いわ!という話です。 問題 高橋君は友達とチャットで逆さしりとりをすることにしました。 逆さしりとりとは、前の人が言った単語の頭文字で終わる単語を答えるしりとりです。 しかし、高橋君は英単語に自信がないのでこっそり「さかさま…

AOJ 0009

問題 n 以下の素数の数を数えろ、という内容。 Prime Number Write a program which reads an integer n (n ≤ 999999) and prints the number of prime numbers which are less than or equal to n. A prime number is a natural number which has exactly t…

AOJ 0003

問題 GCD and LCM Write a program which computes the greatest common divisor (GCD) and the least common multiple (LCM) of given a and b (0 Input Input consists of several data sets. Each data set contains a and b separated by a single space…

AOJ 0003

初心者のC++ 問題 Is it a Right Triangle? Write a program which judges wheather given length of three side form a right triangle. Print "YES" if the given sides (integers) form a right triangle, "NO" if not so. Input Input consists of sever…