Eu quero fazer uma regressão logística ordinal em R sem a suposição de probabilidades de proporcionalidade. Eu sei que isso pode ser feito diretamente usando a vglm()função R, definindo parallel=FALSE. Mas meu problema é como corrigir um conjunto específico de coeficientes nessa configuração de...
Pedimos a 60 pessoas que listassem o maior número possível de franquias de restaurantes em Atlanta. A lista geral incluía mais de 70 restaurantes, mas eliminamos os mencionados por menos de 10% das pessoas, deixando-nos 45. Para esses 45, calculamos a proporção de informantes que listaram a...
Estou tentando aplicar uma série temporal aos dados amostrados trimestralmente (biomassa animal) durante um período de 10 anos com 3 repetições por trimestre. Então, 40 datas, mas 120 observações totais. Eu li até SARIMA'a na análise de séries temporais de Shumway e Stoffer e suas aplicações, bem...
Quero visualizar os resultados de um agrupamento (produzido com protoclust{protoclust}) criando gráficos de dispersão para cada par de variáveis usadas para classificar meus dados, colorindo por classes e sobrepondo as elipses para o intervalo de confiança de 95% para cada uma das classes (para...
Olá, estou fazendo um curso de pós-graduação em estatística e abordamos estatísticas de testes e outros conceitos. No entanto, muitas vezes sou capaz de aplicar as fórmulas e desenvolver uma espécie de intuição sobre como as coisas funcionam, mas muitas vezes fico com a sensação de que talvez se...
Gostaria de testar em que regressão se encaixa melhor nos meus dados. Minha variável dependente é uma contagem e possui muitos zeros. E eu precisaria de ajuda para determinar qual modelo e família usar (poisson ou quasipoisson ou regressão de poisson inflada a zero) e como testar as...
Eu gostaria de usar lme4para ajustar uma regressão de efeitos mistos e multcompcalcular as comparações aos pares. Eu tenho um conjunto de dados complexo com vários preditores contínuos e categóricos, mas minha pergunta pode ser demonstrada usando o ChickWeightconjunto de dados interno como...
A documentação do R para qualquer um deles não esclarece muita coisa. Tudo o que posso obter neste link é que o uso de qualquer um deles deve estar bem. O que eu não entendo é por que eles não são iguais. Fato: A função de regressão gradual em R é step()usada extractAIC(). Curiosamente, a...
Para analisar a contagem de pássaros inflada a zero, eu gostaria de aplicar modelos de contagem inflada a zero usando o pacote R pscl . No entanto, observando o exemplo fornecido na documentação para uma das principais funções ( ? Zeroinfl ), começo a duvidar de qual é a real vantagem desses...
A implementação randomForest não permite amostragem além do número de observações, mesmo quando amostragem com substituição. Por que é isso? Funciona bem: rf <- randomForest(Species ~ ., iris, sampsize=c(1, 1, 1), replace=TRUE) rf <- randomForest(Species ~ ., iris, sampsize=3,...
Em relação à ponderação do escore de propensão (IPTW) ao fazer a modelagem de risco proporcional de Cox dos dados de sobrevivência de tempo para evento: Tenho dados prospectivos de registro nos quais estamos interessados em analisar o efeito do tratamento de um medicamento que, na maioria dos...
No artigo seminal "Rao-Blackwellised Particle Filtering for Dynamic Bayesian Networks" de A. Doucet et. al. é proposto um filtro seqüencial de Monte Carlo (filtro de partículas), que utiliza uma subestrutura linear em um processo markov x k = ( x L k , x N k ) . Pela marginalização dessa estrutura...
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....
Eu tenho um conjunto de dados que contém o número de ações executadas por indivíduos ao longo de sete dias. A ação específica não deve ser relevante para esta pergunta. Aqui estão algumas estatísticas descritivas para o conjunto de dados: AlcanceSignificarVariaçãoNúmero de observações0 -...
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
Considere o seguinte código e saída: par(mfrow=c(3,2)) # generate random data from weibull distribution x = rweibull(20, 8, 2) # Quantile-Quantile Plot for different distributions qqPlot(x, "log-normal") qqPlot(x, "normal") qqPlot(x, "exponential", DB = TRUE) qqPlot(x, "cauchy") qqPlot(x,...
Estou tendo um problema ao usar o 2l.normmétodo de imputação em vários níveis no Windows mice. Infelizmente, não posso postar um exemplo reproduzível por causa do tamanho dos meus dados - quando reduzo o tamanho, o problema desaparece. Para uma variável específica, miceproduz os seguintes erros e...
Eu tenho um conjunto de dados de 20 anos de uma contagem anual de abundância de espécies para um conjunto de polígonos (~ 200 polígonos de forma irregular e contínua). Eu tenho usado a análise de regressão para inferir tendências (mudança na contagem por ano) para cada polígono, bem como agregações...
Me deparei com esta imagem em um post aqui . Fiquei desapontado que a leitura da declaração não provocou a mesma expressão facial para mim como para esse cara. Então, o que se quer dizer com a afirmação de que a hipótese nula é como os freqüentadores expressam um prior não informativo? É mesmo...
Estou tentando obter o prior de Jeffreys para uma distribuição binomial negativa. Não consigo ver onde errei, por isso, se alguém puder ajudar a apontar isso, isso seria apreciado. Ok, então a situação é a seguinte: devo comparar as distribuições anteriores obtidas usando um binômio e um binômio...