带估计的R中经验分布函数图

给出了一个几何分布公式，其中X可以是0,1,2。。。0和1之间的P是未知的。我需要使用扭矩法找到P的估计值，计算数据的估计值，并将其添加到经验分布函数图中

我加载这样一个数据库-

> a<- read.delim(file.choose("data.txt"))
> 
> dput(a)
structure(list(x = c(4L, 1L, 1L, 0L, 1L, 0L, 0L, 1L, 0L, 0L, 
0L, 0L, 1L, 1L, 1L, 4L, 1L, 2L, 0L, 1L, 1L, 0L, 0L, 0L, 0L, 0L, 
0L, 0L, 0L, 1L, 1L, 3L, 0L, 5L, 2L, 2L, 1L, 0L, 0L, 2L, 0L, 0L, 
0L, 1L, 3L, 3L, 0L, 0L, 0L, 2L, 0L, 2L, 1L, 1L, 4L, 4L, 0L, 1L, 
3L, 1L, 0L, 2L, 1L, 2L, 0L, 0L, 0L, 1L, 0L, 1L, 6L, 0L, 2L, 2L, 
0L, 1L, 1L, 2L, 1L, 0L, 1L, 0L, 3L, 0L, 3L, 0L, 4L, 3L, 2L, 2L, 
2L, 1L, 3L, 0L, 3L, 2L, 0L, 1L, 2L, 1L)), class = "data.frame", row.names = c(NA, 
-100L))

m<-mean(x)
p<- m/(m+1)
geom<-pgeom(q=x,prob=p)
lines(x, sort(geom), type="s", lty=2, col="blue")

>a
>dput（a）
结构（列表x=c（4L、1L、1L、0L、1L、0L、0L、1L、0L、0L、，
0L，0L，1L，1L，1L，4L，1L，2L，0L，1L，1L，0L，0L，0L，0L，0L，0L，0L，
0L、0L、0L、1L、1L、3L、0L、5L、2L、2L、1L、0L、0L、2L、0L、0L、0L、0L、，
0L、1L、3L、3L、0L、0L、2L、0L、2L、1L、1L、4L、4L、0L、1L、，
3L，1L，0L，2L，1L，2L，0L，0L，0L，1L，0L，1L，6L，0L，2L，2L，
0L、1L、1L、2L、1L、0L、1L、0L、3L、0L、3L、0L、4L、3L、2L、2L、，
2L，1L，3L，0L，3L，2L，0L，1L，2L，1L），class=“data.frame”，row.names=c（NA，
-（100升）

然后，我对这件事做了一个计划-

d<-sort(a$x)
plot(d, sort(ecdf(d)(d)),type="s", lty=3,col="red", ylab= "P(X<=x)",ylim= 0:1)

---

dany能帮助别人吗？