分析者が独自に考えた非線形モデルを推定するには、まず、計算式を含んだ列を作成する必要があります。この計算式で、パラメータとその初期値も設定します。この方法は、操作手順が少し複雑になりますが、どのような非線形モデルでも指定できます。また、損失関数をユーザ自身で定義したり、反復計算に関する詳細なオプションを選択したりすることもできます。
この節では、モデルの計算式を含む列を作成し、「非線形回帰」プラットフォームを実行する例を述べます。使用するデータは「US Population.jmp」データテーブルです。応答変数は、米国の人口(単位は百万人)で、説明変数は年です。
1. [ヘルプ]>[サンプルデータライブラリ]を選択し、「Nonlinear Examples¥US Population.jmp」を開きます。
2. 新しい列を作成し、その列名を「モデル」とします。
3. 「モデル」列を右クリックして、[列プロパティ]>[計算式]を選択します。
計算式エディタが表示されます。
4. 列のリストの下にあるリストから、[パラメータ]を選択します。
5. [パラメータの新規作成]を選択します。
6. 「名前」に「B0」と入力します。
7. 「値」に「3.9」を入力します。これがパラメータの初期推定値です。
8. [OK]をクリックします。
9. [パラメータの新規作成]を選択します。
10. 「名前」に「B1」と入力し、「値」に「0.022」を入力します。
11. [OK]をクリックします。
12. 計算式エディタの関数、「年」列、およびパラメータを使用してモデル式を入力します。
図15.2 入力が完了したモデル式
ヒント: 指数(Exp)コマンドを見つけるには、「超越関数」の隣のグレーの三角をクリックします。
13. [OK]をクリックします。
14. [分析]>[発展的なモデル]>[非線形回帰]を選択します。
15. 「モデル」を[X,予測式列]に指定します。
16. 「人口」を[Y, 応答変数]に指定します。
17. [OK]をクリックします。
18. 設定パネルの[実行]ボタンをクリックします。
図15.3 「プロット」と「解」レポート
「解」レポートには、最終的なパラメータ推定値とその他の適合度統計量が表示されます。プロットには、あてはめたモデルが表示されます。
計算式エディタでパラメータを追加するときに、[選択された列をカテゴリに展開する]というチェックボックスを使用できます。このオプションは、カテゴリカル変数の水準ごとのパラメータを一度に作成するものです(複数のパラメータが一度に追加されます)。このオプションを選択すると、列を選択するためのダイアログボックスが開きます。列の選択が完了すると、パラメータのリストに「D_列名」という名前の新しいパラメータが表示されます。「D」の部分は、パラメータの名前です。計算式にこのパラメータを使用すると、グループ変数の水準ごとに個別のパラメータを含んだMatch式が挿入されます。