> mdata <-read.table("scrap.csv", header=T, sep=",")
> year <-mdata[,1]
> id   <-mdata[,2]
> y    <-mdata[,3]
> dt   <-mdata[,4]
> d2   <-mdata[,5]
> dtd2 <-mdata[,6]
> library(plm)
> mydata <-data.frame(year,id,y,dt,d2,dtd2)
> fixed <-plm(y~d2+dtd2, data=mydata, index=c("id","year"), model="within")
> summary(fixed)
Oneway (individual) effect Within Model

Call:
plm(formula = y ~ d2 + dtd2, data = mydata, model = "within", 
    index = c("id", "year"))

Balanced Panel: n = 40, T = 2, N = 80

Residuals:
    Min.  1st Qu.   Median  3rd Qu.     Max. 
-1.15074 -0.13889  0.00000  0.13889  1.15074 

Coefficients:
      Estimate Std. Error t-value Pr(>|t|)  
d2    0.096524   0.135352  0.7131   0.4801  
dtd2 -0.470945   0.196390 -2.3980   0.0215 *
---
Signif. codes:  0 �e***�f 0.001 �e**�f 0.01 �e*�f 0.05 �e.�f 0.1 �e �f 1

Total Sum of Squares:    8.7394
Residual Sum of Squares: 7.3098
R-Squared:      0.16359
Adj. R-Squared: -0.73886
F-statistic: 3.71601 on 2 and 38 DF, p-value: 0.033575
>