🚀 go-pugleaf

西村＠秋田県立脳血管研究センターです。

In article <chaa-0202961353300001@news.t.u-tokyo.ac.jp>,
chaa@hgc.ims.u-tokyo.ac.jp (chaa) wrote:
> 土田＠東大医科研と申します
> 次の様なMacintosh 用 soft を探しています
> 
> 実験データを表計算風に入力し，グラフを書き，その面積を計算できるようなソ
フ
> トは
> ないでしょうか．　Excel や普通のグラフ用ソフトでは面積計算ができないよう
な
> のです

　http://www.akita-noken.go.jp にｍｃ(matrix calculator) というフリーの
ソフトがあります。データを表計算ソフトで入力し、データをテキストファイル
にとして保存するかクリップボード経由で行えます。例えばデータが x,y のｍ
行２列のデータであるとすると、このデータを表計算ソフトで入力し処理したい
領域をクリップボードにコピーします。あとは次の命令をｍｃで実行すると望み
の処理ができます。

a = getClipData(0);             ~ クリップボードのデータを変数 a に代入
parametricPlot(a[#,1],a[#,2]);  ~ 第１列をｘ軸に、第２列をｙ軸として作図
menseki(a[#,1],a[#,2]);         ~ ｘとｙで閉じられた面積を計算します。
　　　　　　　　~ 注意、ｘ軸からの面積であれば最初と最後にｙ＝０が必要

　もしサンプル間隔が広く滑らかな曲線として補間（スプライン）してから計算
したい場合には

b = hokan(a,40,2);              ~ 40点に補間
parametricPlot(b[#,1],b[#,2]);
menseki(b[#,1],b[#,2]);

　で計算できます。またデータにアーチファクト等が含まれデータの編集が必要
なら

c = grEdit(a,1,2);

　でデータをマウスで選択できます。また選択条件（大小だけではなく関数でも）

があるなら selectIf 文でできます。

> 
> 2つの重なりのあるピークを別々の式で，別個に面積計算ができるとよいのです
が
> ．

　これはモデル式があるならガウズ・ニュートン法ででき、専用の関数があります
。
ただしパラメータ数が多いと発散しやすいのでシンプレックス法であたりをつけて

からガウズ・ニュートン法を使えばできます。このシンプレックス法はモジュール

の形で含まれています。
　たぶん重なりのあるピークを分離する妥当性の検討（統計計算）もする必要があ

ると思われるので結局はｍｃのようなプログラム可能なソフトが良いのではないか

と思います。

　以上参考まで

-- 
秋田県立脳血管研究センター
脳神経外科学研究部門
西村弘美 nishi@akita-noken.go.jp