大神们，泪求请帮我看看数据框取子集的问题怎么搞？-经管之家官网！

请帮我看看数据框取子集的问题怎么搞？
data4是包含三列的数据框，我想知道这个数据框中x和y值都一样的行数是哪几行？（此时不考虑father值）例如七行和九行值一样，就得到七和九的行数。这里应该有三组是一样的，应返回六个下标。
> data4
yx father
217 15 1
418 14 1
819 15 1
21 17 16 2
81 19 16 2
22 17 17 3
62 18 18 3
82 19 17 3
23 18 18 4
43 19 17 4
83 20 18 4
24 18 19 5
84 20 19 5
25 18 20 6
85 20 20 6
26 18 21 7
66 19 22 7
86 20 21 7
27 19 22 8
47 20 21 8
87 21 22 8
28 19 23 9
88 21 23 9
29 19 24 10
69 20 25 10
89 21 24 10
代码是
x<-c(0,40)
y<-c(0,40)
par(bg="gray")
plot(x,y,type="n",xlab="",ylab="")
for(x in 0:40){for(y in 0:40)
rect(x-.5,y-.5,x+.5,y+.5)}
points(c(15,24),c(18,20),pch=21,bg="yellow")
x[1]=15
y[1]=18
x2=24
y2=20
m=max(abs(y2-y[1]),abs(x2-x[1]))
n=m+1
addy=(y2-y[1])/m
addx=(x2-x[1])/m
i=2
for(i in 2:n){x=x[1]+(i-1)*addx
y=y[1]+(i-1)*addy
y=round(y)}
data<-data.frame(y=y,x=x)
data2=data.frame(x=NULL,Y=NULL,father=NULL)
for(i in 1:10){x<-rep(c((data[i,]$x-1):(data[i,]$x+1)),3)
y<-rep(c((data[i,]$y-1):(data[i,]$y+1)),each=3)
data1<-data.frame(y=y,x=x,father=rep(NA,9))
for(j in 1:9){if (sqrt((data1[j,]$x-data[i,]$x)^2+(data1[j,]$y-data[i,]$y)^2)<=1)
{rect(data1[j,]$x-.5,data1[j,]$y-.5,data1[j,]$x+.5,data1[j,]$y+.5,col="green")
data1[j,]$father=i
data2=data.frame(rbind(data2,data1[j,]))
}
data3=unique(data2)
}}
b=NULL
for(i in 1:dim(data)[1]){
for(j in 1:dim(data3)[1]){
if(data3[j,1]==data[i,1]&data3[j,2]==data[i,2])
b=c(b,j)
}
}
data4=data3[-b, ]

扫码或添加微信号：坛友素质互助

「经管之家」APP：经管人学习、答疑、交友，就上经管之家！
免流量费下载资料----在经管之家app可以下载论坛上的所有资源，并且不额外收取下载高峰期的论坛币。
涵盖所有经管领域的优秀内容----覆盖经济、管理、金融投资、计量统计、数据分析、国贸、财会等专业的学习宝库，各类资料应有尽有。
来自五湖四海的经管达人----已经有上千万的经管人来到这里，你可以找到任何学科方向、有共同话题的朋友。
经管之家（原人大经济论坛），跨越高校的围墙，带你走进经管知识的新世界。
扫描下方二维码下载并注册APP