R2采用plm封装

Question

使用这种语法使用plm包估算合并模型时，得到的R2接近0.8。

// Attach it like the following......after getting your image
PhotoViewAttacher yourAttacher = new PhotoViewAttacher(yourImageView);

使用完全相同的规范并再次使用plm但是现在这个语法我得到的R2小于0.1且R2调整为负。

library(plm)
data <- employmentsez
data$lfirms2 <- data$lfirms*data$lfirms
data$sezname <- as.factor(data$sezname)
data <- plm.data(data, index=c("code", "year"))
fit1 <- plm(lemployment ~ lfirms + lfirms2 + lfirmsfor + lwages +  nuts51 + lgovgrants + leusubs + agrishare, data=data, model="pooling")
    R-Squared:      0.7972
    Adj. R-Squared: 0.79679

我认为差异太大了，在Stata中，使用完全相同的规格，内部模型的R2会更高。当使用完全相同的规范并添加固定效果时，内部模型还应该不提供更高的R2？

有谁知道如何解决这个问题，是否有另一种方法可以使用plm包获得正确的R2？