# 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 # # Low beta scenario (beta = 0.2) # 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 <- "lowbetaLargenMC" comparisonFile <- "highbetaLargenMC" resultDescript <- "beta=0.2 N=1000" compDescript <- "beta=1.0 N=1000" compBeta <- 1 modelNumber <- 2 # Monte Carlo simulation parameters sims <- 100 N <- 1000 beta <- 0.2 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)