Estou tentando replicar o que a função dfbetas()faz R . dfbeta() não é um problema ... Aqui está um conjunto de vetores: x <- c(0.512, 0.166, -0.142, -0.614, 12.72) y <- c(0.545, -0.02, -0.137, -0.751, 1.344) Se eu encaixar dois modelos de regressão da seguinte maneira: fit1 <- lm(y ~...
Eu tenho uma pergunta sobre modelagem de séries temporais em R. meus dados consistem na seguinte matriz: 1 0.03333333 0.01111111 0.9555556 2 0.03810624 0.02309469 0.9387991 3 0.00000000 0.03846154 0.9615385 4 0.03776683 0.03119869 0.9310345 5 0.06606607 0.01201201 0.9219219 6 0.03900325 0.02058505...
Posso criar e implementar modelos clássicos de ML em conjuntos tradicionais de treinamento / teste em R, mas e se um parceiro quiser obter esse modelo para implementar seu próprio (qualquer tipo de) sistema? Salvar e enviar a estrutura do modelo R não ajuda, é claro; e descobrir o mecanismo de...
Se o meu prior é modelado como uma distribuição de probabilidade contínua, digamos, uma distribuição beta distorcida para refletir meu viés em relação a determinados modelos, como posso calcular a probabilidade posterior? O desafio para mim é calcular a probabilidade de um determinado modelo, uma...
[Em questões recentes, eu estava pensando em gerar vetores aleatórios em R e queria compartilhar essa "pesquisa" como uma sessão de perguntas e respostas independente em um ponto específico.] A geração de dados aleatórios com correlação pode ser feita usando a decomposição de Cholesky da matriz de...
Eu tenho um conjunto de dados de 17 pessoas, classificando 77 declarações. Eu quero extrair componentes principais em uma matriz de correlação transposta de correlações entre pessoas (como variáveis) entre instruções (como casos). Eu sei, é estranho, chama-se Q Methodology . Quero ilustrar como o...
Usando o teste T de Student, o T-Critical é calculado através de: t=X¯−μ0s/n√t=X¯−μ0s/nt = \frac{\bar{X} - \mu_{0}}{s / \sqrt{n}} Olhando para o artigo da Wikipedia sobre a Estimação imparcial do desvio padrão, há a seção Resultado para a distribuição normal, que menciona um fator de correção...
Eu tenho duas amostras bem pequenas ( para cada) e recebo uma potência de teste t de duas amostras para elas 0,49, que é relativamente pequena (calculada por R ). No entanto, o teste t Welch Two Sample fornece um valor p de aqui e, portanto, a hipótese nula pode ser rejeitada.0,032n =...
Eu tenho a seguinte tabela de valores: 25 75 38 162 A razão de chances é 0,7037 e log (OR) é -0,3514. Para uma tabela de contingência com os valores a, b, c e d, a variação do log (OR) é dada por (1/a + 1/b + 1/c + 1/d) Como posso calcular o valor p do log (OR) desses dados em R (se é...
Como calcular média, variância, mediana, desvio padrão e modo de distribuição? Se eu gerar aleatoriamente números que formam a distribuição normal, especifiquei a média como m=24.2desvio padrão como sd=2.2: > dist = rnorm(n=1000, m=24.2, sd=2.2) Então eu posso fazer o...
A assimetria e curtose são definidas como: ζ4=E[(X-μ)4]ζ3= E[ ( X- μ )3]E[ ( X- μ )2]3 / 2= μ3σ3ζ3=E[(X-μ)3]E[(X-μ)2]3/2=μ3σ3\zeta_3 = \frac{E[(X-\mu)^3]}{E[(X-\mu)^2]^{3/2}} = \frac{\mu_3}{\sigma^3} ζ4= E[ ( X- μ )4]E[ ( X- μ )2]2= μ4σ4ζ4=E[(X-μ)4]E[(X-μ)2]2=μ4σ4\zeta_4 =...
Estou executando comparações post-hoc após um teste de Kruskal-Wallis. Estou usando o pacote PMCMR . > posthoc.kruskal.nemenyi.test( preference ~ instrument) Pairwise comparisons using Tukey and Kramer (Nemenyi) test with Tukey-Dist approximation for independent samples data: preference...
Eu tenho duas séries temporais (parâmetros de um modelo para homens e mulheres) e pretendo identificar um modelo ARIMA apropriado para fazer previsões. Minha série cronológica se parece com: A plotagem e o ACF mostram não estacionário (os picos do ACF cortam muito lentamente). Assim, uso...
Estou tentando determinar qual alfa usar na minha glmnetfunção, mas o arquivo de ajuda diz-me: Observe que o cv.glmnet NÃO procura valores para alfa. Um valor específico deve ser fornecido, caso contrário, alfa = 1 é assumido por padrão. Se os usuários também desejam validar cruzadamente o alfa,...
Suponha que eu tenha um conjunto de dados no qual a incerteza nas medições (provenientes da propagação de erros sistemáticos do aparato de medição) seja diferente para cada ponto. Se eu fizer uma regressão linear no conjunto de dados, como calculo a incerteza na inclinação? Eu gostaria de um...
Minha pergunta lida com árvores de classificação . Considere o seguinte exemplo do conjunto de dados Iris: Desejo selecionar manualmente o melhor preditor para a primeira divisão. De acordo com o algoritmo CART, o melhor recurso para fazer uma divisão é aquele que maximiza a diminuição da...
Estou tentando entender o Alpha de Cronbach intuitivamente. Qual é a ideia geral por trás dessa construção? Quais propriedades eles estavam tentando
Estou executando algumas otimizações com a implementação do BFGS do optim. A função objetivo é na verdade um algoritmo computacional, não apenas matemática. Descobri que quando adiciono uma penalidade de L1, as coisas diminuem um pouco. Por que isso pode ser? Existe algo no L1 que atrasa as coisas?...
A arimaxfunção no TSApacote é, até onde sei, o único Rpacote que se encaixará em uma função de transferência para modelos de intervenção. Falta uma função de previsão, que às vezes é necessária. A seguir, uma solução alternativa para esse problema, aproveitando o excelente forecastpacote? Os...
Eu tenho um conjunto de dados x, y que estou usando para criar uma floresta aleatória. Os dados x são um vetor de valores que inclui algumas NAs. Então, eu uso rfImputepara manipular os dados ausentes e criar uma floresta aleatória. Agora eu tenho uma nova observação invisível x (com um NA) e quero...