RjpWiki はオープンソースの統計解析システム R に関する情報交換を目的とした Wiki ですgraph [vignette: y; demo: n] グラフデータ構造を取り扱うパッケージ. † ↑ RBGL [vignette: y; demo: n] Interface to boost C++ graph library (based on graph objects from the graph package). † ↑ gRain [vignette: y; demo: y] An alternative interface to graphs. Uses the graph and RBGL and implements additional graph operations. † ↑ dynamicGraph [vignette: n; demo: y] Provid
Twitterでフォローしている人同士を線で結んだグラフ(フォローグラフ)を、一部で流行り(?)のR言語で描いてみます。 手順。 1. Twitterのフォロー関係を取得する。 今回は自分がフォローしている230人ほどのTwitter-erについて、それぞれフォローしている人をTwitterのWebページをスクレイピングして取得しました。 curlを使ったシェルスクリプトとperlの1-linerで (for x in ha_ma sudori syou6162 ; do echo $x ; done) | perl -ne 'print "*$_\n"; chomp; $x = $_; $flag = 1; for($i = 1; $flag; $i++) { $get = `curl -s -u USER:PASSWORD http://twitter.com/$x/followers
17:19 08/11/27 TopCoder Code Jam の練習に……と思ってしばらく前から TopCoder のSRMに参加してたのですが、 せっかくなので cafelier@SRM に記録をつけることにしました。 どういう試行錯誤をしながら提出した時のコードにいたったのかを、 できるだけ詳細にメモろうと思っています。 426以前のは記憶から掘り起こして書いたのでちょい大ざっぱですが。 これまで何回かここで書いたような整然とした考え方を本当に自分がしているかいうと、 してないよなー、と薄々思ってしまっているので、じゃあどういう風にやっているんだろうかと。 自分のふり見て我がふり直す。 20:26 08/11/24 論文 PLAN-X 2009 通ったみたいです。ばんざい。 ただでさえD論まったく間に合う気がしないのに、camera ready版なんて作ってる時間が… オートマトン
今さら感はあるが、MCMC (Markov Chain Monte Carlo; マルコフ連鎖モンテカルロ)を使えるようになろうと、まずは簡単な例から試してみた。 手始めに、正規乱数から生成した標本の平均と標準偏差を推定してみる。 やはりRを使用。MCMCpackパッケージを あらかじめインストールしておいて、呼び出す。MCMCpack中のMCMCmetrop1R()関数を利用して、メトロポリス法によるMCMC推定をおこなう。 library(MCMCpack) 乱数系列を初期化。 set.seed(1) 平均10、標準偏差3の乱数を1000個生成して、xに入れる。 m <- 10 s <- 3 x <- rnorm(1000, m, s) MCMC推定に使用する関数を用意する。betaは要素数2のベクトル。beta[1]が平均、beta[2]が標準偏差で、betaを推定する。関数の返り値
リリース、障害情報などのサービスのお知らせ
最新の人気エントリーの配信
処理を実行中です
j次のブックマーク
k前のブックマーク
lあとで読む
eコメント一覧を開く
oページを開く