tatsu___kun回答ポイント 95ptウォッチ 6

f=194*X^-1.094*Y^-0.278

という関数のグラフを作成したいと考えています。


この数式は、先日日本腎臓学会が発表した、日本人の糸球体濾過量(GFR)推定式で、
変数Xが血清クレアチニン量、変数Yが年齢になります。


http://www.jsn.or.jp/
にノモグラムはあるのですが、2次元で等高線がfの値となるようなグラフ*を作成するソフト(商用ソフト不可)はあるでしょうか？


*例として、このようなものです。(これはBMIの早見表です)
http://www.mmjp.or.jp/stethoscope/bmi.jpeg

• はてなアカウントでログイン
• 無料ユーザー登録

みんなの回答

• るびきち2008/06/05 04:44:56

  満足5pt

  詳しくないですが、gnuplotやgraphvizはどうでしょうか？どちらもフリーソフトウェアです。

  http://d.hatena.ne.jp/rubikitch

  スター

  • tatsu___kun 2008/06/05 07:23:28

    申し訳ございませんが、あなた様より詳しくありませんので、具体的かつ明瞭な例示がないと、回答いただいた内容が有益か否か判定できません。

    それと、URIと回答内容の関係がよくわかりません。URIの参照先には有益と判断されるような内容は見つけられませんでした。

noir_k2008/06/05 09:03:27

満足40pt

gnuplotの3次元プロットしたものを、真上から見れば、等高線になります。

http://t16web.lanl.gov/Kawano/gnuplot/plot3d.html




具体的には、以下のようにすればよいでしょう。




gnuplot> set ticslevel 0

gnuplot> splot 194*(x**(-1.094))*(y**(-0.278))

gnuplot> set contour

gnuplot> set cntrparam levels 10

gnuplot> set contour base

gnuplot> unset surface

gnuplot> set view 0,0

gnuplot> replot

スター

• tatsu___kun 2008/06/05 20:34:18

  ありがとうございました。

  ・・・何故か真上から見た図にならなかったです。

Mook2008/06/05 11:46:19

満足50pt

ツールでも作成はできますが、EXCEL の３D グラフでも一応それらしいグラフができます。

下記は、数式を使用してセルにデータの表を作成し（A1:S17）、挿入から３Dグラフを選んで作成するやり方をマクロで実行した例です。




シートタブをクリックしコードの表示を選択して、そこに下記をコピーしてください。

シートに戻ってAlt＋F8を押し、makeGraph を実行してみてください。

Sub makeGraph()
    Dim i&, j&, age&, creatinine#
    Dim sSheet As Worksheet
    Set gsheet = ActiveSheet
    
    i = 1
    For age = 5 To 80 Step 5
        i = i + 1
        gsheet.Cells(i, "A").Value = age
    Next
    
    j = 1
    For creatinine = 0.6 To 4.01 Step 0.2
        j = j + 1
        gsheet.Cells(1, j).Value = creatinine
    Next

    gsheet.Range("B2").Formula = "=194*B$1^(-1.094)*$A2^(-0.278)"
    gsheet.Range("B2").Copy
    gsheet.Range(gsheet.Cells(2, "B"), gsheet.Cells(i, j)).Select
    gsheet.Paste

    With Charts.Add
        .SetSourceData gsheet.Range("A1").CurrentRegion, PlotBy:=xlRows
        .Location Where:=xlLocationAsObject, Name:=gsheet.Name
    End With
    With ActiveChart
        .ChartType = xlSurface
        .HasTitle = True
        .ChartTitle.Characters.Text = "日本人の糸球体濾過量"
        .Axes(xlCategory).HasTitle = True
        .Axes(xlCategory).AxisTitle.Characters.Text = "クリアチニン"
        .Axes(xlSeries).HasTitle = True
        .Axes(xlSeries).AxisTitle.Characters.Text = "年齢"
        .Axes(xlValue).HasTitle = True
        .Axes(xlValue).AxisTitle.Characters.Text = "GFR"
        
        .Rotation = 120
    End With
End Sub

グラフ上を右クリックして、３D-グラフを選択すると、向きなど細かい調整が可能です。

http://www.relief.jp/itnote/archives/001713.php

スター

• tatsu___kun 2008/06/05 20:34:41

  ありがとうございました。

  マクロで頑張れば出来るものですね。

• その他の回答を読む

この質問へのコメント

• ootatmt2008/06/05 08:36:23
  商用ソフト不可ということは、エクセルもだめなんですね？
• Mook2008/06/05 09:03:03
  graphviz は方向性が異なる気がしますが、gnuplotのサンプルは下記にありました。
  http://t16web.lanl.gov/Kawano/gnuplot/plot3d.html#6.1
  
  直接試してはいないので、コメントで失礼しますが、
  http://search.vector.co.jp/search?query=%82R%8E%9F%8C%B3%83O%83%89%83t
  でも、いろいろとありますので、一度試してみてはどうでしょうか。
  
  参考までに伺いたいのですがグラフを書く場合、血清クレアチニン量の定義域はどの程度でしょうか。
• tatsu___kun2008/06/05 10:02:39
  エクセルは特別なアドオンがなければOKです。クレアチニンの範囲は0.6から4.0を想定しています。
• ita2008/06/05 12:35:06
  こんなかんじ？
  http://f.hatena.ne.jp/ita/20080605120950
  
  gnuplot のスクリプト
  
  
  set xlabel "血清クレアチニン量"
  set xrange [0.6:4.0]
  
  set ylabel "年齢"
  set yrange [1:100]
  
  f(x,a)= (a/194.0*(x**1.094))**(-1.0/0.278)
  
  plot \
  f(x,25) t "25" w l,\
  f(x,50) t "50" w l,\
  f(x,75) t "75" w l,\
  f(x,100) t "100" w l,\
  f(x,125) t "125" w l,\
  f(x,150) t "150" w l
  
  
  pause -1
  
• ita2008/06/05 13:29:11
  ちなみにF(x,y)=a をyについて解いた y=f(x,a)という関数をプロットしてます。
  
• tatsu___kun2008/06/05 20:50:20
  結局id:ita様の回答が最も有用でした。
  ita様には200p送付します。
  皆様有り難うございました。
  gnuplotって面白いですね。
• Mook2008/06/05 21:52:36
  今回はマクロで書きましたが、手作業でできる手順をマクロにしただけですので、
  表ができれば、後はグラフの作成ボタンを押すだけで、今回の結果と同じになります。
  
  でもやはり、グラフは専用のアプリケーションでやったほうがきれいにできそうですね。
  
  

質問の情報

登録日時

2008-06-05 03:36:44

終了日時

2008-06-05 20:48:44

回答条件

回答にURL必須 1人2回まで

この質問のカテゴリ

• コンピュータ
• 科学・統計資料

この質問に含まれるキーワード

人気の質問

• この英文どうなってますか？ The boy seemed to be waiting for his mother to come.（その少年は母親が来…匿名質問者112日前
• 不定詞の形容詞的用法らしいんですけど、副詞的用法と書いたら✖️ですか？ I want to buy a book to read …匿名質問者112日前
• 日向陽葵という○V女優が好きなのですが似ている○V女優はいますか匿名質問者20時間前
• 共通テストで使われるプログラミング言語は共通テスト特有の言語と聞いたのですが、受験勉強で身につけたプ…匿名質問者120時間前
• 内定してほしいポケモンはいますか？匿名質問者20時間前

メニュー

• 質問一覧
• カテゴリ一覧
• 無料ユーザー登録
• 人力検索はてなトップ

PC版

• お問い合わせ
• ヘルプ
• お知らせ
• はてなトップ