如何用R做面板数据模型单位根检验？

如题，谢谢！
附数据：
stata.rar (10.58 KB) 本附件包括：

• stata.csv

2012-3-23 17:47:45 上传

R代码：
xiancheng=read.csv("stata.csv")
library(plm)
attach(xiancheng)
purtest(y1)

提示：
错误于purtest(y1) : the individual dimension is undefined
请问我该如何做？
谢谢！

rm(list=ls())
ls()
library(plm)
data=read.csv("C:\\Users\\acer\\Desktop\\stata.csv", header = TRUE);data#读入数据
levels(data$y1)#y1是因子
data$y1[1]/10#不能运算
xiancheng=data[1:204,]#去掉缺省值
attach(xiancheng)#数据绑定
y=as.vector(y1)
y=as.numeric(y)#把数据转化为可运算的数值型，不能一步到位，否则数据将被更改
y <- data.frame(split(y, xian))#把y按xian分组
purtest(y ,pmax =0,  exo = "none", test = "madwu")#单位根检验
?purtest#查看帮助

library(plm)
purtest()

补充下详细的例子：http://rgm2.lab.nig.ac.jp/RGM2/func.php?rd_id=plm:purtest

谢谢两位朋友，我试了，但提示错误，请看主贴。

使者，你这里有缺省数据啊
xiancheng=read.csv("C:\\Users\\acer\\Desktop\\stata.csv")
library(plm)
attach(xiancheng)
y1 <- data.frame(split(y1, xian))
purtest(y1,pmax = 4, exo = "intercept", test = "madwu")


错误于lm.fit(X, y) : 外接函数调用时不能有NA/NaN/Inf(arg1)
此外: 警告信息：
In Ops.factor(object[2:length(object)], object[1:(length(object) -  :
  - not meaningful for factors



你可以?purtest查看一下帮助的例子，你就懂了

吉林小王子 发表于 2012-3-23 18:13
使者，你这里有缺省数据啊
xiancheng=read.csv("C:\\Users\\acer\\Desktop\\stata.csv")
library(plm)

查了啊，仍不明就里。我也怀疑是我的数据有问题，但不知问题格式出在了哪。
唉

耕耘使者 发表于 2012-3-23 22:43
查了啊，仍不明就里。我也怀疑是我的数据有问题，但不知问题格式出在了哪。
唉

里面不能有缺省值啊，你数据里面有太多的NA了

吉林小王子 发表于 2012-3-23 22:47
里面不能有缺省值啊，你数据里面有太多的NA了

可能不是，我把其它变量删除了，只剩下y1，缺失值也随便填上了，但还是不行。

rm(list=ls())
library(plm)
data=read.csv("C:\\Users\\acer\\Desktop\\stata.csv")
xiancheng=data[1:204,]
attach(xiancheng)
y <- data.frame(split(as.numeric(y1), xian))
purtest(y ,pmax = 4, exo = "intercept", test = "madwu")


       Maddala-Wu Unit-Root Test (ex. var. : Individual Intercepts )

data:  y
chisq = Inf, df = 34, p-value < 2.2e-16
alternative hypothesis: stationarity

我知道了，在读进数据的时候，数据的格式是因子的格式，因此不能参与数值运算，而且你那里面的缺省值应该不能计算，没有意义吧