Eu estou olhando para gerar 450 pontos de dados em R. Existem três conjuntos distintos 150 de cada um distribuídos em uma banda circular com raios diferentes (em 1, 2,8 e 5). Em particular, pretendo reproduzir o primeiro gráfico da p546 de Os elementos do aprendizado estatístico. Ficaria muito...
Eu calculei os valores de ICC para dois grupos e agora gostaria de comparar os valores de ICC para determinar se os grupos diferem em sua repetibilidade. Na literatura, as pessoas simplesmente usaram testes t para comparar a repetibilidade, mas não está claro para mim como fazer isso. Por exemplo,...
Existe alguma maneira de executar o LASSO com regressão binomial negativa em R? Estou executando uma regressão binomial negativa no meu conjunto de dados porque os dados são muito dispersos para impor a regressão de poisson. Enquanto isso, também estou enfrentando algum problema de...
Eu tenho a seguinte tabela em R df <- structure(list(x = structure(c(12458, 12633, 12692, 12830, 13369, 13455, 13458, 13515), class = "Date"), y = c(6080, 6949, 7076, 7818, 0, 0, 10765, 11153)), .Names = c("x", "y"), row.names = c("1", "2", "3", "4", "5", "6", "8", "9"), class =...
Eu tenho o seguinte modelo linear generalizado. O objeto glmDVé modelado como uma proporção de sucessos em relação ao total de tentativas. Os objetos x_isão variáveis contínuas. Como isso se parece na notação matemática? winp.glm = glm(glmDV ~ x1 + x2 + x3 + x4 + x5 + x6 + x7, data=myData,...
Eu tenho um conjunto de dados no qual estou tentando fazer regressão e falhando. A situação: Milhares de operadores de robôs de batalha estão travando batalhas entre si usando robôs de batalha. Alguns robôs de batalha são fortes e poderosos, e outros são fracos; os fortes vencem com mais...
Meu projeto atual pode exigir que eu construa um modelo para prever o comportamento de um determinado grupo de pessoas. o conjunto de dados de treinamento contém apenas 6 variáveis (id é apenas para fins de identificação): id, age, income, gender, job category, monthly spend em que monthly...
Eu tenho ~ 400 indivíduos e> 10k pontos no tempo cada (resultados da simulação) que gostaria de poder monitorar à medida que eles mudam ao longo do tempo. Traçar todos os indivíduos é muito confuso, traçar média + -sd, min / max ou quantis é pouca informação para o meu gosto. Estou imaginando o...
Estou usando o método arima do pacote de estatísticas R com minha série temporal de 17376 elementos. Meu objetivo é obter o valor do critério AIC, observei no meu primeiro teste isso: ts <- arima(serie[,1], order = c(2,1,1), seasonal = list(order=c(2,0,1),period = 24), method = "CSS",...
Quais são algumas abordagens típicas para lidar com dados incompletos no filtro Kalman? Eu estou falando sobre a situação em que alguns elementos do vetor observado ytyty_t estão faltando, distinto do caso em que um inteiro observado vector ytyty_t é perdido. Outra maneira de pensar sobre isso...
Na análise de dados bayesiana , capítulo 13, página 317, segundo parágrafo completo, nas aproximações modal e distributiva, Gelman et al. Escreva: Se o plano é resumir a inferência pelo modo posterior de [o parâmetro de correlação em uma distribuição normal bivariada], substituiríamos a...
Eu tenho um problema com a interpretação de interações de 2 e 3 vias no lmer. Meu DV é a altura, que é uma variável contínua. Todos os IVs são variáveis categóricas. O primeiro fator é animal, seja rato ou leão. O segundo fator é o sexo, masculino ou feminino. O terceiro fator é a cor: vermelho,...
Estou usando lmerTestpara executar modelos mistos lineares (LMM) para obter os valores-p. No entanto, nos artigos que escrevi, eles mostrarão ambos lme4e lmerTest. Então, não tenho certeza se também citei lme4quando uso apenas lmerTestpara executar modelos LMM. voce tem alguma ideia sobre...
Existem vários métodos para fazer previsões de séries temporais equidistantes (por exemplo, Holt-Winters, ARIMA, ...). No entanto, atualmente estou trabalhando no seguinte conjunto de dados espaçados irregulares, que possui uma quantidade variável de pontos de dados por ano e sem intervalos...
Depois de fazer algumas pesquisas sobre o tópico, notei um déficit surpreendente de pacotes de inferência e bibliotecas que dependem de métodos de passagem de mensagens ou otimização para Python e R. Que eu saiba, esses métodos são extremamente úteis. Por exemplo, apenas para uma propagação de...
Estou trabalhando em modelos de dados de contagem inflada a zero usando o psclpacote. Estou apenas me perguntando por que não há desenvolvimento de modelos para modelos de dados de contagem com um inflado! Também por que não há desenvolvimento de modelos de dados bimodais, digamos com zero e 2...
R em ação (Kabacoff, 2011) sugere a seguinte rotina para testar sobredispersão em uma regressão logística: Ajuste a regressão logística usando distribuição binomial: model_binom <- glm(Species=="versicolor" ~ Sepal.Width, family=binomial(), data=iris) Ajuste a regressão logística usando...
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á 2 anos . minha série temporal é obviamente...
Sou um pouco novato em R e seleção de recursos e tentei o pacote Boruta para selecionar (diminuir) meu número de variáveis (n = 40). Eu pensei que esse método também levasse em conta a possível correlação entre variáveis, no entanto, duas (das 20 variáveis selecionadas) são altamente...
Suponha que eu tenha uma função como: f <- function(x){ exp(x) / (1 + exp(x)) } ele deve funcionar com qualquer valor real de x, mas na verdade ele retorna NaN quando x é 710 ou maior. Gostaria de saber qual é a maneira correta de lidar com esse problema. Sei que é fácil fazê-lo retornar...