Digamos que eu gere a probabilidade de um resultado com base em um determinado fator e trace a curva desse resultado. Existe alguma maneira de extrair a equação para essa curva de R? > mod = glm(winner~our_bid, data=mydat, family=binomial(link="logit")) > summary(mod) Call: glm(formula =...
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 de 0...
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...
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...
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...
Existe um princípio geral sobre se alguém deve calcular a correlação de pearson para duas variáveis aleatórias X e Y antes de realizar a transformação logarítmica ou depois? Existe um procedimento para testar o que é mais apropriado? Eles produzem valores semelhantes, mas diferentes, pois a...
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...
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...
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...
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...
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...
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...
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
Temos um conjunto de dados com duas covariáveis e uma variável de agrupamento categórica e queremos saber se existem diferenças significativas entre a inclinação ou a interceptação entre as covariáveis associadas às diferentes variáveis de agrupamento. Usamos anova () e lm () para comparar os...
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...
Eu tenho usado o Root Mean Squared Error(RMSE) para medir a precisão dos valores previstos usando um modelo. Entendo que o valor retornado está usando as unidades das minhas medidas (em vez de uma porcentagem). No entanto, gostaria de citar meus valores como uma porcentagem. A abordagem adotada é...
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} = logit^{-1}(\delta_j - \delta_i)peu...
Meus dados são uma série temporal da população empregada, L, e o período, ano. n.auto=auto.arima(log(L),xreg=year) summary(n.auto) Series: log(L) ARIMA(2,0,2) with non-zero mean Coefficients: ar1 ar2 ma1 ma2 intercept year 1.9122 -0.9567 -0.3082 0.0254 -3.5904 0.0074 s.e. NaN NaN NaN NaN...
Estou lendo A. Agresti (2007), Introdução à análise de dados categóricos , 2ª. edição e não tenho certeza se entendi corretamente este parágrafo (p.106, 4.2.1) (embora deva ser fácil): Na Tabela 3.1 sobre ronco e doenças cardíacas no capítulo anterior, 254 indivíduos relataram ronco todas as...
Fiquei me perguntando se alguém tinha experiência usando a função de ratos, como descrito em ratos: Imputação multivariada por equações encadeadas em R (JSS 2011 45 (3))? Eu tenho um conjunto de dados com um número de variáveis, cada uma com diferentes graus de dados ausentes. Minha pergunta...