Uma pergunta provavelmente muito básica sobre ANOVA multifatorial. Suponha um design bidirecional no qual testamos os efeitos principais A, B e a interação A: B. Ao testar o efeito principal de A com SS tipo I, o efeito SS é calculado como a diferença , em que é a soma dos quadrados dos erros...
Eu tenho tentado aprender quais distribuições usar nos GLMs e estou um pouco confuso sobre quando usar a distribuição normal. Em uma parte do meu livro, ele diz que uma distribuição normal pode ser boa para modelar as notas dos exames. Na próxima parte, ele pergunta qual distribuição seria...
Os gráficos abaixo são gráficos de dispersão residual de um teste de regressão para os quais as suposições "normalidade", "homoscedasticidade" e "independência" já foram atendidas com certeza! Para testar a suposição "linearidade" , embora, olhando os gráficos, possa-se adivinhar que o...
Se uma variável fator (por exemplo, gênero com níveis M e F) for usada na fórmula glm, as variáveis fictícias serão criadas e poderão ser encontradas no resumo do modelo glm juntamente com seus coeficientes associados (por exemplo, gêneroM) Se, em vez de confiar em R para dividir o fator dessa...
Estou procurando conselhos e comentários que tratem da análise de proporções e taxas. No campo em que trabalho, a análise de proporções em particular é generalizada, mas li alguns artigos que sugerem que isso pode ser problemático, estou pensando em: Kronmal, Richard A. 1993. Correlação espúria...
Estive discutindo com vários membros do laboratório sobre esse assunto e fomos a várias fontes, mas ainda não temos a resposta: Quando dizemos que um GLM tem uma família de poisson , digamos que estamos falando sobre a distribuição dos resíduos ou a variável de resposta? Pontos de discórdia Ao...
Até onde eu sei, a diferença entre o modelo logístico e o modelo de resposta fracionária (frm) é que a variável dependente (Y) na qual frm é [0,1], mas a logística é {0, 1}. Além disso, a frm usa o estimador de quase-verossimilhança para determinar seus parâmetros. Normalmente, podemos usar...
Minhas previsões provenientes de um modelo de regressão logística (glm em R) não são delimitadas entre 0 e 1 como eu esperava. Meu entendimento da regressão logística é que seus parâmetros de entrada e modelo são combinados linearmente e a resposta é transformada em probabilidade usando a função de...
Sei que essa pode ser uma questão potencialmente ampla, mas fiquei imaginando se existem suposições generalizáveis que indicam o uso de um GAM (modelo aditivo generalizado) sobre um GLM (modelo linear generalizado)? Alguém me disse recentemente que os GAMs devem ser usados apenas quando eu...
Eu sei como executar uma regressão linear em um conjunto de pontos. Ou seja, eu sei como ajustar um polinômio de minha escolha a um determinado conjunto de dados (no sentido LSE). No entanto, o que não sei é como forçar minha solução a passar por alguns pontos específicos de minha escolha. Já vi...
Aqui está o meu contexto para esta pergunta: Pelo que sei, não podemos executar uma regressão ordinária de mínimos quadrados em R ao usar dados ponderados e o surveypacote. Aqui, temos que usar svyglm(), que em vez disso executa um modelo linear generalizado (que pode ser a mesma coisa? Estou...
Vamos considerar este conjunto de dados hipotéticos: set.seed(12345) num.subjects <- 10 dose <- rep(c(1,10,50,100), num.subjects) subject <- rep(1:num.subjects, each=4) group <- rep(1:2, each=num.subjects/2*4) response <- dose*dose/10 * group + rnorm(length(dose), 50, 30) df...
Atualmente, tenho um problema ao entender a sintaxe do R para ajustar um GLM usando a distribuição Gamma. Eu tenho um conjunto de dados, em que cada linha contém 3 co-variáveis ( X1,X2,X3X1,X2,X3X_1, X_2, X_3 ), uma variável de resposta ( YYY ) e um parâmetro de forma ( KKK ). Quero modelar a...
Eu tenho um conjunto de valores xxx e yyy que são teoricamente relacionados exponencialmente: y=axby=axby = ax^b Uma maneira de obter os coeficientes é aplicando logaritmos naturais em ambos os lados e ajustando um modelo linear: > fit <- lm(log(y)~log(x)) > a <-
Não tenho formação em matemática, mas entendo como o Perceptron simples funciona e acho que entendo o conceito de um hiperplano (imagino-o geometricamente como um plano no espaço 3D que separa duas nuvens de pontos, assim como uma linha separa nuvens de dois pontos no espaço 2D). Mas não entendo...
Eu conduzi um glm.nb por glm1<-glm.nb(x~factor(group)) sendo o grupo uma categoria c e x uma variável métrica. Quando tento obter o resumo dos resultados, obtenho resultados ligeiramente diferentes, dependendo se uso summary()ou nãosummary.glm . summary(glm1)me dá ... Coefficients:...
Uma pergunta embaraçosamente simples - mas parece que nunca foi feita no Cross Validated antes: Qual é a definição de um modelo de regressão? Também uma pergunta de suporte, O que não é um modelo de regressão? No que diz respeito a este último, estou interessado em exemplos complicados em...
Eu costumava diagnosticar meus dados multivariados usando o PCA (dados omicos com centenas de milhares de variáveis e dezenas ou centenas de amostras). Os dados geralmente vêm de experimentos com várias variáveis independentes categóricas que definem alguns grupos, e muitas vezes tenho que...
Eu li algumas explicações sobre as propriedades de modelos lineares versus não lineares, mas ainda não tenho certeza se um modelo disponível é linear ou não linear. Por exemplo, o seguinte modelo é linear ou não linear? yt=β0+β1B(L;θ)Xt+εtyt=β0+β1B(L;θ)Xt+εty_t=\beta_0 +...
Estou executando um modelo linear generalizado, no qual tenho que especificar uma família diferente da normal. Qual é a distribuição esperada de resíduos? Por exemplo, os resíduos devem ser distribuídos