怎么把这个矩阵保存到excel中

epoh 发表于 2012-11-6 14:26
有缺省值,就会形成Levels

> data2$x10

splm_1.0-00.zip (145.31 KB)

2012-11-6 21:13:35 上传

epoh老师，您好！
     这个splm包中的数据Produc,可以通过fix(Produc)看到，
您看有没有什么办法可以把Produc下载下来，然后保存为
excel格式，以供练习之用。
非常感谢！

zhangtao 发表于 2012-11-6 21:15
epoh老师，您好！
     这个splm包中的数据Produc,可以通过fix(Produc)看到，
您看有没有什么办法可以 ...

data(Produc)
write.csv(Produc,"Produc.csv")
Produc.csv

epoh 发表于 2012-11-6 21:52
data(Produc)
write.csv(Produc,"Produc.csv")
Produc.csv

载入需要的程辑包：lmtest
载入需要的程辑包：Ecdat
警告信息：
1: 程辑包‘splm’是用R版本2.15.2 来建造的
2: 程辑包‘spam’是用R版本2.15.2 来建造的
3: 程辑包‘ibdreg’是用R版本2.15.2 来建造的
4: 程辑包‘car’是用R版本2.15.2 来建造的
5: 程辑包‘lmtest’是用R版本2.15.2 来建造的
> data(Produc)
>  write.csv(Produc,"Produc.csv")
>  Produc.csv
错误: 找不到对象'Produc.csv'
> fix(Produc)
>

zhangtao 发表于 2012-11-7 08:54
载入需要的程辑包：lmtest
载入需要的程辑包：Ecdat
警告信息：

哈哈!Produc.csv是我要上传的文件
只是昨日上传失败
####
data(Produc)
write.csv(Produc,"Produc.csv")
以上两句已经把数据写入Produc.csv
要读出只要
Produc=read.csv("Produc.csv")
Produc

epoh 发表于 2012-11-7 11:03
哈哈!Produc.csv是我要上传的文件
只是昨日上传失败
####

example.rar (28.17 KB)

2012-11-8 08:20:57 上传

Sexample=read.csv("D:\Backup\Documents\example\Sexample.csv")

Uexample=read.csv("D:\Backup\Documents\example\Uexample.csv")

data(Sexample, package = "Ecdat")data(Uexample)fm <- log(gsp) ~ log(pcap) + log(pc) + log(emp) + unemp## the two standard specifications (SEM and SAR) one with FE## and the other with RE:## fixed effects panel with spatial errorsfespaterr <- spml(fm, data = Sexample, listw = mat2listw(Uexample), model="within", spatial.error="b")summary(fespaterr)## random effects panel with spatial lagrespatlag <- spml(fm, data = Sexample, listw = mat2listw(Uexample), model="random", spatial.error="none", lag=TRUE)summary(respatlag) epoh老师，您好！  根据已经传给您的splm包，运行以上程序，出现以下问题，您看如何解决？非常感谢！本附件 中是数据。> Sexample=read.csv("D:\Backup\Documents\example\Sexample.csv")
错误: 由"D:\B"开头的字符串中存在'\B'，但没有这种逸出号
> Uexample=read.csv("D:\Backup\Documents\example\Uexample.csv")
错误: 由"D:\B"开头的字符串中存在'\B'，但没有这种逸出号
> data(Sexample, package = "Ecdat")
警告信息：
In data(Sexample, package = "Ecdat") : 没有‘Sexample’这个数据集
> data(Uexample)
警告信息：
In data(Uexample) : 没有‘Uexample’这个数据集
> fm <- log(gsp) ~ log(pcap) + log(pc) + log(emp) + unemp
> ## the two standard specifications (SEM and SAR) one with FE
> ## and the other with RE:
> ## fixed effects panel with spatial errors
> fespaterr <- spml(fm, data = Sexample, listw = mat2listw(Uexample), model="within", spatial.error="b")
错误于mat2listw(Uexample) : x is not a matrix
> summary(fespaterr)
Spatial panel fixed effects error modelCall:
spml(formula = fm, data = Produc, listw = mat2listw(usaww), model = "within",
    spatial.error = "b")Residuals:
   Min. 1st Qu.  Median 3rd Qu.    Max.
-0.1250 -0.0238 -0.0035  0.0171  0.1880 Coefficients:
            Estimate Std. Error t-value Pr(>|t|)   
rho        0.5574013  0.0330749 16.8527  < 2e-16 ***
log(pcap)  0.0051438  0.0250109  0.2057  0.83705   
log(pc)    0.2053026  0.0231427  8.8712  < 2e-16 ***
log(emp)   0.7822540  0.0278057 28.1328  < 2e-16 ***
unemp     -0.0022317  0.0010709 -2.0839  0.03717 *  
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 > ## random effects panel with spatial lag
> respatlag <- spml(fm, data = Sexample, listw = mat2listw(Uexample), model="random", spatial.error="none", lag=TRUE)
错误于mat2listw(Uexample) : x is not a matrix
> summary(respatlag)
Spatial panel random effects ML modelCall:
spml(formula = fm, data = Produc, listw = mat2listw(usaww), model = "random",
    lag = TRUE, spatial.error = "none")Residuals:
   Min. 1st Qu.  Median    Mean 3rd Qu.    Max.
   1.38    1.57    1.70    1.70    1.80    2.13 Error variance parameters:
    Estimate Std. Error t-value Pr(>|t|)  
phi  21.3175     8.3017  2.5678  0.01023 *Spatial autoregressive coefficient:
       Estimate Std. Error t-value  Pr(>|t|)   
lambda 0.161615   0.029099   5.554 2.793e-08 ***Coefficients:
               Estimate  Std. Error t-value  Pr(>|t|)   
(Intercept)  1.65814995  0.15071855 11.0016 < 2.2e-16 ***
log(pcap)    0.01294505  0.02493997  0.5190    0.6037   
log(pc)      0.22555376  0.02163422 10.4258 < 2.2e-16 ***
log(emp)     0.67081075  0.02642113 25.3892 < 2.2e-16 ***
unemp       -0.00579716  0.00089175 -6.5009 7.984e-11 ***
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 > Sexample
      X          state year      pcap      hwy    water     util        pc    gsp     emp unemp
1     1        ALABAMA 1970  15032.67  7325.80  1655.68  6051.20  35793.80  28418  1010.5   4.7
2     2        ALABAMA 1971  15501.94  7525.94  1721.02  6254.98  37299.91  29375  1021.9   5.2
3     3        ALABAMA 1972  15972.41  7765.42  1764.75  6442.23  38670.30  31303  1072.3   4.7

zhangtao 发表于 2012-11-8 08:23
Sexample=read.csv("D:\Backup\Documents\example\Sexample.csv")Uexample=read.csv("D:\Backup\Docume ...

你给的Uexample必须是square matrix
请看function mat2listw的用法
     http://127.0.0.1:13372/library/spdep/html/mat2listw.html
  
mat2listw(x, row.names = NULL, style="M")

x : A square non-negative matrix with no NAs representing spatial weights; may be a matrix of class “sparseMatrix”

epoh 发表于 2012-11-8 11:30
你给的Uexample必须是square matrix
请看function mat2listw的用法
     http://127.0.0.1:13372/libra ...

epoh老师，您好！
我是按您的方法，从splm包中下载了produc和usaww两个数据，然后存为
sexample.csv和uexample.csv，然后用read.csv读入进行估计，现在的问题是
用

Sexample=read.csv("Sexample.csv")

Uexample=read.csv("Uexample.csv")

这两个命令无法读入数据，您运行一下看是什么问题？

Sexample=read.csv("Sexample.csv")

Uexample=read.csv("Uexample.csv")

data(Sexample, package = "Ecdat")data(Uexample)fm <- log(gsp) ~ log(pcap) + log(pc) + log(emp) + unemp## the two standard specifications (SEM and SAR) one with FE## and the other with RE:## fixed effects panel with spatial errorsfespaterr <- spml(fm, data = Sexample, listw = mat2listw(Uexample), model="within", spatial.error="b")summary(fespaterr)## random effects panel with spatial lagrespatlag <- spml(fm, data = Sexample, listw = mat2listw(Uexample), model="random", spatial.error="none", lag=TRUE)summary(respatlag)

> Sexample=read.csv("Sexample.csv")
> Uexample=read.csv("Uexample.csv")
> data(Sexample, package = "Ecdat")
警告信息：
In data(Sexample, package = "Ecdat") : 没有‘Sexample’这个数据集
> data(Uexample)
警告信息：
In data(Uexample) : 没有‘Uexample’这个数据集
> fm <- log(gsp) ~ log(pcap) + log(pc) + log(emp) + unemp
> ## the two standard specifications (SEM and SAR) one with FE
> ## and the other with RE:
> ## fixed effects panel with spatial errors
> fespaterr <- spml(fm, data = Sexample, listw = mat2listw(Uexample), model="within", spatial.error="b")
错误于mat2listw(Uexample) : x is not a matrix
> summary(fespaterr)
错误于summary(fespaterr) :
  在为'summary'函数选择方法时评估'object'参数出了错: 错误: 找不到对象'fespaterr'

> ## random effects panel with spatial lag
> respatlag <- spml(fm, data = Sexample, listw = mat2listw(Uexample), model="random", spatial.error="none", lag=TRUE)
错误于mat2listw(Uexample) : x is not a matrix
> summary(respatlag)
错误于summary(respatlag) :
  在为'summary'函数选择方法时评估'object'参数出了错: 错误: 找不到对象'respatlag'

>
> Example
错误: 找不到对象'Example'
> Sexample
      X          state year      pcap      hwy    water     util        pc    gsp     emp unemp
1     1        ALABAMA 1970  15032.67  7325.80  1655.68  6051.20  35793.80  28418  1010.5   4.7
2     2        ALABAMA 1971  15501.94  7525.94  1721.02  6254.98  37299.91  29375  1021.9   5.2
3     3        ALABAMA 1972  15972.41  7765.42  1764.75  6442.23  38670.30  31303  1072.3   4.7
4     4        ALABAMA 1973  16406.26  7907.66  1742.41  6756.19  40084.01  33430  1135.5   3.9
5     5        ALABAMA 1974  16762.67  8025.52  1734.85  7002.29  42057.31  33749  1169.8   5.5
6     6        ALABAMA 1975  17316.26  8158.23  1752.27  7405.76  43971.71  33604  1155.4   7.7
7     7        ALABAMA 1976  17732.86  8228.19  1799.74  7704.93  50221.57  35764  1207.0   6.8
8     8        ALABAMA 1977  18111.93  8365.67  1845.11  7901.15  51084.99  37463  1269.2   7.4
9     9        ALABAMA 1978  18479.74  8510.64  1960.51  8008.59  52604
本文来自: 人大经济论坛 S-Plus&R专版 版，详细出处参考： https://bbs.pinggu.org/forum.php? ... 5&from^^uid=11232

> Uexample
                X ALABAMA   ARIZONA  ARKANSAS CALIFORNIA  COLORADO CONNECTICUT  DELAWARE FLORIDA GEORGIA     IDAHO  ILLINOIS   INDIANA      IOWA    KANSAS KENTUCKY LOUISIANA     MAINE  MARYLAND MASSACHUSETTS MICHIGAN MINNESOTA
1         ALABAMA   0.000 0.0000000 0.0000000       0.00 0.0000000         0.0 0.0000000    0.25   0.250 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000    0.000 0.0000000 0.0000000 0.0000000     0.0000000     0.00 0.0000000
2         ARIZONA   0.000 0.0000000 0.0000000       0.20 0.2000000         0.0 0.0000000    0.00   0.000 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000    0.000 0.0000000 0.0000000 0.0000000     0.0000000     0.00 0.0000000
3        ARKANSAS   0.000 0.0000000 0.0000000       0.00 0.0000000         0.0 0.0000000    0.00   0.000 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000    0.000 0.1666667 0.0000000 0.0000000     0.0000000     0.00 0.0000000
4      CALIFORNIA   0.000 0.3333333 0.0000000       0.00 0.0000000         0.0 0.0000000    0.00   0.000 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000    0.000 0.0000000 0.0000000 0.0000000     0.0000000     0.00 0.0000000
5        COLORADO   0.000 0.1428571 0.0000000       0.00 0.0000000         0.0 0.0000000    0.00   0.000 0.0000000 0.0000000 0.0000000 0.0000000 0.1428571    0.000 0.0000000 0.0000000 0.0000000     0.0000000     0.00 0.0000000
6     CONNECTICUT   0.000 0.0000000 0.0000000       0.00 0.0000000         0.0 0.0000000    0.00   0.000 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000    0.000 0.0000000 0.0000000 0.0000000     0.3333333     0.00 0.0000000
7        DELAWARE   0.000 0.0000000 0.0000000       0.00 0.0000000         0.0 0.0000000    0.00   0.000 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000    0.000 0.0000000 0.0000000 0.3333333     0.0000000     0.00 0.0000000
8         FLORIDA   0.500 0.0000000 0.0000000       0.00 0.0000000         0.0 0.0000000    0.00   0.500 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000    0.000 0.0000000 0.0000000 0.0000000     0.0000000     0.00 0.0000000
9         GEORGIA   0.200 0.0000000 0.0000000       0.00 0.0000000         0.0 0.0000000    0.20   0.000 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000    0.000 0.0000000 0.0000000 0.0000000     0.0000000     0.00 0.0000000
10          IDAHO   0.000 0.0000000 0.0000000       0.00 0.0000000         0.0 0.0000000    0.00   0.000 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000    0.000 0.0000000 0.0000000 0.0000000     0.0000000     0.00 0.0000000
11       ILLINOIS   0.000 0.0000000 0.0000000       0.00 0.0000000         0.0 0.0000000    0.00   0.000 0.0000000 0.0000000 0.2000000 0.2000000 0.0000000    0.200 0.0000000 0.0000000 0.0000000     0.0000000     0.00 0.0000000
12        INDIANA   0.000 0.0000000 0.0000000       0.00 0.0000000         0.0 0.0000000    0.00   0.000 0.0000000 0.2500000 0.0000000 0.0000000 0.0000000    0.250 0.0000000 0.0000000 0.0000000     0.0000000     0.25 0.0000000
13           IOWA   0.000 0.0000000 0.0000000       0.00 0.0000000         0.0 0.0000000    0.00   0.000 0.0000000 0.1666667 0.0000000 0.0000000 0.0000000    0.000 0.0000000 0.0000000 0.0000000     0.0000000     0.00 0.1666667
14         KANSAS   0.000 0.0000000 0.0000000       0.00 0.2500000         0.0 0.0000000    0.00   0.000 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000    0.000 0.0000000 0.0000000 0.0000000     0.0000000     0.00 0.0000000
15       KENTUCKY   0.000 0.0000000 0.0000000       0.00 0.0000000         0.0 0.0000000    0.00   0.000 0.0000000 0.1428571 0.1428571 0.0000000 0.0000000    0.000 0.0000000 0.0000000 0.0000000     0.0000000     0.00 0.0000000
16      LOUISIANA   0.000 0.0000000 0.3333333       0.00 0.0000000         0.0 0.0000000    0.00   0.000 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000    0.000 0.0000000 0.0000000 0.0000000     0.0000000     0.00 0.0000000
17          MAINE   0.000 0.0000000 0.0000000       0.00 0.0000000         0.0 0.0000000    0.00   0.000 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000    0.000 0.0000000 0.0000000 0.0000000     0.0000000     0.00 0.0000000
18       MARYLAND   0.000 0.0000000 0.0000000       0.00 0.0000000         0.0 0.2500000    0.00   0.000 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000    0.000 0.0000000 0.0000000 0.0000000     0.0000000     0.00 0.0000000
19  MASSACHUSETTS   0.000 0.0000000 0.0000000       0.00 0.0000000         0.2 0.0000000    0.00   0.000 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000    0.000 0.0000000 0.0000000 0.0000000     0.0000000     0.00 0.0000000
20       MICHIGAN   0.000 0.0000000 0.0000000       0.00 0.0000000         0.0 0.0000000    0.00   0.000 0.0000000 0.0000000 0.3333333 0.0000000 0.0000000    0.000 0.0000000 0.0000000 0.0000000     0.0000000     0.00 0.0000000
21      MINNESOTA   0.000 0.0000000 0.0000000       0.00 0

> Sexample=read.csv("Sexample.csv")
> Uexample=read.csv("Uexample.csv")
> data(Sexample, package = "Ecdat")
警告信息：
In data(Sexample, package = "Ecdat") : 没有‘Sexample’这个数据集
> data(Uexample)
警告信息：
In data(Uexample) : 没有‘Uexample’这个数据集
epoh老师，您好！
明明uexample and sexample两个数据集读进去了，为什么会提示
两个数据集找不到呢？