Perguntas com a marcação «r»

Eu tenho uma pergunta sobre a interpretação da chamada tsboot em R. Verifiquei a documentação do Kendall e do pacote de inicialização, mas não sou mais inteligente do que antes. Quando executo um bootstrap usando, por exemplo, o exemplo no pacote Kendall, onde a estatística de teste é a tau …

11 r  bootstrap 

Usei o seguinte código R para ajustar um modelo probit: p1 <- glm(natijeh ~ ., family=binomial(probit), data=data1) stepwise(p1, direction='backward/forward', criterion='BIC') Eu quero saber o que faz stepwisee backward/forwardfaz exatamente e como selecionar as variáveis?

11 r  probit  stepwise-regression 

Usando Amelia em R, obtive vários conjuntos de dados imputados. Depois disso, realizei um teste de medidas repetidas no SPSS. Agora, quero reunir os resultados dos testes. Eu sei que posso usar as regras do Rubin (implementadas por meio de qualquer pacote de imputação múltipla em R) para agrupar meios …

11 r  spss  p-value  multiple-imputation  pooling 

Existe uma estratégia para escolher o número de árvores em um GBM? Especificamente, o ntreesargumento em R's gbmfunção. Não vejo por que você não deve definir ntreeso valor mais alto possível. Percebi que um número maior de árvores reduz claramente a variabilidade dos resultados de vários GBMs. Eu não acho …

11 r  machine-learning  gbm 

Estou tentando usar o neuralnetpacote de R (documentação aqui ) para previsão. Aqui o que estou tentando fazer: library(neuralnet) x <- cbind(runif(50, min=1, max=500), runif(50, min=1, max=500)) y <- x[, 1] * x[, 2] train <- data.frame(x, y) n <- names(train) f <- as.formula(paste('y ~', paste(n[!n %in% 'y'], collapse = …

11 r  neural-networks  prediction  error 

Eu tenho um conjunto de dados de 20 anos de uma contagem anual de abundância de espécies para um conjunto de polígonos (~ 200 polígonos de forma irregular e contínua). Eu tenho usado a análise de regressão para inferir tendências (mudança na contagem por ano) para cada polígono, bem como …

11 r  time-series  correlation  spatial  spatio-temporal 

Considere o seguinte código e saída: par(mfrow=c(3,2)) # generate random data from weibull distribution x = rweibull(20, 8, 2) # Quantile-Quantile Plot for different distributions qqPlot(x, "log-normal") qqPlot(x, "normal") qqPlot(x, "exponential", DB = TRUE) qqPlot(x, "cauchy") qqPlot(x, "weibull") qqPlot(x, "logistic") Parece que o gráfico de QQ para log-normal é quase …

11 r  data-visualization  interpretation  qq-plot 

Para dados de contagem que eu coletei, uso a regressão de Poisson para construir modelos. Eu faço isso usando a glmfunção em R, onde eu uso family = "poisson". Para avaliar possíveis modelos (tenho vários preditores), uso o AIC. Por enquanto, tudo bem. Agora eu quero executar a validação cruzada. …

11 r  generalized-linear-model  poisson-distribution 

Eu tenho um conjunto de dados que contém o número de ações executadas por indivíduos ao longo de sete dias. A ação específica não deve ser relevante para esta pergunta. Aqui estão algumas estatísticas descritivas para o conjunto de dados: AlcanceSignificarVariaçãoNúmero de observações0 - 77218,22791696Range0−772Mean18.2Variance2791Number of observations696 \begin{array}{|c|c|} \hline \text{Range} …

11 r  distributions  poisson-distribution  mean  sample 

A implementação randomForest não permite amostragem além do número de observações, mesmo quando amostragem com substituição. Por que é isso? Funciona bem: rf <- randomForest(Species ~ ., iris, sampsize=c(1, 1, 1), replace=TRUE) rf <- randomForest(Species ~ ., iris, sampsize=3, replace=TRUE) O que eu quero fazer: rf <- randomForest(Species ~ ., …

11 r  sampling  random-forest  stratification  oversampling 

Para analisar a contagem de pássaros inflada a zero, eu gostaria de aplicar modelos de contagem inflada a zero usando o pacote R pscl . No entanto, observando o exemplo fornecido na documentação para uma das principais funções ( ? Zeroinfl ), começo a duvidar de qual é a real …

11 r  poisson-distribution  zero-inflation 

Eu gostaria de usar lme4para ajustar uma regressão de efeitos mistos e multcompcalcular as comparações aos pares. Eu tenho um conjunto de dados complexo com vários preditores contínuos e categóricos, mas minha pergunta pode ser demonstrada usando o ChickWeightconjunto de dados interno como exemplo: m <- lmer(weight ~ Time * …

11 r  mixed-model  multiple-comparisons 

Gostaria de testar em que regressão se encaixa melhor nos meus dados. Minha variável dependente é uma contagem e possui muitos zeros. E eu precisaria de ajuda para determinar qual modelo e família usar (poisson ou quasipoisson ou regressão de poisson inflada a zero) e como testar as suposições. Regressão …

11 r  count-data  zero-inflation  poisson-regression 

Olá, estou fazendo um curso de pós-graduação em estatística e abordamos estatísticas de testes e outros conceitos. No entanto, muitas vezes sou capaz de aplicar as fórmulas e desenvolver uma espécie de intuição sobre como as coisas funcionam, mas muitas vezes fico com a sensação de que talvez se eu …

11 r  hypothesis-testing  sas  simulation  computational-statistics 

Quero visualizar os resultados de um agrupamento (produzido com protoclust{protoclust}) criando gráficos de dispersão para cada par de variáveis ​​usadas para classificar meus dados, colorindo por classes e sobrepondo as elipses para o intervalo de confiança de 95% para cada uma das classes (para inspecionar quais classes de elipses se …

11 r  confidence-interval  ggplot2  scatterplot