2013-01-01から1年間の記事一覧

時系列の移動平均による季節調整,放送大学の教材が役立った.「身近な統計」

形容詞から季節成分の除去がうまく行かない. 7年分あるので,差分に直した後,365日の重みで除算すると,極端になってパルスが多く発生. 区間平滑化とか過去の移動平均にした時系列で同じことをするか. n日単位に直してからの方がよいのか. また,形容詞…

金曜からPapers2 (2.7.2)が起動しなくなった.クラッシュ.端末変えたり,再起動したり,再インストールしたりしても同じ. それで窓口に依頼しても,Papers3が出たからそちらを試してよ,とのこと. って前からのアップグレードでも29.5ドルとかかかるんで…

プログラムのソースコードをそのままのせる. \usepackage{ascmac} \begin{itembox}[l]{Rを使った偏相関係数の計算} \begin{verbatim} x <- c(8, 5) n <- x(35, 30) prop.test(x, n, correct=FALSE) #連続性の補正なし \end{verbatim} \end{itembox}とする.…

SWMW2013 in Ehime on 2013.10.19[1] 鳥海さん ネット選挙とイージーマイノリティ 2012年衆議院選挙と2013年参議院選挙 サイレントマジョリティの逆.少数の人間が騒いで大勢が騒いでいるように見える → ツイートしたユーザがどれだけ偏っているか? 指標を…

Gephiで,最初の画面にグラフが現れない場合,メニューバーから, ウィンドウ>グラフ,と選択すれば表示される. Betweenness centralityを計算したい場合は,右の統計メニューから,平均経路長を選ぶと同時に計算してくれる. コミュニティ別に分割して色…

Gphiはコミュニティ表示が分かりやすく,きれいなので使えそう. データをインポートするには,CSVとかPajek形式とかいろいろあるが,一番分かりやすそうなのが,GML形式. データの記述は以下のようにする. ノードの大きさはwとhで指定. graph [ node [ i…

まとめてファイル処理. Corr_1.txt,Corr_3.txt,Corr_2.txtとかに対して行数や,削除処理をしたい場合. wc -l Corr_* 1 Corr_1.txt 1 Corr_2.txt 1 Corr_3.txt 3 total rm Corr_* (Corr_1.txt,Corr_3.txt,Corr_2.txtが削除される.)

ワードでインデントの位置がそろわない場合,まずインデントの数字を右クリックして,段落を選択, その後,左のインデント幅などをそろえるとよい. 意外と,できそうでできなかった.

R R

日経225*1のデータをRからゲット. install.packages(quantmod) library(quantmod) getSymbols('^N225') head(N225) #ファイル確認 sink("day.tmp") N225$N225.Open sink() write.table(N225,file="prices.tmp",append=FALSE,quote=FALSE,sep=" ")どういうわ…

日本語でshift-jisで書き込んだはずのファイル,読み込むと文字化け. 前はちゃんと日本語も読めていたのになぜかと思ったら,macのシステム文字に依存して,読み込みができないことがわかった. macのシステム環境設定で,言語とテキストの優先する言語の順…

行列でランダムシャッフル

R

applyの引数2つ目が「1」か「2」かで行(=1)か列(=2)かが決まる. > tmp <- matrix(1:10,ncol=2,nrow=5) > tmp [,1] [,2] [1,] 1 6 [2,] 2 7 [3,] 3 8 [4,] 4 9 [5,] 5 10 > apply(tmp,1,sample) #行方向のシャッフル [,1] [,2] [,3] [,4] [,5] [1,] 1 2 8…

SSDに取り替えたらnkfがインストールされなかったので,再インストール. sudo port install nkf でnkfをインストールしようと思ったら,そもそもmac portがインストールされていないので,そこから. このサイトの該当(lion)のpkgをダウンロードして,イン…

去年の2月頃に購入していたMac Pro用のSSDをやっと取り付けた. まず専用のマウンタが必要で,それに取り付けるのに一苦労し,どうもねじの大きさが合わなくてピッタリくっつけられなかった部分は,接着剤+セロテープで代用.ベイ1に取り付けると,まずは…

課題図書3+冊.[博士として]田崎晴明 統計力学(上・下)増田直紀,今野紀雄 複雑ネットワーク 基礎から応用までC. M. ビショップ パターン認識と機械学習(上・下)東京大学教養学部統計学教室 統計学入門洋書 George Kingsley Zipf Human Behavior and the Pri…

FIT2013メモ(情報科学技術フォーラム)2013.09.05 [1] まつもとゆきひろ 講演 言語マニア,現存する唯一の松江市名誉市民,初代は小泉八雲 鶏口牛後,代替不可能な人材となるべし すでにあるものを使う,巨人の肩の上に乗る 仕事が楽しくて何が悪い 動いてい…

Dynamical Classes of Collective Attention on Twitter J. Lehmann, B. Gonçalves, J. J. Ramasco, and C. Cattuto WWW '12 Pages 251-260 Twitterのハッシュタグの増減時系列で,(1)ピークと前に盛り上がったもの,(2)ピークと後に盛り上がったもの,(3)ピ…

Macでラジオ音源を録音. RAdikaインストール (NRF = NHK FM) SoundFlowerインストール サウンド設定で,出力をSoundFlower(2ch)に変更 Garageband起動 Garagebandの環境設定で入力をSoundFlower(2ch)に変更 Garagebandに適当なプロジェクトを作って,録音…

「風立ちぬ」という宮崎映画を観た. 零戦を設計した堀越二郎さんの人生を,堀辰雄の風立ちぬに被せた話. 空や雲の書き方が美しくて,それはそれで満足だった. でも,ストーリーは,というと一言で言うと,「妙齢のおっさんのファンタジー」. 二郎にも,…

博士最終試験も終わった.まだまだ,考察が足りない部分とか,パッと数式とけないとか, 筋道立てて話せないとかいろいろ改善点は多い.これまらも気を抜かず精進しなさいということだろう. 私の人生はいつだってそうだ.はっきりと区切りは付けられなくて…

bibtexで,出てきた順で引用文献を出力するには, \bibliographystyle{unsrt}を使う.で,unsrt.bstというファイルが,書式を決定するプログラムになっている. そこで,自分専用のmy_unsrt.bstを作業ローカルフォルダにおいて, \bibliographystyle{my_unsr…

Lean in/Sheryl SandbergKindleにて初めて読破した洋書.不明な英語をぱっと検索できるので,読み遂げられた. でも,左下に読み終わるまでの時間「7時間10分」などと表示されるのだが,読めば読むほど,読み終えるまでにかかる時間が延びて,ちょっと鬱に…

Papers2を使って選択した論文のみ,bibtex形式でコピーしたい場合, Edit -> Copy As -> BibTex \cite Commandを選択すれば, あとはペーストで,bibtex形式のコピーができる.

パワポで作った概念図などを,TeXに取り込む方法 for mac. パワポ1枚で絵を描き,「プリント」ボタンを押す. 左下の「PDF」のコマンド一覧より「PostScriptとして保存」を選択する. ps形式で出力されるので,そのままTeXで表示できる.

博論書いていて,細かいことをチマチマ調べながら書くのが楽しい.すごく. 天文少女だった小学生のころ,独自の天文ノートを作っていたことを思い出した. 天文ノートのページは,水星,金星,地球,火星,木製...と決めて,そこに惑星の直径,主な成分…

失敗の本質 戦場のリーダーシップ篇 /野中 郁次郎 (著) 複数の学者(主に大学教授)が,大東亜戦争時のリーダーシップについて研究してまとめた内容(論文?)をまとめた本. 2012年に出版されていて,前書きには2011年東日本大震災時の民主党リーダーシップに…

博士論文β版提出完了,当日出張につき不在の審査員の一人対する発表会終了で小休止. 5月中から国際学会の準備,博士論文の仕上げなどで息つく間も無く,という感じだった.机の周りも汚いし,やらなければいけないことで忘れていることは無いか,読もうと…

awkで時間をシフトして,任意の時間形式で出力(strftime). BEGIN{ #Shift time #2009-07-18 12:23:38 a = 5 #5 hours shift } { gsub(/-/," ",$1) gsub(/:/," ",$2) str=$1" "$2 #raw time t = mktime(str) #shift time t2 = mktime(str)-(60*60*5) #remove…