Perguntas com a marcação «r»

10
Ajustando a distribuição log-normal em R vs. SciPy

Eu ajustei um modelo lognormal usando R com um conjunto de dados. Os parâmetros resultantes foram: meanlog = 4.2991610 sdlog = 0.5511349 Gostaria de transferir esse modelo para o Scipy, que nunca usei antes. Usando o Scipy, consegui obter uma forma e escala de 1 e 3,1626716539637488e + 90 -...

10
Fatores de Bayes com antecedentes impróprios

Eu tenho uma pergunta sobre a comparação de modelos usando fatores Bayes. Em muitos casos, os estatísticos estão interessados ​​em usar uma abordagem bayesiana com anteriores impróprios (por exemplo, alguns anteriores de Jeffreys e anteriores de referência). Minha pergunta é: nos casos em que a...

10
Por que Anova () e drop1 () forneceram respostas diferentes para os GLMMs?

Eu tenho um GLMM do formulário: lmer(present? ~ factor1 + factor2 + continuous + factor1*continuous + (1 | factor3), family=binomial) Quando uso drop1(model, test="Chi"), obtenho resultados diferentes dos que utilizo Anova(model, type="III")na embalagem do carro ou summary(model). Estes dois...

10
REML vs ML stepAIC

Sinto-me impressionado depois de tentar pesquisar na literatura sobre como executar minha análise de modelo misto, seguindo-a com o uso da AIC para selecionar o melhor modelo ou modelos. Não acho que meus dados sejam tão complicados, mas estou procurando uma confirmação de que o que fiz foi correto...

10
Diferenças R e EViews nas estimativas de AR (1)

O principal problema é: não consigo obter estimativas de parâmetros semelhantes com EViews e R. Por razões que eu não conheço, preciso estimar parâmetros para determinados dados usando EViews. Isso é feito escolhendo a opção NLS (mínimos quadrados não lineares) e usando a seguinte...

10
Escolha do número de componentes principais a serem retidos

Um método que me foi sugerido é examinar um gráfico de seixos e verificar "cotovelo" para determinar o número correto de PCs a serem usados. Mas se o gráfico não estiver claro, R tem um cálculo para determinar o número? fit <- princomp(mydata,