Perguntas com a marcação «r»

O modelo de Bradley-Terry-Luce (BTL) afirma que , em que é a probabilidade de o objeto ser considerado "melhor", mais pesado, etc, do que o objecto i , e \ delta_i , e \ delta_j são parâmetros.p i j j ipj i= l o geu t- 1( δj- δEu)pjEu=euogEut-1 1(δj-δEu)p_{ji} …

9 r  logistic  categorical-data  paired-data  bradley-terry-model 

Sei que essa é provavelmente uma pergunta básica ... Mas não acho a resposta. Estou ajustando um GLM com uma família Poisson e tentei dar uma olhada nas previsões, no entanto, o deslocamento parece ser levado em consideração: model_glm=glm(cases~rhs(data$year,2003)+lhs(data$year,2003), offset=(log(population)), data=data, subset=28:36, family=poisson()) predict (model_glm, type="response") Recebo casos, não taxas …

9 r  generalized-linear-model  prediction  offset 

O @whuber demonstrou como simular resultados multivariados ( , e ) por um ponto no tempo.y 2y1 1y1y_1y2y2y_2y3y3y_3 Como sabemos, dados longitudinais geralmente ocorrem em estudos médicos. Minha pergunta é como simular resultados multivariados de medidas repetidas em R? Por exemplo, medimos repetidamente , e em 5 pontos no tempo …

9 r  repeated-measures  simulation 

Estou precisando de um pacote que possa me dar a equação para um modelo SVM linear. Atualmente estou usando o e1071 assim: library(e1071) m = svm(data, labels, type='C', kernel='linear', cost=cost, probability=FALSE, scale=scale) w = t(m$coefs) %*% data[m$index,] #Weight vector b = -model$rho #Offset No entanto, não tenho certeza de como …

9 r  svm  e1071 

Tenho problemas ao usar as funções cor()e cor.test(). Eu só tenho duas matrizes (apenas valores numéricos e o mesmo número de linhas e colunas) e quero ter o número de correlação e o valor p correspondente. Quando uso cor(matrix1, matrix2), obtenho os coeficientes de correlação para todas as células. Eu …

9 r  correlation 

Estou usando o pacote kernlab no R para criar um SVM para classificar alguns dados. O SVM está funcionando bem, pois fornece 'previsões' de uma precisão decente, no entanto, minha lista de variáveis ​​de entrada é maior do que eu gostaria e não tenho certeza quanto à importância relativa das …

9 r  machine-learning  svm  genetic-algorithms 

Digamos que você tenha dados de sobrevivência como este: obs <- data.frame( time = c(floor(runif(100) * 30), floor((runif(100)^2) * 30)), status = c(rbinom(100, 1, 0.2), rbinom(100, 1, 0.7)), group = gl(2,100) ) Para executar um teste padrão de classificação de log, pode-se usar survdiff(Surv(time, status) ~ group, data = obs, …

9 r  survival  wilcoxon-signed-rank 

Como um passatempo paralelo, tenho explorado séries temporais de previsão (em particular, usando R). Para os meus dados, tenho o número de visitas por dia, para todos os dias, voltando quase quatro anos. Nestes dados, existem alguns padrões distintos: De segunda a sexta tem muitas visitas (a mais alta de …

9 r  time-series  seasonality  forecasting 

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á 3 anos . Estou muito interessado em aprender a criar relatórios recorrentes a partir do meu código R …

9 r 

Estou procurando um exemplo de código simples de como executar um Filtro de Partículas em R. O pacote pomp parece oferecer suporte ao bit matemático do espaço de estado, mas os exemplos são um pouco complicados de seguir programaticamente para um desenvolvedor OO simples como eu, particularmente como carregar os …

9 r 

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á 8 meses . O problema mencionado nesta pergunta foi corrigido na versão 1.7.3 do pacote R glmnet. Estou …

9 r  multinomial  glmnet 

Criei minha própria versão ligeiramente aprimorada do termplot que utilizo neste exemplo, você pode encontrá-la aqui . Eu publiquei anteriormente no SO, mas quanto mais penso nisso, acredito que isso provavelmente esteja mais relacionado à interpretação do modelo de riscos proporcionais de Cox do que com a codificação real. O …

9 r  survival  cox-model 

Eu tenho alguns dados de séries temporais em que a variável medida é números inteiros positivos discretos (contagens). Quero testar se há uma tendência ascendente ao longo do tempo (ou não). A variável independente (x) está no intervalo de 0 a 500 e a variável dependente (y) está no intervalo …

9 r  regression  python 

Se eu tenho uma variável com 4 níveis, em teoria eu preciso usar 3 variáveis ​​fictícias. Na prática, como isso é realmente realizado? Uso 0-3, uso 1-3 e deixo o 4 em branco? Alguma sugestão? NOTA: Eu vou trabalhar na R. ATUALIZAÇÃO: O que aconteceria se eu apenas usasse uma …

9 r  regression  categorical-data  categorical-encoding 

Estou tentando visualizar alguns dados do consumidor, que tem quatro categorias. Os usuários podem alternar entre diferentes categorias. Gostaria de visualizar os últimos três ou quatro comutadores para cada indivíduo. Então começamos com um gráfico com uma coluna com 4 proporções empilhadas. Depois disso, teríamos 16, conforme cada categoria se …

9 r  data-visualization  sankey-diagram