geom.mlc<-function(x,censor=F){ num<-sum(x)-x[length(x)] inx<-seq(length(x)) if (censor) denom<-sum(inx*x)-inx[length(x)]*x[length(x)]+inx[length(x)-1]*x[length(x)] else denom<-sum(inx*x) return(num/denom) } geom.plotllc<-function(x) { p<-seq(0,1,,500) y<-log(1-p)*sum((seq(length(x))-1)*x)+log(p)*sum(sum(x)-x[length(x)]) plot(p,y,xlab="p",ylab="llike(p)") }