In [127]:
data <- read.csv("adinda.clean2.csv")
In [132]:
summary(data[4:26])
Out[132]:
In [133]:
head(data)
Out[133]:
In [134]:
tail(data)
Out[134]:
In [135]:
names(data)
Out[135]:
In [137]:
# ncol(data)
dataMean <- data[4:26]
dataMean
Out[137]:
In [112]:
library(stargazer)
In [113]:
summary(data)
Out[113]:
In [114]:
library("ggplot2")
In [115]:
qplot( y = data$TOBINQ, x = data$SIZE)
In [138]:
library("plm")
In [139]:
data_panel <- pdata.frame(data, index = c("Firm", "Year"), drop.index = TRUE, row.names = TRUE)
In [140]:
head(data_panel)
Out[140]:
In [141]:
tabel2_reg1 <- plm(CEI ~ BRD_INDP + BRD_MEET + BRD_SIZE + BRD_COMPT, data = data_panel, model = "within" )
tabel2_reg2 <- plm(CEI ~ BRD_INDP + BRD_MEET + BRD_SIZE + BRD_COMPT + LOSS + TOT_ASSETS + SLACK + ROE , data = data_panel, model = "within" )
tabel2_reg3 <- plm(BRIB_CORR ~ BRD_INDP + BRD_MEET + BRD_SIZE + BRD_COMPT + LOSS + TOT_ASSETS + SLACK + ROE , data = data_panel, model = "within" )
tabel2_reg4 <- plm(BUSS_ETH ~ BRD_INDP + BRD_MEET + BRD_SIZE + BRD_COMPT + LOSS + TOT_ASSETS + SLACK + ROE , data = data_panel, model = "within" )
tabel2_reg5 <- plm(FAIR_COMP ~ BRD_INDP + BRD_MEET + BRD_SIZE + BRD_COMPT + LOSS + TOT_ASSETS + SLACK + ROE , data = data_panel, model = "within" )
tabel2_reg6 <- plm(POL_CONTR ~ BRD_INDP + BRD_MEET + BRD_SIZE + BRD_COMPT + LOSS + TOT_ASSETS + SLACK + ROE , data = data_panel, model = "within" )
tabel2_reg7 <- plm(INDIG_PPL ~ BRD_INDP + BRD_MEET + BRD_SIZE + BRD_COMPT + LOSS + TOT_ASSETS + SLACK + ROE , data = data_panel, model = "within" )
tabel2_reg8 <- plm(IND_EC_IMP ~ BRD_INDP + BRD_MEET + BRD_SIZE + BRD_COMPT + LOSS + TOT_ASSETS + SLACK + ROE , data = data_panel, model = "within" )
tabel2_reg9 <- plm(X0TH_ENG ~ BRD_INDP + BRD_MEET + BRD_SIZE + BRD_COMPT + LOSS + TOT_ASSETS + SLACK + ROE , data = data_panel, model = "within" )
In [ ]:
In [142]:
library(lmtest)
In [143]:
tabel3_reg2 <- plm(CEI ~ BRD_INDP + BRD_MEET + BRD_SIZE + BRD_COMPT + DIVIDEND + LOSS + TOT_ASSETS + SLACK + ROE , data = data_panel, model = "within" )
tabel3_reg2 <- plm(CEI ~ BRD_INDP + BRD_MEET + BRD_SIZE + BRD_COMPT + DIVIDEND + BRD_INDP_DIV + LOSS + TOT_ASSETS + SLACK + ROE , data = data_panel, model = "within" )
tabel3_reg3 <- plm(CEI ~ BRD_INDP + BRD_MEET + BRD_SIZE + BRD_COMPT + DIVIDEND + BRD_MEET_DIV +LOSS + TOT_ASSETS + SLACK + ROE , data = data_panel, model = "within" )
tabel3_reg4 <- plm(CEI ~ BRD_INDP + BRD_MEET + BRD_SIZE + BRD_COMPT + DIVIDEND + BRD_SIZE_DIV +LOSS + TOT_ASSETS + SLACK + ROE , data = data_panel, model = "within" )
tabel3_reg5 <- plm(CEI ~ BRD_INDP + BRD_MEET + BRD_SIZE + BRD_COMPT + DIVIDEND + BRD_COMPT_DIV +LOSS + TOT_ASSETS + SLACK + ROE , data = data_panel, model = "within" )
In [144]:
summary(tabel3_reg5)
Out[144]:
In [145]:
names(data)
Out[145]:
In [146]:
tabel4_reg1 <- plm(CEI ~ BRD_EFFC + DIVIDEND + LOSS + TOT_ASSETS + SLACK + ROE, data = data_panel, model = "within" )
tabel4_reg2 <- plm(ALT_CEI ~ BRD_EFFC + DIVIDEND + LOSS + TOT_ASSETS + SLACK + ROE, data = data_panel, model = "within" )
tabel4_reg3 <- plm(CEI ~ BRD_EFFC + DIVIDEND + BRD_EFFC_DIV + LOSS + TOT_ASSETS + SLACK + ROE, data = data_panel, model = "within" )
tabel4_reg4 <- plm(ALT_CEI ~ BRD_EFFC + DIVIDEND + BRD_EFFC_DIV + LOSS + TOT_ASSETS + SLACK + ROE, data = data_panel, model = "within" )
tabel4_reg5 <- plm(BRD_EFFC ~ ALT_CEI + LOSS + TOT_ASSETS + SLACK + ROE, data = data_panel, model = "within" )
tabel4_reg6 <- plm(ALT_CEI ~ BRD_EFFC +LOSS + TOT_ASSETS + SLACK + ROE, data = data_panel, model = "within" )
In [147]:
robust111 <- coeftest(tabel4_reg1, vcov=function(x) vcovHC(x, method="arellano",
type="HC1", cluster="group"))
robust222 <- coeftest(tabel4_reg2, vcov=function(x) vcovHC(x, method="arellano",
type="HC1", cluster="group"))
robust333 <- coeftest(tabel4_reg3, vcov=function(x) vcovHC(x, method="arellano",
type="HC1", cluster="group"))
robust445 <- coeftest(tabel4_reg4, vcov=function(x) vcovHC(x, method="arellano",
type="HC1", cluster="group"))
robust555 <- coeftest(tabel4_reg5, vcov=function(x) vcovHC(x, method="arellano",
type="HC1", cluster="group"))
robust666 <- coeftest(tabel4_reg5, vcov=function(x) vcovHC(x, method="arellano",
type="HC1", cluster="group"))
In [76]:
skripsi.fe <- plm(TOBINQ ~ PDEBT + LogAssets + LEV + LIQ + INTCOV + PROF + DR + GROWTH + FIRMAGE, data = data_panel, model = "within")
skripsi.re <- plm(TOBINQ ~ PDEBT + LogAssets + LEV + LIQ + INTCOV + PROF + DR + GROWTH + FIRMAGE, data = data_panel, model = "random")
In [77]:
summary(skripsi.fe)
In [78]:
summary(skripsi.re)
In [79]:
phtest( skripsi.fe, skripsi.re)
In [80]:
pooltest.billy <- pooltest(TOBINQ ~ PDEBT + FIRMAGE + LIQ, data = data_panel, model="within")
In [81]:
pooltest.billy
In [82]:
pool_rergression <- plm(TOBINQ ~ PDEBT + SIZE + LEV + LIQ + INTCOV + PROF + DR + GROWTH + FIRMAGE, data = data_panel, model="pooling")
bpagan_test <- plmtest(pool_rergression,effect="twoways",type="bp")
In [83]:
bpagan_test
In [84]:
summary(skripsi.fe)
In [85]:
summary(data_panel)
Out[85]:
In [86]:
qplot(y = data$LogTobinsQ, x = data$LIQ, geom = c("point", "smooth"), method = lm)
In [87]:
var(data$LIQ)
In [88]:
var(data$TOBINQ)
In [89]:
qplot(data$TOBINQ)
In [90]:
skripsi.try1 <- plm(log(TOBINQ) ~ PDEBT + log(SIZE) + log(LEV) + log(LIQ) + log(INTCOV) + log(PROF) + DR + log(GROWTH) + FIRMAGE + y1 + y2 + y3 , data = data_panel, model = "within")
In [91]:
summary(skripsi.try1)
In [92]:
install.packages("R2wd")
In [ ]:
In [ ]:
In [ ]:
In [ ]:
In [ ]:
In [ ]: