Perguntas com a marcação «r»

Vou ver os vídeos do curso de aprendizado de máquina on - line gratuito de Andrew Ng em Stanford. Ele discute o Gradient Descent como um algoritmo para resolver funções de regressão linear e escrita no Octave para executá-lo. Presumivelmente, eu poderia reescrever essas funções em R, mas minha pergunta …

14 r  regression  machine-learning  gradient-descent 

Alguém pode me dizer como R estimar o ponto de interrupção em um modelo linear por partes (como um parâmetro fixo ou aleatório), quando eu também precisar estimar outros efeitos aleatórios? Incluí um exemplo de brinquedo abaixo que se encaixa em uma regressão de taco de hóquei / taco quebrado …

14 r  mixed-model  lme4-nlme  change-point  piecewise-linear 

Fechadas. Esta questão está fora de tópico . No momento, não está aceitando respostas. Deseja melhorar esta pergunta? Atualize a pergunta para que ela esteja no tópico de Validação cruzada. Fechado há 7 meses . É fácil encontrar uma área de cálculo de pacotes no ROC, mas existe um pacote …

14 r  precision-recall 

Estou prestes a mergulhar na aprendizagem R e meu projeto de aprendizagem envolverá a aplicação de regressão de efeitos mistos ou aleatórios a um conjunto de dados para desenvolver uma equação preditiva. Partilho da preocupação do escritor neste post. Como escolher a biblioteca nlme ou lme4 R para modelos de …

14 r  regression  random-effects-model 

Um erro padrão estimado é o desvio padrão σ ( θ ) de um estimador θ para um parâmetro θ .σ^(θ^)σ^(θ^)\hat \sigma(\hat\theta)θ^θ^\hat\thetaθθ\theta Por que o desvio padrão estimado dos resíduos é chamado de "erro padrão residual" (por exemplo, na saída da summary.lmfunção de R ) e não "desvio padrão residual"? …

14 r  standard-error  residuals  terminology 

Eu gostaria de realizar uma previsão com base em um modelo ARIMA de várias séries temporais com várias variáveis ​​exógenas. Como não sou tão habilidoso em relação a estatísticas nem ao RI, quero mantê-lo o mais simples possível (a previsão de tendências por três meses é suficiente). Eu tenho 1 …

14 r  time-series  arima 

Estou tentando entender como calcular o ponto de corte ideal para uma curva ROC (o valor no qual a sensibilidade e a especificidade são maximizadas). Estou usando o conjunto aSAHde dados do pacotepROC . A outcomevariável pode ser explicada por duas variáveis ​​independentes: s100be ndka. Usando a sintaxe do Epipacote, …

14 r  roc  sensitivity-analysis  sensitivity-specificity 

Eu sou novo na linguagem R. Eu gostaria de saber como simular a partir de um modelo de regressão linear múltipla que atenda às quatro suposições da regressão. ok .. obrigado. Digamos que eu queira simular os dados com base neste conjunto de dados: y<-c(18.73,14.52,17.43,14.54,13.44,24.39,13.34,22.71,12.68,19.32,30.16,27.09,25.40,26.05,33.49,35.62,26.07,36.78,34.95,43.67) x1<-c(610,950,720,840,980,530,680,540,890,730,670,770,880,1000,760,590,910,650,810,500) x2<-c(1,1,3,2,1,1,3,3,2,2,1,3,3,2,2,2,3,3,1,2) fit<-lm(y~x1+x2) summary(fit) então …

14 r  multiple-regression  simulation 

Preciso agrupar unidades em clusters para minimizar a soma de quadrados dentro do grupo (WSS), mas preciso garantir que cada um dos clusters contenha pelo menos m unidades. Alguma idéia se alguma das funções de cluster de R permitir agrupar em k clusters sujeitos a uma restrição mínima de tamanho …

14 r  clustering 

Eu tenho tentado discernir o que exatamente a saída "coef" e "(exp) coef" do coxph significa. Parece que o "(exp) coef" são comparações da primeira variável no modelo de acordo com o grupo designado no comando. Como a função coxph chega aos valores de "coef" e "(exp) coef"? Além disso, …

14 r  survival  interpretation 

Eu tenho um vetor de números que eu enviei aqui (... / code / MyData.Rdata) usando dput. Gostaria de obter o bca ci, então eu escrevi este código: my.mean <- function(dat, idx){ return (mean(dat[idx], na.rm = TRUE)) } boot.out<-boot(data=my.data, statistic = my.mean, R=1000) Mas quando executo o seguinte, recebo o …

14 r  bootstrap 

Com base em uma amostra de tempos de sobrevivência, gostaria de estimar a probabilidade de tempo de sobrevivência , para alguns t específicos , usando o estimador de Kaplan-Meier. É possível fazer isso ? Observe que t não é necessariamente um horário do evento.nnnttttttRttt

14 r  kaplan-meier 

Acabei de executar um GLM binomial negativo e esta é a saída: Call: glm.nb(formula = small ~ method + site + depth, data = size.dat, init.theta = 1.080668549, link = log) Deviance Residuals: Min 1Q Median 3Q Max -2.2452 -0.9973 -0.3028 0.3864 1.8727 Coefficients: Estimate Std. Error z value Pr(>|z|) …

14 r  categorical-data  generalized-linear-model  interpretation  negative-binomial 

Eu estou usando um script. É para registros principais. Eu tenho um quadro de dados que mostra as diferentes composições elementares nas colunas em uma determinada profundidade (na primeira coluna). Quero executar um PCA com ele e estou confuso sobre o método de padronização que tenho que escolher. Alguém de …

13 r  pca  normalization  compositional-data 

Douglas Bates afirma que os seguintes modelos são equivalentes "se a matriz de variância-covariância para efeitos aleatórios com valor vetorial tiver uma forma especial, chamada simetria composta" ( slide 91 desta apresentação ): m1 <- lmer(y ~ factor + (0 + factor|group), data) m2 <- lmer(y ~ factor + (1|group) …

13 r  anova  mixed-model  repeated-measures  lme4-nlme