R程序改写为matlab程序

ywh19860616

18698

101

收藏 2011-10-30

程序1-R程序
x<-c(1,2,3)
n<-length(x)
xibar <- (rep(sum(x), n) - x)/(n - 1)
si2 <- ((rep(sum(x^2), n) - x^2)/(n - 1)) - xibar^2
Wi <- sapply(listw$weights, sum)
Wi2=Wi^2

S1i <- sapply(listw$weights, function(x) sum(x^2))

lx <- lag.listw(listw, x, zero.policy = zero.policy)

res <- (lx - Wi * xibar)

res <- res/sqrt(si2 * (((n - 1) * S1i - Wi^2)/(n - 2)))

程序2-matlab程序

x=input('x=')
weights=input('weights=')
n=length(x)
xibar=(sum(x)*ones(n,1)-x)/(n-1)
si2=((sum(x.^2).*ones(n,1)-x.^2)/(n-1))-xibar.^2
Wi=weights(:)' %这句对应Wi <- sapply(listw$weights, sum)是否正确？
Wi2=Wi.^2
S1i=sum(weights.^2)
lx=Wi*x %这句如何从lx <- lag.listw(listw, x, zero.policy = zero.policy)改写

res=lx-Wi*xibar

res =res./sqrt(si2 .* (((n - 1) * S1i - Wi.^2)/(n - 2)))

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

全部回复

ywh19860616

2011-10-30 12:38:05

这个程序哪里不对

%输入一个列向量x（n行1列）和矩阵（n行n列）
x=input('x=')
weights=input('weights=')
n=length(x)
xibar=(sum(x)*ones(n,1)-x)/(n-1)
si2=((sum(x.^2).*ones(n,1)-x.^2)/(n-1))-xibar.^2
%Wi=weights(:)'
Wi=sum(weights')
Wi2=Wi.^2
S1i=sum(weights.^2')
lx=Wi*x
res=lx-Wi*xibar
res =res./sqrt(si2 .* (((n - 1).*S1i-Wi2)./(n - 2)))

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

epoh

2011-10-30 21:17:12

library(spdep)
example(columbus)
col.listw <- nb2listw(col.gal.nb)
x=seq(1:49)
n<-length(x)
xibar <- (rep(sum(x), n) - x)/(n - 1)
si2 <- ((rep(sum(x^2), n) - x^2)/(n - 1)) - xibar^2
Wi <- sapply(col.listw$weights, sum)
#[1] 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
#[37] 1 1 1 1 1 1 1 1 1 1 1 1 1
Wi2=Wi^2
S1i <- sapply(col.listw$weights, function(x) sum(x^2))
# [1] 0.5000000 0.3333333 0.2500000 0.2500000 0.1428571 0.5000000 0.2500000
# [8] 0.1666667 0.1250000 0.2500000 0.2000000 0.1666667 0.2500000 0.1666667
#[15] 0.1666667 0.1428571 0.3333333 0.2500000 0.3333333 0.1000000 0.3333333
#[22] 0.1666667 0.3333333 0.1428571 0.1428571 0.1666667 0.2500000 0.1111111
#[29] 0.1428571 0.2500000 0.5000000 0.2500000 0.2500000 0.2500000 0.1428571
#[36] 0.2000000 0.1666667 0.1666667 0.5000000 0.2000000 0.3333333 0.5000000
#[43] 0.1666667 0.2000000 0.2500000 0.5000000 0.5000000 0.2500000 0.3333333
lx <- lag.listw(col.listw, x)
res <- (lx - Wi * xibar)
res <- res/sqrt(si2 * (((n - 1) * S1i - Wi^2)/(n - 2)))
res
[1]  -2.37332110 -2.90731303 -3.33404263 -3.10096296 -3.52712088 -1.87601018
[7]  -2.00517335 -3.07281364 -2.00272695 -1.21686333 -2.32942618 -2.51955856
[13] -2.18641859 -2.01890149 -1.80019026 -1.63006479 -0.93223336 -1.00013383
[19] -0.80158055 -0.07496224  0.52694104 -0.91849373 -0.25298002 -0.48562024
[25] -0.79300074 -0.66119098  0.11454945  1.21862840  0.95192246  0.41343356
[31]  1.01485374  0.89135794  0.38706277  1.11682531  2.06289781  2.24404180
[37]  1.96425487  2.06215004  1.64678694  1.70734475  1.88352479  1.05182914
[43]  3.05170115  3.00125143  2.90624770  1.32511796  1.63837602  3.08496360
[49]  2.70457746
############

但是function lag.listw()

要转成matlab比较麻烦

因为lag.listw()是调用C function

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

ywh19860616

2011-10-30 21:33:00

epoh老师，我是想根据附件中的4个公式计算Gi（d），其中已知x和w（为一个矩阵）。
第（1）个公式是第（2）个公式的标准化
呵呵，我看这个很好实现的，但是程序总是运行不出来

附件列表

原图尺寸 13.9 KB

原图尺寸 11.21 KB

原图尺寸 10.53 KB

原图尺寸 3.92 KB

原图尺寸 4 KB

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

ywh19860616

2011-10-30 21:34:55

我写的程序：
%输入一个列向量x（n行1列）和矩阵（n行n列）
x=input('x=')
weights=input('weights=')
n=length(x)
xibar=(sum(x)*ones(n,1)-x)/(n-1)
si2=((sum(x.^2).*ones(n,1)-x.^2)/(n-1))-xibar.^2
%Wi=weights(:)'
Wi=sum(weights')
Wi2=Wi.^2
S1i=sum(weights.^2')
lx=Wi*x
res=lx-Wi*xibar
res =res./sqrt(si2 .* (((n - 1).*S1i-Wi2)./(n - 2)))

这里我暂时假定了G与d无关，比较好编写。而且假定了w的对角线全为0
后面可能还要与d相关，所以如果用R比较难实现，因此想转换为matlab

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

epoh

2011-10-30 22:07:54

未命名3.jpg & 未命名2.jpg

有没矛盾?两个都是 j not equal to i

Wij(d)是甚么关系?要知道

如果公式对,就有方向了

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

点击查看更多内容…

ywh19860616

2011-10-30 23:11:22

epoh 发表于 2011-10-30 22:07
未命名3.jpg & 未命名2.jpg 有没矛盾?两个都是 j not equal to i Wij(d)是甚么关系?要知道如果公式对,就有方 ...

epoh老师，未命名（3）与（2）没有矛盾，（3）是（2）的标准化形式。
Wij（d）是一个矩阵，为已知的。
若x为n*1的列向量，Wij（d）就为n*n的矩阵，且x也为已知数。
附件中为参考文章

reference.rar
大小:(1.11 MB)

马上下载

本附件包括：

reference.pdf

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

epoh

2011-10-31 09:21:44

function [GI,ZG]=Getis(Wij,Xij)

% 计算G统计值：G值是不包含自己的G统计值

% 输入：Wij—空间权值矩阵

% 输入：Xij—研究区域的空间属性数据

% 输出：GI— 空间局部自相关指标-G统计值

% 输出：ZG— 对于GIJ的检验Z值的计算结果

这里有现成程序,请先看看.

局部G统计计算的Matlab程序.pdf

局部G统计计算的Matlab程序.pdf
大小:(256.93 KB)

马上下载

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

ywh19860616

2011-10-31 14:07:24

epoh老师，谢谢您，我先用这个程序试下
如需修改，而自己不会，再找您

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

epoh

2011-10-31 14:47:41

ct_data.mat,distance_wm.m,Getis.m

Getis.rar
大小:(2.92 KB)

马上下载

distance_wm.m from
  http://community.wvu.edu/~djl041/matlab.html
%%%%%%run script
load ct_data % Load Connecticut county latitude and longitude coordinates
                     % 8 counties in all
                     % http://en.wikipedia.org/wiki/List_of_counties_in_Connecticut
result = distance_wm(yc,xc); % Make sure that you enter latitude and longitude
                                             % in that order
Xij=result.dista                      % This martrix contains the pairwise distances
Wij=result.dw                         % This matrix contains the connectivity structure
                                             % 1 = neighbor, 0 = not a neighbor
[GI,ZG]=Getis(Wij,Xij)
%%%%%%%%%%
Getis_revised(2011.10.31  21:03)

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

ywh19860616

2011-10-31 20:15:30

epoh老师，我想显示出在Getis程序中的E(G(i))，但是总是显示最后一个，而前面的都为0。
如
form=1:N; % compute E]Gi(d)^
for n=1:N
  if n==m
   WI=WI;
  else
   WI=WI+Wij(m,n);
  end
end %end for
EG(m)=WI/(N-1);

我把最后面的一句改为如下循环，也出不来想要结果
for m=1:N
EG(m)=WI/(N-1);
end
EG

Getis
GI =
0.0847
0.1061
0.1861
0.2221
0.2074
0.1125
0.1555
0.1262

ZG =
5.3371
6.6847
11.7238
13.9902
13.0668
7.0897
9.7949
-1.0495

EG =
      0
      0
      0
      0
      0
      0
      0
0.1429

EG不应该前面都为0的，我不知道如何取出前面的循环结果

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

epoh

2011-10-31 21:00:12

哈哈!抱歉!抱歉!

不是你的错

是我的错!

copy 过来的时候,产生的错误

Getis_revised.rar
大小:(531 Bytes)

马上下载

本附件包括：

Getis_revised.m

%%%%%%%%result
Xij =

      0 26.3976 46.9654 64.7889 62.6534 36.7131 43.1493 30.9348
26.3976       0 20.6691 38.3926 39.3141 50.1330 25.8122 29.2114
46.9654 20.6691       0 18.0368 28.1406 67.6762 28.9354 43.4492
64.7889 38.3926 18.0368       0 24.4042 82.0487 37.9407 56.7230
62.6534 39.3141 28.1406 24.4042       0 68.3261 22.4538 43.2943
36.7131 50.1330 67.6762 82.0487 68.3261       0 46.2011 25.5506
43.1493 25.8122 28.9354 37.9407 22.4538 46.2011       0 20.8906
30.9348 29.2114 43.4492 56.7230 43.2943 25.5506 20.8906       0

Wij =

   0    1    0    0    0    0    0    0
   1    0    1    0    0    0    1    0
   0    1    0    1    0    0    0    0
   0    0    1    0    1    0    0    0
   0    0    0    1    0    0    1    0
   0    0    0    0    0    0    0    1
   0    1    0    0    1    0    0    1
   0    0    0    0    0    1    1    0

GI =

0.0847
0.3160
0.3446
0.4441
0.4336
0.1125
0.4470
0.2845

EG =

0.1429
0.4286
0.2857
0.2857
0.2857
0.1429
0.4286
0.2857

VarG =

1.0e-003 *

0.2520
0.5039
0.4199
0.4199
0.4199
0.2520
0.5039
0.4199

ZG =

-3.6629
-5.0163
2.8722
7.7309
7.2155
-1.9103
0.8188
-0.0571

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

ywh19860616

2011-10-31 21:09:40

呵呵，谢谢epoh老师了，又帮忙解决了一个难题
我只看到了最后部分漏了一些语句，没有注意中间了
运行也没有出错，也就没有注意找了

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

ywh19860616

2012-5-9 14:16:00

epoh 发表于 2011-10-31 21:00
哈哈!抱歉!抱歉!不是你的错是我的错!copy 过来的时候,产生的错误 %%%%%%%%result
Xij = 0 26.3 ...

epoh老师，您好
我想请教您下面的程序应该如何传递参数的？
在您帮忙编写的Getis_revised：
function [GI,ZG,ZZ]=Getis(Wij,Xij)
在您基础上，我加了最后一句计算ZZ的ZZ=sum(abs(ZG))
根据这个函数function [GI,ZG,ZZ]=Getis(Wij,Xij)，我需要输入Wij和Xij。现在我想实现如下操作，即这里的Wij是可以有多种取值的，然后根据哪个Wij得出的ZZ值最大，选得最终的Wij。
比如我先用下面程序计算多个Wij：
load ct_data
dista = GeogDistance(yc,xc);
for nn=1:28
Wij=exp(-nn*dista);
end
Wij

那么如何传递到function [GI,ZG,ZZ]=Getis(Wij,Xij)？

% GeogDistance(yc,xc) 这个函数是您上传的distance_wm.m里面一部分。

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

epoh

2012-5-9 15:41:57

ywh19860616 发表于 2012-5-9 14:16
epoh老师，您好
我想请教您下面的程序应该如何传递参数的？
在您帮忙编写的Getis_revised：

load ct_data % Load Connecticut county latitude and longitude coordinates
                  % 8 counties in all
                  % http://en.wikipedia.org/wiki/List_of_counties_in_Connecticut

result = distance_wm(yc,xc); % Make sure that you enter latitude and longitude
                                             % in that order
Xij=result.dista                % This martrix contains the pairwise distances
Wij=result.dw                % This matrix contains the connectivity structure
                                    % 1 = neighbor, 0 = not a neighbor
[GI,ZG]=Getis(Wij,Xij)

%%%%%%
dista = GeogDistance(yc,xc);
for nn=1:28
   WWij=exp(-nn*dista);
end
WWij

[GGI,ZZG,ZZ]=Getis1(WWij,Xij)  %ZZ=sum(abs(ZG))

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

ywh19860616

2012-5-9 16:06:42

epoh 发表于 2012-5-9 15:41
load ct_data % Load Connecticut county latitude and longitude coordinates
...

epoh老师，您上面的[GGI,ZZG,ZZ]=Getis1(WWij,Xij)作了哪些修改吗？
我运行下面程序只得到一个ZZ值

load ct_data
result = distance_wm(yc,xc);
Xij=result.dista;
dista = GeogDistance(yc,xc);
for nn=1:28
WWij=exp(-nn*dista);
end
WWij
[GGI,ZZG,ZZ]=Getis_revised(WWij,Xij)

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

epoh

2012-5-9 16:17:02

ywh19860616 发表于 2012-5-9 16:06
epoh老师，您上面的[GGI,ZZG,ZZ]=Getis1(WWij,Xij)作了哪些修改吗？
我运行下面程序只得到一个ZZ值

呵呵,补充一下
Getis_revised我更名为Getis

Getis1就是在Getis_revised最后加一句ZZ=sum(abs(ZG))

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

ywh19860616

2012-5-9 16:24:06

epoh 发表于 2012-5-9 16:17
呵呵,补充一下
Getis_revised我更名为Getis

恩，谢谢epoh老师。不知道为何我运行16楼的程序只得出对应一个Wij的结果。

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

epoh

2012-5-9 16:43:56

ywh19860616 发表于 2012-5-9 16:24
恩，谢谢epoh老师。不知道为何我运行16楼的程序只得出对应一个Wij的结果。

Getis_ZZ.rar
大小:(4.12 KB)

马上下载

load ct_data

result = distance_wm(yc,xc);
Xij=result.dista
Wij=result.dw

[GI,ZG]=Getis(Wij,Xij)

%%%%%%
dista = GeogDistance(yc,xc);
for nn=1:28
   WWij=exp(-nn*dista);
end
WWij

[GGI,ZZG,ZZ]=Getis1(WWij,Xij)  %ZZ=sum(abs(ZG))

%%%%%%%%%%%
Xij =

      0 26.3976 46.9654 64.7889 62.6534 36.7131 43.1493 30.9348
26.3976       0 20.6691 38.3926 39.3141 50.1330 25.8122 29.2114
46.9654 20.6691       0 18.0368 28.1406 67.6762 28.9354 43.4492
64.7889 38.3926 18.0368       0 24.4042 82.0487 37.9407 56.7230
62.6534 39.3141 28.1406 24.4042       0 68.3261 22.4538 43.2943
36.7131 50.1330 67.6762 82.0487 68.3261       0 46.2011 25.5506
43.1493 25.8122 28.9354 37.9407 22.4538 46.2011       0 20.8906
30.9348 29.2114 43.4492 56.7230 43.2943 25.5506 20.8906       0

Wij =

   0    1    0    0    0    0    0    0
   1    0    1    0    0    0    1    0
   0    1    0    1    0    0    0    0
   0    0    1    0    1    0    0    0
   0    0    0    1    0    0    1    0
   0    0    0    0    0    0    0    1
   0    1    0    0    1    0    0    1
   0    0    0    0    0    1    1    0

GI =

0.0847
0.3160
0.3446
0.4441
0.4336
0.1125
0.4470
0.2845

ZG =

-3.6629
-5.0163
2.8722
7.7309
7.2155
-1.9103
0.8188
-0.0571

WWij =

1.0000 0.0000       0       0       0       0       0       0
0.0000 1.0000 0.0000       0       0       0 0.0000       0
      0 0.0000 1.0000 0.0000       0       0       0       0
      0       0 0.0000 1.0000 0.0000       0       0       0
      0       0       0 0.0000 1.0000       0 0.0000       0
      0       0       0       0       0 1.0000       0 0.0000
      0 0.0000       0       0 0.0000       0 1.0000 0.0000
      0       0       0       0       0 0.0000 0.0000 1.0000

GGI =

  1.0e-219 *

0.0000
0.0000
0.1141
0.0887
0.0000
0.0000
0.0000
0.0000

ZZG =

  1.0e-108 *

   -Inf
0.0000
0.1284
0.0597
0.0000
-0.0000
-0.0000
0.0000

ZZ =

Inf

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

ywh19860616

2012-5-9 18:34:52

epoh 发表于 2012-5-9 16:43
load ct_data

result = distance_wm(yc,xc);

呵呵，epoh老师，非常抱歉，可能我没有说清楚问题。
您上面的程序得出的还只是一个ZZ值，而我想要的其实是根据多个WWij(i,j)计算出多个ZZ值，然后再取sum(abs(ZZ))最大值所对应的那个WWij(i,j)。

而我WWij的生成公式为：
WWij(i,j)=exp(-k*dista(i,j)) %其中dista(i,j)就是dista = GeogDistance(yc,xc); 而我这里变化k，所以WWij(i,j)应该是多个的。

计算出WWij(i,j)后再传递到函数Getis_revised(WWij(jj),Xij)，这样每个WWij(i,j)都可以得到ZZ值的。

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

epoh

2012-5-9 19:45:41

ywh19860616 发表于 2012-5-9 18:34
呵呵，epoh老师，非常抱歉，可能我没有说清楚问题。
您上面的程序得出的还只是一个ZZ值，而我想要的其 ...

load ct_data
result = distance_wm(yc,xc);
Xij=result.dista ;
Wij=result.dw ;
dista = GeogDistance(yc,xc);
k=[0.6,0.7,0.4,0.5,0.8,0.9,1];
s=length(k);
ZZmat=zeros(1,length(k));
[n,m]=size(dista);
WW=zeros(n,m,s);
for s =1:length(k)
WWij=zeros(n,m);
for i = 1:n
for j = 1:m
   WWij(i,j)=exp(-k(s)*dista(i,j));
end
   WW(:,:,s)=WWij;
end %end j
[GGI,ZZG,ZZ]=Getis1(WWij,Xij);
ZZmat(s)=ZZ;
end %end s
ZZmat
index=find(ZZmat==max(ZZmat))
WW
WW(:,:,index)
%%%%%%%%%%%%%%%%%%%%%%
ZZmat =

0.0447 0.0175 0.3017 0.1153 0.0070 0.0028 0.0011

index =

   3

WW(:,:,1) =

1.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000
0.0000 1.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000
0.0000 0.0000 1.0000 0.0000 0.0000 0.0000 0.0000 0.0000
0.0000 0.0000 0.0000 1.0000 0.0000 0.0000 0.0000 0.0000
0.0000 0.0000 0.0000 0.0000 1.0000 0.0000 0.0000 0.0000
0.0000 0.0000 0.0000 0.0000 0.0000 1.0000 0.0000 0.0000
0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 1.0000 0.0000
0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 1.0000

WW(:,:,2) =

1.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000
0.0000 1.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000
0.0000 0.0000 1.0000 0.0000 0.0000 0.0000 0.0000 0.0000
0.0000 0.0000 0.0000 1.0000 0.0000 0.0000 0.0000 0.0000
0.0000 0.0000 0.0000 0.0000 1.0000 0.0000 0.0000 0.0000
0.0000 0.0000 0.0000 0.0000 0.0000 1.0000 0.0000 0.0000
0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 1.0000 0.0000
0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 1.0000

WW(:,:,3) =

1.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000
0.0000 1.0000 0.0003 0.0000 0.0000 0.0000 0.0000 0.0000
0.0000 0.0003 1.0000 0.0007 0.0000 0.0000 0.0000 0.0000
0.0000 0.0000 0.0007 1.0000 0.0001 0.0000 0.0000 0.0000
0.0000 0.0000 0.0000 0.0001 1.0000 0.0000 0.0001 0.0000
0.0000 0.0000 0.0000 0.0000 0.0000 1.0000 0.0000 0.0000
0.0000 0.0000 0.0000 0.0000 0.0001 0.0000 1.0000 0.0002
0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0002 1.0000

WW(:,:,4) =

1.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000
0.0000 1.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000
0.0000 0.0000 1.0000 0.0001 0.0000 0.0000 0.0000 0.0000
0.0000 0.0000 0.0001 1.0000 0.0000 0.0000 0.0000 0.0000
0.0000 0.0000 0.0000 0.0000 1.0000 0.0000 0.0000 0.0000
0.0000 0.0000 0.0000 0.0000 0.0000 1.0000 0.0000 0.0000
0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 1.0000 0.0000
0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 1.0000

WW(:,:,5) =

1.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000
0.0000 1.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000
0.0000 0.0000 1.0000 0.0000 0.0000 0.0000 0.0000 0.0000
0.0000 0.0000 0.0000 1.0000 0.0000 0.0000 0.0000 0.0000
0.0000 0.0000 0.0000 0.0000 1.0000 0.0000 0.0000 0.0000
0.0000 0.0000 0.0000 0.0000 0.0000 1.0000 0.0000 0.0000
0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 1.0000 0.0000
0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 1.0000

WW(:,:,6) =

1.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000
0.0000 1.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000
0.0000 0.0000 1.0000 0.0000 0.0000 0.0000 0.0000 0.0000
0.0000 0.0000 0.0000 1.0000 0.0000 0.0000 0.0000 0.0000
0.0000 0.0000 0.0000 0.0000 1.0000 0.0000 0.0000 0.0000
0.0000 0.0000 0.0000 0.0000 0.0000 1.0000 0.0000 0.0000
0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 1.0000 0.0000
0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 1.0000

WW(:,:,7) =

1.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000
0.0000 1.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000
0.0000 0.0000 1.0000 0.0000 0.0000 0.0000 0.0000 0.0000
0.0000 0.0000 0.0000 1.0000 0.0000 0.0000 0.0000 0.0000
0.0000 0.0000 0.0000 0.0000 1.0000 0.0000 0.0000 0.0000
0.0000 0.0000 0.0000 0.0000 0.0000 1.0000 0.0000 0.0000
0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 1.0000 0.0000
0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 1.0000

ans =

1.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000
0.0000 1.0000 0.0003 0.0000 0.0000 0.0000 0.0000 0.0000
0.0000 0.0003 1.0000 0.0007 0.0000 0.0000 0.0000 0.0000
0.0000 0.0000 0.0007 1.0000 0.0001 0.0000 0.0000 0.0000
0.0000 0.0000 0.0000 0.0001 1.0000 0.0000 0.0001 0.0000
0.0000 0.0000 0.0000 0.0000 0.0000 1.0000 0.0000 0.0000
0.0000 0.0000 0.0000 0.0000 0.0001 0.0000 1.0000 0.0002
0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0002 1.0000

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

ywh19860616

2012-5-9 20:11:18

epoh 发表于 2012-5-9 19:45
load ct_data
result = distance_wm(yc,xc);
Xij=result.dista ;

谢谢epoh老师，呵呵，循环有点复杂
我原来把s那个循环和i，j循环放在一起，总是出来一个结果。
现在程序运行很好。

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

ywh19860616

2012-5-10 18:55:09

epoh 发表于 2012-5-9 19:45
load ct_data
result = distance_wm(yc,xc);
Xij=result.dista ;

epoh老师，您好。今天利用这个程序试算了一下，感觉有点问题，但是我找不出哪里。
您这里取 k=[0.6,0.7,0.4,0.5,0.8,0.9,1]; 所以运行出来index=3，即最后选择了第三个k值，k=0.4
也是这组k值中最小的。

我试了好多k值，比如k=0.01:0.1:1，得出的结果永远是k值最小的那个，按照理论应该不可能这样的，不知道问题出在哪？

起初我想 Xij应该是一个变量，数据应该是一个行向量或者列向量，所以我试了自己找的数据，但是也是一样的结果。
就是实现这篇文章第5部分提到的：Spatial Filtering with EViews and MATLAB http://www.stat.tugraz.at/AJS/ausg071/071Ferstl.pdf

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

epoh

2012-5-10 19:54:09

ywh19860616 发表于 2012-5-10 18:55
epoh老师，您好。今天利用这个程序试算了一下，感觉有点问题，但是我找不出哪里。
您这里取 k=[0.6,0. ...

load ct_data
result = distance_wm(yc,xc);
Xij=result.dista ;
Wij=result.dw ;
dista = GeogDistance(yc,xc);
k=0.01:0.01:0.09
s=length(k);
ZZmat=zeros(1,length(k));
[n,m]=size(dista);
WW=zeros(n,m,s);
for s =1:length(k)
WWij=zeros(n,m);
for i = 1:n
for j = 1:m
   if (i~=j)
      WWij(i,j)=exp(-k(s)*dista(i,j));
   elseif (i==j)
      WWij(i,j)=0;
   end %end if
end  %end i
   WW(:,:,s)=WWij;
end %end j
[GGI,ZZG,ZZ]=Getis1(WWij,Xij);
ZZmat(s)=ZZ;
end %end s
ZZmat
index=find(ZZmat==max(ZZmat))
WW
WW(:,:,index)

%%%%%%%%%%%%%%%%%%%%%
k =

0.0100 0.0200 0.0300 0.0400 0.0500 0.0600 0.0700 0.0800 0.0900

ZZmat =

10.5321 13.1748 14.1048 14.1081 13.5731 12.7366 11.7508 10.7131 9.6844

index =

   4

WW(:,:,1) =

      0 0.7680 0.6252 0.5231 0.5344 0.6927 0.6495 0.7339
0.7680       0 0.8133 0.6812 0.6749 0.6057 0.7725 0.7467
0.6252 0.8133       0 0.8350 0.7547 0.5083 0.7487 0.6476
0.5231 0.6812 0.8350       0 0.7835 0.4402 0.6843 0.5671
0.5344 0.6749 0.7547 0.7835       0 0.5050 0.7989 0.6486
0.6927 0.6057 0.5083 0.4402 0.5050       0 0.6300 0.7745
0.6495 0.7725 0.7487 0.6843 0.7989 0.6300       0 0.8115
0.7339 0.7467 0.6476 0.5671 0.6486 0.7745 0.8115       0

WW(:,:,2) =

      0 0.5898 0.3909 0.2737 0.2856 0.4799 0.4219 0.5386
0.5898       0 0.6614 0.4640 0.4555 0.3669 0.5968 0.5575
0.3909 0.6614       0 0.6972 0.5696 0.2583 0.5606 0.4194
0.2737 0.4640 0.6972       0 0.6138 0.1938 0.4682 0.3216
0.2856 0.4555 0.5696 0.6138       0 0.2550 0.6382 0.4207
0.4799 0.3669 0.2583 0.1938 0.2550       0 0.3969 0.5999
0.4219 0.5968 0.5606 0.4682 0.6382 0.3969       0 0.6585
0.5386 0.5575 0.4194 0.3216 0.4207 0.5999 0.6585       0

WW(:,:,3) =

      0 0.4530 0.2444 0.1432 0.1527 0.3324 0.2740 0.3953
0.4530       0 0.5379 0.3161 0.3075 0.2222 0.4610 0.4163
0.2444 0.5379       0 0.5821 0.4299 0.1313 0.4198 0.2716
0.1432 0.3161 0.5821       0 0.4809 0.0853 0.3204 0.1824
0.1527 0.3075 0.4299 0.4809       0 0.1288 0.5099 0.2729
0.3324 0.2222 0.1313 0.0853 0.1288       0 0.2501 0.4646
0.2740 0.4610 0.4198 0.3204 0.5099 0.2501       0 0.5343
0.3953 0.4163 0.2716 0.1824 0.2729 0.4646 0.5343       0

WW(:,:,4) =

      0 0.3479 0.1528 0.0749 0.0816 0.2303 0.1780 0.2901
0.3479       0 0.4375 0.2153 0.2075 0.1346 0.3561 0.3108
0.1528 0.4375       0 0.4860 0.3244 0.0667 0.3143 0.1759
0.0749 0.2153 0.4860       0 0.3768 0.0376 0.2192 0.1034
0.0816 0.2075 0.3244 0.3768       0 0.0650 0.4073 0.1770
0.2303 0.1346 0.0667 0.0376 0.0650       0 0.1575 0.3599
0.1780 0.3561 0.3143 0.2192 0.4073 0.1575       0 0.4336
0.2901 0.3108 0.1759 0.1034 0.1770 0.3599 0.4336       0

WW(:,:,5) =

      0 0.2672 0.0955 0.0392 0.0436 0.1595 0.1156 0.2129
0.2672       0 0.3558 0.1467 0.1401 0.0815 0.2751 0.2321
0.0955 0.3558       0 0.4058 0.2449 0.0339 0.2353 0.1139
0.0392 0.1467 0.4058       0 0.2952 0.0165 0.1500 0.0587
0.0436 0.1401 0.2449 0.2952       0 0.0328 0.3254 0.1148
0.1595 0.0815 0.0339 0.0165 0.0328       0 0.0993 0.2787
0.1156 0.2751 0.2353 0.1500 0.3254 0.0993       0 0.3519
0.2129 0.2321 0.1139 0.0587 0.1148 0.2787 0.3519       0

WW(:,:,6) =

      0 0.2052 0.0597 0.0205 0.0233 0.1105 0.0751 0.1563
0.2052       0 0.2893 0.0999 0.0945 0.0494 0.2125 0.1733
0.0597 0.2893       0 0.3388 0.1848 0.0172 0.1762 0.0738
0.0205 0.0999 0.3388       0 0.2313 0.0073 0.1026 0.0333
0.0233 0.0945 0.1848 0.2313       0 0.0166 0.2600 0.0744
0.1105 0.0494 0.0172 0.0073 0.0166       0 0.0625 0.2159
0.0751 0.2125 0.1762 0.1026 0.2600 0.0625       0 0.2855
0.1563 0.1733 0.0738 0.0333 0.0744 0.2159 0.2855       0

WW(:,:,7) =

      0 0.1576 0.0373 0.0107 0.0125 0.0765 0.0488 0.1147
0.1576       0 0.2353 0.0681 0.0638 0.0299 0.1642 0.1294
0.0373 0.2353       0 0.2829 0.1395 0.0088 0.1319 0.0478
0.0107 0.0681 0.2829       0 0.1812 0.0032 0.0702 0.0189
0.0125 0.0638 0.1395 0.1812       0 0.0084 0.2077 0.0483
0.0765 0.0299 0.0088 0.0032 0.0084       0 0.0394 0.1672
0.0488 0.1642 0.1319 0.0702 0.2077 0.0394       0 0.2317
0.1147 0.1294 0.0478 0.0189 0.0483 0.1672 0.2317       0

WW(:,:,8) =

      0 0.1210 0.0233 0.0056 0.0067 0.0530 0.0317 0.0842
0.1210       0 0.1914 0.0464 0.0431 0.0181 0.1268 0.0966
0.0233 0.1914       0 0.2362 0.1053 0.0045 0.0988 0.0309
0.0056 0.0464 0.2362       0 0.1419 0.0014 0.0481 0.0107
0.0067 0.0431 0.1053 0.1419       0 0.0042 0.1659 0.0313
0.0530 0.0181 0.0045 0.0014 0.0042       0 0.0248 0.1295
0.0317 0.1268 0.0988 0.0481 0.1659 0.0248       0 0.1880
0.0842 0.0966 0.0309 0.0107 0.0313 0.1295 0.1880       0

WW(:,:,9) =

      0 0.0929 0.0146 0.0029 0.0036 0.0367 0.0206 0.0618
0.0929       0 0.1556 0.0316 0.0291 0.0110 0.0980 0.0721
0.0146 0.1556       0 0.1972 0.0794 0.0023 0.0740 0.0200
0.0029 0.0316 0.1972       0 0.1112 0.0006 0.0329 0.0061
0.0036 0.0291 0.0794 0.1112       0 0.0021 0.1325 0.0203
0.0367 0.0110 0.0023 0.0006 0.0021       0 0.0156 0.1003
0.0206 0.0980 0.0740 0.0329 0.1325 0.0156       0 0.1526
0.0618 0.0721 0.0200 0.0061 0.0203 0.1003 0.1526       0

ans =

      0 0.3479 0.1528 0.0749 0.0816 0.2303 0.1780 0.2901
0.3479       0 0.4375 0.2153 0.2075 0.1346 0.3561 0.3108
0.1528 0.4375       0 0.4860 0.3244 0.0667 0.3143 0.1759
0.0749 0.2153 0.4860       0 0.3768 0.0376 0.2192 0.1034
0.0816 0.2075 0.3244 0.3768       0 0.0650 0.4073 0.1770
0.2303 0.1346 0.0667 0.0376 0.0650       0 0.1575 0.3599
0.1780 0.3561 0.3143 0.2192 0.4073 0.1575       0 0.4336
0.2901 0.3108 0.1759 0.1034 0.1770 0.3599 0.4336       0

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

ywh19860616

2012-5-10 20:27:50

epoh 发表于 2012-5-10 19:54
load ct_data
result = distance_wm(yc,xc);
Xij=result.dista ;

epoh老师，谢谢
还有一个问题再请教您，上面程序的Xij其实为8*8的矩阵。
这两组命令得出的结果是一样的：
.....
Xij=result.dista;
.....
与
.....
Xij1=result.dista ;
Xij=Xij1(1,:);
......
也就是说其实现在的命令也只取了Xij中的第一行数据运行。如果我想对每一行数据都执行同样操作，应该在哪里加入循环？
比如循环操作：
for kk=1:8 % 8即为Xij的行数。
[GGI,ZZG,ZZ]=Getis_revised(WWij,Xij(kk,:));
end

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

ywh19860616

2012-5-10 20:44:32

epoh老师，您看我下面执行结果正常吗？
我的想法是先对求出的Wij执行标准化，命令更改为：
......
for i = 1:n
for j = 1:m
   if (i~=j)
      WWij(i,j)=exp(-k(s)*dista(i,j));
   elseif (i==j)
      WWij(i,j)=0;
   end %end if
end  %end i
   WW(:,:,s)=WWij;
end %end j
WW1ij=normw(WWij);  %normw是jvpl7工具箱的一个函数
[GGI,ZZG,ZZ]=Getis1(WW1ij,Xij);
ZZmat(s)=ZZ;
end %end s
......
k =
0.0100 0.0200 0.0300 0.0400 0.0500 0.0600 0.0700 0.0800 0.0900

ZZmat =
4.2811 8.1618 11.5776 14.5126 16.9869 19.0427 20.7317 22.1072 23.2194

index =
   9
我尝试了几组数据，结果显示ZZmat总是从小到大排列，感觉不对。

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

epoh

2012-5-10 21:57:19

ywh19860616 发表于 2012-5-10 20:44
epoh老师，您看我下面执行结果正常吗？
我的想法是先对求出的Wij执行标准化，命令更改为：
......

distance_wm.m已经标准化
对求出的Wij执行标准化
恐怕重复

ZZmat =

10.5321 13.1748 14.1048 14.1081 13.5731 12.7366 11.7508 10.7131 9.6844
7.0585 8.2455 8.2969 7.8566 7.2101 6.5026 5.8074 5.2754 4.7729
14.8975 18.3559 19.4375 19.3021 18.4963 17.3349 16.0100 14.6396 13.2945
16.2223 19.8581   20.8892 20.6055 19.6146 18.2636 16.7615 15.2342 13.7550
13.6478 16.6916 17.5135 17.2049 16.2868 15.0609 13.7105 12.3468 11.0344
13.8836 16.9895 17.8799 17.6583 16.8396 15.7157 14.4612 13.1810 11.9363
7.5528 8.9430 9.1190 8.7473 8.1289 7.4203 6.7042 6.0220 5.3923
10.1812 13.0251 14.2785 14.6337 14.4292 13.8759 13.1149 12.2428 11.3248

WW(:,:,C)

ans =

      0 0.4530 0.2444 0.1432 0.1527 0.3324 0.2740 0.3953
0.4530       0 0.5379 0.3161 0.3075 0.2222 0.4610 0.4163
0.2444 0.5379       0 0.5821 0.4299 0.1313 0.4198 0.2716
0.1432 0.3161 0.5821       0 0.4809 0.0853 0.3204 0.1824
0.1527 0.3075 0.4299 0.4809       0 0.1288 0.5099 0.2729
0.3324 0.2222 0.1313 0.0853 0.1288       0 0.2501 0.4646
0.2740 0.4610 0.4198 0.3204 0.5099 0.2501       0 0.5343
0.3953 0.4163 0.2716 0.1824 0.2729 0.4646 0.5343       0

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

ywh19860616

2012-5-10 22:36:51

epoh 发表于 2012-5-10 21:57
distance_wm.m已经标准化
对求出的Wij执行标准化
恐怕重复

epoh老师，谢谢。能否根据ZZmat的每行都分别找出最大值的位置，而不是对ZZmat矩阵整体找最大值？
其实我是想对Xij的每行分别运行一次的，但是这样操作量太大，因此才让您帮忙修改的，用循环实现。

是的，我明白您说的，程序中的确提到：
results.W    = (n x n) sparse row-standardized weight matrix
例如，我把results.W和normw函数得出的结果对比了，不一致：
load ct_data
result = distance_wm(yc,xc);
Xij=result.dista
Xij =
      0 26.3976 46.9654 64.7889 62.6534 36.7131 43.1493 30.9348
26.3976       0 20.6691 38.3926 39.3141 50.1330 25.8122 29.2114
46.9654 20.6691       0 18.0368 28.1406 67.6762 28.9354 43.4492
64.7889 38.3926 18.0368       0 24.4042 82.0487 37.9407 56.7230
62.6534 39.3141 28.1406 24.4042       0 68.3261 22.4538 43.2943
36.7131 50.1330 67.6762 82.0487 68.3261       0 46.2011 25.5506
43.1493 25.8122 28.9354 37.9407 22.4538 46.2011       0 20.8906
30.9348 29.2114 43.4492 56.7230 43.2943 25.5506 20.8906       0
Xij_1=result.W
Xij_1=
(2,1)    0.3333
(1,2)    1.0000
(3,2)    0.5000
(7,2)    0.3333
(2,3)    0.3333
(4,3)    0.5000
(3,4)    0.5000
(5,4)    0.5000
(4,5)    0.5000
(7,5)    0.3333
(8,6)    0.5000
(2,7)    0.3333
(5,7)    0.5000
(8,7)    0.5000
(6,8)    1.0000
(7,8)    0.3333

Xij_row=normw(Xij)

Xij_row =
      0 0.0847 0.1507 0.2079 0.2011 0.1178 0.1385 0.0993
0.1148       0 0.0899 0.1670 0.1710 0.2180 0.1123 0.1270
0.1850 0.0814       0 0.0710 0.1108 0.2666 0.1140 0.1711
0.2010 0.1191 0.0560       0 0.0757 0.2545 0.1177 0.1760
0.2171 0.1362 0.0975 0.0846       0 0.2368 0.0778 0.1500
0.0975 0.1331 0.1797 0.2178 0.1814       0 0.1227 0.0678
0.1914 0.1145 0.1284 0.1683 0.0996 0.2050       0 0.0927
0.1237 0.1168 0.1738 0.2268 0.1731 0.1022 0.0835       0

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

epoh

2012-5-10 22:46:53

ywh19860616 发表于 2012-5-10 22:36
epoh老师，谢谢。能否根据ZZmat的每行都分别找出最大值的位置，而不是对ZZmat矩阵整体找最大值？
其实 ...

我没看懂底下怎么得来
Xij_1=result.W
Xij_1=
(2,1)    0.3333
(1,2)    1.0000
(3,2)    0.5000
(7,2)    0.3333
(2,3)    0.3333
(4,3)    0.5000
(3,4)    0.5000
(5,4)    0.5000
(4,5)    0.5000
(7,5)    0.3333
(8,6)    0.5000
(2,7)    0.3333
(5,7)    0.5000
(8,7)    0.5000
(6,8)    1.0000
(7,8)    0.3333

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

ywh19860616

2012-5-10 23:18:33

epoh 发表于 2012-5-10 22:46
我没看懂底下怎么得来
Xij_1=result.W
Xij_1=

epoh老师，就是distance_wm函数的返回结果的results.W
%       results.dista = (n x n) matrix of pairwide distances
%       results.W    = (n x n) sparse row-standardized weight matrix
%       results.dw = (n x n) matrix of neighbors (0,1)

您指的原程序已经标准化是指W = sparse(normw(dw)) 吗？我以为是指results.W.
我一开始没有看到这个，呵呵，那我可以直接在distance_wm里取出Wij=normw(results.dista)就行。

另外，epoh老师，能否根据ZZmat的每行都分别找出一个最大值及其对应位置，而不是对ZZmat矩阵整体找最大值？load ct_data
result = distance_wm(yc,xc);
Xij=result.dista ;
[n1,m1]=size(Xij)
Wij=result.Wn ;
dista = GeogDistance(yc,xc);
k=0.01:0.01:0.09;
s=length(k);
ZZmat=zeros(n1,length(k));
[n,m]=size(dista);
WW=zeros(n,m,s);
for s =1:length(k)
WWij=zeros(n,m);
for i = 1:n
for j = 1:m
   if (i~=j)
      WWij(i,j)=exp(-k(s)*dista(i,j));
   elseif (i==j)
      WWij(i,j)=0;
   end %end if
end  %end i
   WW(:,:,s)=WWij;
end %end j
for p = 1:n1
   [GGI,ZZG,ZZ]=Getis_revised(WWij,Xij(p,:));
   ZZmat(p,s)=ZZ;
end % end p
end %end s
ZZmat

for q=1:n1
%[value, location] = max(ZZmat(:))
[value, location] = max(ZZmat(q,:))
%[R,C] = ind2sub(size(ZZmat),location)
end
这样可以取出每一行最大值所在location，但是如何同时根据s（表示变量）和location 取出WW？
比如对应s=1，即数据的第一行，因为k=0.01:0.01:0.09，所以就算出来有9个WW，根据ZZmat第一行最大值，可以取出一个矩阵。
  对应s=2，同时可以取出一个矩阵。
....
      s=9，也取出一个矩阵。

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

分享

扫码加好友，拉您进群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群