Quero ajustar um GLMM multinível com uma distribuição Poisson (com excesso de dispersão) usando R. No momento, estou usando o lme4, mas notei que recentemente a quasipoissonfamília foi removida. Já vi em outro lugar que você pode modelar a super dispersão aditiva para distribuições binomiais...
Digamos que temos as estatísticas abaixo gender mean sd n f 1.666667 0.5773503 3 m 4.500000 0.5773503 4 Como você executa um teste t de duas amostras (para ver se há uma diferença significativa entre as médias de homens e mulheres em alguma variável) usando estatísticas como essa em vez de dados...
Você pode ter dados em formato amplo ou em formato longo. Isso é muito importante, pois os métodos utilizáveis são diferentes, dependendo do formato. Eu sei que você tem que trabalhar com melt()e a cast()partir do pacote de reformulação, mas parece que algumas coisas eu não entendo. Alguém pode...
Imagine Você executa uma regressão linear com quatro preditores numéricos (IV1, ..., IV4) Quando apenas IV1 é incluído como preditor, o beta padronizado é +.20 Quando você também inclui IV2 a IV4, o sinal do coeficiente de regressão padronizado de IV1 muda para -.25(isto é, tornou-se...
Eu gostaria de determinar a importância relativa de conjuntos de variáveis para um randomForestmodelo de classificação em R. A importancefunção fornece a MeanDecreaseGinimétrica para cada preditor individual - é tão simples quanto resumir isso em cada preditor de um conjunto? Por exemplo: #...
Muito foi escrito sobre as opções de cores compatíveis com daltônicos para mapas, polígonos e regiões sombreadas em geral (veja, por exemplo, http://colorbrewer2.org ). Não consegui encontrar recomendações para cores de linha e espessura de linha variável para gráficos de linha. Os objetivos...
É possível calcular valores de AIC ou BIC para modelos de regressão de laço e outros modelos regularizados em que os parâmetros estão inserindo apenas parcialmente a equação. Como se determina os graus de liberdade? Estou usando R para ajustar modelos de regressão de laço com a glmnet()função do...
Quais são algumas diretrizes úteis para testar parâmetros (por exemplo, profundidade de interação, filho pequeno, taxa de amostragem etc.) usando o GBM? Digamos que tenho 70-100 recursos, uma população de 200.000 e pretendo testar a profundidade de interação de 3 e 4. Claramente, preciso fazer...
Entendo o que é um posterior, mas não tenho certeza do que o último significa? Como os 2 são diferentes? Kevin P Murphy indicou em seu livro, Aprendizado de Máquina: uma Perspectiva Probabilística , que é "um estado de crença interno". O que isso realmente significa? Fiquei com a impressão de que...
Contexto: Em resposta a uma pergunta anterior sobre pesquisa reproduzível, Jake escreveu Um problema que descobrimos ao criar nosso arquivo JASA foi que as versões e padrões dos pacotes CRAN foram alterados. Portanto, nesse arquivo, também incluímos as versões dos pacotes que usamos. O sistema...
Quero gerar o gráfico descrito no livro ElemStatLearn "Os elementos do aprendizado estatístico: mineração de dados, inferência e previsão. Segunda edição" de Trevor Hastie e Robert Tibshirani e Jerome Friedman. O enredo é: Gostaria de saber como posso produzir esse gráfico exato R,...
É possível visualizar a saída da Análise de componentes principais de maneiras que fornecem mais informações do que apenas tabelas de resumo? É possível fazer isso quando o número de observações é grande, digamos ~ 1e4? E é possível fazer isso em R [outros ambientes são...
Estou tentando analisar alguns dados usando um modelo de efeito misto. Os dados que eu coletei representam o peso de alguns animais jovens de diferentes genótipos ao longo do tempo. Estou usando a abordagem proposta aqui:
Eu encontrei algumas distribuições para as quais BUGS e R têm diferentes parametrizações: Normal, log-Normal e Weibull. Para cada uma delas, entendo que o segundo parâmetro usado por R precisa ser transformado inversamente (1 / parâmetro) antes de ser usado no BUGS (ou JAGS no meu caso). Alguém...
/ editar: Acompanhamento adicional agora você pode usar o irlba :: prcomp_irlba / edit: acompanhando meu próprio post. irlbaagora possui argumentos de "centro" e "escala", que permitem usá-lo para calcular componentes principais, por exemplo: pc <- M %*% irlba(M, nv=5, nu=0,...
Estou tentando representar graficamente o número de ações dos usuários (neste caso, "curtidas") ao longo do tempo. Portanto, tenho "Número de ações" como meu eixo y, meu eixo x é o tempo (semanas) e cada linha representa um usuário. Meu problema é que quero analisar esses dados para um conjunto...
Sou muito novo nas estatísticas bayesianas e isso pode ser uma pergunta boba. Mesmo assim: Considere um intervalo credível com um prior que especifique uma distribuição uniforme. Por exemplo, de 0 a 1, em que 0 a 1 representa toda a gama de valores possíveis de um efeito. Nesse caso, um intervalo...
Acho que R pode levar muito tempo para gerar gráficos quando milhões de pontos estão presentes - o que não é surpreendente, uma vez que os pontos são plotados individualmente. Além disso, essas parcelas costumam ser muito confusas e densas para serem úteis. Muitos dos pontos se sobrepõem e formam...
Estou tentando interpretar fatores de inflação de variância usando a viffunção no pacote de R car. A função imprime um generalizado e também . De acordo com o arquivo de ajuda , esse último valorVIFVIF\text{VIF}GVIF1 / ( 2 ⋅ df )GVIF1/(2⋅df)\text{GVIF}^{1/(2\cdot\text{df})} Para ajustar a...
Estou tentando acelerar a resposta com R. Eu finalmente quero usar as bibliotecas R para fazer a classificação de texto. Fiquei me perguntando quais são as experiências das pessoas com relação à escalabilidade de R quando se trata de fazer a classificação de texto. É provável que eu ocorra com...