devianceTest

越轨行为的分析广义线性回归模型

描述

TBL= devianceTest(MDL返回广义线性回归模型的偏差表中的分析MDL。桌子TBL给出了测试的确定是否该模型的结果MDL配合显著优于一常数模型。

例子

全部收缩

执行对广义线性回归模型中的偏差测试。

使用泊松随机数具有两个潜在的预测因子生成的样本数据X(:,1)X(:,2)

RNG('默认'%用于重现rndvars = randn(100,2);X = [2个+ rndvars(:,1),rndvars(:,2)];亩= EXP(1 + X * [1; 2]);Y = poissrnd(亩);

创建泊松数据的广义线性回归模型。

MDL = fitglm(X,Y,'Y〜X1 + X2''分配'“泊松
MDL =广义线性回归模型:日志(y)的〜1 + X1 + X2分布=泊松估计系数:估计SE TSTAT p值________ _________ ______(截距)1.0405 0.022122 47.034 0 X1 0.9968 0.003362 296.49 0×2 1.987 0.0063433 313.24 0 100观察,97个错误自由度分散体:1驰^ 2统计量与常数模型:2.95e + 05,p值= 0

测试是否从一个恒定的统计显著方式的模型不同。

TBL = devianceTest(MDL)
TBL =2×4表越轨DFE chi2Stat p值__________ ___ __________ ______日志(Y)〜1个2.9544e + 05 99日志(y)的〜1 + X1 + X2 107.4 97 2.9533e + 05 0

p- 值表明,该模型显著从常数不同。请注意,模型显示MDL包括表中的第二行中示出的统计信息。

输入参数

全部收缩

广义线性回归模型,指定为GeneralizedLinearModel使用对象的创建fitglm要么stepwiseglmCompactGeneralizedLinearModel使用对象的创建紧凑

输出参数

全部收缩

越轨汇总统计分析,返回表。

TBL包含两个恒定的模型和模型偏差的统计分析MDL。该表包括这些列每个模型。

描述
越轨

偏差是(相应的模型的loglikelihoods之间两倍的差异MDL或恒定)和饱和模式。欲了解更多信息,请参阅越轨

DFE

自由的误差(残差)的度,等于ñ-p,其中ñ是观测值的数量,并且p是估计系数的数目

chi2Stat

Ft-统计或卡方统计,这取决于分散体是否被估计(Ft-统计)或未被(卡方统计)

  • Ft-统计是常数模型的偏差和完整模型的偏差,通过所估计的分散划分之间的差。

  • 卡方统计是常数模型的偏差和完整模型的偏差之间的差异。

p值

p卡方统计与: - 值与测试相关联p- 1自由度,或F与t-统计p- 1分子自由度和DFE分母自由度,其中p是估计系数的数目

更多关于

全部收缩

越轨

偏差模型的中号1是两次模型的对数似然之间的差中号1和饱和模型中号小号。饱和模型是,你可以进行参数估计的最大数量的模型。

举例来说,如果你有ñ意见(ÿ一世一世= 1,2,...,ñ)具有用于潜在地不同的值X一世Ťβ,那么你可以定义与饱和模型ñ参数。让大号bÿ)表示似然函数的最大值与所述参数的模型b。然后,该模型的偏差中号1

- 2 日志 大号 b 1 ÿ - 日志 大号 b 小号 ÿ

哪里b1b小号包含估计的参数模型中号1和分别饱和模型。该偏差与卡方分布ñ-p自由度,其中ñ是的参数在饱和模型的数目和p是参数模型中的数中号1

假设你有两个不同的广义线性回归模型中号1中号2中号1有条款的一个子集中号2。您可以通过比较deviances评估模型的拟合d1d2这两款车型。该deviances的差别

d = d 2 - d 1 = - 2 日志 大号 b 2 ÿ - 日志 大号 b 小号 ÿ + 2 日志 大号 b 1 ÿ - 日志 大号 b 小号 ÿ = - 2 日志 大号 b 2 ÿ - 日志 大号 b 1 ÿ

渐近,差异d具有自由度的卡方分布v等于差在在估计的参数数量中号1中号2。你可以得到p通过使用-VALUE此测试1 - chi2cdf(d,v)的

通常情况下,你检查d使用模型中号2与常数项,没有预测。因此,d与卡方分布p- 1自由程度。如果分散估计,差值除以估计色散具有F与分配p- 1分子自由度和ñ-p分母自由度。

介绍了在R2012a