Описание тега glm

Описание тега Вопросы с тегом
По вопросам, связанным с обобщенными линейными моделями. Для математической библиотеки GLM см. [Tag:glm-math].
1 ответ

GLMER: ошибка: (maxstephalfit) разделение на шаги PIRLS не смогло уменьшить отклонение в pwrssUpdate

Я изучаю влияние различных характеристик на решение суда по конкретным преступлениям. Набор данных довольно большой (28928 наблюдений с 86 единицами уровня 2). Я смотрю на решение, заключать ли кого-либо в тюрьму или нет (= бинарная переменная резул…
24 июн '16 в 14:19
1 ответ

Подгонка и прогнозирование модели на уровне факторов в r

Я хотел бы подобрать модель на уровне фактора и использовать название подходящей модели на лету для прогнозирования новых данных на таком уровне фактора соответствия. Я ошибаюсь в предсказании в этой логике, может кто-то руководить этим, рассматрива…
02 июн '17 в 07:43
1 ответ

glmulti с дробной моделью логита

Я хочу установить модель дробного логита с помощью команды: glmfit <- glm(tr1 ~ period + male + stib+ income, family = quasibinomial(link = "logit"), data=mydata) куда tr1 переменная, которая лежит между нулем и единицей (включая некоторые нули).…
02 мар '16 в 18:27
1 ответ

В логистической регрессии, как я могу найти способ упорядочить коэффициенты по их Pr(>|z|)?

Я хотел бы упорядочить переменные из сводной таблицы функции glm по p-значению. Есть идеи как это сделать? вот мой пример фрейма данных: dat <- read.table(text = " TargetVar Var1 Var2 Var3 0 0 0 7 0 0 1 1 0 1 0 3 0 1 1 7 1 0 0 5 1 0 1 1 1 1 0 0 1…
11 май '14 в 07:23
1 ответ

Прогнозируемые значения для glmer AICcmodavg

Я знаю, что с помощью AICcmodavg можно получить прогнозируемые значения (в исходном масштабе ~ вероятность) и их SE для фиксированных эффектов, но я пытаюсь безуспешно... Может ли кто-нибудь помочь мне с этим? заранее спасибо library(lme4) (gm1 <…
22 сен '16 в 12:37
1 ответ

R: Ошибка в остановке if (any(y < 0)) ("недопустимые отрицательные значения для семейства Пуассона")

Я пытался использовать glm для оценки сильных сторон футбольных команд. # data is dataframe (structure on bottom). model &lt;- glm(Goals ~ Home + Team + Opponent, family=poisson(link=log), data=data) но получите ошибку: Error in if (any(y &lt; 0)) s…
30 апр '18 в 11:11
0 ответов

Как обработать отсутствующие особенности данных в выводе R GLM

У меня есть набор данных с рядом наблюдений, каждое из которых имеет внешние данные, добавленные к ним из разных источников. Некоторые из этих наблюдений (учетных записей) имеют данные для одного источника данных, но отсутствуют в другом источнике д…
19 окт '18 в 12:44
1 ответ

Логит биномиальная регрессия с кластерными стандартными ошибками

Я пытаюсь повторить оценку glm из stata: sysuse auto logit foreign weight mpg, cluster(rep78) Logistic regression Number of obs = 69 Wald chi2(2) = 31.57 Prob &gt; chi2 = 0.0000 Log pseudolikelihood = -22.677963 Pseudo R2 = 0.4652 (Std. Err. adjuste…
25 ноя '15 в 22:48
1 ответ

Как можно отладить / получить логи для сбоев Java-сервера SparkR?

Я одурманен No status is returned. Java SparkR backend might have failed. ошибка при подгонке glm с помощью Spark. Похоже, что задание фактически выполняется до конца на основе веб-интерфейса Spark, но в какой-то момент во время подбора модели (оно …
15 сен '16 в 08:22
0 ответов

AIC по умолчанию glm() отличается от lm ()

Я выполняю линейную регрессию на данных, которая объединяет непрерывные и категориальные предикторы. Мой ответ, а также остатки обычно распределяются. Я сравнивал результаты, заданные glm() с настройками по умолчанию, с параметрами lm(). Насколько я…
31 окт '18 в 17:41
2 ответа

Как извлечь фиксированные эффекты с помощью наблюдения?

У меня есть объект Ime, построенный на основе данных повторного измерения потребления питательных веществ (два 24-часовых периода потребления для каждого RespondentID): Male.lme2 &lt;- lmer(BoxCoxXY ~ -1 + AgeFactor + IntakeDay + (1|RespondentID), d…
30 дек '11 в 18:50
1 ответ

partykit - Как построить glmtree без перекрытия терминальных узлов?

Я хотел бы построить дерево, полученное из glmtree (пакет partykit). К сожалению, терминальные узлы перекрываются, и метки графиков больше не отображаются должным образом. Вот код, который генерирует дерево: library(partykit) library(aVirtualTwins) …
09 ноя '18 в 14:31
1 ответ

R: как передать ссылку на переменную в glm или lm?

Так скажем, у меня есть названный вектор: sorted = c(1,2,3) names(sorted) = c("A","B","C") и это будет выглядеть следующим образом: &gt; sorted A B C 1 2 3 так что это вектор с именем A,B,C и имеет значение 1,2,3 соответственно. и у меня также есть …
07 мар '16 в 01:45
0 ответов

Полная парная сводная таблица для GLM в R

Я пытаюсь запустить логистическую модель, используя 4 различных переменных предиктора: Глубина (непрерывная), DPM (2 уровня), Близость (3 уровня) и Площадь (2 уровня). Я знаю, что многие взаимодействия усложняют ситуацию, но мне было интересно, возм…
24 ноя '15 в 19:09
2 ответа

Полиномиальные данные и R's glm()

Как вы можете заставить R glm () соответствовать полиномиальным данным? Я пробовал несколько итераций 'family=AAA(link="BBB") ", но мне кажется, что тривиальные прогнозы не совпадают. Например, пожалуйста, помогите с R GLM, чтобы соответствовать пол…
12 фев '14 в 19:26
1 ответ

glm `Foregnet()` Ошибка: нет применимого метода для "предиката" применяется к объекту класса "список"

У меня проблема с управлением типами объектов, поступающими в функцию предсказания. Вот моя упрощенная функция, которая генерирует glm объект. fitOneSample &lt;- function(x,data,sampleSet) { #how big of a set are we going to analyze? Pick a number b…
02 сен '16 в 19:22
1 ответ

[r]: интерпретация результатов блеска, ретрансформирующие оценки

РЕДАКТИРОВАТЬ: В настоящее время я пишу свою магистерскую диссертацию о влиянии определенного инсектицида на колонии шмелей. Я, например, проверял, были ли поврежденные / больные появляющиеся пчелы более распространенными в колониях, которые подверг…
05 авг '15 в 08:59
1 ответ

Как игнорировать линейно коррелированные переменные, введенные с помощью кодирования эталонной ячейки

Предположим, у меня есть набор данных, содержащий две категориальные предикторные переменные (a, b) и двоичную целевую переменную (y). &gt; df &lt;- data.frame( &gt; a = factor(c("cat1","cat2","cat3","cat1","cat2")), &gt; b = factor(c("cat1","cat1",…
17 ноя '16 в 13:35
2 ответа

Регрессия в цикле, но получим "Ошибка в коэффиценте (сводка (подгонка))[2, 4]: индекс за пределами"

При запуске следующего кода в наборе данных большого размера я получаю сообщение об ошибке Ошибка в coef(итоговая (подходящая))[2, 4]: индекс за пределами Вектор beta для которого p-значения моделей логистической регрессии сохраняются, имеет длину 1…
22 июл '18 в 15:38
1 ответ

Как пройтись по GLM, тестируя разные модели

У меня есть вектор a который содержит некоторые комбинации моих переменных. Чтобы получить наилучшую возможную модель (наименьший BIC), я надеялся, что смогу перебрать glm функция - ну, это не сработало. Данные выглядят примерно так: de.dat &lt;- da…
03 окт '14 в 11:09