2012-03-01から1ヶ月間の記事一覧

doMC

RでdoMCを使ったお手軽並列計算から、doMCの使い方。 Mac、Linux、Windowsで並列計算用Rパッケージが使える、使えない、があるらしい。とりあえず今はMacを使っているのでそこらへんの検証は後回し。 install.packages("doMC") library(doMC) #計算コア数ベ…

foreachとdoSNOW

foreach+doSNOWパッケージを使って、並列処理をやってみたから、foreachとdoSNOWの使い方その1。 install.packages("foreach") install.packages("doSNOW") library(foreach) library(doSNOW) #基本系(返却値リスト) foreach(i = 1:3) %do% {sqrt(i)} #.c…

使うCPUはいくつがいいのか

こんな質問があったので、これを改変してみる。 aaa2<-function(x){ mean(rnorm(100)) } elapsed.time <- matrix(0, 8, 2) colnames(elapsed.time) <- c("single", "multi") for(i in 1:8){ cl<-makeCluster(i, type="SOCK") clusterExport(cl, "aaa2") elap…

snow

Rで並列処理 (snowライブラリ)から、パッケージsnowの使い方その1。 install.packages("snow") library(snow) #localhostに2個のノードを割り当てられます。 cl <- makeSOCKcluster(c("localhost", "localhost")) #試しにgetwd()でそれぞれの環境のカレント…

並列計算をやる

とにかく、並列計算をしたいけど、スクリプトをどう書けばいいかわからない人(自分のことだが)用に、いくつかのサイトからスクリプトを集めてきた。勝手に転載してます。何かあったら言ってください。 環境 Mac OS X 10.6.8 2.4 GHz Intel Core 2 Duo 8 GB…

どうすればRで並列計算できるのか〜そもそも並列計算って何?

という訳でまず先生の記事から。snowパッケージというものがあるらしい。うん、なるほどわからん。 こちらにもまとめてある。うん、なるほどわk(ry。 毎度おなじみのR wiki先生とグーグル先生に聞いてみる。 そもそも、並列計算とは何か、という話なのだが、…

並列計算したい

計算量が多くなってきた。 最近、並列計算というのが流行っているらしい。 今までは、Rコンソールとターミナル(コマンドプロンプト)でRを複数起動させて、loopを適度に分割して計算していた、という超アナログ手法だった。 そんなことしなくてもRにも実装…

Rで並列計算

並列計算したい どうすればRで並列計算できるのか〜そもそも並列計算って何? 並列計算をやる snow 使うCPUはいくつがいいのか foreachとdoSNOW doMC

Isotonic regressionの高速化

Isotonic regressionを扱っている。 制約条件行列の作成で、行列が大きくなると時間がかかるので、for loopをなんとか回避して関数の高速化を行った。 library(quadprog) #for loopの関数 f1 <- function(Mat, multiple=3){ N <- nrow(Mat) # 観測用量点の数…

randomForest関数

randomForest Random Forests不徹底入門 NA not permitted in predictors Can not handle categorical predictors with more than The response has five or fewer unique values

3次元空間に拡張する

2次元平面から3次元平面でやってみたいと思った。 ただ、微分方程式的に解くのは難しそうだったので、シミュレーション的に計算してしまおう。 飛行機は初期位置から、速度で移動する。 ミサイルは初期位置から発射され、速度で飛行機を追尾する。 微小時間…

6.4化学反応速度論

という反応式があり、濃度をそれぞれと表すと、質量作用の法則から となる。以降、簡単に という反応を考える。がの濃度、がの濃度とすると、 が成り立つ。を消去して という、非線型2階常微分方程式が得られる。を用いると、なので、上式は ただし、とした…

6.3追跡曲線

MGSでスネークがスティンガーミサイルを打って、リキッドやらソリダスやらメタルギアと闘うのだが、発射したミサイルがどうやって目標物にあたるか、その弾道と捕捉時間を求めてみよう。 においてミサイルは原点から、点にいる飛行機(カサッカ)に向かって…

6.2惑星の運動

惑星の運動を考える。導入する値は以下の通り。 質量: 位置ベクトル: 万有引力定数: 太陽の質量: 方向の単位ベクトル: こられを用いると、万有引力の法則より、惑星の運動方程式は となる。と書き、極座標を持つ平面を考えると、であるから、上式は と…

6.1はじめに

前回、線型第2階微分方程式を解いたが、今回は役に立たないらしい。そもそも、一般的にほとんどのものが解析的に解けないらしい。 2階常微分方程式の一般形は、 と書ける。 がのなかに現れない場合、を用いて、1階微分方程式 に帰着できる。 例題として、微…

非線型2階微分方程式

数学モデルで微分方程式を作ろうの第6章。 6.1はじめに 6.2惑星の運動 6.3追跡曲線 6.4化学反応速度論

意識障害(担当MO)

68F 意識障害 救急車で搬送 ABCに異常なし JCS2点 意識障害は定常状態よりもむしろ増悪していくことが多い。 定点比較することが重要。 意識障害の鑑別:AIUEOTIPS A:アルコール I:インスリン U:尿毒症 E:脳症、ホルモン、電解質 O:酸素、麻薬 T:外傷…

初音ミク 大感謝祭

初音ミク大感謝祭 泣いた。 泣いた。 1万人きたって下田さんと咲ちゃんが言ってたよ。 ロイター通信も報じているぞ! 物販?行くのが遅すぎてお目当てのグッズが買えませんでしたけど何か?

初音ミク ライブパーティー ミクパ2012

初音ミク大感謝祭ミクパ 実況の様子 まずは前哨戦として、ファミリーマートの戦争に参加。 こんな状況に。やべえぇぇぇぇぇ。 なんとか、お目当てのクリアファイルはゲットできた。 ミクさんマジ天使。さすがオレの嫁。 あとはお待ちかねのライブだあぁぁぁ…

初音ミク大感謝祭ライブ前日

緊張するな! オフ会の抽選で当たった。うれしい。

ライフサイエンス統合データベースセンター訪問

ライフサイエンス統合データベースセンター(HP)を訪問した。 声優統計は面白そうだ。オレもやってみたい。 合宿の様子。英知の無駄遣いはこのようなものをさすのだろうか...

5.5糖尿病の検査

5.4電気回路網

5.3個人の消費動向

4.2広告に対する売上げ反応で、広告効果をモデル化したが、実際に消費者が買う決心をする過程をモデル化する。 :時刻 :ブランドXの購入レベル :ブランドX志向の動機 :ブランドXの宣伝などによる情報レベル とし、次の微分方程式を仮定する。 これらから…

5.2力学的振動

自然長から力を加えて長さだけバネを伸ばしたとき、質量の質点の運動は となる。ここで、簡単に抵抗なし、外力なしの場合 を解くと、 となる。 抵抗を考えると、 となるが、これはとすると を解くことになる。であるから、判別式による。 a <- 10 b <- 0.000…

5.1はじめに

と表される微分方程式を。線型2階微分方程式というらしい。 これを解くには、特解をひとつ求めることがコツらしい。一般解と特解の差 とすると、とは上式を満たすから、 となる。従って、関数は、斉次微分方程式 を満たす。 一般解は、 の形になる。 、が定…

線型2階微分方程式

微分方程式で数学モデルを作ろうの第5章。 5.1はじめに 5.2力学的振動 5.3個人の消費動向 5.4電気回路網 5.5糖尿病の検査

息切れ(担当EI)

主訴:息切れ 62F 息切れ、喘鳴随伴症状がなにがあるかも考えながら鑑別を挙げる。 喘息:日内変動 肺水腫:チアノーゼ COPD:喫煙 貧血:眼瞼結膜 アレルギー:皮疹 MI:胸痛 狭心症:冷汗(MI、緊張性気胸などでもある) 気道内異物:喘鳴 頸部腫瘍:リン…