Eu uso o pacote "boot" para calcular um valor p aproximado de inicialização de dois lados, mas o resultado está muito longe do valor p de usar t.test. Não consigo descobrir o que fiz de errado no meu código R. Alguém pode me dar uma dica para isso time = c(14,18,11,13,18,17,21,9,16,17,14,15,...
A maioria das distribuições padrão em R tem uma família de comandos - pdf / pmf, cdf / cmf, quantil, desvios aleatórios (por exemplo: dnorm, pnorm, qnorm, rnorm). Sei que é fácil o suficiente usar alguns comandos padrão para reproduzir essas funções para distribuições uniformes discretas, mas já...
Encontrei uma fórmula para o pseudo no livro Estendendo o modelo linear com R, Julian J. Faraway (p. 59).R2R2R^2 1 - ResidualDevianceNullDeviance1−ResidualDevianceNullDeviance1-\frac{\text{ResidualDeviance}}{\text{NullDeviance}} . Essa é uma fórmula comum para o pseudo para...
É útil estudar o código de análise de dados de especialistas. Recentemente, estive lendo o github e há várias pessoas compartilhando código de análise de dados por lá. Isso inclui alguns pacotes R (que obviamente estão disponíveis diretamente no CRAN), mas também vários exemplos de pesquisa...
Estou usando os modelos R (3.1.1) e ARIMA para previsão. Gostaria de saber qual deve ser o parâmetro "frequency", atribuído na ts()função , se estiver usando dados de séries temporais que sejam: separados por minutos e estão espalhados por 180 dias (1440 minutos / dia) separados por segundos e...
Olhe para essa foto: Se extrairmos uma amostra da densidade vermelha, espera-se que alguns valores sejam menores que 0,25, ao passo que é impossível gerar essa amostra a partir da distribuição azul. Como conseqüência, a distância Kullback-Leibler da densidade vermelha à densidade azul é infinito....
Acabei de me deparar com este artigo , que descreve como calcular a repetibilidade (também conhecida como confiabilidade, também conhecida como correlação intraclasse) de uma medição via modelagem de efeitos mistos. O código R seria: #fit the model fit = lmer(dv~(1|unit),data=my_data) #obtain the...
Ao ajustar um modelo de regressão, o que acontece se as suposições das saídas não forem atendidas, especificamente: O que acontece se os resíduos não forem homocedásticos? Se os resíduos mostrarem um padrão crescente ou decrescente na plotagem Residuais vs. Ajustados. O que acontece se os...
Entendo que essa questão é bastante ampla, mas me pergunto quais devem ser os pontos decisivos na decisão de criar (ou não) um novo pacote para R. Para ser mais específico, acrescentaria que a questão não é sobre os motivos para use R em si mesmo, mais sobre a decisão de compilar vários scripts e...
Como exemplo, considere o ChickWeightconjunto de dados em R. A variação obviamente aumenta com o tempo, portanto, se eu usar uma regressão linear simples como: m <- lm(weight ~ Time*Diet, data=ChickWeight) Minhas perguntas: Quais aspectos do modelo serão questionáveis? Os problemas estão...
Preciso usar variáveis binárias (valores 0 e 1) em médias k. Mas k-means funciona apenas com variáveis contínuas. Eu sei que algumas pessoas ainda usam essas variáveis binárias em k-means, ignorando o fato de que k-means é projetado apenas para variáveis contínuas. Isso é inaceitável para...
Talvez essa pergunta seja ingênua, mas: Se a regressão linear está intimamente relacionada ao coeficiente de correlação de Pearson, existem técnicas de regressão intimamente relacionadas aos coeficientes de correlação de Kendall e
Bloqueado . Esta pergunta e suas respostas estão bloqueadas porque a questão está fora do tópico, mas tem um significado histórico. No momento, não está aceitando novas respostas ou interações. Eu queria saber se é possível fazer computação simbólica em R? Por exemplo,...
Eu tenho uma matriz de números de ponto flutuante 336x256 (336 genomas bacterianos (colunas) x 256 frequências tetranucleotídicas normalizadas (linhas), por exemplo, todas as colunas somam 1). Eu obtenho bons resultados quando executo minha análise usando a análise de componentes principais....
Esta é uma pergunta sobre terminologia. Um "anterior vago" é o mesmo que um anterior não informativo ou existe alguma diferença entre os dois? Minha impressão é de que são iguais (de procurar vagas e não informativas juntas), mas não tenho
Considere um Jeffreys anterior onde p(θ)∝|i(θ)|−−−−√p(θ)∝|i(θ)|p(\theta) \propto \sqrt{|i(\theta)|}, ondeiiié a informação de Fisher. Eu continuo vendo esse prior sendo mencionado como um prior não informativo, mas nunca vi um argumento sobre por que ele não é informativo. Afinal, não é um...
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 no ano passado . Eu tenho um gráfico de dispersão. Como...
Estou surpreso que isso não tenha sido feito antes, mas não consigo encontrar a pergunta em stats.stackexchange. Esta é a fórmula para calcular a variação de uma amostra normalmente distribuída: ∑(X−X¯)2n−1∑(X−X¯)2n−1\frac{\sum(X - \bar{X}) ^2}{n-1} Esta é a fórmula para calcular o erro...
Estou refletindo sobre a discussão em torno desta questão e, em particular, o comentário de Frank Harrell de que a estimativa de variação em um modelo reduzido (ou seja, do qual várias variáveis explicativas foram testadas e rejeitadas) deve usar os graus de liberdade generalizados de Ye . O...
Costumo ler sobre uma função ser 'altamente não linear'. No meu entendimento, há "linear" e "não linear", então o que é isso 'altamente'? Existe uma diferença formal de não linear? Como é