Perguntas com a marcação «r»

Qual é a diferença entre (1|DNA.concentration/mouse.id)e (DNA.concentration|mouse.id)? O que significam os símbolos |e /significam dentro da sintaxe para o efeito aleatório?

8 r  lme4-nlme 

Estou usando a função Deming fornecida por Terry T. neste segmento de ajuda r arquivado . Estou comparando dois métodos, então tenho dados parecidos com este: y x stdy stdx 1 1.2 0.23 0.67 2 1.8 0.05 0.89 4 7.5 1.13 0.44 ... ... ... ... Fiz minha regressão de …

8 r  r-squared  deming-regression  total-least-squares 

Eu tenho uma regressão logística binária com apenas um preditor de fator fixo binário. A razão de eu não fazer isso como um quadrado de Chi ou o teste exato de Fisher é que também tenho vários fatores aleatórios (existem vários pontos de dados por indivíduo e os indivíduos estão …

8 r  logistic  generalized-linear-model  odds-ratio  relative-risk 

Fechadas. Esta questão está fora de tópico . No momento, não está aceitando respostas. Deseja melhorar esta pergunta? Atualize a pergunta para que ela esteja no tópico de Validação cruzada. Fechado há 4 anos . Estou trabalhando nos benchmarks de algumas redes neurais e eles levam uma quantidade significativa de …

8 r  neural-networks 

Estou trabalhando com o conjunto de dados "geyser" do pacote MASS e comparando as estimativas de densidade de kernel do pacote np. Meu problema é entender a estimativa de densidade usando a validação cruzada de mínimos quadrados e o kernel Epanechnikov: blep<-npudensbw(~geyser$waiting,bwmethod="cv.ls",ckertype="epanechnikov") plot(npudens(bws=blep)) Para o kernel gaussiano, parece estar bem: …

8 r  nonparametric  kernel-smoothing 

Estou usando R para executar regressão linear. Vi maneiras de calcular intervalos de previsão, mas estes dependem de dados homoscedásticos. Existe uma maneira de calcular intervalos de previsão com dados heterocedásticos?

8 r  regression  least-squares  heteroscedasticity  prediction-interval 

Eu tenho um conjunto de dados x, y que estou usando para criar uma floresta aleatória. Os dados x são um vetor de valores que inclui algumas NAs. Então, eu uso rfImputepara manipular os dados ausentes e criar uma floresta aleatória. Agora eu tenho uma nova observação invisível x (com …

8 r  random-forest  missing-data 

A arimaxfunção no TSApacote é, até onde sei, o único Rpacote que se encaixará em uma função de transferência para modelos de intervenção. Falta uma função de previsão, que às vezes é necessária. A seguir, uma solução alternativa para esse problema, aproveitando o excelente forecastpacote? Os intervalos preditivos estarão corretos? …

8 r  time-series  forecasting  arima  intervention-analysis 

Estou executando algumas otimizações com a implementação do BFGS do optim. A função objetivo é na verdade um algoritmo computacional, não apenas matemática. Descobri que quando adiciono uma penalidade de L1, as coisas diminuem um pouco. Por que isso pode ser? Existe algo no L1 que atrasa as coisas? Então, …

8 r  optimization  lasso 

Minha pergunta lida com árvores de classificação . Considere o seguinte exemplo do conjunto de dados Iris: Desejo selecionar manualmente o melhor preditor para a primeira divisão. De acordo com o algoritmo CART, o melhor recurso para fazer uma divisão é aquele que maximiza a diminuição da impureza da partição, …

8 r  machine-learning  classification  data-mining  cart 

Estou tentando determinar qual alfa usar na minha glmnetfunção, mas o arquivo de ajuda diz-me: Observe que o cv.glmnet NÃO procura valores para alfa. Um valor específico deve ser fornecido, caso contrário, alfa = 1 é assumido por padrão. Se os usuários também desejam validar cruzadamente o alfa, eles devem …

8 r  lasso  ridge-regression  glmnet 

Eu tenho duas séries temporais (parâmetros de um modelo para homens e mulheres) e pretendo identificar um modelo ARIMA apropriado para fazer previsões. Minha série cronológica se parece com: A plotagem e o ACF mostram não estacionário (os picos do ACF cortam muito lentamente). Assim, uso diferenciação e obtenho: Esse …

8 r  time-series  modeling  arima  autocorrelation 

Estou executando comparações post-hoc após um teste de Kruskal-Wallis. Estou usando o pacote PMCMR . > posthoc.kruskal.nemenyi.test( preference ~ instrument) Pairwise comparisons using Tukey and Kramer (Nemenyi) test with Tukey-Dist approximation for independent samples data: preference by instrument Cello Drums Guitar Drums 0.157 - - Guitar 0.400 0.953 - Harp …

8 r  ties 

Eu tenho a seguinte tabela de valores: 25 75 38 162 A razão de chances é 0,7037 e log (OR) é -0,3514. Para uma tabela de contingência com os valores a, b, c e d, a variação do log (OR) é dada por (1/a + 1/b + 1/c + 1/d) …

8 r  variance 

Como calcular média, variância, mediana, desvio padrão e modo de distribuição? Se eu gerar aleatoriamente números que formam a distribuição normal, especifiquei a média como m=24.2desvio padrão como sd=2.2: > dist = rnorm(n=1000, m=24.2, sd=2.2) Então eu posso fazer o seguinte: Significar: > mean(dist) [1] 24.17485 Variação: > var(dist) [1] …

8 r  distributions  mean