Perguntas com a marcação «r»

Eu gostaria de testar a diferença na resposta de duas variáveis ​​para um preditor. Aqui está um exemplo mínimo reproduzível. library(nlme) ## gls is used in the application; lm would suffice for this example m.set <- gls(Sepal.Length ~ Petal.Width, data = iris, subset = Species == "setosa") m.vir <- gls(Sepal.Length …

29 r  data-visualization  multivariate-analysis  hypothesis-testing 

Eu uso a função auto.arima () no pacote de previsão para ajustar os modelos ARMAX a uma variedade de covariáveis. No entanto, muitas vezes tenho um grande número de variáveis ​​para selecionar e geralmente termino com um modelo final que funciona com um subconjunto delas. Não gosto de técnicas ad-hoc …

29 r  time-series  lasso  regularization  elastic-net 

Sou novato em estatísticas e preciso de sua ajuda. Eu tenho uma pequena amostra, da seguinte maneira: H4U 0.269 0.357 0.2 0.221 0.275 0.277 0.253 0.127 0.246 Eu executei o teste Shapiro-Wilk usando R: shapiro.test(precisionH4U$H4U) e obtive o seguinte resultado: W = 0.9502, p-value = 0.6921 Agora, se eu assumir …

29 r  distributions  interpretation  goodness-of-fit  normality-assumption 

Eu sempre uso lm()em R para executar regressão linear de yyy em . Essa função retorna um coeficiente tal quexxxββ\betay=βx.y=βx.y = \beta x. Hoje eu aprendi sobre o total de mínimos quadrados e essa princomp()função (análise de componentes principais, PCA) pode ser usada para realizá-lo. Deve ser bom para mim …

29 r  pca  least-squares  deming-regression  total-least-squares 

Podemos usar lm()para prever um valor, mas ainda precisamos da equação da fórmula do resultado em alguns casos. Por exemplo, adicione a equação aos gráficos.

29 r  regression  lm 

Estou usando o sinal de intercalação para executar uma floresta aleatória validada cruzada em um conjunto de dados. A variável Y é um fator. Não há NaN, Inf ou NA no meu conjunto de dados. No entanto, ao executar a floresta aleatória, recebo Error in randomForest.default(m, y, ...) : NA/NaN/Inf …

29 r  random-forest  caret  regression  prediction  fitting  social-science  poisson-distribution  distributions  characteristic-function  bayesian  prior  regression  normal-distribution  interaction  nonparametric  skewness  svm  standard-deviation  standard-error  regression-coefficients  igraph  natural-language  word2vec  word-embeddings  regression  machine-learning  sampling  r  regression  machine-learning  random-forest  ensemble  sampling  unbiased-estimator  proof  estimators  mse  probability  conditional-probability  bayes  anova  missing-data  neural-networks  recommender-system  r  confidence-interval  sample  multiple-imputation  r  time-series  forecasting  mase 

Estou usando o pacote robustbase para executar uma estimativa de glm. No entanto, quando faço isso, recebo o seguinte erro: Error in solve.default(crossprod(X, DiagB * X)/nobs, EEq) : system is computationally singular: reciprocal condition number = 1.66807e-16 O que isso significa / indica? E como posso depurá-lo? PS. Se você …

29 r  generalized-linear-model  robust 

Estou usando os modelos R (3.1.1) e ARIMA para previsão. Gostaria de saber qual deve ser o parâmetro "frequency", atribuído na ts()função , se estiver usando dados de séries temporais que sejam: separados por minutos e estão espalhados por 180 dias (1440 minutos / dia) separados por segundos e estão …

28 r  time-series  arima  multiple-seasonalities  mape 

[O título inicial "Medição de similaridade para árvores de cluster hierárquicas" foi posteriormente alterado por @ttnphns para refletir melhor o tópico] Estou realizando várias análises hierárquicas de cluster em um quadro de dados de registros de pacientes (por exemplo, semelhante a http://www.biomedcentral.com/1471-2105/5/126/figure/F1?highres=y ) Estou experimentando diferentes medidas de distância , …

28 r  clustering  distance-functions  similarities  dendrogram 

Gostaria de obter um valor-p e um tamanho de efeito de uma variável categórica independente (com vários níveis) - que é "geral" e não para cada nível separadamente, como é o resultado normal de lme4R. É exatamente como a coisa que as pessoas relatam ao executar uma ANOVA. Como posso …

28 r  hypothesis-testing  anova  mixed-model  lme4-nlme 

Entendo que essa questão é bastante ampla, mas me pergunto quais devem ser os pontos decisivos na decisão de criar (ou não) um novo pacote para R. Para ser mais específico, acrescentaria que a questão não é sobre os motivos para use R em si mesmo, mais sobre a decisão …

28 r  software 

Estou usando o glmnet para calcular estimativas de regressão de crista. Eu obtive alguns resultados que me fizeram suspeitar que o glmnet está realmente fazendo o que eu acho que ele faz. Para verificar isso, escrevi um script R simples, onde comparo o resultado da regressão de cume feita por …

28 r  ridge-regression  glmnet 

Espero que todos não se importem com essa pergunta, mas preciso de ajuda para interpretar a saída de um modelo linear de efeitos mistos que tenho tentado aprender na R. Sou novo na análise de dados longitudinal e na regressão linear de efeitos mistos. Eu tenho um modelo que me …

28 r  mixed-model  interpretation  panel-data 

Como exemplo, considere o ChickWeightconjunto de dados em R. A variação obviamente aumenta com o tempo, portanto, se eu usar uma regressão linear simples como: m <- lm(weight ~ Time*Diet, data=ChickWeight) Minhas perguntas: Quais aspectos do modelo serão questionáveis? Os problemas estão limitados à extrapolação fora do Timeintervalo? Quão tolerante …

28 r  regression  heteroscedasticity  assumptions 

É útil estudar o código de análise de dados de especialistas. Recentemente, estive lendo o github e há várias pessoas compartilhando código de análise de dados por lá. Isso inclui alguns pacotes R (que obviamente estão disponíveis diretamente no CRAN), mas também vários exemplos de pesquisa reproduzível, principalmente usando o …

28 r  reproducible-research