As funções empíricas do CDF são geralmente estimadas por uma função de etapa. Existe uma razão para que isso seja feito dessa maneira e não usando uma interpolação linear? A função step possui propriedades teóricas interessantes que nos fazem preferir? Aqui está um exemplo dos dois: ecdf2 <-...
Eu gostaria de obter intervalos de confiança de 95% nas previsões de um nlmemodelo misto não linear . Como nada é fornecido para fazer isso dentro de nlmemim, eu me perguntava se é correto usar o método de "intervalos de previsão da população", conforme descrito no capítulo do livro de Ben Bolker...
Estou lendo o artigo Propagação de erros pelo método de Monte Carlo em cálculos geoquímicos, Anderson (1976) e há algo que não entendo direito. Considere alguns dados medidos e um programa que os processe e retorne um determinado valor. No artigo, este programa é usado para obter primeiro o melhor...
Percebi que, ao criar modelos aleatórios de regressão florestal, pelo menos em R, o valor previsto nunca excede o valor máximo da variável de destino vista nos dados de treinamento. Como um exemplo, veja o código abaixo. Estou construindo um modelo de regressão para prever com mpgbase nos...
Eu realmente tenho uma falta de entendimento de como o filtro de autoinicialização funciona. Conheço os conceitos aproximadamente, mas não consigo apreender certos detalhes. Esta pergunta é para eu esclarecer a desordem. Aqui vou usar esse algoritmo de filtro popular a partir de uma referência por...
Depois de um ano na pós-graduação, meu entendimento de "mínimos quadrados ponderados" é o seguinte: seja , seja algumas \ n p vezes matriz de design \ boldsymbol \ beta \ in \ mathbb {R} ^ p seja um vetor de parâmetro, \ boldsymbol \ epsilon \ in \ mathbb {R} ^ n seja um vetor de erro tal que \...
Eu tenho um conjunto de dados composto de proporções que medem o "nível de atividade" de girinos individuais, tornando os valores vinculados entre 0 e 1. Esses dados foram coletados contando o número de vezes que o indivíduo se moveu dentro de um determinado intervalo de tempo (1 para movimento, 0...
Polinômios ortogonais em um conjunto univariado de pontos são polinômios que produzem valores nesses pontos de maneira que seu produto escalar e a correlação pareada sejam zero. R pode produzir polinômios ortogonais com a função poli . A mesma função possui um polímero variante que produz...
Aqui está um pequeno exemplo: MyDf<-data.frame(x=c(1,2,3,4), y=c(1.2, .7, -.5, -3)) Agora com o base::lm: > lm(y~x, data=MyDf) %>% summary Call: lm(formula = y ~ x, data = MyDf) Residuals: 1 2 3 4 -0.47 0.41 0.59 -0.53 Coefficients: Estimate Std. Error t value Pr(>|t|)...
Eu tenho lido sobre o método Nyström para aproximação de kernel de baixo escalão. Este método é implementado no scikit-learn [1] como um método para projetar amostras de dados para uma aproximação de baixo escalão do mapeamento de recursos do kernel. Para o melhor de meu conhecimento, dado um...
Eu tenho esse enorme conjunto de dados com 2500 variáveis e 142 observações. Eu quero executar uma correlação entre a variável X e o resto das variáveis. Mas para muitas colunas, há entradas ausentes. Tentei fazer isso no R usando o argumento "pairwise-complete" ( use=pairwise.complete.obs) e...
No R, tenho uma amostra de 348 medidas e quero saber se posso assumir que ela é normalmente distribuída para testes futuros. Essencialmente, seguindo outra resposta da pilha , estou analisando o gráfico de densidade e o gráfico QQ
Eu tenho uma trajetória de um objeto em um espaço 2D (uma superfície). A trajetória é dada como uma sequência de (x,y)coordenadas. Sei que minhas medições são barulhentas e às vezes tenho discrepâncias óbvias. Então, eu quero filtrar minhas observações. Tanto quanto eu entendi o filtro Kalman, ele...
Fechado . Esta questão é baseada em opiniões . No momento, não está aceitando respostas. Deseja melhorar esta pergunta? Atualize a pergunta para que ela possa ser respondida com fatos e citações editando esta postagem . Fechado há 2 anos . Atualmente,...
Eu estive explorando várias ferramentas de previsão e constatei que os Modelos Aditivos Generalizados (GAMs) têm o maior potencial para essa finalidade. GAMs são ótimos! Eles permitem que modelos complexos sejam especificados de maneira muito sucinta. No entanto, essa mesma sucessão está me...
O seguinte problema foi publicado na página do Mensa International no Facebook: \quad\quad\quad\quad\quad\quad\quad\quad O post em si recebeu mais de 1000 comentários, mas não vou entrar em detalhes sobre o debate, pois sei que esse é o paradoxo da caixa de Bertrand e a resposta é . O que me...
Contexto Um grupo de cientistas sociais e estatísticos ( Benjamin et al., 2017 ) sugeriu recentemente que a taxa de falso positivo típica ( αα\alpha = 0,05) usada como limiar para determinar a "significância estatística" precisa ser ajustada para um limiar mais conservador ( αα\alpha = 0,005). Um...
Num artigo recente, Norton et al. (2018) [ 1 ] afirmam que[ 1 ][1]^{[1]} Odds ratio diferentes do mesmo estudo não podem ser comparados quando os modelos estatísticos que resultam em estimativas de odds ratio possuem variáveis explicativas diferentes, pois cada modelo possui um fator de escala...
Na análise bayesiana conjugada da distribuição gaussiana de Kevin Murphy , ele escreve que a distribuição preditiva posterior é p(x∣D)=∫p(x∣θ)p(θ∣D)dθp(x∣D)=∫p(x∣θ)p(θ∣D)dθ p(x \mid D) = \int p(x \mid \theta) p(\theta \mid D) d \theta onde são os dados nos quais o modelo está ajustado e são...
Suponha que tenhamos um conjunto (mensurável e adequadamente comportado) S⊆B⊂RnS⊆B⊂RnS\subseteq B\subset\mathbb R^n , onde BBB é compacto. Além disso, suponha que possamos extrair amostras da distribuição uniforme sobre BBB na medida de Lebesgue λ(⋅)λ(⋅)\lambda(\cdot) e que conhecemos a medida...