Mathematica 9 is now available
THIS IS DOCUMENTATION FOR AN OBSOLETE PRODUCT.
SEE THE DOCUMENTATION CENTER FOR THE LATEST INFORMATION.
Mathematica > データの操作 > 数値データ > 統計モデル分析 >
Mathematica > データの操作 > 統計 > 統計モデル分析 >
Mathematica > 数学とアルゴリズム > 統計 > 統計モデル分析 >

GeneralizedLinearModelFit

GeneralizedLinearModelFit[{y1, y2, ...}, {f1, f2, ...}, x]
連続する変数 x の値1, 2, ...について yi をフィットする g-1(Beta0+Beta1 f1+Beta2 f2+...)の形式の一般化された線形モデルを構築する.
GeneralizedLinearModelFit[{{x11, x12, ..., y1}, {x21, x22, ..., y2}, ...}, {f1, f2, ...}, {x1, x2, ...}]
fi が変数 xk に依存する g-1(Beta0+Beta1 f1+Beta2 f2+...)という形式の一般化された線形モデルを構築する.
GeneralizedLinearModelFit[{m, v}]
計画行列 m と応答ベクトル v から一般化された線形モデルを構築する.
  • GeneralizedLinearModelFitは,自身が構築した一般化された線形モデルを表す記号的なFittedModelオブジェクトを返す.モデルの特性と診断は model["property"]で得ることができる.
  • 特定の点 x1, ...におけるGeneralizedLinearModelFitからの最もよくフィットした関数の値は model[x1, ...]で得ることができる.
  • {{x_(11),x_(12),... ,y_(1)},{x_(21),x_(22),... ,y_(2)},...}の形式のデータでは,座標 xi1, xi2, ... の数は変数 xi の数と同じでなければならない.
  • {y1, y2, ...}の形式のデータは{{1, y1}, {2, y2}, ...}の形式のデータに等しい.
  • GeneralizedLinearModelFitは,もとの yi が平均が の指数型分布に従う独立の観察であり関数 g が可逆なリンク関数であるという仮定の下に,という形式の一般化された非線形モデルを作成する.
AccuracyGoalAutomatic目標確度
ConfidenceLevel95/100パラメータと予測の信頼性のレベル
CovarianceEstimatorFunction"ExpectedInformation"パラメータ共分散行列のための推定法
DispersionEstimatorFunctionAutomatic分散パラメータを推定する関数
ExponentialFamilyAutomaticy の指数型分布
IncludeConstantBasisTrue定数基底関数を入れるかどうか
LinearOffsetFunctionNone線形予測量における既知のオフセット
LinkFunctionAutomaticモデルのリンク関数
MaxIterationsAutomatic使用する最大反復回数
NominalVariablesNone名義的とみなされる関数
PrecisionGoalAutomatic目標とする精度
WeightsAutomaticデータ要素の重み
WorkingPrecisionAutomatic内部計算で使われる精度
  • LinearOffsetFunction->h の設定では,g-1(h[x1, x2, ...]+Beta0+Beta1 f1+Beta2 f2+...)の形式のモデルがフィットされる.
  • ConfidenceLevel->p のとき,確率 p の信頼区間はパラメータと予測区間に基づいて計算される.
  • DispersionEstimatorFunction->f と設定すると,共通分散は で推定される.ただし,y={y1, y2, ...}は観察リスト,は予測値のリスト,w={w1, w2, ...}は測度 yiの重みのリストである.
  • ExponentialFamilyの可能な設定値は,"Gaussian""Binomial""Poisson""Gamma""InverseGaussian""QuasiLikelihood"である.
  • model["property"]を使って得られるデータとフィットされた関数に関する特性:
"BestFit"フィットされた関数
"BestFitParameters"パラメータの推定
"Data"入力情報あるいは計画行列と反応ベクトル
"DesignMatrix"モデルのための計画行列
"Function"最もよくフィットした純関数
"LinearPredictor"フィットされた線形結合
"Response"入力データの応答値
  • 分散とモデルからの逸脱に関連する特性:
"Deviances"尤離度
"DevianceTable"尤離度表
"DevianceTableDegreesOfFreedom"表からの自由度の差
"DevianceTableDeviances"表からの尤離度の差
"DevianceTableEntries"表からのフォーマットされていない値の配列
"DevianceTableResidualDegreesOfFreedom"表からの自由度残差
"DevianceTableResidualDeviances"表からの残差の尤離度
"EstimatedDispersion"予測分散パラメータ
"NullDeviance"ヌルモデルの尤離度
"NullDegreesOfFreedom"ヌルモデルの自由度
"ResidualDeviance"モデル尤離度とヌル尤離度の差
"ResidualDegreesOfFreedom"モデル自由度とヌル自由度の差
  • 残差のタイプ:
"AnscombeResiduals"アンスコム(Anscombe)残差
"DevianceResiduals"尤離度残差
"FitResiduals"実際の応答と予想された応答の差
"LikelihoodResiduals"尤度残差
"PearsonResiduals"ピアソン(Pearson)残差
"StandardizedDevianceResiduals"標準化された尤離度残差
"StandardizedPearsonResiduals"標準化されたピアソン残差
"WorkingResiduals"作業残差
  • パラメータ推定の特性と診断:
"CorrelationMatrix"漸近的なパラメータ相関行列
"CovarianceMatrix"漸近的なパラメータ共分散行列
"ParameterConfidenceIntervals"パラメータ推定の偏り
"ParameterConfidenceIntervalTable"フィットされたパラメータの信頼区間情報の表
"ParameterConfidenceIntervalTableEntries"表からのフォーマットされていない値の配列
"ParameterConfidenceRegion"楕円体パラメータ信頼領域
"ParameterTableEntries"表からのフォーマットされていない値の配列
"ParameterErrors"予測パラメータの標準誤差
"ParameterPValues"パラメータ z 統計の p
"ParameterTable"フィットされたパラメータ情報の表
"ParameterZStatistics"パラメータ推定の z 統計
  • 影響尺度関連特性:
"CookDistances"クック(Cook)の距離のリスト
"HatDiagonal"ハット行列の対角要素
  • 予測値特性:
"PredictedResponse"データのフィットされた値
  • 適合度尺度の特性:
"AdjustedLikelihoodRatioIndex"Ben-AkivaとLermanの修正尤度比指数
"AIC"赤池情報量基準
"BIC"ベイズ(Bayes)情報量基準
"CoxSnellPseudoRSquared"CoxとSnellの擬似 R^2
"CraggUhlerPseudoRSquared"CraggとUhlerの擬似 R^2
"EfronPseudoRSquared"Efronの擬似 R^2
"LikelihoodRatioIndex"McFaddenの尤度比指数
"LikelihoodRatioStatistic"尤度比
"LogLikelihood"フィットされたモデルのための対数尤度
"PearsonChiSquare"ピアソンの chi^2 統計
  • GeneralizedLinearModelFit[m, v]では,計画行列 m は,{{f1, f2, ...}, {f1, f2, ...}, ...}の形式のデータ点における基底関数 fi の値から形成される.応答ベクトル v は応答のリスト{y1, y2, ...}である.
  • 計画行列 m と応答ベクトル v について,モデルは である.ただし,Beta は推定されるパラメータのベクトルである.
  • 計画行列が使われるとき,基底関数 fiGeneralizedLinearModelFit[{m, v}, {f1, f2, ...}]の形式で指定できる.
データセットを定義する:
In[1]:=
Click for copyable input
対数線形ポアソン(Poisson)モデルをデータにフィットする:
In[2]:=
Click for copyable input
Out[2]=
モデルの関数形を見る:
In[3]:=
Click for copyable input
Out[3]=
ある点でモデルを評価する:
In[4]:=
Click for copyable input
Out[4]=
データ点とモデルをプロットする:
In[5]:=
Click for copyable input
Out[5]=
モデルの尤離度残差を計算し,プロットする:
In[6]:=
Click for copyable input
Out[6]=
In[7]:=
Click for copyable input
Out[7]=
バージョン 7 の新機能
Ask a question about this page  |  Suggest an improvement  |  Leave a message for the team