# Monte Carlo of FE estimator # # Set up follows Judson & Owen 1999 # except short-run effects are held fixed # and long-run effects allowed to vary # across scenarios # # High beta scenario (beta = 1) # Large N (N=1000) # # Christopher Adolph faculty.washington.edu/cadolph # 27 June 2015 rm(list=ls()) library(plm) # for FE estimator library(tile) # for graphing tools library(RColorBrewer) # for nice colors source("helperNickell.R") # all MC functions and code # User settings useSaved <- FALSE makeGraphs <- TRUE file <- "highbetaLargenMC" comparisonFile <- NULL resultDescript <- "beta=1.0 N=1000" compDescript <- " " compBeta <- NA modelNumber <- 1 # Monte Carlo simulation parameters sims <- 100 N <- 1000 beta <- 1 sigma <- 1 sigmaX <- 2 phiX <- 0.5 sigmaAlpha <- 1 # Key varied MC parameters T <- c(2, 3, 4, 5, 6, 7, 8, 9, 10, 12, 15, 20, 25, 30, 40, 50, 75, 100, 150, 200) phi <- c(0.2, 0.5, 0.8, 0.95) # Seed comes from digits from runif(1) set.seed(1037494) # Run Monte Carlo results <- runMC(file, useSaved, comparisonFile, resultDescript, compDescript, compBeta, modelNumber, makeGraphs, sims, N, beta, sigma, sigmaX, phiX, sigmaAlpha, T, phi)