雑記

雑記: sklearn.cluster.KMeans で k-means 法

sklearn.cluster.KMeans をつかってみたいときとかあると思います。iris データにつかうと以下のようになります。 # -*- coding: utf-8 -*- import pandas from sklearn import datasets from sklearn.cluster import KMeans from sklearn.metrics import co…

雑記: Keras で自作レイヤー(Dense)

Keras でオリジナルの自作レイヤーを追加したいときとかあると思います。 自作レイヤー自体は以下の記事でつかったことがありますが、これはウェイトをもつレイヤーではなく、最後にかぶせて損失関数のみをカスタマイズするためのレイヤーでした。 Keras で…

雑記: 対象日付が米国サマータイム期間か判定(R, Python)

時刻カラムがあるようなデータのプレ処理で、各レコードが米国サマータイム期間かどうか自動判定したいときとかあると思います。サマータイムかそうでないかで東京とニューヨークの時差が違うからね。 とりわけ、時刻カラムのタイムゾーンは日本時間なのです…

雑記: 交差エントロピーって何

機械学習をしているとよく交差エントロピーを最小化させられると思います。 でも冷静に考えると交差エントロピーが何かよくわかりませんでした。むしろエントロピーがわかりませんでした。 以下の記事を読み、もっと無理がない例で短くまとめたかったのです…

雑記

今日勉強会で、論文を読むのに、内容を1枚のスライドにまとめて数を何十本もこなすという話を聞きました。 同じ話が以下の記事に取り上げられています。論文1本にかける時間を切っておくのもコツみたいです。 lafrenze.hatenablog.com 自分は読みたいと思っ…

雑記: Python statsmodels をつかってみるだけ

以下の記事を全面的に参考にさせていただきました。 data.gunosy.io 上の記事では電力使用量のトレンドを抽出していますが、ここでは(電力使用量とかなり似通っていますが)R の組み込みデータである UKgas(1960年から1986年までの四半期ごとのイギリスの…

雑記: X-means法

X-means アルゴリズムを実装してみます。 参考文献 使用データ やること 手順 スクリプト 参考文献 TLに以下の記事の引用ツイートが流れてきたのでやってみます。 スクリプトを全面的に参考にさせていただきました。 aaaazzzz036.hatenablog.com 使用データ …

雑記: 穴を数える手続き

あなたは以下の図形 にいくつ穴があるか求めるようお客様に要求されたとします。 灰色の部分は面で埋まっているとします。 あなたはこう答えたいです。 だっても何も、図形 の穴は がなす三角形の1個だからです。 しかし、これではお客様へ要件定義の説明責…

雑記: XLConnect の読み込みで数式が NA になる場合

XLConnect でエクセルから R にデータを読み込みたいのに意図通りにいかないときとかあると思います。 今回、読み込み対象のセルが数式になっていて、読み込み結果が目的の数値ではなく NA になってしまいました。 読み込み対象が数式セルであっても大抵は数…

雑記: 読書メモ

数学ガールのシリーズ3作目に「数学ガール/ゲーデルの不完全性定理」がある。 帯の「数学って、不完全だったの?」というコピーがいい。数学ガール/ゲーデルの不完全性定理 (数学ガールシリーズ 3)結城 浩 SBクリエイティブ 2009-10-27売り上げランキング : …

確率論のダイジェスト

「確率論って何?」と訊かれるときとかあると思います。ないけど。そんなときのためにダイジェスト版にしました。 相手のセリフは固定(グレー四角部分)なので協力を仰いでください。 なお、セミナー進度の都合上、以下の本の最初3分の1くらいの内容です。…

雑記: R パッケージのデタッチ

R で library() してしまったけど、やっぱりはがしたいときとかあると思います。 自宅 Windows 上の R 3.1.2では以下でできました。ありがとうございます。 How to unload a package without restarting R? - Stack Overflow > library(ggplot2) > ggplot(ir…

OpenShogiLib をつかってみる

将棋ライブラリをつかってみる(けど将棋の話はほぼない)。 経緯 環境準備 中身を少しみてみる 経緯 強化学習の勉強に、「風が吹く格子世界(教科書の例)」のような模式的な例でなく、経済データの学習のような「これは学習が上手くいっているのか」がよく…

雑記

関連記事: 基礎からのベイズ統計学 Skype読書会(13): 参加メモ - クッキーの日記 ハミルトニアンモンテカルロ法での母数 の探索のようす。初期点除いた3個目の候補点まで。 コマ送りのアニメーションをつくるフリーソフトを探そうと思っていたら GIMP でで…

雑記

意味がないであろうことをやったら意味がなかったけどMOTTAINAI精神で記事を上げておく。 データは以下の日経平均株価(日足)を使用する。 日経平均株価 日足 時系列データ CSVダウンロード データは 2010 年以降を取得しておく。 x_2016 <- read.table("in…

強化学習: ノート10 大反省会

読んでいる本(出典): 強化学習 : Richard S.Sutton, Andrew G.Barto, 三上 貞芳, 皆川 雅章 : 本 : Amazon.co.jp この記事はただのデバッグ記録です。 経緯 前回の記事の、「レンタカー2営業所問題」のスクリプトがバグっていた(現在は修正済)。 調査 …

雑記

マハラノビス距離 参考文献: Mahalanobis distance - Wikipedia, the free encyclopedia 点 と分布の距離であって次のように定義される。ただし、 は分布の平均ベクトルと共分散行列。 意味的には、ある点が、ある分布の平均から標準偏差の何倍離れたところ…

雑記: AlphaGoって何

AlphaGo のニュースでもちきりだけどよくわかっていない自分(囲碁わからない)への自分なりのまとめ。 参考文献(出典) AlphaGo | Google DeepMind Official Google Blog: AlphaGo: using machine learning to master the ancient game of Go Mastering th…

雑記

確率論セミナーの発表準備 テキスト: 確率論 (岩波基礎数学選書) | 伊藤 清 | 本 | Amazon.co.jpひっかかったのが13ページの例題1.2 (ii)確率ベクトル に対して を示せ。 ヒント: 日本語に訳すと「確率ベクトルの平均ベクトルの長さは、長さの平均以下」 …

雑記: シンタックスハイライトの配色変更

シンタックスハイライトの配色を再考して更新したのでまとめ。飽きたらまた変えるかも。 Ex. は自分の過去の記事での例。R と CSS はコードがほぼハイライトされるがハイライトとは何だったのか。 CSS class (property) color Ex. Python Ex. R Ex. CSS .ent…

少しだけマークダウンしたい

要求仕様 少しだけマークダウンしたい 場合によってはHTMLタグも自由に書きたいのでパッケージは使いたくない 頻繁に使用する箇条書きやテーブルだけマークダウンで対応したい 設計 convert.bat に hoge.txt をドロップしたらテンプレートをもとに hoge.html…

弱定常過程

弱定常過程でかつ正規過程(ガウス過程)であったら(強)定常過程である理屈のメモ:(強)定常過程の定義(Wikipedia より) 確率過程 が時刻 にとる値の同時分布の累積確率分布 が任意の 、 に対して次式を満たすならば は(強)定常過程である。 弱定常…

雑記

信用リスクモデルにはどんなものがあるか知りたい このドキュメントはいろいろ載っていそう? http://ms.mcmaster.ca/~grasselli/Credittext2011.pdf最初から読んでみると、 定義 1. zero-coupon bond を「保有していると満期 に1ドルが支払われる債券」とす…

スタイルシートの上書き内容メモ

(後日追記)注意事項 スタイルシートは随時更新しているので、この記事の css をコピペしてもこのブログと完全に同じ外観にはなりません。最新の css はソースを表示して、blog_style が含まれるリンクから最新のスタイルシートをたどって、/*** css を上書…

前受収益とか

前受家賃を負債とするのはなぜか、が過去に質問サイトに寄せられている。 http://detail.chiebukuro.yahoo.co.jp/qa/question_detail/q1121308461これに対して、 「もし契約を解除されたら返さなければならないから。」 「サービスを提供する義務を負債とし…