Temos dados com resultado binário e algumas covariáveis. Eu usei regressão logística para modelar os dados. Apenas uma análise simples, nada de extraordinário. O resultado final deve ser uma curva dose-resposta, onde mostramos como a probabilidade muda para uma covariável específica. Algo...
Y=μ+σGY=μ+σGY=\mu+\sigma Gμμ\muWWWGGGRnRn\mathbb{R}^nFFFH0:{μ∈U}H0:{μ∈U}H_0\colon\{\mu \in U\}U⊂WU⊂WU \subset W Como podemos saber que essa estatística fornece o teste mais poderoso paraH0(talvez depois de descartar casos particulares incomuns)? Isso não deriva do teorema de Neyman-Pearson, porque...
Eu tenho desenvolvido um modelo de regressão logística baseado em dados retrospectivos de um banco de dados nacional de trauma de traumatismo craniano no Reino Unido. O principal resultado é a mortalidade em 30 dias (denominada medida "Sobreviver"). Outras medidas com evidências publicadas de...
Tanto quanto posso dizer, curvilíneo é definido vagamente, mas significa o mesmo que não-linear . Isso está correto? Ou curvilínea tem uma definição
Este artigo utiliza modelos lineares generalizados (distribuições binomiais e negativas de erros binomiais negativos) para analisar dados. Mas então, na seção de análise estatística dos métodos, há esta declaração: ... e segundo modelando os dados de presença usando Modelos de Regressão...
Eu tenho um conjunto de dados muito grande e faltam cerca de 5% de valores aleatórios. Essas variáveis estão correlacionadas entre si. O exemplo a seguir do conjunto de dados R é apenas um exemplo de brinquedo com dados correlatos simulados. set.seed(123) # matrix of X variable xmat <-...
Tentando aprender um pouco de Python e Sklearn, mas para o meu trabalho eu preciso executar regressões que usem distribuições de erro das famílias Poisson, Gamma e, especialmente, Tweedie. Não vejo nada na documentação sobre eles, mas eles estão em várias partes da distribuição R, então eu queria...
Estou tentando modelar as intensidades médias de parasitas que afetam um hospedeiro em R usando um modelo binomial negativo. Eu recebo 50 ou mais avisos que dizem: In dpois(y, mu, log = TRUE) : non-integer x = 251.529000 Como posso lidar com isso? Meu código fica assim: mst.nb =...
Alguém sabe qual é a fórmula para a distância de Cook? A fórmula de distância original de Cook usa resíduos estudados, mas por que R está usando std. Os resíduos de Pearson ao calcular o gráfico de distância de Cook para um GLM. Sei que resíduos estudantis não são definidos para GLMs, mas como é a...
Devo confessar que nunca ouvi falar desse termo em nenhuma das minhas aulas, graduação ou pós-graduação. O que significa uma regressão logística ser bayesiana? Estou procurando uma explicação com uma transição da logística regular para a logística bayesiana semelhante à seguinte: Esta é a equação...
Minha pergunta revela meu fraco entendimento da regressão de Poisson e GLMs em geral. Aqui estão alguns dados falsos para ilustrar minha pergunta: ### some fake data x=c(1:14) y=c(0, 1, 2, 3, 1, 4, 9, 18, 23, 31, 20, 25, 37, 45) Algumas funções personalizadas para retornar psuedo-R2: ###...
Atualmente, estou trabalhando em um projeto envolvendo GLMs (e eventualmente GAMs) de alguns dados de contagem ao longo do tempo. Normalmente eu faria isso no SAS, mas estou tentando mudar para R e tendo ... problemas. Quando encaixo um GLM para contar dados usando o seguinte: cdi_model <-...
Antecedentes Biológicos Com o tempo, algumas espécies de plantas tendem a duplicar todo o seu genoma, ganhando uma cópia adicional de cada gene. Devido à instabilidade dessa configuração, muitos desses genes são excluídos e o genoma se rearranja e se estabiliza, pronto para duplicar novamente....
Para dados de contagem que eu coletei, uso a regressão de Poisson para construir modelos. Eu faço isso usando a glmfunção em R, onde eu uso family = "poisson". Para avaliar possíveis modelos (tenho vários preditores), uso o AIC. Por enquanto, tudo bem. Agora eu quero executar a validação cruzada....
Qual é o significado t valuee Pr(>|t|)quando se utiliza a summary()função no modelo de regressão linear em R? Coefficients: Estimate Std. Error t value Pr(>|t|) (Intercept) 10.1595 1.3603 7.469 1.11e-13 *** log(var) 0.3422 0.1597 2.143 0.0322
Medimos duas variáveis, e o gráfico de dispersão parece sugerir vários modelos "lineares". Existe uma maneira de tentar destilar esses modelos? Identificar outras variáveis independentes acabou sendo difícil. Ambas as variáveis são fortemente inclinadas para a esquerda (em direção aos...
No GLM, assumindo um escalar e para a distribuição subjacente com pdf Pode ser mostrado que . Se a função de link g (\ cdot) satisfizer o seguinte, g (\ mu) = \ theta = X '\ beta em que X' \ beta é o preditor linear, então g (\ cdot) será chamada de função de link canônico para este...
Eu tenho 100.000 observações (9 variáveis indicadoras fictícias) com 1000 positivos. A regressão logística deve funcionar bem neste caso, mas a probabilidade de corte me intriga. Na literatura comum, escolhemos um ponto de corte de 50% para prever 1s e 0s. Não posso fazer isso, pois meu modelo...
O mgcvpacote para Rpossui duas funções para ajustar as interações do produto tensorial: te()e ti(). Entendo a divisão básica do trabalho entre os dois (ajustando uma interação não linear versus decompondo essa interação em efeitos principais e uma interação). O que não entendo é o porquê te(x1,...
Eu identifiquei vários lugares nos livros onde o GLM é descrito com 5 distribuições (a saber, Gama, Gaussiana, Binomial, Gaussiana Inversa e Poisson). Isso também é exemplificado na função familiar em R. Ocasionalmente, encontro referências ao GLM, onde distribuições adicionais estão incluídas (...