Eu tenho uma amostra de dados gerados Rpor rnorm(50,0,1), então os dados obviamente assumem uma distribuição normal. No entanto, Rnão "conhece" essas informações distributivas sobre os dados. Existe um método para Restimar que tipo de distribuição minha amostra vem? Caso contrário, usarei a...
Fiz a seguinte pergunta como parte dos trabalhos de casa: Projete e implemente um estudo de simulação para examinar o desempenho do bootstrap para obter intervalos de confiança de 95% na média de uma amostra univariada de dados. Sua implementação pode estar em R ou SAS. Os aspectos de...
Eu tenho vários conjuntos de dados relacionados. As correlações de Pearson entre pares deles são tipicamente definitivamente maiores que as correlações de lanceiros. Isso sugere que qualquer correlação é linear, mas pode-se esperar que, mesmo que pearson e spearman sejam os mesmos. O que significa...
Estou tentando analisar dados de medidas repetidas e estou lutando para fazê-lo funcionar R. Meus dados são essencialmente os seguintes, eu tenho dois grupos de tratamento. Cada sujeito de cada grupo é testado todos os dias e recebe uma pontuação (a porcentagem correta em um teste). Os dados estão...
θθ\thetaθ ∈ [0,1]θ∈[0,1] \theta \in [0,1]θ ∼U( 0 , 1 )θ∼U(0 0,1)\theta \sim
Em um post anterior, eu me perguntei como lidar com as pontuações EQ-5D . Recentemente, me deparei com a regressão logística quantílica sugerida por Bottai e McKeown, que apresenta uma maneira elegante de lidar com resultados limitados. A fórmula é simples: l o geu t ( y) = l o g( y- ym i nym a...
Nota: esta pergunta é um repost, pois minha pergunta anterior teve que ser excluída por razões legais. Ao comparar o PROC MIXED do SAS com a função lmedo nlmepacote no R, deparei-me com algumas diferenças bastante confusas. Mais especificamente, os graus de liberdade nos diferentes testes...
Eu usei o rjags para executar o MCMC em um modelo, especificado na linguagem JAGS. Existe uma boa maneira de extrair esse modelo e executar previsões com ele (usando as distribuições posteriores dos meus parâmetros)? Posso re-especificar o modelo em R e conectar os modos dos meus parâmetros...
Dado um conjunto de dados: x <-
Estou tentando testar várias abordagens funcionais de análise de dados. Idealmente, eu gostaria de testar o painel de abordagens que tenho em dados funcionais simulados. Eu tentei gerar FD simulado usando uma abordagem baseada em um somar ruídos gaussianos (código abaixo), mas as curvas resultantes...
Estou tendo problemas de desempenho usando o MCMCglmmpacote no R para executar um modelo de efeitos mistos. O código fica assim: MC1<-MCMCglmm(bull~1,random=~school,data=dt,family="categorical" , prior=list(R=list(V=1,fix=1), G=list(G1=list(V=1, nu=0))) , slice=T, nitt=iter, ,burnin=burn,...
Eu tenho a seguinte função de probabilidade: Prob=11+e−zProb=11+e−z\text{Prob} = \frac{1}{1 + e^{-z}} Onde z=B0+B1X1+⋯+BnXn.z=B0+B1X1+⋯+BnXn.z = B_0 + B_1X_1 + \dots + B_nX_n. Meu modelo parece Pr(Y=1)=11+exp(−[−3.92+0.014×(bid)])Pr(Y=1)=11+exp(−[−3.92+0.014×(bid)])\Pr(Y=1) = \frac{1}{1 +...
É melhor diferenciar uma série (supondo que seja necessária) antes de usar um Arima OU melhor para usar o parâmetro d no Arima? Fiquei surpreso com a diferença entre os valores ajustados, dependendo de qual rota é tomada com o mesmo modelo e dados. Ou estou fazendo algo
Os pesquisadores costumam usar duas medidas que possuem itens muito semelhantes e argumentam que medem coisas diferentes (por exemplo, "sempre me preocupo quando estou perto de carros"; "tenho medo de carros"). Vamos chamar as medidas hipotéticas de Medição e Ansiedade do Medo dos Carros na Escala...
Digamos que tenho os seguintes dados e estou executando um modelo de regressão: df=data.frame(income=c(5,3,47,8,6,5), won=c(0,0,1,1,1,0), age=c(18,18,23,50,19,39), home=c(0,0,1,0,0,1)) Por um lado, corro um modelo linear para prever a renda: md1 = lm(income ~ age + home + home,...
Estou usando o pacote randomForest no R e usando os dados da íris, a floresta aleatória gerada é uma classificação, mas quando uso um conjunto de dados com cerca de 700 recursos (os recursos são cada pixel em uma imagem de 28 x 28 pixels) e a coluna do rótulo é denominada label, o...
Estou tentando decompor uma matriz de covariância com base em um conjunto de dados esparsos / gappy. Estou percebendo que a soma do lambda (variação explicada), conforme calculada svd, está sendo amplificada com dados cada vez mais escassos. Sem lacunas, svde eigenproduz os mesmos resultados. Isso...
Estou tentando detectar valores anômalos em uma série temporal de dados climáticos com algumas observações ausentes. Pesquisando na web, encontrei muitas abordagens disponíveis. Dessas, a decomposição do stl parece atraente, no sentido de remover componentes de tendência e sazonais e estudar o...
Não consegui descobrir como executar a regressão linear em R in para um design de medida repetida. Em uma pergunta anterior (ainda sem resposta), me foi sugerido não usar, lmmas usar modelos mistos. Eu usei lmda seguinte maneira: lm.velocity_vs_Velocity_response <-
Eu estou tentando interpretar a saída de nls (). Eu li este post, mas ainda não entendo como escolher o melhor ajuste. Dos meus ajustes, tenho duas saídas: > summary(m) Formula: y ~ I(a * x^b) Parameters: Estimate Std. Error t value Pr(>|t|) a 479.92903 62.96371 7.622 0.000618 *** b...