Perguntas com a marcação «r»

Tenho gráficos de caixa de 13 grupos que mostro em um gráfico. Os grupos têm populações desequilibradas e não são normalmente distribuídas. Quero mostrar quais pares são estatisticamente semelhantes (ou seja, têm valor de p kruskal.test <0,05) colocando a, b, c, etc. em cima das caixas correspondentes. Aqui está um …

9 r  statistical-significance  boxplot 

Em R, como especifico o modelo lmer sem efeito fixo global? Por exemplo, se eu disser algo como lmer(y ~ (1 | group) + (0 + x | group), data = my_df) o modelo ajustado será yij=a+αi+βixijyij=a+αi+βixijy_{ij} = a + \alpha_{i} + \beta_{i} x_{ij} Como encaixo o modelo yij=αi+βixijyij=αi+βixijy_{ij} = …

9 r  mixed-model  lme4-nlme 

Eu tenho o seguinte tipo de dados. Avaliei 10 indivíduos, cada um repetido 10 vezes. Tenho matriz de relação 10x10 (relação entre todas as combinações dos indivíduos). set.seed(1234) mydata <- data.frame (gen = factor(rep(1:10, each = 10)), repl = factor(rep(1:10, 10)), yld = rnorm(10, 5, 0.5)) Essa geração é de …

9 r  mixed-model 

Estou tentando ajustar uma função de decaimento exponencial aos valores y que se tornam negativos em valores x altos, mas não consigo configurar minha nlsfunção corretamente. Alvo Estou interessado na inclinação da função decay ( acordo com algumas fontes ). Como obtenho essa inclinação não é importante, mas o modelo …

9 r  model  exponential 

Estou procurando neste site um workshop sobre o GAM em R: http://qcbs.ca/wiki/r_workshop8 No final da seção, 2. Multiple smooth termseles mostram um exemplo, onde eles usam anovapara comparar três modelos diferentes para determinar o melhor modelo de ajuste. A saída é Analysis of Deviance Table Model 1: y ~ x0 …

9 r  anova  gam  mgcv 

Comecei a trabalhar com o GAM em R e adquiri o excelente livro de Simon Wood sobre o tópico ( "Modelos aditivos generalizados, uma introdução ao R" ). Com base em um de seus exemplos, estou analisando o seguinte: library(mgcv) data(trees) ct1<-gam(log(Volume) ~ Height + s(Girth), data=trees) Eu tenho duas …

9 r  regression  machine-learning  gam  mgcv 

Eu tenho ~ 1 milhão de pontos de dados. Aqui está o link para o arquivo data.txt. Cada um deles pode ter um valor entre 0 e 145. É um conjunto de dados discreto. Abaixo está o histograma do conjunto de dados. No eixo x é a contagem (0-145) e …

9 r  statistical-significance  goodness-of-fit  negative-binomial 

Antecedentes : Israel (e o Oriente Médio em geral) está repleto de protestos. Estou curioso, quando recebi uma foto, para estimar quantas pessoas estão nessa foto (geralmente uma foto de uma grande multidão). Que modelagem pode oferecer alguma solução para esse problema? (e, é claro, isso pode ser feito com …

8 r  image-processing  pattern-recognition 

Ajudo um colega a inicializar um modelo de efeitos mistos de meta-análise usando a estrutura do pacote metafor R da autoria de @Wolfgang. Interessante e preocupante, para um dos coeficientes do modelo, recebo uma distribuição bimodal durante a inicialização (consulte o painel inferior direito da figura abaixo). Eu acho que …

8 r  mixed-model  bootstrap  meta-analysis  bimodal 

A pergunta a seguir se baseia na discussão encontrada nesta página . Dada uma variável de resposta y, uma variável explicativa contínua xe um fator fac, é possível definir um Modelo Aditivo Geral (GAM) com uma interação entre xe facusando o argumento by=. De acordo com o arquivo de ajuda …

8 r  interaction  gam  mgcv 

Estou tendo dificuldades para acessar a mesma página que meu supervisor quando se trata de validar meu modelo. Analisei os resíduos (observados em relação aos valores ajustados) e usei isso como argumento para discutir os resultados obtidos pelo meu modelo, no entanto, meu supervisor insiste que a única maneira de …

8 r  references  glmm  validation  glmmtmb 

Estimei algumas medidas repetidas modelos de efeitos fixos, com um componente de erro aninhado, com base em variáveis ​​de agrupamento, ou seja, modelos não aninhados, usando plm . Agora estou interessado em teste se os modelos completos são significativamente diferentes, ou seja, onde é o modelo completo para e é …

8 r  hypothesis-testing  repeated-measures  nested-data  plm 

Eu tenho um conjunto de dados com 5.818.446 linhas e 51 colunas, onde 50 delas são preditores. Minha resposta é quantitativa, por isso estou interessado em um modelo de regressão. Estou tentando ajustar uma floresta aleatória aos meus dados usando o pacote de intercalação. No entanto, não tenho RAM suficiente …

8 r  random-forest  large-data 

Estou vendo como a distância euclidiana mínima esperada entre pontos aleatoriamente uniformes e a origem muda à medida que aumentamos a densidade de pontos aleatórios ( pontos por unidade quadrada ) ao redor da origem. Eu consegui chegar a um relacionamento entre os dois descritos como tal: Distância mínima esperada …

8 r  expected-value  monte-carlo  uniform  minimum 

Eu tenho um conjunto de dados com 338 preditores e 570 instâncias (infelizmente não é possível carregar) no qual estou usando o Lasso para executar a seleção de recursos. Em particular, estou usando a cv.glmnetfunção da glmnetseguinte maneira, onde mydata_matrixé uma matriz binária de 570 x 339 e a saída …

8 r  classification  lasso  glmnet