超幾何級数を扱うhypergeoパッケージ

世の中には興奮することっていっぱいあるけども、 一番興奮するのはやっぱり超幾何級数を扱うときですよね。超幾何級数についてのWikipediaはこちら。 活用例などはこちらがおすすめです。 ちなみに僕は、非心ベータ分布の平均を計算するために用いました。R…

Shinyで作った確率分布を動かせるページを, NVD3.jsでヌルヌルでインタラクティブなグラフにしました。そしてShinyでD3.jsを使う方法3つ。

Shinyで確率分布を学ぶアプリ作った。あるいは、ShinyでD3.jsを使う方法 from ksmzn 今回のブログの内容を、TokyoR#45で発表しました。 上は、その時の発表資料です。 先日、Shinyを使った確率分布を動かせるページを作りました。 Shinyで確率分布を動かして…

Fedora20に、Python3.4によるvirtualenvwrapperでデータ分析環境を構築する。

Fedora20にPython3.4をインストールし、virtualenvwrapper環境を構築します。 その後、virtualenv環境に分析系ライブラリを入れてみます。 メモし忘れたので、覚えている範囲で。 あくまで自分が試したことであり、ベストプラクティスではないのであしからず…

ggplot2の色をTeXで使うための設定を書いてみた

TeX R

Keynoteもパワポも持っていないので、 スライド作りには現在 TeX のBeamerを使っているわけです。 ggplot2で書いたグラフを張り付ける際には、その他の図形の色もなるべくggplot2の色の範囲におさめたくなります [R] ggplot2の色をプレゼンでも使いたい なる…

Shinyで確率分布を動かして遊べるページ作った

確率分布ってたくさんあるので、統計を学びたてのうちはよくわからないですよね。 パラメータによって確率分布のカタチも変わってしまうし、 似たような種類が多くて混乱してしまいます。 そこで、R言語のWebフレームワークであるShinyを使って、 確率分布の…

latexmkとquickrun.vimを使ってVimでLaTeXを書く。

先日TeXShopの記事(TeXShopのSolarize_darkカラースキーマ - Dimension Planet Adventure 最終章 最終話『栄光なる未来』)を書きましたが、 やはりVimを使いたいということで、 やはりVimしかないということで、 VimでTeXを書く環境を整えましょう。 Vim+Tex…

TeXShopのSolarize_darkカラースキーマ

論文やスライドはTeXで書いてます。 TeXShop、そのままでは味気ないのでカラーを変更しましょう。 世に知られている多くのカラースキーマのうち、特に有名なSolarize_darkに変更するスクリプトを書いたので紹介します。 なお、動作確認はMacOSX mavericksの…

Python3 で MeCab を使う

他のブログにも情報が載ってるけどメモ代わりに。 友人の研究を手伝う可能性が出てきたので、OSXにMeCabを入れましょう。 Pythonのバージョンは以下です。 $ python Python 3.4.2 (default, Oct 19 2014, 17:52:17) [GCC 4.2.1 Compatible Apple LLVM 6.0 (c…

pandasとggplotで回帰分析

Python版ggplotの練習です。 コードと元データはこちらを参考にしました。 ITエンジニアのためのデータサイエンティスト養成講座(6):回帰分析I:回帰分析って何? から、最小二乗法、モデル評価、妥当性検討の実際まで (1/3) - @IT import pandas as pd …

Python版ggplotで日本語を使う

このエントリで書いた例は対症療法的な解決方法です。もっとスマートに書くやり方があったら教えて下さい。 オヒサシブリデス。 Rのggplot2のような綺麗なグラフをPythonでも使いたい! ということで、Python移植のggplotパッケージを先日インストールしまし…

Python3.4にPython版ggplotをインストールしたい

最近環境構築しかしてない。Marvericksにしたせいである。 このまま環境構築し続けて死ぬ。 前回の続き。pandasを忘れていました。 ついでに、Python版ggplotもvirtualenv環境にインストールしてしまいましょう。 pip install pandas OK。 pip install ggplo…

Mac OSX Marvericks にPythonの環境を整える

こん○○わ。 Marvericks にアップデートしたせいで環境構築のやり直しを迫られている今日このごろです。 ムカシ作ったVirtualenv環境が使えなくなっていたのでぶっ壊して、作りなおしました。 Python のインストール $ brew install python python3 さくっと…

Julia 。をインストール

OS X Mavericks でJuliaを使ってみる。 せっかくなのでJulia Studio をインストールする。 ありがたいことにcaskが用意されてるので使ってしまおうそうしよう。 インストール $ brew tap caskroom/cask $ brew cask install julia-studio らくちん! 設定 さ…

Fedora 20 KDE をインストールした際にやったこと

fedora kde config Japanese

gzファイルを解凍せずに調べる

中身を見るだけじゃなくて、いろいろ調べる 中身を見る ファイルが小さいとき # zcat test.gz ファイルが大きいとき # gzip -dc access_log.1.gz | less参考 http://open-groove.net/linux/gzcat-and-gzip-dc/ ファイルサイズと、解凍後のファイルサイズを調…

Pythonで自然数の分割(要素最大値に制限があり・順序を考慮)を求める。

若干研究にも関係する話。自然数の分割とはどういうものか。 例えばWikipediaの例をそのまま借りると、 4の分割は、 4, 3 + 1, 2 + 2, 2 + 1 + 1, 1 + 1 + 1 + 1 である。このとき、普通は順序を考慮しないので 1 + 3 は 3 + 1 に含まれる。 今回、順序を考…

Numpyの固有値と固有ベクトル

勘違いしがちなのでメモNumpyで固有値と固有ベクトルを求める。 In [2]: import numpy as np In [4]: A = np.array([[2,3],[1,4]]) In [6]: la, v = np.linalg.eig(A) In [8]: la, v Out[8]: (array([ 1., 5.]), array([[-0.9486833 , -0.70710678], [ 0.316…

要素が1つのタプルの和

タプル同士の和を求めるとき、要素がひとつしかないとint型とみなされてうまくいかない。Python 3.3.2 In [1]: (1) + (2,3) --------------------------------------------------------------------------- TypeError Traceback (most recent call last) <ipython-input-30-608273b9461c> in <module></module></ipython-input-30-608273b9461c>…

MacVim のLeaderの記述変更

Vim

ちょっと気分転換で、GUIのMacVimでなく、iTerm2からMacVimを起動させるようにした。 ら、Vimのが効かない。いろいろ調べると、.vimrcの let mapleader='¥'を、普通に map ¥ <leader>と変更したら効くようになった。ちなみに、 "<Leader><Leader>で変更があれば保存 noremap <Leader><Leader> :up<CR>とし</cr></leader></leader></leader></leader></leader>…

PyQuery + Requests で画像取得

あえて BeautifulSoup ではなく。 この程度、別にライブラリ使わなくて良さそうだけど。 ホントは何の画像を保存したのかはナイショだ!!

Facebookでも批判書く人はいる

人ははてブを開設すると、無駄に語りたくなる生き物なのです。大学生に見る若者の「承認」欲求の現在(難波 功士) - 個人 - Yahoo!ニュース FBコメントでボロクソに言われていて引いた! FBって批判とか少ないものだと思ってたけど、するひとはするか。 批判…

Pythonで要素がすべて数値のリストを、カンマ区切りの文字列にする

書いたものを載っけるシリーズPython 2.7.5 ','.join(map(str, [1,2,3,4,5])) >> '1,2,3,4,5' これをcsvとかで書き出せばよし。

Pythonで配列の要素をすべて1加算する

書いたものを載っけるシリーズ Gistに載せろや感があるPython 2.7.5 map(lambda x: x+1, [1,2,3,4,5]) >> [2, 3, 4, 5, 6]

Web上でデータの可視化ができるPlotlyをRで試してみる

R

この記事は『R Advent Calendar 2013』12日目の記事です。Plotlyはオンライン上でデータの可視化できるWebサービスです。 作成したグラフをブラウザ上で確認したり、ズームイン・ズームアウトが可能です。有名ドコロですと、ワシントン・ポストが採用したこ…

続 uniq -c 後のスペースをタブ区切りに変換する

教えていただいた。 普通にもう一度awkにつなげればよかったんだ。 なんで気が付かなかったんだろう。。。 $ sort hoge.txt | uniq -c | awk -F"," '{print $1"¥t"$2}' > result.txt

複数のcsvファイルの2行目以降の値を合計してカンマ区切りで出力する.

awk

書いたコマンドを消してしまったので、 覚えているうちに書く。。。複数のcsvファイルの2行目以降の値を合計してカンマ区切りで出力する $ cat *.csv | awk -v OFS=, -F"," 'NR>1 {sum1+=$1; sum2+=$2; num++} END{print sum1,sum2,num} ' > result.txt

uniq -c 後のスペースをタブ区切りに変換する

ソートして重複した行の個数を数えようと、uniq -c をすると、 値が行頭に配置され、しかもスペース区切りになってしまう。 行頭に置くだけならまだしも、文字数を揃えようとしてスペースを勝手に加えてきやがるんですねぇ。 いろいろとめんどいことになるの…

仮にもほどがあるタイトルだなぁ

ブログ作っちゃ消してブログ作っちゃ消してブログ作っちゃ消して まだしっくり来るブログタイトルを思いついた事がないなにか考えついたら変えます。

はてなブログつくった

意識高いやつ