Доверительные интервалы оценок коэффициентов обобщенной линейной регрессионой модели
Найдите доверительные интервалы для коэффициентов подобранной обобщенной линейной регрессионой модели.
Сгенерируйте выборочные данные с помощью случайных чисел Пуассона с двумя базовыми предикторами X(:,1)
и X(:,2)
.
rng('default') % For reproducibility rndvars = randn(100,2); X = [2 + rndvars(:,1),rndvars(:,2)]; mu = exp(1 + X*[1;2]); y = poissrnd(mu);
Создайте обобщенную линейную регрессионую модель данных Пуассона.
mdl = fitglm(X,y,'y ~ x1 + x2','Distribution','poisson')
mdl = Generalized linear regression model: log(y) ~ 1 + x1 + x2 Distribution = Poisson Estimated Coefficients: Estimate SE tStat pValue ________ _________ ______ ______ (Intercept) 1.0405 0.022122 47.034 0 x1 0.9968 0.003362 296.49 0 x2 1.987 0.0063433 313.24 0 100 observations, 97 error degrees of freedom Dispersion: 1 Chi^2-statistic vs. constant model: 2.95e+05, p-value = 0
Найдите 95% (по умолчанию) доверительных интервалов для коэффициентов модели.
ci = coefCI(mdl)
ci = 3×2
0.9966 1.0844
0.9901 1.0035
1.9744 1.9996
Найдите 99% доверительные интервалы для коэффициентов.
alpha = 0.01; ci = coefCI(mdl,alpha)
ci = 3×2
0.9824 1.0986
0.9880 1.0056
1.9703 2.0036
mdl
- Обобщенная линейная регрессионая модельGeneralizedLinearModel
| объекта CompactGeneralizedLinearModel
объектОбобщенная линейная регрессионая модель, заданная как GeneralizedLinearModel
объект, созданный с помощью fitglm
или stepwiseglm
, или CompactGeneralizedLinearModel
объект, созданный с помощью compact
.
alpha
- Уровень значимостиУровень значимости для интервала доверия, заданный как числовое значение в область значений [0,1]. Уровень доверия ci
равно 100 ( 1 - alpha
)%. alpha
- вероятность того, что доверительный интервал не содержит истинного значения.
Пример: 0.01
Типы данных: single
| double
ci
- Доверительные интервалыДоверительные интервалы, возвращенные как k -на-2 числовая матрица, где k - количество коэффициентов. j строка ci
- доверительный интервал j-го коэффициента mdl
. Имя j коэффициентов сохранено в CoefficientNames
свойство mdl
.
Типы данных: single
| double
Доверительные интервалы коэффициента обеспечивают меру точности для оценок коэффициента регрессии.
100 ( 1 - α)% доверительный интервал дает область значений, в котором соответствующий коэффициент регрессии будет со 100 (1 - α)% доверительным, что означает, что 100 ( 1 - α)% интервалов, полученных в результате повторных экспериментов, будут содержать истинное значение коэффициента.
Программа находит доверительные интервалы с помощью метода Уолда. 100 * (1 - α)% доверительных интервалов для коэффициентов регрессии
где b i - оценка коэффициента, SE (b i) - стандартная ошибка оценки коэффициента, а t (1-α/2, n - p) - 100 (1 - α/2) процентиль t - распределение с n - p степенями свободы. n - количество наблюдений, а p - количество коэффициентов регрессии.
Указания и ограничения по применению:
Эта функция поддерживает объекты модели, оснащенные входными параметрами массива GPU.
Для получения дополнительной информации смотрите Запуск функций MATLAB на графическом процессоре (Parallel Computing Toolbox).
coefTest
| CompactGeneralizedLinearModel
| devianceTest
| GeneralizedLinearModel
У вас есть измененная версия этого примера. Вы хотите открыть этот пример с вашими правками?
Вы щелкнули по ссылке, которая соответствует команде MATLAB:
Выполните эту команду, введя её в командном окне MATLAB.
1. Если смысл перевода понятен, то лучше оставьте как есть и не придирайтесь к словам, синонимам и тому подобному. О вкусах не спорим.
2. Не дополняйте перевод комментариями “от себя”. В исправлении не должно появляться дополнительных смыслов и комментариев, отсутствующих в оригинале. Такие правки не получится интегрировать в алгоритме автоматического перевода.
3. Сохраняйте структуру оригинального текста - например, не разбивайте одно предложение на два.
4. Не имеет смысла однотипное исправление перевода какого-то термина во всех предложениях. Исправляйте только в одном месте. Когда Вашу правку одобрят, это исправление будет алгоритмически распространено и на другие части документации.
5. По иным вопросам, например если надо исправить заблокированное для перевода слово, обратитесь к редакторам через форму технической поддержки.