Osaka.Stan #5 で LT してきた

『StanとRでベイズ統計モデリング』読書会(Osaka.Stan#5) : ATND speakerdeck.com これ↓についてLTで発表してきました。 kujira16.hateblo.jp 社会心理学や教育学が専門の人が多いコミュニティなのでIT系勉強会とは違った趣があっていいですね。

順列の簡潔データ構造

助教の先生と @yurahuna で簡潔データ構造の輪講をやっています。今週は順列のところを担当したのですが,とてもトリッキーで楽しいデータ構造だと思ったので一部を公開しておきます。これまでの章の内容(定数時間の rank, succ, pred を実現する Bitvector…

報酬が分散未知の正規分布に従うときのThompson Sampling

前回 kujira16.hateblo.jp 問題設定 行ったのはHondaらのAISTATS 2014の論文 [1] の中の実験の追試で,報酬が , の正規分布に従うアームと , の正規分布に従うアームがそれぞれ1本ずつあるという設定の多腕バンディット問題です。 原理 最初にパラメータ を …

報酬がベルヌーイ分布に従うときのThompson Sampling

報酬がベルヌーイ分布に従うときのThompson Samplingの追試をだいぶ前にやったのですが,あまり受けが良くなかったのでブログで供養をしていきます。 問題設定 Thompson Samplingは多腕バンディット問題に対する方策の1つです。多腕バンディット問題について…

StanでAizu Online Judgeの難易度・習熟度を推定したい(3:IRTモデルによる習熟度推定)

シリーズ一覧 kujira16.hateblo.jp kujira16.hateblo.jp はじめに 前回までの記事を公開したところ,Twitterで「問題に取り組んだときの正答確率の部分を項目応答理論でモデリングしないのはなぜか」というコメントをいただきました。 …すいません,項目応答…

StanでAizu Online Judgeの難易度・習熟度を推定したい(2:人工データによる実験)

シリーズ一覧 kujira16.hateblo.jp kujira16.hateblo.jp 目的 あるユーザがある問題に正答したというデータは得られますが,正答していない問題は,取り組んだけれども実力不足で解けなかったのか,そもそも取り組んでいないのかを区別することができません…

StanでAizu Online Judgeの難易度・習熟度を推定したい(1:モデル式)

シリーズ一覧 kujira16.hateblo.jp kujira16.hateblo.jp はじめに Aizu Online Judge (AOJ) という競技プログラミングの練習サイトがあります。 AIZU ONLINE JUDGE: Programming Challenge 解けるか解けないかくらいのちょうど良い難易度の問題に取り組むこ…

Stanの実行時間がchainごとに違う理由

warm-upよりもsamplingのほうが早く終わるのもよく分からない— しょラー (@shora_kujira16) 2017年4月15日 leapfrog法のステップ数LはPyStanでは fit.get_sampler_params()[i]['n_leapfrog__']) で取れるようです— しょラー (@shora_kujira16) 2017年4月17日

Indeed Machine Learning CodeSprint 2017

www.hackerrank.com 問題 求人の文面が与えられるので,以下のタグを付けるべきかどうかそれぞれのタグについて二値分類してください。 アルバイト フルタイム 時給制 月給制 短大卒対象 大卒対象 修士または博士対象 免許が必要 1年の経験が必要 2〜4年の経…

Go言語で競技プログラミングするときに使ってる入出力

Go

入力について 普通は fmt.Scan で問題ないのですが,高速な入力が要求される問題では fmt.Scan だと遅すぎることがあります。fmt.Scan では 3*106 個の整数の入力に10秒以上かかりました。 高速な入力のためには bufio を使います。bufio を使った方法では 3…

RUPC 2017 Day2 I: Tree-Light

http://judge.u-aizu.ac.jp/onlinejudge/cdescription.jsp?cid=RitsCamp17Day2&pid=I 深さ優先探索をしたときの訪問順序によって頂点を並べると,ある頂点を根とする部分木に含まれる頂点が連続した区間に含まれる性質を使います。以下のページが詳しいです…

AtCoder Regular Contest 038 C:茶碗と豆

C: 茶碗と豆 - AtCoder Regular Contest 038 | AtCoder 100点解法 見るからにNimとかGrundy数っぽい雰囲気を感じるので,ググってみると良質な記事がたくさん見つかります。今回は以下の記事を参考にしました。 pekempey.hatenablog.com 最初に,左からi番目…

NAISTから自転車で行ける本屋の紹介

私がNAISTの寮に引っ越してからもうすぐで1年になります。私は原付や自動車を持っていないので,移動はほとんど自転車頼りです。本記事では同じような状況の人のためにNAISTから自転車で行くことができる本屋を紹介したいと思います。 私は情報科学研究科の…

Androidスマートフォンのディスプレイが映らなくなったときの悪あがき

東京への遠征中にスマートフォン(ASUS Zenfone 2 Laser)のディスプレイが突然映らなくなってしまいました。「USBケーブルを接続するとバイブレーションが動く」「カメラを起動してから画面をタッチするとフォーカスを合わせる音がする」といった状況から,…

/usr/bin/xauth: timeout in locking authority file .Xauthority

研究室の計算サーバでX転送を使おうとして ssh -X hoge としてログインしようとすると,なぜかログインに数十秒の時間がかかり,その後タイトルのメッセージが表示されました。 最初はパーミッションあたりに問題があるのかなと思って試行錯誤してみたのです…

Go言語のデバッガをMacのIntelliJ IDEA CEで使う

IntelliJ IDEAのGo言語プラグインをインストールするとdelveというデバッガもインストールされます デバッガのように他のプログラムにアクセスするプログラムは安全なものではないため,delveを実行する際にパスワードの入力を求められます 証明書を作ってデ…

パワーポイントの発表者ツールを無効にする (PowerPoint 2016 for Mac)

学生・研究者のための 使える!PowerPointスライドデザイン 伝わるプレゼン1つの原理と3つの技術作者: 宮野公樹出版社/メーカー: 化学同人発売日: 2009/04/10メディア: 単行本購入: 31人 クリック: 748回この商品を含むブログ (15件) を見る 研究発表のための…

Preztoのforkを探してみる

zsh

github.com Preztoはzshの設定ファイルの詰め合わせで,よく書かれた便利な設定をプラグイン機構によって手軽に導入することが出来ます。 私は数ヶ月前にzshの設定ファイルを作り直したときにPreztoを導入しました。とても便利です。 kujira16.hateblo.jp 残…

Matplotlib 2.0.0rc1を試す

ダサくなくなったという話を聞いたので,試してみました。 私はpyenv+Anacondaの環境を使っています。せっかくなので分離した環境を作って,そこにインストールすることにしましょう。 # `matplotlib-v2` という名前で新しく環境を作る。 # Pythonのバージョ…

セグメント木をあきらめた人のための平方分割

この記事はCompetitive Programming Advent Calendar 2016(その2)の12月15日の記事です。 www.adventar.org はじめに 基本事項 1点に対する変更クエリ・区間に対する質問クエリ Range Sum Query Range Minimum Query 区間に対する変更クエリ・1点に対する…

Jupyter Notebookの小ネタ (display, tqdm, SSHポート転送)

この記事は jupyter notebook Advent Calendar 2016 の12月13日の記事です。 qiita.com Jupyter Notebookを使う上で,この機能を知らない人を見つけたら教えてあげたくなるような小ネタをいくつか紹介します。 IPython.display.display Jupyter Notebook上で…

ACM-ICPC 2016アジア地区つくば大会 参加記

結果 5完26位でした。 Results and Standings | ACM-ICPC 2016 Asia Tsukuba Regional メンバー @shora_kujira16 @menphim M2の先輩 @odashi_t(コーチ) 会場に行くまで 奈良に新幹線なんてあるわけないじゃん! 移動時間の見積もりで,1時間に2本のバスに…

末尾がトップレベルドメインになっている英単語を探してみた

del.icio.us や honokak.osaka のように接尾辞がトップレベルドメインになっている単語のドメインを買いたくなったので探してみました。 英単語のリストは以下のページから拝借 www.eigo-duke.com トップレベルドメインのリストはWikipediaから トップレベル…

リクルートホールディングスのサマーインターンシップ 2016に参加しました

r.recruit-jinji.jp 9月1日(木)~ 9月16日(金)の期間でデータアナリストコースに参加していました 私のスキルセットならエンジニアコースでは?と思う方もいらっしゃるかもしれませんが,学部の卒業研究では情報検索に取り組んでいましたし,リクルート…

サイボウズのエンジニアサマーインターン 2016に参加しました

blog.cybozu.io 8月22日(月) 〜 26日(金) の回に参加していました。 私が参加したのはWebアプリケーション開発コースで,kintoneに新機能のプロトタイプを実装するというのが課題でした。 1日目 初日から台風に見舞われました… 私が出勤した頃には雨はそれほ…

ISUCON6予選にチーム「試運転」で参加しました

isucon.net 私と@menphim, @yurahunaで予選2日目に参加しました。 事前にやったこと 使用する言語の選択では,メンバー全員が研究で多少は使っていることから,Pythonを選択することに決めました。過去のISUCONでの本選出場者数から考えるとPythonはやや分が…

FXで有り金全部溶かして「ぬ」と「ね」の区別がつかなくなったときのための手書き文字データセットを作った

経緯 パターン認識の講義で高次元データを用意して次元削減を適用する課題が出たので,FXで有り金全部溶かした人のためのデータセット作ってる pic.twitter.com/K6emS0YVXj— しょラー (@shora_kujira16) 2016年7月14日 . @shora_kujira16 左:「ぬ」と「ね」…

ICFPC 2016にチーム「試運転」で参加しました

ICFP Programming Contest 2016: Official Site: Task Description ICFPC初参加です。 チームメンバー @yazaten @yurahuna @menphim @shora_kujira16 チーム名「試運転」でNAISTの学生とNAISTの学生予備軍で参加しました。 pic.twitter.com/zGE84UeIT1— ゆら…

TypeScript + power-assert + Mocha + Karma + webpack を使ってブラウザ上でJavaScriptをテストする(2016/12/19追記)

github.com 2016/12/19追記 現在では様々な状況が変わっており,この記事の内容だけでは動作しません。コメント欄に追加の情報がありますのでそちらを参照してください。 やりたいこと JavaScriptを書くのがつらい → TypeScriptを使って型のあるプログラミン…

JavaScriptで作った微妙なツールの紹介

単独の記事で紹介するほどではないかな,と思ってブログには書いていなかったのですが,少し前からよくJavaScriptで微妙なツールを作ってGitHub Pagesにアップロードしています。作ったツールの中から,多少は使い道がありそうなものをいくつか紹介します。 …

ACM ICPC 2016 国内予選参加記

進学したので新しいチームで参加しました。 チーム構成 学部生がいない大学院に進学してしまったので,M1と早生まれ or 飛び入学のM2だけでチームを作る必要がありましたが,なんとかなりました。 大学の近所でこういう看板を見つけたので,チーム名はchasen…

100 numpy exercisesで見つけたNumPyの便利そうな機能

github.com 今まで知らなかった関数をいくつか見つけたので,それについて調べました。 numpy.ndarray.flags numpy.ndarray.flags にはいろいろなフラグが設定できる。たとえば writeable を False に設定した行列に対して書き換えを行うと例外が発生する。 …

hyperoptでハイパーパラメータをチューニングしてみた

Kaggleとかやってる人に人気のあるアレです。 github.com インストール ここ数ヶ月の間,開発が停滞していて,Python3の互換性の問題を修正するプルリクエストが取り込まれていない状況です。 不幸にもPython3に関する問題に追突してしまった場合には,修正…

Wordの数式ショートカット (Word 2016 for Mac)

Mac

学生・研究者のための 使える!PowerPointスライドデザイン 伝わるプレゼン1つの原理と3つの技術作者: 宮野公樹出版社/メーカー: 化学同人発売日: 2009/04/10メディア: 単行本購入: 31人 クリック: 748回この商品を含むブログ (15件) を見る 設定方法 「ツー…

NAIST周辺の生活環境(2017年版)

基本的な情報については他のサイトでも見つかるので,私は最近更新されたことについて書きたいと思います。 (2017年版)最寄りのローソンが閉店した NAISTから徒歩10分(寮からなら5分)のところにあるローソン高山サイエンスタウン店が2017年7月末に閉店し…

Preztoを使ってzshrcをシンプルに作りなおした

追記 ここ10ヶ月ほど本家リポジトリのmasterブランチの更新が無く,先行きが不透明なのでforkを探してみました。 kujira16.hateblo.jp はじめに 研究室にzshのヘビーユーザーがいたので,それに触発されてzshの設定を作りなおしました。 私がzshを使い始めた…

サークルでPython勉強会を開いてみた

speakerdeck.com speakerdeck.com . @scaled_wurm 「今インストールしている公式処理系を消してanacondaで入れなおそう」の意味でした。スライドだけ見ると意味不明ですね— しょラー (@shora_kujira16) March 28, 2016 サークルの後輩が「今までC言語しか触…

kd木のライブラリjtsiomb/kdtreeで最近傍探索を高速化してみた

次のような問題を解く必要がありました 3次元空間中に点群 と がある。ユークリッド距離が最小となるペア を報告せよ(決定的に振る舞うなら近似でもOK) もっと良い方法もあるのかもしれませんが,今回はKD木を使って くらいの計算量で解くことにします。調…

CMakeとNuGetの相性が怪しい?

C++

C++のビルドツールとして、先月くらいからCMakeを使い始めました。理由はVisual Studio用のソリューションファイル生成機能を使いたかったからです。 このプロジェクトではBoostも使用しているのですが、せっかくなのでNuGetを使ってBoostを導入してみました…

私がIIJmioを契約したときの罠の話

今月 (2016年1月) がドコモの2年縛りの更新月だったので,紆余曲折ありましたがIIJmioに乗り換えました。 紆余曲折の具体的な内容 ( ^o^)デビットカードでIIJmioを契約するぞ! ( ˘⊖˘) 。o(ん,待てよ?先月調べた時には使えたはずだぞ?) IIJ┗(☋` )┓三 IIJ…

Okapi BM25をスパース行列のまま計算するPythonライブラリを作った

情報検索で使われる単語の重み付け方法のひとつにOkapi BM25というものがあります。文献によって細かな違いはありますが,今回は Wikipediaに載ってるやつ を使うことにします。 それぞれの文書におけるそれぞれの語の重みを表す行列を計算するクラスを作り…

MS Wordで論文を書くことを強いられているんだ! Vol.2

前回 kujira16.hateblo.jp 数式ツール Word 2007以降かなり使いやすくなっているのだが,インライン数式を多用するような場合にはさすがにLaTeXのほうが使いやすいと思う。同意できない人は以下の文をWordで入力してみてほしい。 $adj(v)$は頂点$v$に隣接す…

FF13 FF13-2 LR-FF13 (Steam PC版) レビュー

無印から順番に数ヶ月前からプレイしていたのですが,正月休み中に一通りクリアできました。PS3を持っていないので,全て初見プレイです。 コントローラー ゲーム自体にはキーコンフィグを変更する機能が付いていないので,コントローラーによってはまともに…

くまのプーさんのホームランダービー!

本年の打ち初めはオウルに惨敗という結果に終わりました。 応援ありがとうございました #打ち初め pic.twitter.com/So53e3FF9S— しょラーbot (@shora_kujira16) 2015, 12月 31 pic.twitter.com/KmzdBuo606— しょラーbot (@shora_kujira16) 2015, 12月 31 pic…

お誕生日コンテスト X - この問題はほんとうにひどい問題であるため,できれば先に他の問題のほうをお楽しみいただければと思っておりまして,ですので他の問題を通し終えて暇になり,かつその暇を

問題 birthday0410.contest.atcoder.jp 解法 343点までしか取れていませんが,満点を取ろうとするとキツいので勘弁して下さい。 まずはノイズを取るためにメディアンフィルタをかけます。注目している画素と,そのマスの周囲8画素の合計9画素の中で黒と白の…

MS Wordで論文を書くことを強いられているんだ! Vol.1

若干はみ出ても良いから1行に収めたい 画面上方(文書とリボンの間)にあるやつを使って右インデントを右側にずらす. 脚注 参考資料タブの中にあり,普通の使い方をする分には使いやすいのだが,「左下に置きたい」,「右下に置きたい」,「線の太さを変え…

Matplotlibでお絵かきする練習

# coding: UTF-8 from __future__ import absolute_import, division, print_function, unicode_literals import matplotlib.pyplot as plt from matplotlib.patches import Circle, Wedge, Polygon, Ellipse plt.figure(figsize=(8, 8)) ax = plt.gca() # …

PythonでEMアルゴリズムを可視化してみた(混合正規分布)

可視化してみました。 解説 以下の動画が一番分かりやすいです。 www.youtube.com Q関数を最大化する, , の求め方は,実際に解いてみると,本当にラグランジュの未定乗数法を使うだけという感じでした。ただし という制約を組み込むのを忘れずに… について求…

楽しいGNU Parallel

「スクリプト言語で書いた処理を合計で5ギガバイトくらいあるHTMLファイル群に適用したい。数時間なら待っても良いが,数日かかるのは困る」という時に役立ちました。 find input -name "*.html" | parallel --dry-run --progress --joblog /tmp/parallel.lo…

ICPC アジア地区予選 2015 C問題 Sibling Rivalry

問題 http://judge.u-aizu.ac.jp/onlinejudge/contest/ICPCOOC2015/all_problems.pdf 解説 ゲームっぽい問題でもゲーム木みたいなものを書けば良いとは限らない。特に,状態遷移に閉路が生じるときはダメ。ベルマンフォードっぽい方法で,ゴールから決めてい…