Let\'s consider data following :
library(plm) data("EmplUK", package="plm") df1<-EmplUK df1 <- cbind(df1,"Trend" = as.nume