glm
JAGS У меня есть логистическая модель только для перехвата в JAGS, определенная следующим образом: model{ for(i in 1:Ny){ y[i] ~ dbern(mu[s[i]]) } for(j in 1:Ns){ mu[j] <- ilogit(b0[j]) b0[j] ~ dnorm(0, sigma) } sigma ~ dunif(0, 100) } Когда я строю заднее распределение b0 коллапса по вс...
Я занимаюсь логистической регрессией в R. может кто-нибудь объяснить, в чем разница между этими двумя линиями? 1. glm(Response ~ Temperature, data=temp, family = binomial(link="logit")) 2. glm(cbind(Response, n - Response) ~ Temperature, data=temp, family =binomial, ...
Мой вопрос имеет отношение к оценке темпов роста численности населения в модель роста Мальтуса. В качестве примера игрушки рассмотрим набор данных игрушек df: structure(list(x= c(0L, 24L, 48L, 72L, 96L, 120L, 144L, 168L ), y = c(10000, 18744.0760659189, 35134.0387564953, 65855.509495469, 123440.067934292, ...
При выполнении логит-регрессии в R можно получить коэффициенты после того, как алгоритм оптимизации сошелся (или не сошелся) с функцией coefficients(): library(MASS) data(menarche) glm.out = glm(cbind(Menarche, Total-Menarche) ~ Age, family=binomial(logit), data=menarche) coefficients(glm.out...
Я пытаюсь создать модель, используя пакет MCMCglmm в R. Данные структурированы следующим образом, где dyad, focal, other-все случайные эффекты, predict1 - 2-переменные предиктора, а response 1-5-переменные результата, которые фиксируют # наблюдаемых поведений различных подтипов: dyad focal other r pre...
У меня есть файл rda, в котором есть модель логистической регрессии. Если я только хочу знать, какие типы данных для независимых переменных используются в glm, Как я могу их проверить? Например, если мой Y~ X1+X2+X3 как я узнаю, какие типы данныхX1, X2 а X3 есть? На данный момент у меня нет фрейма данных для ...
Пакет "prais" содержит функцию prais.winsten для запуска регрессионных моделей с оценкой Prais Winsten. Однако ни stargazer, ни broom пакеты, по-видимому, не работают с результатами функции prais.winsten. После подгонки модели вида pw<- prais.winsten(speed ~ dist, cars) Когда я пытаюсь сообщить резуль...
Может ли кто-нибудь помочь мне воспроизвести эти расчеты относительного риска (и его доверительный интервал) в R? Аналогичная процедура, используемая в Stata, описаназдесь . Может ли кто-нибудь показать мне, как это сделать в R (мои данные имеют кластеры и слои, но я взял более простой пример)? Я уже пробова...
Предположим, у меня есть переменная ответа и данные, содержащие три ковариаты (как пример игрушки): y = c(1,4,6) d = data.frame(x1 = c(4,-1,3), x2 = c(3,9,8), x3 = c(4,-4,-2)) Я хочу подогнать линейную регрессию к данным: fit = lm(y ~ d$x1 + d$x2 + d$y2) есть ли способ, чтобы написать формулу, так что я...
в R предсказать.lm вычисляет Прогнозы на основе результатов линейной регрессии, а также предлагает вычислить доверительные интервалы для этих прогнозов. Согласно руководству, эти интервалы основаны на дисперсии ошибок подгонки, но не на интервалах ошибок коэффициента. с другой стороны предсказать.glm, котор...