Como eu usaria a saída de uma análise de componentes principais (PCA) em um modelo linear generalizado (GLM), assumindo que o PCA é usado para seleção de variáveis para o GLM? Esclarecimento: Eu quero usar o PCA para evitar o uso de variáveis correlacionadas no GLM. No entanto, o PCA me...
Sei que essa é provavelmente uma pergunta básica ... Mas não acho a resposta. Estou ajustando um GLM com uma família Poisson e tentei dar uma olhada nas previsões, no entanto, o deslocamento parece ser levado em
Estou lendo A. Agresti (2007), Introdução à análise de dados categóricos , 2ª. edição e não tenho certeza se entendi corretamente este parágrafo (p.106, 4.2.1) (embora deva ser fácil): Na Tabela 3.1 sobre ronco e doenças cardíacas no capítulo anterior, 254 indivíduos relataram ronco todas as...
Atualmente, estou trabalhando em uma série de modelos de séries temporais de Poisson tentando estimar o efeito de uma alteração na forma como as contagens foram obtidas (alternando de um teste de diagnóstico para outro) enquanto controlo outras tendências ao longo do tempo (digamos, um aumento...
Fechadas. Esta questão está fora de tópico . No momento, não está aceitando respostas. Deseja melhorar esta pergunta? Atualize a pergunta para que ela esteja no tópico de Validação cruzada. Fechado há 2 anos . O código a seguir PredictNew <- predict...
Com dados contínuos, uma regressão linear Y=β1+β2X2+uY=β1 1+β2X2+vocêY=\beta_1+\beta_2X_2+u assume que o termo de erro está distribuído N (0, σ2σ2\sigma^2 ) 1) Assumimos que Var (Y | x) é igualmente ~ N (0, σ2σ2\sigma^2 )? 2) Qual é essa distribuição de erro na regressão logística? Quando os...
Eu tenho uma pergunta sobre a interpretação dos coeficientes de uma interação entre variável contínua e variável categórica. aqui está o meu modelo: model_glm3=glm(cog~lg_hag+race+pdg+sex+as.factor(educa)+(lg_hag:as.factor(educa)), data=base_708) Coefficients: Estimate Std. Error t value...
Na página 232 de "Um companheiro R para regressão aplicada", Fox e Weisberg observam Somente a família gaussiana tem variação constante e, em todos os outros GLMs, a variação condicional de y em depende dexx\bf{x}μ ( x )μ(x)\mu(x) Anteriormente, eles observam que a variação condicional do...
Por padrão, quando usamos uma glmfunção em R, utiliza a mínimos quadrados iterativo reponderadas método (IWLS) para encontrar a estimativa da probabilidade máxima dos parâmetros. Agora eu tenho duas perguntas. As estimativas da IWLS garantem o máximo global da função de probabilidade? Com base no...
Eu tenho esses dados: set.seed(1) predictor <- rnorm(20) set.seed(1) counts <- c(sample(1:1000, 20)) df <- data.frame(counts, predictor) Fiz uma regressão de Poisson poisson_counts <- glm(counts ~ predictor, data = df, family = "poisson") E uma regressão binomial...
Eu tenho um conjunto de dados de cerca de 5000 recursos. Para esses dados, usei o teste Chi Square para seleção de recursos; depois disso, obtive cerca de 1500 variáveis que mostraram relação de significância com a variável resposta. Agora eu preciso ajustar a regressão logística nisso. Estou...
Quando usar um modelo linear generalizado sobre um modelo linear? Eu sei que o modelo linear generalizado permite, por exemplo, que os erros tenham alguma outra distribuição além da normal, mas por que alguém se preocupa com as distribuições dos erros? Como por que diferentes distribuições de erro...
Estou ajustando vários modelos do formulário .. glm(DV ~ I(1/IV), family = Gamma(link = "log") .. e estou procurando maneiras de comparar os modelos obtidos para diferentes variáveis. Gostaria de saber se o valor alfa é de alguma utilidade prática? Para as três parcelas abaixo, os valores alfa...
Tenho dados do seguinte desenho experimental: minhas observações são contagens do número de sucessos ( K) do número correspondente de tentativas ( N), medidos para dois grupos, cada um composto por Iindivíduos, a partir de Ttratamentos, nos quais em cada combinação de fatores existem Rrepetições ....
Estou tentando ajudar um aluno de um colega. O aluno observou e contou o comportamento das aves (número de chamadas) em uma configuração experimental. O número de chamadas atribuíveis a uma ave observada específica durante cada experimento não pôde ser determinado, mas foi possível contar o número...
Estou simulando ensaios de Bernoulli com um aleatório entre grupos e, em seguida, encaixo o modelo correspondente em o pacote:logitθ∼N(logitθ0,12)logitθ∼N(logitθ0,12)\text{logit}\, \theta \sim {\cal N}(\text{logit}\, \theta_0, 1^2)lme4 library(lme4) library(data.table) I <- 30 # number of...
O seguinte modelo logístico multinível com uma variável explicativa no nível 1 (nível individual) e uma variável explicativa no nível 2 (nível do grupo): π 0 j = γ 00 + γ 01 z j + u 0 j … ( 2 ) π 1 j = γ 10 + γ 11 z j + u 1 j … ( 3
No meu campo, a maneira usual de plotar dados emparelhados é como uma série de segmentos finos de linhas inclinadas, sobrepondo-os à mediana e ao IC da mediana para os dois grupos: No entanto, esse tipo de gráfico se torna muito mais difícil de ler à medida que o número de pontos de dados se...
Eu estou usando usando um modelo linear generalizado gama (GLM) com um link de identidade. A variável independente é a compensação de um grupo específico. O resumo do statsmodels do Python está me dando um aviso sobre a função de vínculo de identidade ( "DomainWarning: a função de vínculo de...
Estou estudando a diferença entre regularização na regressão RKHS e regressão linear, mas tenho dificuldade em entender a diferença crucial entre as duas. Dados os pares entrada-saída , quero estimar uma função seguinte forma que é uma função do kernel. Os coeficientes podem ser encontrados...