老兄你可能是需要

row number & column number

来修改你的数据

底下是一个300 x 300 matrix,

ywh.dat

  

你可以由matlab读入,修改,存档.

再由R读入.

epoh老师， 为什么我运行您2楼给出的程序会出现错误
> w1=nb2listw(COL.nb)   
> W <- as(as_dgRMatrix_listw(nb2listw(COL.nb)), "CsparseMatrix")                                                               
> x=as.matrix(W)   
> write.csv(x,"Wmatrix.csv",row.names = F,col.names = F)
警告信息：
In write.csv(x, "Wmatrix.csv", row.names = F, col.names = F) :
  不能修改'col.names'

而，您运行又可以呢？
我不想把列名和行名读入的，而现在只有把col.names=F删除才能运行

epoh老师，我的数据矩阵是261*261的，这个用您的命令，通过R可以保存到excel，但是excel打不开那么打的数据
我还要用这个 261*261数据矩阵做运算的

附件是我的数据矩阵，如何把第一行的行名称“V1”那些删除？

第二个数据才是，第一次不是，传错了，不能删除

头脑坏了,
应该在s-plus修改数据才方便
####in R
x=matrix(rnorm(261*261),261,261)
dim(x)    #[1] 261 261
write.csv(x,"ywh.csv",row.names = F)
#####in S-plus
Data\Select Data\Import file\......ywh.csv
#after modified
File\Export Data\To File

######
click "1",then delete row 1.
   

谢谢epoh老师了，软件功能就是比excel强大一些
您按这种方式导出的csv，哈哈，在excel也是看不全的，只有256列

epoh老师，您抽空能否看下看下这个winbugs程序，
模型和数据确认都可以通过，但是获取初始值
时就出错。

老兄,你的朋友似乎是要做Multivariate Stochastic Volatility Models
中的Model 3 (MSV with Granger causality or GC-MSV).
只是不知他的公式依据什么而来:
  thmean[i,1]<-mu1+phi11*(th[i-1,1]-mu1)+phi12*(th[i-1,2]+mu2);
  thmean[i,2]<-mu2+phi21*(th[i-1,1]-mu1)+phi22*(th[i-1,2]+mu2);

若要依他的程序运行也可以
只要加上initial values
list(phi1star=0.99, phi2star=0.99,mu1=0,mu2=0,itaua2=100,itaub2=100,rhoep=0.7,phi21=0,phi12=0)
ywh_gc_msv.txt   

是的，epoh老师，就是Multivariate Stochastic Volatility Models。
非常感谢epoh老师提醒，那个公式我让他再确认下，很有可能是
他写错了，他是根据别人的程序自己修改的，自己不是很懂。
epoh老师，这个模型您有相应的程序吗。

MULTIVARIATE STOCHASTIC VOLATILITY MODELS:
BAYESIAN ESTIMATION AND MODEL COMPARISON.pdf
  http://harrisd.net/papers/ARCHSV ... s/YuMeyer2006ER.pdf

谢谢epoh老师，记得他让我下载过这篇文献，那应该是
原来的程序存在问题，非常谢谢。

% Read data
A=xlsread('X:\freret\Data1.xls');
W=wk1read('X:\freret\weight.wk1');
请问题Dtata1.xls和weight.wk1这两个数据文件epoh老师您有吗？要是有，请帮忙提供，非常感谢！
本文来自: 人大经济论坛 Matlab及其他计量软件专版 版，详细出处参考： https://bbs.pinggu.org/forum.php? ... 1&from^^uid=11232

帮你看了,似乎Paul Elhorst的网站也没提供DEMO数据
ttp://www.regroningen.nl/elhorst/software.shtml   

To get started, select MATLAB Help or Demos from the Help menu.
??? Error using ==> wk1read
File not found.
Error in ==> demoLMsarsem_panel at 1
A=wk1read('E:\lotus\cigarette.wk1',1,0); % data set with T=30
??? Error using ==> wk1read
File not found.
Error in ==> demopanelscompare at 38
A=wk1read('E:\lotus\cigarette.wk1',1,0);
>>
epoh老师，您好！
以上问题如何解决？
非常感谢！

clear all;

A=wk1read('cigardemo.wk1',1,0);
W1=wk1read('Spat-Sym-US.wk1');
为什么这个wk1read又没有问题？可以读'cigardemo.wk1'
% Dataset downloaded from www.wiley.co.uk/baltagi/
% Spatial weights matrix constructed by Elhorst
% -------------------------------------------------------------------------
% written by: J.Paul Elhorst 5/2008
% University of Groningen
% Faculty of Economics and Business
% 9700AV Groningen
% the Netherlands
% j.p.elhorst@rug.nl
%
% REFERENCES:
% Elhorst J.P. (2008) Serial and spatial autocorrelation. Economics Letters
% http://dx.doi.org/10.1016/j.econlet.2008.03.009
% -------------------------------------------------------------------------
% dimensions of the problem
T=6; % number of time periods
N=46; % number of regions
% row-normalize W
W=normw(W1); % function of LeSage
y=A(:,[3]); % column number in the data matrix that corresponds to the dependent variable
x=A(:,[4,5,6]); % column numbers in the data matrix that correspond to the independent variables
xconstant=ones(N*T,1);

% serial and spatial error correlation, including logpn variable
% one constraint on the paramaters active
vnames=strvcat('logcit','constant','logp','logpn','logy');
serialspatial(y,[xconstant x],W,N,T,vnames);

% serial and spatial error correlation, without logpn variable
% no constraint active, numerical problems with one of the t-values
x=A(:,[4,6]); % column numbers in the data matrix that correspond to the independent variables
vnames=strvcat('logcit','constant','logp','logy');
serialspatial(y,[xconstant x],W,N,T,vnames);



Optimization terminated: Search direction less than 2*options.TolX
and maximum constraint violation is less than options.TolCon.
Active inequalities (to within options.TolCon = 1e-006):
  lower      upper     ineqlin   ineqnonlin
                          3           

model with both serial and spatial autocorrelation
Dependent Variable =         logcit   
R-squared          =    0.9921   
Rbar-squared       =    0.9920   
sigma^2            =    0.0021   
log-likelihood     =       -380.51003  
***************************************************************
Variable       Coefficient  Asymptot t-stat    z-probability
constant          3.184573        12.127214         0.000000
logp             -0.480758        -9.024798         0.000000
logpn             0.226449         3.151678         0.001623
logy              0.368082         6.157506         0.000000
ser.aut           0.963577       116.546888         0.000000
spat.aut.         0.036423        97.200039         0.000000

Optimization terminated: Search direction less than 2*options.TolX
and maximum constraint violation is less than options.TolCon.
Active inequalities (to within options.TolCon = 1e-006):
  lower      upper     ineqlin   ineqnonlin
                          3           

model with both serial and spatial autocorrelation
Dependent Variable =         logcit   
R-squared          =    0.9918   
Rbar-squared       =    0.9917   
sigma^2            =    0.0021   
log-likelihood     =       -375.63073  
***************************************************************
Variable       Coefficient  Asymptot t-stat    z-probability
constant          3.355565        12.831304         0.000000
logp             -0.438336        -8.355384         0.000000
logy              0.322546         5.463177         0.000000
ser.aut           0.963577       110.933671         0.000000
spat.aut.         0.036423        94.006837         0.000000

>>

zhangzhao兄，你的数据是放在E:\lotus\cigarette.wk1吗？
应该不是命令错误，只是提示找不到数据文件，自己设置过
路径就行了。

% PURPOSE: An example of two types of sampling for rho in the semip_g function
% Gibbs sampling spatial Probit model with individual effects
% using 1996 presidential election data set
%               
%---------------------------------------------------
% USAGE: semip_gd2
%---------------------------------------------------

clear all;
load semip.mat;

vnames = strvcat('y','highs','college','grad','non-white', ...
'income','urban');

% x-matrix contains 3,110 x 6 matrix with:
%   col1 = high school graduates as a percent of population     
%   col2 = college percent  
%   col3 = graduate school
%   col4 = non-white
%   col5 = median income     
%   col6 = urban  

% z = 0,1 with 0 = Dole wins, 1 = Clinton wins, 3,110 counties
% W = a 48x48 spatial weight matrix (standardized)
% nregions = 48
% regionobs = a 48 x 1 vector with the # of counties in each state
% states organized alphabetically

[n k] = size(x);


ndraw = 2500;
nomit = 500;
prior.rval = 4;
prior.dflag = 1; % metropolis-hastings for rho
prior.rmin = 0;
prior.rmax = 1;

% use M-H sampling for rho
result = semip_g(z,x,W,nregions,regionobs,ndraw,nomit,prior);
prt(result,vnames);

prior.dflag = 0; % use alternative approach to sample for rho

result2 = semip_g(z,x,W,nregions,regionobs,ndraw,nomit,prior);
prt(result2,vnames);

[h1,f1,y1] = pltdens(result.pdraw);
[h2,f2,y2] = pltdens(result2.pdraw);

plot(y1,f1,'.',y1,f2,'o');
legend('Met-H','Inversion');
在以上程序中pltdens和legend这两个函数哪儿？为什么在原函数包中看不到?

epoh老师，我还想请教您一下，winbugs提示这个错误
sorry something wrong in procedure Stack.Value in module Graphstack
通常是因为数据量太大吗？ 程序赋初始值后，把迭代次数放大，也会出现
这个错误。

lenged是matlab自带的命令，就是给图形注释标签。

%%%%%  plotdens
function [h,f,xx] = plotdens(x,h,positive,kernel)
% PURPOSE: Draw a nonparametric density estimate.
%---------------------------------------------------
% USAGE: [h f y] = pltdens(x,h,p,kernel)
%        or pltdens(x) which uses gaussian kernel default
% where:
%        x is a vector
%        h is the kernel bandwidth
%          default=1.06 * std(x) * n^(-1/5); Silverman page 45
%        p is 1 if the density is 0 for negative values
%        k is the kernel type:
%          =1 Gaussian (default)
%        =2 Epanechnikov
%        =3 Biweight
%        =4 Triangular
%   A jittered plot of the
%   observations is shown below the density.
%---------------------------------------------------
% RETURNS:
%        h = the interval used
%        f = the density
%        y = the domain of support
%        plot(y,f) will produce a plot of the density
% --------------------------------------------------
% SEE ALSO hist, histo
%---------------------------------------------------

%       Anders Holtsberg, 18-11-93
%       Copyright (c) Anders Holtsberg

x = x(:); n = length(x);
if nargin < 4, kernel = 1; end
if nargin < 3, positive = 0; end
if nargin < 2, h = []; end
if isempty(h)
   h = 1.06 * std(x) * n^(-1/5);  % Silverman page 45
end
if positive & any(x < 0)
   error('There is a negative element in X')
end
mn1 = min(x); mx1 = max(x);
mn = mn1 - (mx1-mn1)/3;
mx = mx1 + (mx1-mn1)/3;
gridsize = 256;
xx = linspace(mn,mx,gridsize)';
d = xx(2) - xx(1);
xh = zeros(size(xx));
xa = (x-mn)/(mx-mn)*gridsize;
for i=1:n
   il = floor(xa(i));
   a  = xa(i) - il;
   xh(il+[1 2]) = xh(il+[1 2])+[1-a, a]';
end
xk = [-gridsize:gridsize-1]'*d;
if kernel == 1
   K = exp(-0.5*(xk/h).^2);
elseif kernel == 2
   K = max(0,1-(xk/h).^2/5);
elseif kernel == 3
   c = sqrt(1/7);
   K = (1-(xk/h*c).^2).^2 .* ((1-abs(xk/h*c)) > 0);
elseif kernel == 4
   c = sqrt(1/6);
   K = max(0,1-abs(xk/h*c));
end
K = K / (sum(K)*d*n);
f = ifft(fft(fftshift(K)).*fft([xh ;zeros(size(xh))]));
f = real(f(1:gridsize));
if positive
   m = sum(xx<0);
   f(m+(1:m)) = f(m+(1:m)) + f(m:-1:1);
   f(1:m) = zeros(size(f(1:m)));
   xx(m+[0 1]) = [0 0];
end
if nargout == 0
plot(xx,f)
if ~ishold
   hold on
   d = diff(get(get(gcf,'CurrentAxes'),'Ylim'))/100;
   plot(x,(-rand(size(x))*6-1)*d,'.')
   plot([mn mx],[0 0])
   axis([mn mx -0.2*max(f) max(f)*1.2])
   hold off
end
end;

ywh19860616学兄，那这个plotdens是哪儿的函数？
非常感谢！

http://rgm2.lab.nig.ac.jp/RGM2/f ... Smodelling:plotDens
上面链接有比较详细的介绍。
如果你的电脑上没有这个函数，那就把这个函数下载保存到你电脑上，
应该就可以用了。

1.cigardemo.wk1,Spat-Sym-US.wk1两个文件都在
   jplv7\spatial\panel
  所以可以读到

2.cigarette.wk1,你只要放在jplv7\spatial\panel,就可读到

3.pltdens.m 在
   jplv7\graphs
  Elhorst JP (2011) MATLAB Software for Spatial Panels
  很多都要用到 LeSage's Spatial Econometrics Toolbox

老兄这个问题,你可参考
  https://www.jiscmail.ac.uk/cgi-bin/webadmin?A2=bugs;f81d94fc.1205
我不清楚程序,数据,不便回答

好的，我看下您给的链接，谢谢epoh老师。

??? Error using ==> wk1read
Too many output arguments.

Error in ==> demoLMsarsem_panel at 1
A=wk1read('E:\lotus\cigarette.wk1',1,0); % data set with T=30

>>
现在又出现以上问题，epoh都是，您看如何解决？