Não sei se isso já foi perguntado antes, mas não encontrei nada sobre isso. Minha pergunta é se alguém pode fornecer uma boa referência para aprender como obter a proporção de variância explicada por cada um dos fatores fixos e aleatórios em um modelo de efeitos
Realizei uma avaliação baseada em computador de diferentes métodos de ajuste de um tipo específico de modelo usado nas ciências paleo. Eu tinha um conjunto de treinamento amplo e, por isso, aleatoriamente (amostragem aleatória estratificada), anotei um conjunto de testes. Ajustei métodos diferentes...
Diga que você está na biblioteca do seu departamento de estatística e se depara com um livro com a figura a seguir na primeira página. Você provavelmente pensará que este é um livro sobre coisas de regressão linear. Qual seria a imagem que faria você pensar em modelos lineares...
Estou realmente impressionado com o fato de o Poisson GLM aceitar números não inteiros! Veja: Dados (conteúdo de data.txt): 1 2001 0.25 1 1 2002 0.5 1 1 2003 1 1 2 2001 0.25 1 2 2002 0.5 1 2 2003 1 1 Script R: t <- read.table("data.txt") names(t) <- c('site', 'year', 'count', 'weight') tm...
Treinei um modelo de regressão linear, usando um conjunto de variáveis / recursos. E o modelo tem um bom desempenho. No entanto, percebi que não há variável com uma boa correlação com a variável prevista. Como isso é
Suponha que você tenha uma única seção de dados em que os indivíduos estejam localizados dentro de grupos (por exemplo, estudantes nas escolas) e deseje estimar um modelo da forma em Y_i = a + B*X_ique Xseja um vetor de características de nível individual e auma constante. Nesse caso, suponha que...
Estou seguindo AW van der Vaart, estatística assintótica (1998). Ele fala de experimentos estatísticos, alegando que eles são diferentes de um modelo estatístico, mas ele também não os define. Minha pergunta: O que é (1) um experimento estatístico, (2) um modelo estatístico e (3) qual é o...
Estamos trabalhando com algumas regressões logísticas e percebemos que a probabilidade média estimada sempre é igual à proporção de uma na amostra; isto é, a média dos valores ajustados é igual à média da amostra. Alguém pode me explicar o motivo ou me dar uma referência para encontrar essa...
Atualização : desculpe por outra atualização, mas encontrei algumas soluções possíveis com polinômios fracionários e o pacote de risco concorrente com o qual preciso de ajuda. O problema Não consigo encontrar uma maneira fácil de fazer uma análise de coeficiente dependente do tempo em R. Quero...
Pergunta simples: como especificar uma distribuição lognormal no argumento da família GLM em R? Não consegui encontrar como isso pode ser alcançado. Por que lognormal (ou exponencial) não é uma opção no argumento da família? Em algum lugar nos R-Archives, li que é preciso simplesmente usar o link...
A regressão beta (ou seja, GLM com distribuição beta e geralmente a função de link de logit) é frequentemente recomendada para lidar com a resposta, também conhecida como variável dependente, recebendo valores entre 0 e 1, como frações, proporções ou probabilidades: Regressão para um resultado...
Pode-se executar uma regressão logit em R usando esse código: > library(MASS) > data(menarche) > glm.out = glm(cbind(Menarche, Total-Menarche) ~ Age, + family=binomial(logit), data=menarche) > coefficients(glm.out) (Intercept) Age -21.226395 1.631968 Parece que o algoritmo de...
Depois de ajustar um coxmodel, é possível fazer previsões e recuperar o risco relativo de novos dados. O que não entendo é como o risco relativo é calculado para um indivíduo e o que é relativo (ou seja, a média da população)? Alguma recomendação de recursos para ajudar a entender (eu não sou muito...
Eu tenho seis variáveis dependentes (dados de contagem) e várias variáveis independentes, vejo que em um MMR o script é assim: my.model <- lm(cbind(DV1,DV2,DV3,DV4,DV5,DV6) ~ IV1 + IV2 + ... + IVn) Mas, como meus dados são contados, quero usar um modelo linear generalizado e tentei o...
Eu tenho dados de densidade de peixes que estou tentando comparar entre várias técnicas de coleta diferentes, os dados têm muitos zeros e o histograma parece vaugley apropriado para uma distribuição de poisson, exceto que, como densidades, não são dados inteiros. Sou relativamente novo nos GLMs e...
Eu tenho aprendido sobre o modelo de risco proporcional de Cox. Eu tenho muita experiência montagem modelos de regressão logística e modelos para construir intuição Eu estive comparando ajuste usando coxpha partir do R "sobrevivência", com modelos de regressão logística ajuste usando glmcom...
Existe um pacote que preveja a estimativa do modelo de efeitos mistos binomiais negativos inflados a zero em R? Com isso quero dizer: Inflação zero, onde você pode especificar o modelo binomial para inflação zero, como na função zeroinfl no pacote pscl: zeroinfl (y ~ X | Z, dist = "negbin") onde...
Estou tendo problemas para entender a saída do meu lmer()modelo. É um modelo simples de uma variável de resultado (Suporte) com interceptações de estado variáveis / efeitos aleatórios de estado: mlm1 <- lmer(Support ~ (1 | State)) Os resultados de summary(mlm1)são: Linear mixed model fit by...
Considere o modelo linear simples: yy=X′ββ+ϵyy=X′ββ+ϵ\pmb{y}=X'\pmb{\beta}+\epsilon onde ϵi∼i.i.d.N(0,σ2)ϵi∼i.i.d.N(0,σ2)\epsilon_i\sim\mathrm{i.i.d.}\;\mathcal{N}(0,\sigma^2) e X∈Rn×pX∈Rn×pX\in\mathbb{R}^{n\times p} ,p≥2p≥2p\geq2 eXXX contém uma coluna de constantes. Meu questão é, dado...
Eu tenho usado o MCMCglmmpacote recentemente. Estou confuso com o que é referido na documentação como estrutura R e estrutura G. Eles parecem estar relacionados aos efeitos aleatórios - em particular especificando os parâmetros para a distribuição anterior sobre eles, mas a discussão na...