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...
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...
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...
A padronização de uma variável dependente dentro do grupo de identificação faz sentido? O documento de trabalho a seguir (desaceleração do desmatamento na Amazônia Legal; Preços ou políticas ?, pdf ) usa uma variável dependente padronizada para analisar o efeito da mudança de política geral no...
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...
Preciso estimar a função de risco de linha de base em um modelo Cox dependente do tempoλ0(t)λ0(t)\lambda_0(t) λ(t)=λ0(t)exp(Z(t)′β)λ(t)=λ0(t)exp(Z(t)′β)\lambda(t) = \lambda_0(t) \exp(Z(t)'\beta) Enquanto eu fazia o curso de Sobrevivência, lembro que a derivada direta da função de risco...
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:...
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...
Estou trabalhando em algo como o seguinte problema. Eu tenho vários usuários e N livros. Cada usuário cria uma classificação ordenada de todos os livros que leu (o que provavelmente é um subconjunto dos N livros), por exemplo, Livro 1> Livro 40> Livro 25. Agora, desejo transformar essas...
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 estudando modelos gráficos probabilísticos , um livro para auto-estudo. As arestas em um gráfico acíclico direcionado (DAG) representam relações causais? E se eu quiser construir uma rede bayesiana , mas não tiver certeza da direção das setas nela? Todos os dados vão me dizer são as...
Para o modelo linear , podemos ter uma boa interpretação geométrica do modelo estimado via OLS: . é a projeção de y no espaço medido por xe residual é perpendicular a esse espaço medido por x.y= x β+ey=xβ+ey=x\beta+ey^= x β^+ e^y^=xβ^+e^\hat{y}=x\hat{\beta}+\hat{e}y^y^\hat{y}e^e^\hat{e} Agora,...
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...
Estive analisando alguns dados usando a modelagem linear de efeitos mistos em R. Estou planejando fazer um pôster com os resultados e fiquei imaginando se alguém com modelos de efeitos mistos poderia sugerir quais gráficos usar para ilustrar os resultados do modelo. Eu estava pensando em gráficos...
Eu tenho lido sobre o cálculo de valores de R2R2R^2 em modelos mistos e depois de ler as perguntas frequentes do R-sig, outras postagens neste fórum (eu vincularia algumas, mas não tenho reputação suficiente) e várias outras referências que eu entendo que usando valores no contexto de modelos...
Aprendi nas estatísticas elementares que, com um modelo linear geral, para que as inferências sejam válidas, as observações devem ser independentes. Quando ocorre o agrupamento, a independência não pode mais ser mantida, levando a inferência inválida, a menos que isso seja explicado. Uma maneira de...