Perguntas com a marcação «r»

Eu tenho um GLMM do formulário: lmer(present? ~ factor1 + factor2 + continuous + factor1*continuous + (1 | factor3), family=binomial) Quando uso drop1(model, test="Chi"), obtenho resultados diferentes dos que utilizo Anova(model, type="III")na embalagem do carro ou summary(model). Estes dois últimos dão as mesmas respostas. Usando um monte de dados fabricados, …

10 r  anova  glmm  r  mixed-model  bootstrap  sample-size  cross-validation  roc  auc  sampling  stratification  random-allocation  logistic  stata  interpretation  proportion  r  regression  multiple-regression  linear-model  lm  r  cross-validation  cart  rpart  logistic  generalized-linear-model  econometrics  experiment-design  causality  instrumental-variables  random-allocation  predictive-models  data-mining  estimation  contingency-tables  epidemiology  standard-deviation  mean  ancova  psychology  statistical-significance  cross-validation  synthetic-data  poisson-distribution  negative-binomial  bioinformatics  sequence-analysis  distributions  binomial  classification  k-means  distance  unsupervised-learning  euclidean  correlation  chi-squared  spearman-rho  forecasting  excel  exponential-smoothing  binomial  sample-size  r  change-point  wilcoxon-signed-rank  ranks  clustering  matlab  covariance  covariance-matrix  normal-distribution  simulation  random-generation  bivariate  standardization  confounding  z-statistic  forecasting  arima  minitab  poisson-distribution  negative-binomial  poisson-regression  overdispersion  probability  self-study  markov-process  estimation  maximum-likelihood  classification  pca  group-differences  chi-squared  survival  missing-data  contingency-tables  anova  proportion 

Não obtenho a diferença entre the rfobject$importancee importance(rfobject)na coluna MeanDecreaseAccuracy. Exemplo: > data("iris") > fit <- randomForest(Species~., data=iris, importance=TRUE) > fit$importance setosa versicolor virginica MeanDecreaseAccuracy MeanDecreaseGini Sepal.Length 0.027078501 0.019418330 0.040497602 0.02898837 9.173648 Sepal.Width 0.008553449 0.001962036 0.006951771 0.00575489 2.472105 Petal.Length 0.313303381 0.291818815 0.280981959 0.29216790 41.284869 Petal.Width 0.349686983 0.318527008 0.270975757 0.31054451 46.323415 …

10 r  random-forest  importance 

Estrutura de dados > str(data) 'data.frame': 6138 obs. of 10 variables: $ RT : int 484 391 422 516 563 531 406 500 516 578 ... $ ASCORE : num 5.1 4 3.8 2.6 2.7 6.5 4.9 2.9 2.6 7.2 ... $ HSCORE : num 6 2.1 7.9 1 6.9 …

10 r  mixed-model  repeated-measures  model-selection  lme4-nlme 

Eu ajustei um modelo lognormal usando R com um conjunto de dados. Os parâmetros resultantes foram: meanlog = 4.2991610 sdlog = 0.5511349 Gostaria de transferir esse modelo para o Scipy, que nunca usei antes. Usando o Scipy, consegui obter uma forma e escala de 1 e 3,1626716539637488e + 90 - …

10 r  python  numpy  scipy 

Percebi que quando tento encontrar os valores críticos para o Mann-Whitney U usando R, os valores são sempre 1 + valor crítico. Por exemplo, para , o valor crítico (bicaudal) é 8, enquanto para , a (bicaudal) ) o valor crítico é 22 (verifique as tabelas ), mas:α = 0,05 …

10 r  hypothesis-testing  nonparametric  wilcoxon-mann-whitney 

Estou fazendo uma ANOVA de sentido único (por espécie) com contrastes personalizados. [,1] [,2] [,3] [,4] 0.5 -1 0 0 0 5 1 -1 0 0 12.5 0 1 -1 0 25 0 0 1 -1 50 0 0 0 1 onde eu comparo a intensidade 0,5 contra 5, 5 …

10 r  anova  contrasts  generalized-least-squares 

Estou tentando analisar o lead-lag entre as séries temporais de dois preços das ações. Em análises regulares de séries temporais, podemos fazer a correlação cruzada, VECM (Granger Causality). No entanto, como lidar com o mesmo em séries temporais com espaçamento irregular. A hipótese é que um dos instrumentos lidera o …

10 r  time-series  cross-correlation  unevenly-spaced-time-series 

Quero modelar duas variáveis ​​de tempo diferentes, algumas das quais são fortemente colineares nos meus dados (idade + coorte = período). Ao fazer isso, tive alguns problemas com ee lmerinterações poly(), mas provavelmente não está limitado a lmer, obtive os mesmos resultados com o nlmeIIRC. Obviamente, falta entender o que …

10 r  lme4-nlme  polynomial 

Estou tentando replicar o trabalho de um colega e estou movendo a análise de Stata para R. Os modelos que ela emprega invocam a opção "cluster" na função nbreg para agrupar os erros padrão. Veja http://repec.org/usug2007/crse.pdf para obter uma descrição bastante completa do quê e por que dessa opção Minha …

10 r  stata  negative-binomial 

Se eu ajustar meus dados com algo como lm(y~a*b), na sintaxe R, onde aé uma variável binária e buma variável numérica, o a:btermo de interação é a diferença entre a inclinação de y~bat a= 0 e at a= 1. Agora, digamos que a relação entre ye bé curvilínea. Se agora …

10 r  multiple-regression  splines 

Existem muitos métodos para executar a regularização - regularização baseada na norma , L 1 e L 2, por exemplo. De acordo com Friedman Hastie e Tibsharani , o melhor regularizador depende do problema: a natureza da verdadeira função do alvo, a base específica usada, a relação sinal / ruído …

10 r  regression  machine-learning  regularization 

Eu realizei uma ANOVA de três medidas repetidas; que análises post-hoc são válidas? Este é um design totalmente equilibrado (2x2x2) com um dos fatores que possui uma medida repetida dentro dos sujeitos. Estou ciente das abordagens multivariadas para ANOVA de medidas repetidas em R, mas meu primeiro instinto é prosseguir …

10 r  anova  repeated-measures  interaction  post-hoc 

Estou tentando determinar se probabilidades simples funcionarão para o meu problema ou se será melhor usar (e aprender sobre) métodos mais sofisticados, como regressão logística. A variável de resposta nesse problema é uma resposta binária (0, 1). Eu tenho um número de variáveis ​​preditoras que são todas categóricas e não …

10 r  probability  logistic 

Estou modelando alguns dados nos quais acho que tenho dois efeitos aleatórios cruzados. Mas o conjunto de dados não é equilibrado e não tenho certeza do que precisa ser feito para justificá-lo. Meus dados são um conjunto de eventos. Um evento ocorre quando um cliente se reúne com um provedor …

10 r  mixed-model  random-effects-model  logistic 

Estou com muitos problemas com um conjunto de dados ao qual estou tentando aplicar o SEM. Supomos a existência de 5 fatores latentes A, B, C, D, E, com indicadores resp. A1 a A5 (fatores ordenados), B1 a B3 (quantitativo), C1, D1, E1 (todos os três últimos fatores ordenados, com …

10 r  modeling  multiple-regression  sem