Perguntas com a marcação «r»

Eu tenho os seguintes vetores X e Y simples: > X [1] 1.000 0.063 0.031 0.012 0.005 0.000 > Y [1] 1.000 1.000 1.000 0.961 0.884 0.000 > > plot(X,Y) Quero fazer a regressão usando o log do X. Para evitar o log (0), tento colocar +1 ou +0,1 ou …

10 r  regression  lognormal 

Imagine uma situação: temos registros históricos (20 anos) de três minas. A presença de prata aumenta a probabilidade de encontrar ouro no próximo ano? Como testar essa pergunta? Aqui estão dados de exemplo: mine_A <- c("silver","rock","gold","gold","gold","gold","gold", "rock","rock","rock","rock","silver","rock","rock", "rock","rock","rock","silver","rock","rock") mine_B <- c("rock","rock","rock","rock","silver","rock","rock", "silver","gold","gold","gold","gold","gold","rock", "silver","rock","rock","rock","rock","rock") mine_C <- c("rock","rock","silver","rock","rock","rock","rock", "rock","silver","rock","rock","rock","rock","silver", "gold","gold","gold","gold","gold","gold") time <- …

10 r  time-series  hypothesis-testing  stochastic-processes 

Dadas três variáveis, ye x, que são positivas contínuas e z, que são categóricas, tenho dois modelos candidatos dados por: fit.me <- lmer( y ~ 1 + x + ( 1 + x | factor(z) ) ) e fit.fe <- lm( y ~ 1 + x ) Espero comparar esses …

10 r  hypothesis-testing  mixed-model  lme4-nlme 

Estou procurando agrupar / mesclar nós em um gráfico usando o agrupamento de gráficos em 'r'. Aqui está uma variação impressionante do meu problema. Existem dois "clusters" Existe uma "ponte" conectando os clusters Aqui está uma rede candidata: Quando olho para a distância da conexão, a "contagem de hop", se …

10 r  clustering  data-visualization  numerics 

Eu tenho algumas variáveis ​​e estou interessado em encontrar relações não lineares entre elas. Então, decidi encaixar um spline ou loess e imprimir bons gráficos (veja o código abaixo). Mas também quero ter algumas estatísticas que me dêem uma idéia do quanto é provável que o relacionamento seja uma questão …

10 r  regression  splines  loess 

Fiz um bootstrap com um modelo misto (várias variáveis ​​com interação e uma variável aleatória). Eu obtive este resultado (apenas parcial): > boot_out ORDINARY NONPARAMETRIC BOOTSTRAP Call: boot(data = a001a1, statistic = bootReg, R = 1000) Bootstrap Statistics : original bias std. error t1* 4.887383e+01 -1.677061e+00 4.362948e-01 t2* 3.066825e+01 1.264024e+00 …

10 r  confidence-interval  bootstrap  bias 

Estou executando uma regressão OLS em pool usando o pacote plm em R. No entanto, minha pergunta é mais sobre estatísticas básicas, então tento publicá-la aqui primeiro;) Como meus resultados de regressão produzem resíduos heterocedásticos, eu gostaria de tentar usar erros padrão robustos de heterocedasticidade. Como resultado coeftest(mod, vcov.=vcovHC(mod, type="HC0")), …

10 r  anova  multiple-regression  heteroscedasticity  robust-standard-error 

Estou executando glms em R (modelos lineares generalizados). Eu pensei que conhecia pvalues ​​- até que vi que chamar um resumo para um glm não lhe desse um pvalor predominante representativo do modelo como um todo - pelo menos não no lugar onde os modelos lineares o fazem. Gostaria de …

10 r  statistical-significance  generalized-linear-model  p-value  descriptive-statistics 

Estou brincando com o conjunto de dados de câncer de mama e criei um gráfico de dispersão de todos os atributos para ter uma idéia de quais têm mais efeito na previsão da classe malignant(azul) de benign(vermelho). Entendo que a linha representa o eixo x e a coluna representa o …

10 r  data-visualization  interpretation  scatterplot 

Toda a literatura de Modelagem de Distribuição de Espécies sugere que, ao prever a presença / ausência de uma espécie usando um modelo que gera probabilidades (por exemplo, RandomForests), é importante escolher a probabilidade de limiar pela qual classificar uma espécie como presença ou ausência é importante e deve-se nem …

10 r  classification  random-forest  threshold 

Estou tendo problemas para entender como a varImpfunção funciona para um modelo randomForest com o caretpacote. No exemplo abaixo, o recurso var3 recebe zero importância usando a varImpfunção de cursor , mas o modelo final randomForest subjacente tem importância diferente de zero para o recurso var3. Por que esse é …

10 r  caret  random-forest 

Estou trabalhando com alguns dados do mundo real e os modelos de regressão estão produzindo alguns resultados contra-intuitivos. Normalmente confio nas estatísticas, mas, na realidade, algumas dessas coisas não podem ser verdadeiras. O principal problema que estou vendo é que um aumento em uma variável está causando um aumento na …

10 r  regression  regression-coefficients  computational-statistics 

Minha situação: tamanho pequeno da amostra: 116 variável de resultado binário lista longa de variáveis ​​explicativas: 44 variáveis ​​explicativas não vieram do alto da minha cabeça; sua escolha foi baseada na literatura. a maioria dos casos na amostra e a maioria das variáveis ​​têm valores ausentes. Abordagem da seleção de …

10 r  spss  lasso  data-imputation 

Sou relativamente novo em estatística e R. Gostaria de conhecer o processo para determinar os parâmetros ARIMA para meu conjunto de dados. Você pode me ajudar a descobrir o mesmo usando R e teoricamente (se possível)? Os dados variam de 12 de janeiro a 14 de março e retratam as …

10 r  arima  box-jenkins 

Eu tenho algumas experiências com modelagem de séries temporais, na forma de modelos ARIMA simples e assim por diante. Agora eu tenho alguns dados que exibem clustering de volatilidade e gostaria de tentar começar ajustando um modelo GARCH (1,1) nos dados. Eu tenho uma série de dados e acho que …

10 r  regression  garch