Tentando replicar os resultados do artigo publicado recentemente, Aghion, Philippe, John Van Reenen e Luigi Zingales. 2013. "Inovação e Propriedade Institucional". American Economic Review, 103 (1): 277-304. (O código de dados e dados está disponível em
Sabemos que se um estimador é um estimador imparcial de teta e se sua variância tende a 0 como n tende ao infinito, é um estimador consistente para teta. Mas essa é uma condição suficiente e não necessária. Estou procurando um exemplo de um estimador que seja consistente, mas cuja variação não...
Estou tentando ajustar um modelo log-linear a um grande número de variáveis dos dados da pesquisa. Existem alguns motivos pelos quais pode ser preferível ajustar as regressões logísticas a esses dados. Várias autoridades sugerem que estes são equivalentes. No entanto, tenho alguns motivos para...
Eu sei que R não é particularmente útil para analisar grandes conjuntos de dados, uma vez que R carrega todos os dados na memória, enquanto algo como o SAS faz análise sequencial. Dito isto, existem pacotes como o bigmemory que permitem aos usuários realizar análises de grandes dados (análise...
Usando o excelente pacote de previsões de Rob Hyndman, deparei-me com a necessidade de não apenas ter intervalos de previsão, mas simular vários caminhos futuros, considerando observações passadas de uma série temporal com sazonalidades complexas. Há algo para séries temporais menos complexas com...
Eu tenho um conjunto variável de respostas que são expressas como um intervalo, como a amostra abaixo. > head(left) [1] 860 516 430 1118 860 602 > head(right) [1] 946 602 516 1204 946 688 onde esquerda é o limite inferior e direita é o limite superior da resposta. Quero estimar os...
Estou tentando calcular esta distribuição posterior: ( θ | - ) = ∏ni = 1pyEuEu( 1 - pEu)1 - yEu∑todosθ , pEu| θ∏ni = 1pyEuEu( 1 - pEu)1 -
Na previsão de séries temporais usando vários modelos como AR, MA, ARMA, etc, geralmente focamos na modelagem dos dados na mudança de tempo. Mas quando temos duas séries temporais em que o coeficiente de correlação de Pearson mostra que elas são altamente correlacionadas, é possível modelar a...
Preciso calcular o CDF Dirichlet , mas só consigo encontrar implementações do PDF . Vocês conhecem alguma biblioteca (de preferência a R)
Eu tenho alguns conjuntos de dados EEG que estou testando em duas classes. Posso obter uma taxa de erro decente do LDA (as distribuições condicionais de classe não são gaussianas, mas têm caudas semelhantes e separação suficientemente boa) e, portanto, quero plotar o ROC do preditor LDA contra...
Esta é a pergunta de um iniciante em um exercício na “Computação Bayesiana com R” de Jim Albert. Observe que, embora isso possa ser um dever de casa, no meu caso não é, pois estou aprendendo métodos bayesianos em R porque acho que posso usá-lo em minhas análises futuras. De qualquer forma, embora...
Escrevo esta pergunta com referência a um exemplo na p138-142 do seguinte documento: ftp://ftp.software.ibm.com/software/analytics/spss/documentation/amos/20.0/en/Manuals/IBM_SPSS_Amos_User_Guide.pdf . Aqui estão figuras ilustrativas e uma tabela: Entendo que a variável latente não possui...
Quando executo este código: require(nlme) a <- matrix(c(1,3,5,7,4,5,6,4,7,8,9)) b <- matrix(c(3,5,6,2,4,6,7,8,7,8,9)) res <- lm(a ~ b) print(summary(res)) res_gls <- gls(a ~ b) print(summary(res_gls)) Eu obtenho os mesmos coeficientes e a mesma significância estatística nos...
Eu gostaria de fazer análises de poder para testes de hipóteses de (não) igualdade de proporções nas quais as proporções são muito pequenas. Eu gostaria de fazer isso sem usar aproximações normais (ou Poisson) da distribuição binomial. Existem vários tipos gerais de questões de poder que eu...
Por favor, permita-me fazer uma pergunta básica. Entendo a mecânica de Naive Bayes para variáveis discretas e posso refazer os cálculos "à mão". (código de HouseVotes84 todo o caminho abaixo). No entanto - estou lutando para ver como a mecânica funciona para variáveis contínuas (código de...
Eu tenho um parâmetro que fica entre . Digamos que eu possa executar um experimento e obter , em que é um gaussiano padrão. O que eu preciso é de uma estimativa de que seja 1) imparcial 2) quase certamente limitada. O requisito (2) é crucial para mim.θθ\thetaθ = θ + w w
Estou tentando executar uma regressão de Cox em um exemplo de conjunto de dados de 2.000.000 de linhas da seguinte maneira, usando apenas R. Esta é uma tradução direta de um PHREG no SAS. A amostra é representativa da estrutura do conjunto de dados original. ## library(survival) ### Replace...
Estou trabalhando em um modelo que conta com uma função parametrizada feia que atua como uma função de calibração em uma parte do modelo. Usando uma configuração bayesiana, preciso obter informações preliminares não informativas para os parâmetros que descrevem minha função. Eu sei que, idealmente,...
Eu tenho 2 estudos analisando a resposta do paciente ao mesmo medicamento. O estudo 1 encontrou 10.000 genes expressos acima do plano de fundo e 500 deles são diferencialmente expressos e referidos como a assinatura da resposta ao medicamento. O estudo 2 encontrou 1.000 genes representando a...
Eu preciso de um único conjunto de dados imputados (por exemplo, para criar um manequim de grupo de países a partir dos dados de renda per capita do país imputado). O R oferece pacotes de pacotes para criar vários dados imputados (por exemplo, Amelia) e combinar resultados de vários conjuntos de...