★データ解析備忘録★

ゆる〜い技術メモ

Tokyo.R #54の振り返りと補足

R

6/18にTokyo.Rに参加したので、遅ればせながら振り返りと(自分の知識の範囲内での)補足をしたいと思います。 初心者セッション 10分で分かるR言語入門 ver.2.19 質疑 パッケージを手動で入れるには? パッケージはinstalll.packagesや{githubinstall}パッ…

【SAS】PROC SQL では LIMIT 句が使えない

SAS

普通のSQLを書くとき、データが巨大だと一部をちら見するためにLIMIT 1000と書いて一部だけを抽出したりします。 SASのPROC SQLではSQLの構文がほぼそのまま使えて便利なのですが、LIMIT句はどうやら使えないようです。 代わりにOUTOBS=1000で指定すれば良い…

SASの(よく使う)ショートカットキー一覧

SAS

よく使いそうなSASのショートカットキーをまとめました。 製品版でしか試してないのでSAS雲丹でも同じかはわかりません。 他のソフトと同じもの キー 機能 Ctrl+C コピー Ctrl+X 切り取り Ctrl+V ペースト Ctrl+A 全選択 Ctrl+S 上書き保存 Ctrl+Shift+S 新…

MacTeXのbstファイル置き場

メモbstファイルは /usr/local/texlive/2015/texmf-dist/pbibtexstyファイルは /usr/local/texlive/2015/texmf-dist/tex/latexにある。usrフォルダをファインダーで開くにはターミナルで open /usr

Treasure Dataの日付範囲指定はイコールを含むのか

Treasure Dataでクエリを書くとき、where句で日付指定をする場合が多いわけですが、範囲はイコールを含むのでしょうか? 以前書いた記事でちょろっと触れてそこに答えもあるのですが、イコールの有無に焦点を当てたものも欲しかったのでメモ。 なのか なのか…

統計の誤解と濫用や「p値至上主義」を憂慮しp値の6原則を発表したASAの声明に対する統計学徒の素人意見

アメリカ統計学会がp値に関して以下のような声明を3月7日(現地時間)に発表したということで注目を集めているようです。AMERICAN STATISTICAL ASSOCIATION RELEASES STATEMENT ON STATISTICAL SIGNIFICANCE AND P-VALUES https://www.amstat.org/newsroom/p…

Matlabで変数の削除

新しい言語、特に対話型のをやると必ず気になるのが変数の除去について。Matlabの場合、は以下のようにやります。変数をすべて除去 clearvars 特定の変数を削除 変数名をvariablesで指定 clearvars variables 特定の変数を残してあとは削除 keepVariablesで…

多項分布とディリクレ分布のまとめと可視化

多項分布とその共役事前分布について、可視化をしながら整理してみたいと思います。 どちらかというと、可視化をしてパラメーターで分布の形がどう変わるのかを見ることがメインです。 多項分布とは 二項分布の一般化と考えればよいです。 「コインを投げた…

『R言語徹底解説』(原著 "Advanced R")を読んだので感想を書く

待望の『R言語徹底解説』が届きました。(2/10) R言語徹底解説作者: Hadley Wickham,石田基広,市川太祐,高柳慎一,福島真太朗出版社/メーカー: 共立出版発売日: 2016/02/10メディア: 単行本この商品を含むブログ (29件) を見る本来は昨年のクリスマスあたりに…

sas7bdatファイルの作成と各種ソフトでの読み込み

sas7bdatとは何か sas7bdatというファイル形式があります。これは一言で言えばSASで作ったデータのファイル形式です。 このファイル形式はSASのいろいろなアウトプットを格納できて、例えば回帰分析後のパラメーター推定値や予測値、ベイズ推定のMCMCであれ…

RMeCabの顔文字の誤解析を後処理で消す

顔文字を含んだ文の解析 {RMeCab}、日本語の形態素解析に便利なのですが顔文字の分類は苦手です。 例えば、 words ふう。お腹いっぱい( ´•౪•`) 暇だな(☝ ՞ਊ ՞)☝ 試験頑張るぞ\\\\ ٩( ‘ω’ )و //// というデータフレームを解析して単語の頻度を数えるとし…

SASでライブラリ内のデータセットを全て削除する

SAS

タイトルの通りSASでライブラリ内のデータセットを全て削除する方法です。 Rでいうと rm(list=ls()) にあたるものですが、SASだとライブラリという概念がある点で異なります。 特にライブラリ名を指定してない時はWORKライブラリが割り当てられています。 例…

【SAS】PROC SQLが強かった件

SAS

SASでデータ集計をやらなきゃいけなくなって、SAS初心者*1の僕がいろいろ調べた結果、SQLプロシージャがよさそう。使ってみた結果、とても便利だということが分かったので使い方を自分用にメモ。 SQL文がそのまま使える MySQLなり、SQLiteなり、HiveQLやPres…

『ベイズ法の基礎と応用』を読みました

間瀬茂先生による『ベイズ法の基礎と応用-条件付き確率分布による統計モデリングとMCMC法を用いたデータ解析-』が昨日に発売されたようなので早速買ってみました。ベイズ法の基礎と応用 条件付き分布による統計モデリングとMCMC法を用いたデータ解析作者: 間…

IMLの行列とSASデータセットを行き来する

SAS

IMLで作った行列をSASデータセットとして使いたいとき。 SASデータセットをIMLに取り込んで行列にしたいとき。 IMLで作った行列をSASデータセットにする 必要なコマンドは create A var { 取り出したい変数名 }; /*Aは作りたいデータセット名*/ append; clos…

TreasureData Command Line のインストール

TreasureDataで、以下のように結果が重すぎてブラウザからダウンロードできないとき、ターミナルからコマンドライン(CLI)でダウンロードできます。 その手順を説明します。 手順 環境設定 Rubyを使います。 1.9.0以上じゃないと対応してないみたいです。 ま…

GitHubを初めて使う準備

GitHubを初めて使うときの準備的な。 要点は1個だけですが(笑) 環境はWindows10 GitHub for Windowsのインストール これをインストールしないと始まりません。 desktop.github.com 最新版をインストールします。 インストールすると二つのショートカットが…

Prestoでの日付の扱い方

Prestoで日付の処理に躓いたのでメモ。TreasureData (TD)ってUDF(user defined functions)が多いから、HiveにしろPrestoにしろ結構躓いたりするんですよね。 公式のマニュアルページが手放せません。 docs.treasuredata.com今回はこのなかで初心者やSQLユー…

重回帰のパラメーター推定理論とSAS/IMLによる行列実装

統計学の「基本の「き」」である重回帰を復習しながらSAS/IMLで行列で実装していこうと思います。 理論 重回帰を式で表すと、 となります。 これを行列で表すと 要素を書き下すと となります。

史上最大の素数が本当に素数なのか確かめてみた

史上最大の素数が更新されたらしいです。wired.jpこういうの見つける人ってすごいですよね、、、 でも この数って本当に素数なの? という疑問がわきます。気になります。 どうしましょう。 確かめてみましょう。素数を見つけるのは大変でも、ある数が素数か…

LaTeXで左上付き文字

僕の通う大学の数学の教授で、転置行列の t を右上ではなく左上に書く人がいます。 これをTexでどう書くのかなーと思って調べたら判明したので自分用にメモしておきます。結論から言うと超簡単でした。 %数式内で ^{t}A 出力 化学の論文とかで質量数12の元素…

Python3.5でワードクラウドを描く

Pythonでワードクラウドを描く方法は、すでに結構な数の記事があるので、基本はそこを見ればいけます。 なので、今更僕が解説するようなことはしません。 例えば qiita.com ライブラリの入れ方も丁寧に説明してくださってます。 今回もこのページをほぼその…

LaTeXで「1」を白抜きにする

インディケーター関数とかで「1」を白抜きにしたいとき。 これは元々のtexには入ってないので自分で定義します。 1とl(Lの小文字)を組み合わせて作るといった感じです。 %プリアンブルで \newcommand{\1}{\mbox{1}\hspace{-0.25em}\mbox{l}} %数式内で \1 (y…

jsonを楽に扱えるjqコマンドを導入(Mac, Windows)

データ解析の一連のプロセスの中で、前処理の占める割合はたいへん大きいです。 その前処理のなかで、json形式のものを楽に扱えるjqコマンドというものがあります。json形式はデータベースなんかではよく見る形式で、自分でも個人的に必要になったのでその導…

WindowsのPython3.5でlxmlを入れようとしたら詰んだ

Webスクレイピングなどで使われるlxmlパッケージですが、WindowsのPython3.5で導入にてこずりました。 pip install lxml でいけるかなあと思ったのですが、普通にエラーになりました。

SASの中でRからSASにデータを渡す

R SAS

前回の記事でSAS/STATでもRを使えることを書きました。 y-mattu.hatenablog.com今回はSAS上で動かしたRのデータをSASに渡す方法です。簡単です。 /*proc iml内でendsubmit;した後に*/ run ImportDataSetFromR("work.A","B"); /*Aは後でSASデータとして使いた…

SAS/STATでもPROC IMLしてRを使う

前回の記事でSAS/IML StudioでRを走らせる方法を書きました。 y-mattu.hatenablog.comで、その最後にSAS/STATでproc imlの中でRを使う方法を書いている方がいたのでそれを紹介しました。SAS上でRを走らせよう SAS/IML | うずまき(4th Gen) | Jun-Systemsその…

SAS/IML StudioでRを使う

SAS R

自然言語処理をやるとき、僕はSAS/Minerを使える環境にないのでなんとかならないかなあと思っていたのですが、 SAS/IML Studio *1 を使うとSAS上でRを動かせるそうです。 ということでやってみます。 support.sas.com 手始めということでMeCabで簡単な形態素…

rstanで個人のパラメーターを推定した話(JapanR2015のLT補足)

昨年の12/5に開催されたJapanRでLTをしました。 rstanで個人のパラメーターを推定した話 from Yuya Matsumura www.slideshare.net 5分しかなかったのであんまりちゃんと説明できなかったのですが要約すると以下のような感じです。 rstanで階層ベイズモデルを…

RとPythonでMeCabを使う環境構築(Windows, Mac)

自然言語処理まわりはずっとRでやっていたのですがそろそろPythonでもやろうと思い環境構築をまとめてみます。 ついでなのでRのほうもまとめてみます。 PythonについてはWindowsとMacで違うので一応そちらも触れたいと思います。Windows10とYosemiteです。 …