Estou interessado em saber como os Filtros Kalman podem ser usados para atribuir valores ausentes nos Dados de Séries Temporais. Também é aplicável se faltarem alguns pontos consecutivos no tempo? Não consigo encontrar muito sobre esse assunto. Quaisquer explicações, comentários e links são...
Eu construí um modelo de glm em R e o testei usando um grupo de teste e treinamento, por isso estou confiante de que funciona bem. Os resultados de R são: Coefficients: Estimate Std. Error t value Pr(>|t|) (Intercept) -2.781e+00 1.677e-02 -165.789 < 2e-16 *** Coeff_A 1.663e-05 5.438e-06...
Estou tentando especificar uma fórmula para um modelo linear de efeitos mistos (com lme4) para o meu design experimental, mas não tenho certeza se estou fazendo certo. O design: basicamente estou medindo um parâmetro de resposta nas plantas. Eu tenho 4 níveis de tratamento e 2 níveis de irrigação....
Estou usando o RandomForestpacote R e estou confuso sobre como interpretar os valores do eixo Y em seus gráficos de dependência parcial. Os documentos da Ajuda afirmam que o gráfico é uma "representação gráfica do efeito marginal de uma variável na probabilidade de classe". No entanto, ainda estou...
Estou procurando o método ideal de binning (discretização) de uma variável contínua em relação a uma variável binária de resposta (destino) e com o número máximo de intervalos como parâmetro. exemplo: eu tenho um conjunto de observações de pessoas com variáveis "height" (numeral contínuo) e...
Eu tenho um quadro de dados que contém valores em 4 colunas: Por exemplo: ID, price, click count,rating O que eu gostaria de fazer é "dividir" esse quadro de dados em N grupos diferentes, onde cada grupo terá um número igual de linhas com a mesma distribuição de preço, contagem de cliques e...
Estou lendo o artigo da teoria de Doug Bates no pacote lme4 de R para entender melhor os detalhes dos modelos mistos e me deparei com um resultado intrigante que eu gostaria de entender melhor sobre o uso da máxima verossimilhança restrita (REML) para estimar a variação . Na seção 3.3, no critério...
Isso é semelhante ao Bootstrap: a estimativa está fora do intervalo de confiança Eu tenho alguns dados que representam contagens de genótipos em uma população. Quero estimar a diversidade genética usando o índice de Shannon e também gerar um intervalo de confiança usando o bootstrapping. Notei, no...
Estou lendo o livro de Kevin Murphy: Machine Learning - uma perspectiva probabilística. No primeiro capítulo, o autor está explicando a maldição da dimensionalidade e há uma parte que eu não entendo. Como exemplo, o autor declara: Considere que as entradas são distribuídas uniformemente ao longo...
Gostaria de saber como posso simular uma amostra de n vidas úteis de distribuição Weibull que incluem observações tipo I censuradas à direita. Por exemplo, vamos ter n = 3, forma = 3, escala = 1 e a taxa de censura = 0,15, e o tempo de censura = 0,88. Sei como gerar uma amostra Weibull, mas não sei...
Estou examinando a seção LAB §6.6 sobre Regressão de Ridge / Lasso no livro 'Uma Introdução à Aprendizagem Estatística com Aplicações em R' de James, Witten, Hastie, Tibshirani (2013). Mais especificamente, estou tentando aplicar o Ridgemodelo scikit-learn ao conjunto de dados 'Hitters' do pacote...
Eu identifiquei vários lugares nos livros onde o GLM é descrito com 5 distribuições (a saber, Gama, Gaussiana, Binomial, Gaussiana Inversa e Poisson). Isso também é exemplificado na função familiar em R. Ocasionalmente, encontro referências ao GLM, onde distribuições adicionais estão incluídas (...
Gostaria de simular um processo de excursão browniana (um movimento browniano condicionado sempre será positivo quando a em ). Como um processo de excursão browniana é uma ponte browniana condicionada a ser sempre positiva, eu esperava simular o movimento de uma excursão browniana usando uma ponte...
Estou tentando modelar as intensidades médias de parasitas que afetam um hospedeiro em R usando um modelo binomial negativo. Eu recebo 50 ou mais avisos que dizem: In dpois(y, mu, log = TRUE) : non-integer x = 251.529000 Como posso lidar com isso? Meu código fica assim: mst.nb =...
Quando soube dos antecedentes de Jeffrey na minha aula de inferência estatística, meus professores fizeram parecer que isso era interessante principalmente por razões históricas, e não porque alguém jamais as usaria. Então, quando fiz a análise de dados bayesiana, nunca nos pediram para usar os...
Estou lendo um artigo cujo método é totalmente baseado no teste da razão de verossimilhança. O autor diz que o teste LR contra alternativas unilaterais é o UMP. Ele procede alegando que "... mesmo quando não é possível demonstrar que [o teste LR] é uniformemente mais poderoso, o teste LR...
Estou analisando uma tabela 2x2 de um pequeno conjunto de dados de 30 pacientes. Estamos tentando retrospectivamente encontrar algumas variáveis que dão uma dica sobre qual tratamento escolher. As variáveis (obs normal / estranho) e decisão de tratamento (A / B) são de interesse especial e,...
A seguir, um trecho da Introdução às Estatísticas Bayesianas de Bolstad . Para todos os especialistas por aí, isso pode ser trivial, mas não entendo como o autor conclui que não precisamos fazer nenhuma integração para calcular a probabilidade posterior de algum valor de . Entendo a segunda...
Considere um modelo AR ( ) (assumindo média zero por simplicidade):ppp xt= φ1 1xt - 1+ … + Φpxt - p+ εtxt=φ1 1xt-1 1+…+φpxt-p+εt x_t = \varphi_1 x_{t-1} + \dotsc + \varphi_p x_{t-p} + \varepsilon_t O estimador OLS (equivalente ao estimador condicional de verossimilhança máxima) para é conhecido...
Minhas perguntas reais estão nos dois últimos parágrafos, mas para motivá-las: Se estou tentando estimar a média de uma variável aleatória que segue uma distribuição Normal com uma variação conhecida, li que colocar um uniforme antes da média resulta em uma distribuição posterior proporcional à...