data(WhittenWilliams)
# expect warning urging to use 'complete.case=FALSE'
out <- panelAR(milex_gdp~lag_milex_gdp+GOV_rl+gthreat+GOV_min+GOV_npty+election_yr+
lag_real_GDP_gr+cinclag+lag_alliance+lag_cinc_ratio+lag_us_change_milex_gdp,
data=WhittenWilliams, panelVar="ccode", timeVar="year", autoCorr="psar1",
panelCorrMethod="pcse", complete.case=TRUE)
summary(out)
summary(out)$rho # psar1 coefficients
summary(out)$Sigma # panel covariances
summary(out)$wald # results of Wald test
Run the code above in your browser using DataLab