Perguntas com a marcação «r»

Tenho uma série temporal um tanto barulhenta que paira em diferentes níveis. Por exemplo, os seguintes dados: Eu tenho os dados de linha sólida disponíveis e gostaria de obter uma estimativa para a linha tracejada. Deve ser constante por partes. Quais algoritmos são apropriados para testar aqui? Até agora, minhas …

7 r  time-series  smoothing 

Em R, os parâmetros estimados por glm()e seus erros padrão são fornecidos na escala de links, como alguém recentemente me esclareceu aqui . Faz sentido fornecer os parâmetros e seus erros padrão na mesma escala, mas por que não exibi-los na escala original dos dados? Eu imagino que a maioria …

7 r  generalized-linear-model  standard-error 

Se meu qqplot é linear, mas o gradiente não é o mesmo da linha de 45 graus, o que isso sugere? Estou tentando examinar o ajuste da distribuição laplace aos meus dados de amostra, por isso gerei aleatoriamente observações distribuídas de laplace (com parâmetros estimados a partir da minha amostra) …

7 r  qq-plot 

Considere um planejamento fatorial dentro do sujeito e dentro do item, onde a variável de tratamento experimental possui dois níveis (condições). Seja m1o modelo máximo e m2o modelo sem correlações aleatórias. m1: y ~ condition + (condition|subject) + (condition|item) m2: y ~ condition + (1|subject) + (0 + condition|subject) + …

7 r  mixed-model  lme4-nlme  categorical-encoding  machine-learning  pandas  proportion  r  irt  distributions  conditional-probability  kernel-smoothing  r  data-visualization  r  mixed-model  sas  curve-fitting  matplotlib  data-visualization  python  matplotlib  regression  logistic  simulation  sas  jmp  logit  beta-regression  regression  maximum-likelihood  posterior 

library(datasets) library(nlme) n1 <- nlme(circumference ~ phi1 / (1 + exp(-(age - phi2)/phi3)), data = Orange, fixed = list(phi1 ~ 1, phi2 ~ 1, phi3 ~ 1), random = list(Tree = pdDiag(phi1 ~ 1)), start = list(fixed = c(phi1 = 192.6873, phi2 = 728.7547, phi3 = 353.5323))) Eu ajustei um …

7 r  mixed-model  sas 

Eu tenho uma pergunta sobre como obter os erros padrão dos coeficientes no meu modelo GLM. Eu tenho a matriz de informações de fisher que calculei manualmente, mas não é dimensionada. Como posso escalar a matriz de informações do fisher para obter os mesmos erros padrão da função GLM?

7 r  generalized-linear-model  covariance 

O código abaixo gera um conjunto de dados de teste que consiste em uma série de probabilidades de "sinal" com ruído binomial ao seu redor. O código então usa 5000 conjuntos de números aleatórios como séries "explicativas" e calcula o valor p da regressão logística para cada um. Acho que …

7 r  regression  logistic  p-value 

Estou trabalhando em uma distribuição específica cujo cdf inverso não existe na forma fechada. O cdf da distribuição é dado por F( x ; d, m , p , α , β) =1 - ( 1 +xm)- dexp( - βxα)1 - p ( 1 +xm)- dexp( - βxα)F(x;d,m,p,α,β)=1 1-(1 1+xm)-dexp⁡(-βxα)1 …

7 r  distributions  random-generation 

Sou autodidata sobre probabilidade usando R, modelos lineares e cálculos de probabilidade. Atualmente, estou preso em como comparar duas previsões de um modelo. Os dados que estou usando são baixados (grátis) a partir daqui: wmbriggs.com/public/sat.csv df <- read.csv("sat.csv") # Load data lm <- lm(cgpa~hgpa+sat+ltrs,data=df) # model to predict College GPA …

7 r  regression  predictive-models  interpretation 

Eu sou um estudante de estatística de doutorado. Estou trabalhando com um conjunto de dados de contagem. É a contagem de usuários envolvidos em uma conversa de bate-papo em tempo real. O número de usuários varia de 1 a 6 e existem aproximadamente 300 dados no conjunto. Minha motivação inicial …

7 r  distributions  modeling  poisson-distribution 

Antecedentes: Imagine uma pizza cortada em 8 fatias. [ Em cada extremidade reta da fatia, insiro um ímã com polaridades opostas voltadas para fora. Se eu separar esses componentes, evitar sacudi-los e sacudi-los, eles deverão formar uma pizza cheia. Agora, se eu colocar uma fatia adicional (mesmas dimensões, 1/8 de …

7 r  confidence-interval  generalized-linear-model  multivariate-analysis  multinomial 

Para uma apresentação, tenho que visualizar dados tridimensionais. Eu deveria visualizá-los no "estilo de um gráfico de dispersão". As primeiras idéias de maio foram Um gráfico de dispersão tridimensional Uma matriz de dispersão Redução de dimensionalidade (PCA) e, posteriormente, um gráfico de dispersão bidimensional Quais são as alternativas para esses …

7 r  data-visualization  scatterplot 

Qual é a diferença entre normalizando as variáveis ​​e fazendo PCA; usando a scale=TRUEopção (sem normalizar as variáveis) na prcompfunção em R?

7 r  pca 

Eu tenho um modelo de regressão logística binária derivada passo a passo. Eu usei a calibrate(, bw=200, bw=TRUE)função no rmspacote em R para estimar sua calibração futura. A saída é fornecida abaixo e mostra a estimativa da curva de calibração corrigida por overfitting de autoinicialização para o modelo de logística …

7 r  logistic  multiple-regression  regression-strategies  calibration 

Parece um problema comum, mas não consigo encontrar uma solução. Eu tenho um conjunto de observações binárias e dois modelos diferentes, cada um com previsões para cada observação. Eu quero comparar a calibração dos modelos. Existem várias abordagens para comparar a discriminação desses modelos (ou seja, consulte o roc.test no …

7 r  machine-learning  logistic  calibration  scoring-rules