R里面的chisq.test()函数可以直接拿来解决这个问题：
chisq.test (x, y=Null, correct=TRUE,
                    p=rep(1/length (x), lenghth(x)),
                    rescale.p= FALSE,
                    simulate.p.value = FALSE, B=2000)
x_是由观测数据构成的向量或矩阵
y_是数据向量（当x为矩阵时，y无效）
correct_是逻辑变量，表明是否用于连续修正，TRUE（默认值）表示修正，FALSE表示不修正
p_是原假设落在小区间的理论概率，默认均匀分布
rescale.p_是逻辑变量，选择FALSE（默认值）时，要求输入的P满足∑Pi=1，选择TRUE时不要求这一点，程序会重新计算P值
simulate.p.value_是逻辑变量（默认为FALSE），当为TRUE时，将用仿真的方法计算P值，此时，B表示仿真的次数。
PS：来自《R统计建模与R软件》薛毅等编著，上面有具体例子，简单明了~~

求大神，麻烦告诉一下空间误差模型和空间滞后模型的拟合优度的程序？