Twoways effects Random Effect Model
(Nerlove's transformation)
Call:
plm(formula = formula_, data = data_plm, effect = c("twoways"),
    model = c("random"), random.method = "nerlove", index = c("kraj",
        "rok"))
Unbalanced Panel: n = 77, T = 2-27, N = 1859
Effects:
                     var    std.dev share
idiosyncratic 0.0023478 0.0484546 0.898
individual     0.0001121 0.0105867 0.043
time           0.0001549 0.0124477 0.059
theta:
             Min.   1st Qu.     Median     Mean   3rd Qu.      Max.
id     0.04456980 0.3390225 0.3390225 0.3267491 0.3390225 0.3390225
time 0.54487559 0.5624653 0.5756626 0.5758629 0.5922814 0.5944996
total 0.04142385 0.2891521 0.2955809 0.2865027 0.3007343 0.3018645
Residuals:
    Min. 1st Qu.     Median      Mean   3rd Qu.      Max.
-0.69811 -0.01590   0.00333   0.00014   0.02186   0.28729
Coefficients:
               Estimate   Std. Error z-value Pr(>|z|)
(Intercept) 0.01835626    0.13849695 0.1325    0.8946
Xdk          0.56050153   0.79657493 0.7036    0.4817
Xdl          0.00036219   1.39571914 0.0003    0.9998
Xdh         -0.01245029   1.59971294 -0.0078   0.9938
Xyr         -0.00061122   0.00708244 -0.0863   0.9312
Total Sum of Squares:    5.6227
Residual Sum of Squares: 4.8113
R-Squared:      0.14436
Adj. R-Squared: 0.14252
Chisq: 0.524963 on 4 DF, p-value: 0.97103