glmfit
Обобщенная линейная регрессия модели
Синтаксис
b = glmfit(X,y,distr)
b = glmfit(X,y,distr,param1,val1,param2,val2,...)
[b,dev] = glmfit(...)
[b,dev,stats] = glmfit(...)
Описание
b = glmfit(X,y, возвращается (p + 1)-by-1 векторный distr)b из коэффициента оценивает для обобщенной линейной регрессии ответов в y на предикторах в X, использование распределения distrX n-by-p матрица p предикторов при каждом из n наблюдений. distr может быть любое следующее: 'binomial'\Gamma, 'inverse gaussian', 'normal' (значение по умолчанию), и 'poisson'.
В большинстве случаев, y n-1 вектор наблюдаемых ответов. Для биномиального распределения, y может быть бинарный вектор, указывающий на успешность или неуспешность при каждом наблюдении или две матрицы столбца с первым столбцом, указывающим на количество успехов для каждого наблюдения и второго столбца, указывающего на количество испытаний за каждое наблюдение.
Этот синтаксис использует каноническую ссылку (см. ниже) связать распределение с предикторами.
Примечание
По умолчанию, glmfit добавляет первый столбец 1 с к X, соответствие постоянному члену в модели. Не вводите столбец 1 с непосредственно в X. Можно изменить поведение по умолчанию glmfit использование 'constant' параметр, ниже.
glmfit обработки NaNs в любом X или y как отсутствующие значения, и игнорирует их.
b = glmfit(X,y, дополнительно позволяет вам задавать дополнительное название параметра / пары значения, чтобы управлять подгонкой модели. Приемлемые параметры следующие.distr,param1,val1,param2,val2,...)
| Параметр | Значение | Описание |
|---|---|---|
'link' |
| µ = X b |
| регистрируйте (µ) = X b | |
| журнал (µ / (1 – µ)) = X b | |
'probit' |
| |
'comploglog' | журнал (-журнал (1 – µ)) = X b | |
'reciprocal', значение по умолчанию для распределения 'gamma' | 1/µ = X b | |
| журнал (-журнал (µ)) = X b | |
| µ p = X b | |
массив ячеек формы | Пользовательски заданная функция ссылки. Необходимо обеспечить
| |
массив структур, имеющий эти поля:
Значение каждого поля является вектором символов, соответствующим функции, которая находится на пути или указателе на функцию (созданное использование | Пользовательски заданная функция ссылки, ее производная и ее инверсия. | |
'estdisp' | 'on' |
|
|
| |
'offset' | Вектор |
|
'weights' | Вектор предшествующих весов, таких как инверсии относительного отклонения каждого наблюдения | |
'constant' |
|
|
'off' |
|
[b,dev] = glmfit(...) возвращает dev, отклонение подгонки в векторе решения. Отклонение является обобщением остаточной суммы квадратов. Возможно выполнить анализ отклонения, чтобы сравнить несколько моделей, каждый подмножество другого, и протестировать, значительно лучше ли модель с большим количеством условий, чем модель с меньшим количеством условий.
[b,dev,stats] = glmfit(...) возвращает dev и stats.
stats структура со следующими полями:
beta— Коэффициент оцениваетbdfe— Степени свободы для ошибкиsfit— Предполагаемый дисперсионный параметрs— Теоретический или оцененный дисперсионный параметрestdisp— 0, когда'estdisp'значением аргумента пары "имя-значение" является'off'и 1, когда'estdisp'значением аргумента пары "имя-значение" является'on'.covb— Предполагаемая ковариационная матрица для Bse— Вектор стандартных погрешностей коэффициента оцениваетbcoeffcorr— Корреляционная матрица дляbt— статистика t дляbp— p-значения дляbresid— Вектор остаточных значенийresidp— Вектор остаточных значений Пирсонаresidd— Вектор остаточных значений отклоненияresida— Вектор остаточных значений Anscombe
Если вы оцениваете дисперсионный параметр для бинома или распределения Пуассона, то stats.s установлен равный stats.sfit. Кроме того, элементы stats.se отличайтесь факторным stats.s от их теоретических значений.
Примеры
Подбирайте обобщенную линейную модель со ссылкой пробита
Введите выборочные данные.
x = [2100 2300 2500 2700 2900 3100 ...
3300 3500 3700 3900 4100 4300]';
n = [48 42 31 34 31 21 23 23 21 16 17 21]';
y = [1 2 0 3 8 8 14 17 19 15 17 21]';Каждый y значение является количеством успехов в соответствующем количестве испытаний в n, и x содержит значения переменного предиктора.
Подбирайте модель регрессии пробита для y на x.
b = glmfit(x,[y n],'binomial','link','probit');
Вычислите предполагаемое количество успехов. Постройте процент наблюдаемый и оцененный успех процента по сравнению с x значения.
yfit = glmval(b,x,'probit','size',n); plot(x, y./n,'o',x,yfit./n,'-','LineWidth',2)
Используйте пользовательски заданную функцию ссылки
Загрузите выборочные данные.
load fisheririsВектор-столбец, species, состоит из ирисовых цветов трех различных разновидностей, setosa, versicolor, virginica. Двойной матричный meas состоит из четырех типов измерений на цветах, длине и ширине чашелистиков и лепестков в сантиметрах, соответственно.
Задайте переменные отклика и переменные предикторы.
X = meas(51:end,:);
y = strcmp('versicolor',species(51:end));Задайте три указателя на функцию, созданное использование @, это задает ссылку, производную ссылки и обратной ссылки для функции ссылки логита. Сохраните их в массиве ячеек.
link = @(mu) log(mu ./ (1-mu));
derlink = @(mu) 1 ./ (mu .* (1-mu));
invlink = @(resp) 1 ./ (1 + exp(-resp));
F = {link, derlink, invlink};Соответствуйте логистической регрессии с помощью glmfit с функцией ссылки, которую вы задали.
b = glmfit(X,y,'binomial','link',F)
b = 5×1
42.6378
2.4652
6.6809
-9.4294
-18.2861
Подбирайте обобщенную линейную модель при помощи logit соединитесь функционируют и сравнивают результаты.
b = glmfit(X,y,'binomial','link','logit')
b = 5×1
42.6378
2.4652
6.6809
-9.4294
-18.2861
Ссылки
[1] Добсон, A. J. Введение в обобщенные линейные модели. Нью-Йорк: Chapman & Hall, 1990.
[2] Маккуллаг, P. и Дж. А. Нелдер. Обобщенные линейные модели. Нью-Йорк: Chapman & Hall, 1990.
[3] Collett, D. Моделирование двоичных данных. Нью-Йорк: Chapman & Hall, 2002.
Смотрите также
GeneralizedLinearModel | fitglm | glmval | regress | regstats | stepwiseglm