Eu tenho os seguintes dados localizados aqui . Estou tentando calcular o intervalo de confiança de 95% sobre a pureza média quando a porcentagem de hidrocarbonetos é 1,0. Em R, insiro o seguinte. > predict(purity.lm, newdata=list(hydro=1.0), interval="confidence", level=.95) fit lwr upr 1...
Estou analisando um conjunto de dados sobre comunidades entre marés. Os dados são cobertura percentual (de algas, cracas, mexilhões etc.) em quadrantes. Estou acostumado a pensar na análise de correspondência (CA) em termos de contagem de espécies e na análise de componentes principais (PCA) como...
Eu já vi na literatura publicada (e postada aqui) que a eficiência relativa assintótica do teste de posto assinado de Wilcoxon é de pelo menos 0,864 quando comparada ao teste t. Também ouvi dizer que isso se aplica apenas a grandes amostras, embora alguns livros não mencionem isso (o que é...
É possível plotar um gráfico de dependência parcial para exibir a probabilidade da classe e estimar os efeitos de um preditor para um modelo GBM ? Algo semelhante ao partialPlotdo randomForestpacote. De acordo com este artigo , uma plotagem parcial é possível com o gbm. Agradeço antecipadamente...
Em Dixon, Coles ( 1997 ), eles usaram a estimativa de máxima verossimilhança para os dois modelos independentes de Poisson modificados em (4.3) para modelar as pontuações no futebol. Estou tentando usar R para "reproduzir" os parâmetros alfa e beta, bem como os parâmetros de efeito doméstico (pág....
Estou interessado em encontrar um procedimento para simular dados que sejam consistentes com um modelo de mediação especificado. De acordo com a estrutura geral do modelo de equações estruturais lineares para testar modelos de mediação, primeiro descritos por Barron e Kenny (1986) e descritos em...
Estou trabalhando na estimativa da auto-regressão vetorial (VARs) e da função de resposta ao impulso (IRFs) com base nos dados do painel com 33 indivíduos em 77 trimestres. Como esse tipo de situação deve ser analisado? Quais algoritmos existem para esse fim? Eu preferiria realizar essas análises...
Estou tentando estimar o efeito de 2 medicamentos ( drug1, drug2) na probabilidade de um paciente cair ( event). Os pacientes podem cair mais de uma vez e podem ser colocados ou retirados dos medicamentos a qualquer momento. Minha pergunta é como os dados devem ser estruturados com relação ao...
Eu tenho um modelo de regressão logística (adequado via glmnet em R com regularização líquida elástica) e gostaria de maximizar a diferença entre verdadeiros positivos e falsos positivos. Para fazer isso, veio à mente o seguinte procedimento: Ajuste o modelo de regressão logística padrão Usando o...
Quero avaliar a precisão dos testes de normalidade em diferentes tamanhos de amostra em R (percebo que os testes de normalidade podem ser enganosos ). Por exemplo, para examinar o teste de Shapiro-Wilk, estou realizando a seguinte simulação (bem como plotando os resultados) e esperaria que, à...
Eu sou novo em estatística e estou tentando entender a diferença entre ANOVA e regressão linear. Estou usando R para explorar isso. Eu li vários artigos sobre por que a ANOVA e a regressão são diferentes, mas ainda são as mesmas, e como elas podem ser visualizadas, etc. Entendo que o ANOVA compara...
Não encontrei nenhuma literatura sobre a aplicação de florestas aleatórias no MNIST, CIFAR, STL-10, etc., então pensei em experimentá-las com o MNIST invariante, permutador . Em R , tentei: randomForest(train$x, factor(train$y), test$x, factor(test$y), ntree=500) Isso foi executado por 2 horas e...
O índice de estabilidade populacional quantifica a alteração de uma distribuição de uma variável comparando amostras de dados em dois períodos de tempo. É muito comumente usado para medir mudanças nas pontuações. É calculado da seguinte forma: 1) A amostra do período base é discretizada....
Um colega meu me enviou esse problema aparentemente fazendo rondas na internet: If $3 = 18, 4 = 32, 5 = 50, 6 = 72, 7 = 98$, Then, $10 =$ ? A resposta parece ser 200. 3*6 4*8 5*10 6*12 7*14 8*16 9*18 10*20=200 Quando eu faço uma regressão linear em R: data <-
Eu faço uma regressão linear usando a função R lm: x = log(errors) plot(x,y) lm.result = lm(formula = y ~ x) abline(lm.result, col="blue") # showing the "fit" in blue mas não se encaixa bem. Infelizmente, não consigo entender o manual. Alguém pode me apontar na direção certa para ajustar isso...
Estou olhando para uma planilha do Excel que afirma estar calculando o , mas não reconheço essa maneira de fazê-lo, e fiquei pensando se estou perdendo alguma coisa.χ2χ2\chi^2 Aqui estão os dados que está analisando: +------------------+----------+----------+ | Total Population | Observed |...
Quero ajustar um DLM com coeficientes variáveis no tempo, ou seja, uma extensão da regressão linear usual, yt=θ1 1+θ2x2yt=θ1 1+θ2x2y_t = \theta_1 + \theta_2x_2 . Eu tenho um preditor ( ) e uma variável de resposta ( ), capturas anuais de peixes no mar e no interior, respectivamente, de 1950 a...
Tenho muitas séries temporais com períodos: dia, semana ou mês. Com stl()função ou com loess(x ~ y), posso ver como são as tendências de determinadas séries temporais. Preciso detectar se a tendência das séries temporais está aumentando ou diminuindo. Como posso gerenciar isso? Tentei calcular...
O livro Elements of Statistical Learning (disponível em PDF online) discute o viés otimista (7.21, página 229). Ele afirma que o viés de otimismo é a diferença entre o erro de treinamento e o erro dentro da amostra (erro observado se coletarmos novos valores de resultado em cada um dos pontos de...
Estou trabalhando nos exemplos da Doing Bayesian Data Analysis de Kruschke , especificamente na ANOVA exponencial de Poisson no cap. 22, que ele apresenta como uma alternativa aos testes qui-quadrado freqüentes de independência para tabelas de contingência. Eu posso ver como obtemos informações...