http://w.atwiki.jp/pianissimo/ pianissimo @ ウィキ ja 2010-04-30T19:23:14+09:00 1272622994 https://w.atwiki.jp/pianissimo/pages/29.html このページでは、RとMaximaの文法を比較します。 私はどちらのソフトもそれほど頻繁には使用しないのですが、たまに使用するときに文法を混乱して、毎回成書や解説websiteを当たっています。時間の無駄なので、頻繁に使用する文法だけでも早見表にまとめておこうという趣旨です。 参考文献 [[Maxima 入門ノート1.2.1 中川義行氏>http://www.eonet.ne.jp/~kyo-ju/maxima.pdf]] [[Maxima による数式処理 弘前大学 葛西 真寿氏>http://phys.hirosaki-u.ac.jp/wiki.cgi/maxima?page=FrontPage]] *代入と基本演算 ||R|Maxima| |べき乗|^|^| |100÷7の剰余|100 %% 7|mod(100,7)| |100÷7の整数商|100 %/% 7|| |日常使用関数|sqrt exp log sin cos tan abs|sqrt exp log sin cos tan abs| |変数に代入|a=5|a:5| |ベクトルを代入|a=c(1,2,3)|a:[1,2,3]| |ベクトル&bold(){a},&bold(){b}の要素ごとの積|a*b|a*b| |ベクトル&bold(){a},&bold(){b}の内積|a%*%b　または　sum(a*b)|a.b| |ベクトル&bold(){a},&bold(){b}の外積|outer(a,b)|自前で定義が必要？| |行列Aに代入&ref(matrix.PNG)|A=matrix(c(1,2,3,4),2,byrow=T)|A:matrix([1,2],[3,4])| |円周率を浮動小数点表示||float(%pi)| |自然対数の底e||%e| |虚数単位||%i| *関数 ||R|Maxima| |関数を定義|f=function(x) x^2+1|f(x):=x^2+1| |x^2をxで微分||diff(x^2,x)| |x^2をxで2階微分||diff(x^2,x,2)| |x^2をxで不定積分||integrate(x^2,x)| |数値積分||romberg(x^2,x,0,1)| |cos(x)を0～pi/2で積分||int 2010-04-30T19:23:14+09:00 1272622994 https://w.atwiki.jp/pianissimo/pages/28.html [[Gaussian Fitting>http://nuclear.phys.tohoku.ac.jp/~ykoba/latex2html/gaussian-fitting/]] 問題 (x1, y1) (x2, y2) (x3, y3)の3点を通る放物線(2次関数)を求めよ。ただしx1≠x2, x1≠x3, x2≠x3は保証されているとする。 解答1 y=ax^2+bx+c とおいて、ここに3点を代入すれば求められる。 かなり大変な計算になるが、 a = {(y1-y2)/(x1-x2)-(y2-y3)/(x2-x3)}/(x1-x3) b = (y1-y2)/(x1-x2)-(x1+x2)a c = y1-ax1^2-bx1 となる。 解答2 3点を通る放物線は以下の式で書ける。[[ref>http://detail.chiebukuro.yahoo.co.jp/qa/question_detail/q1113582701?fr=rcmd_chie_detail]] y=a(x-x1)(x-x2)+b(x-x2)(x-x3)+c(x-x3)(x-x1) ここへ3点を代入すると、簡単に、 (x1,y1) → b = y1/{(x1-x2)(x1-x3)} (x2,y2) → c = y2/{(x2-x3)(x2-x1)} (x3,y3) → a = y3/{(x3-x1)(x3-x2)} が求まる。 解答3 2次関数に限らず高次になっても応用しやすい最も一般的な解法は、行列を用いて連立方程式を解くことだろう。[[ref1>http://cse.naro.affrc.go.jp/takezawa/r-tips/r/20.html]] [[ref2>http://flashjp.com/column/gaussjordan.php]] 2010-04-21T07:51:41+09:00 1271803901 https://w.atwiki.jp/pianissimo/pages/27.html R memo C言語とC#とR 単数の変数を扱うのであれば、ほとんど同じと思われる。違いはC言語とC#は型にうるさいのに対して、Rは型のことは気にしなくても欲しい結果が得られることだろう。 ベクトル、行列、テーブルを扱うときに、話が変わってくる。C言語でやろうとすると、forループが必須である。C#では、forでもできるが、個数を意識しなくてよいforeachや、私は使いこなせないがLINQを使うことによって、あたかもブロック単位に演算しているようなコードが書ける。Rではもっと気楽に、単数とベクトルと行列を同じ書式で記述できてしまう。 たとえばn個の実数(double)からなるベクトルaの成分をすべて2倍したものをbに代入するには、 C言語 // 配列bはすでに確保されているものとして、、、 for (i=0; i<n; i++){ b[i] = a[i] *2; } C# その1 var b = new double[a.Length]; for (inti i=0; i<a.Length; i++){ b[i] = a[i] * 2; } //これは単に 配列長 n が必要なくなったというだけ。 その2 var b = from i in a select i * 2; //LINQを使った書き方。速度は遅くなるが、今後PLINQが使えるようになれば別かも。 その3 //ベクトルクラスにして演算子*をオーバーロードすれば、 var b=a*2; R b = a*2 # b <- a*2のほうが推奨される書き方らしいが・・・ ベクトル以上のデータを扱うときは、Rのほうが短いコードで済んでしまうであろう。プログラミングにかかる時間もRのほうが短いであろう（慣れれば）。しかし問題は実行時間である。Rはインタープリター言語なので、実行速度が遅い。プログラミングにかかる時間と実行時間をよく天秤にかける必要がある。1度実行すれば目的が達せられるプログラムであれば、仮にRではコーディングに1時間、実行に24時間、C#ではコーディングに3時間、実行に10分だとしても、Rでコードする人が多いかもしれない。しかし、これを100回実行する必要がある場合（たとえば患 2010-08-27T11:54:26+09:00 1282877666 https://w.atwiki.jp/pianissimo/pages/26.html *MedConの使い方 **インストール [[MedCon website>http://xmedcon.sourceforge.net/]] このサイトで、左のコラムにあるDownloadを選択し、win32 versionをDLすればよい。 2009/12/30現在、xmedcon-0.10.5-1-win32.zipというファイルが取得される。 解凍して、setup.exeを実行する。 インストール後に表示されるreadme.txtの要点は、command lineで実行したいなら、\bin\ディレクトリにある実行ファイルにパスを通しておくこと。絶対パスでいうと、 C:\Program Files\XMedCon\bin である。 **コマンドラインでの使用 そのパスで、medcon.exe -hを実行すると、下記のヘルプファイルが手に入る。 全文は最後に載せるとして、基本的な使い方を書いておきたい。 まず、もっとも単純な使い方は、ヘッダ内容の確認である。 例 medcon -f aaaaaa.v aaaaaa.vというecat7ファイルを開いて、そのヘッダを表示する。 dicomヘッダも表示できるのでなかなかよい。 むしろ、MedConの主な使用目的は、ファイルのコンバートであろう。 例 medcon -c dicom -f *.v ワイルドカードで指定した".v"で終わるファイルをすべて、DICOM形式に変換する。 DICOMは1ボリューム1ファイルなので、いつもの感覚と違うが、スライス数は、(0054,0081) US[1] NumberOfSlices: 63 (2 bytes)のタグに入っているようだ。 いつもの1スライス1ファイルにするには、オプションで指定してやればよい。 例 medcon -c dicom -split3d -f *.v 出力形式には以下のものが選べる。 "ascii" = Raw Ascii (.asc) "bin" = Raw Binary (.bin) "acr" = Acr/Nema (.ima) "anlz" = Analyze (.hdr)+(.img) 2009-12-30T20:41:32+09:00 1262173292 https://w.atwiki.jp/pianissimo/pages/16.html - [[Vistaにおけるファイル同期>http://www10.atwiki.jp/pianissimo/archive/20080511/59e8dbcce1c93f1356675b960e608fb2]] (2008-05-11 22:13:45) - [[Vista usersが共有されるバグ>http://www10.atwiki.jp/pianissimo/archive/20080512/7ccfd326e731052f17c312e8370dda01]] (2008-05-12 22:37:13) - [[Excel 時間差を分だけで表示したいとき(75分など)...>http://www10.atwiki.jp/pianissimo/archive/20080516/ecbd5d28874fad50cb58e5bdde6b43a6]] (2008-05-16 08:23:46) - [[deinococcus radiodurans>http://www10.atwiki.jp/pianissimo/archive/20080523/72e88529e3f5a17b4b175984332c670b]] (2008-05-23 13:19:07) - [[ポリゴン内部にあることを判定する方法...>http://www10.atwiki.jp/pianissimo/archive/20080709/7e9a088df172bf3639829d0466e81bc5]] (2008-07-09 05:27:37) - [[Analyze format spec>http://www10.atwiki.jp/pianissimo/archive/20080907/ade948c54682f9d2780b9010c50986b7]] (2008-09-07 13:03:09) - [[dicom c#>http://www10.atwiki.jp/pianissimo/archive/20090825/7be4258595ee8b0b3967a22292582627]] (2009-08-25 23:56:13) - [[bisphosphonate>http://www10.atwiki.jp/pianissimo/archive/20091 2010-03-09T12:54:37+09:00 1268106877 https://w.atwiki.jp/pianissimo/pages/13.html 2008-05-10T18:40:07+09:00 1210412407 https://w.atwiki.jp/pianissimo/pages/12.html <p> </p> <p> </p> <hr /><p> </p> 2010-01-23T09:46:28+09:00 1264207588 https://w.atwiki.jp/pianissimo/pages/11.html 2008-05-06T14:07:37+09:00 1210050457 https://w.atwiki.jp/pianissimo/pages/10.html * コメントプラグイン @wikiのwikiモードでは #comment() と入力することでコメントフォームを簡単に作成することができます。 詳しくはこちらをご覧ください。 ＝＞http://atwiki.jp/guide/17_60_ja.html ----- たとえば、#comment() と入力すると以下のように表示されます。 #comment 2008-05-06T14:02:13+09:00 1210050133 https://w.atwiki.jp/pianissimo/pages/9.html * 関連ブログ @wikiのwikiモードでは #bf(興味のある単語) と入力することで、あるキーワードに関連するブログ一覧を表示することができます 詳しくはこちらをご覧ください。 ＝＞http://atwiki.jp/guide/17_161_ja.html ----- たとえば、#bf(ゲーム)と入力すると以下のように表示されます。 #bf(ゲーム) 2008-05-06T14:02:13+09:00 1210050133