2016-11-23 86 views
0

我幾天前開始使用R studio,並且我正在努力計算一個VIF。這是的情況:在R中執行VIF測試

我有一個面板數據,並運行固定效應和隨機效應迴歸。我有一個從屬變量(New_biz_density)和2個獨立變量(Cost_to_start,Capital_requirements)。我想通過計算它們的方差膨脹因子來檢查我的兩個獨立變量是否表現出多重共線性,包括固定效應模型和隨機效應模型。

我已經安裝了一些包來執行BIF(遠景,汽車),但沒有設法做到這一點。有人知道該怎麼做嗎?

非常感謝!

這裏是我的腳本:

# install.packages("plm") 
library(plm) 

mydata<- read.csv("/Users/juliantabone/Downloads/DATAweakoutliers.csv") 


Y <- cbind(new_biz_density) 
X <- cbind(capital_requirements, cost_to_start) 

# Set data as panel data 
pdata <- plm.data(mydata, index=c("country_code","year")) 

# Descriptive statistics 
summary(Y) 
summary(X) 

# Pooled OLS estimator 
pooling <- plm(Y ~ X, data=pdata, model= "pooling") 
summary(pooling) 

# Between estimator 
between <- plm(Y ~ X, data=pdata, model= "between") 
summary(between) 

# First differences estimator 
firstdiff <- plm(Y ~ X, data=pdata, model= "fd") 
summary(firstdiff) 

# Fixed effects or within estimator 
fixed <- plm(Y ~ X, data=pdata, model= "within") 
summary(fixed) 

# Random effects estimator 
random <- plm(Y ~ X, data=pdata, model= "random") 
summary(random) 

# LM test for random effects versus OLS 
plmtest(pooling) 

# LM test for fixed effects versus OLS 
pFtest(fixed, pooling) 

# Hausman test for fixed versus random effects model 
phtest(random, fixed) 

回答

0

似乎有要計算的VIF(方差膨脹因子,檢測迴歸變量之間的共線性)R中的兩種常用方法:

的VIF()功能在汽車包裝中,其中輸入是模型。這要求您在檢查模型中變量之間的VIF之前首先擬合模型。

corvif()函數,其中的輸入是實際的候選解釋變量(即模型擬合之前的變量列表)。該功能是AED軟件包的一部分(Zuur et al。2009),已停用。這似乎只適用於變量列表,而不適用於擬合迴歸模型。