Perguntas com a marcação «r»

Estou estimando o modelo onde e são parâmetros , é um vetor de parâmetros length, é uma matriz de dados , a variável dependente é binária e é um modelo probit, portanto, a função de distribuição cumulativa da normal normal distribuição. Para derivar a expectativa, assumiu-se que os erros são …

7 r  maximum-likelihood  nonlinear  nls 

Eu executei uma ANOVA simples em R e, em seguida, gerei as seguintes comparações de médias de TukeyHSD (): Eu tenho uma idéia muito boa (eu acho) do que tudo isso significa, exceto o 'p adj'. Se eu estiver correto: A diferença nas pontuações dos testes entre os Juniors e …

7 r  tukey-hsd 

Variável dependente Eu tenho um valor dependente no intervalo de [0,1]. Significado 0 e 1, e todos os valores intermediários estão incluídos. Portanto, esse é um valor proporcional, como, por exemplo, a porcentagem de terra que um agricultor fertiliza. Modelo O modelo no qual estou focando atualmente é um modelo …

7 r  logistic  generalized-linear-model  proportion  logit 

Talvez essa pergunta dependa dos dados fornecidos, mas existe um método "melhor" de inicialização do que os outros? Estou simplesmente usando um conjunto de dados de uma variável (que consiste nas diferenças entre as pontuações de futebol (2 equipes) nas últimas 15 semanas) .. Primeiro, observe a inclinação correta desses …

7 r  confidence-interval  bootstrap 

Tenho muitas medidas para vários indivíduos, mas não sei como explicar essa estrutura de repetição de medidas ao executar um modelo de floresta aleatório. Existe uma maneira de explicar a estrutura de dados subjacente dos dados longitudinais usando um modelo de floresta aleatória? Isso é mesmo necessário? - parece-me que …

7 r  repeated-measures  random-forest  panel-data 

Então, vamos comparar duas distribuições normais Do this x times: runs <- 100000 a.samples <- rnorm(runs, mean = 5) b.samples <- rbeta(runs, mean = 0) mc.p.value <- sum(a.samples > b.samples)/runs Os valores mc.p. abaixo de nosso alfa (0,05) dividido por x dariam a taxa de erro do tipo 1. Nosso …

7 r  hypothesis-testing  monte-carlo 

Eu tenho duas amostras. Observando suas densidades, uma parece simétrica e a outra, de alguma distribuição da direita. Eu gostaria de testar se os dois não têm a mesma assimetria (ignorando questões de inferência seletiva aqui). Meu plano é pegar uma replicação de bootstrap das duas amostras, calcular a diferença …

7 r  hypothesis-testing  bootstrap  skewness 

Estou usando o pacote Causal Impact em R para inferir o efeito causal de uma intervenção em alguns dados altamente correlacionados e sazonais. Especificamente, recebi 17 dias de dados por hora, a intervenção ocorrendo no final do dia 13. Tenho dois conjuntos de dados de controle que não são afetados …

7 r  intervention-analysis  causalimpact 

Estou recebendo o erro "Novos níveis de fatores não presentes nos dados de treinamento". Mas verifiquei os nlevels e a classe de todas as colunas em desenvolvimento, bem como os dados de teste, e eles são os mesmos. Alguma explicação plausível?

7 r  machine-learning  random-forest  many-categories 

Estou realizando uma regressão logística. Criei os seguintes dados de teste (os dois preditores e o critério são variáveis ​​binárias): UV1 UV2 AV 1 1 1 1 2 1 1 1 3 1 1 1 4 1 1 1 5 1 1 1 6 1 1 1 7 1 1 …

7 r  regression  logistic  logit 

Estou trabalhando em um problema no qual tenho vários pares de homens atualmente vivos, icada um com um suposto ancestral paterno nigerações atrás (com base em evidências genealógicas) e onde tenho informações sobre se há uma incompatibilidade no genótipo cromossômico Y (exclusivamente paternalmente herdadoxi= 1 para incompatibilidade, 0 se houver …

7 r  confidence-interval  maximum-likelihood  bootstrap  delta-method 

Para explorar como a LASSOregressão funciona, escrevi um pequeno pedaço de código que deve otimizar a LASSOregressão escolhendo o melhor parâmetro alfa. Não consigo descobrir por que a LASSOregressão está me dando resultados tão instáveis ​​para o parâmetro alfa após a validação cruzada. Aqui está o meu código Python: from …

7 r  regression  python  lasso  penalized 

Ao usar o comando drop1 em R para construção de modelo, é dito que a variável com o menor valor de AIC deve ser descartada. Qual poderia ser a razão para o mesmo? Eu sei que a AIC fala sobre perda de informações e um valor mais baixo da AIC …

7 r  regression  aic 

Suponha que eu tenha uma série temporal xtxt x_t que eu quero ajustar usando um modelo ARIMA (1,1,0) do formulário: Δxt= α Δxt - 1+WtΔxt=αΔxt−1+wt \Delta x_t = \alpha \Delta x_{t-1} + w_t Isso pode ser reescrito como: xt-xt - 1= α (xt - 1-xt - 2) +Wtxt−xt−1=α(xt−1−xt−2)+wt x_t - …

7 r  time-series  arima 

Estou pensando em maneiras de calcular elipses de dados e confiança em torno de uma mediana bivariada. Por exemplo, eu posso calcular facilmente uma elipse de dados ou uma elipse de confiança para a média bivariada dos seguintes dados (aqui apenas mostrando uma elipse de dados) library("car") set.seed(1) df <- …

7 r  confidence-interval  median  bivariate  ellipse