Chapter 16: Minimization estimators

Figure 1

mixnum<-10 
set.seed(1)
M<-runif(10)  
si<-0.01
sig<-rep(si,mixnum)
p0<-1/mixnum
p<-p0*rep(1,mixnum)    
N<-c(200)

pcf<-pcf.func("mixt",N,sig=sig,M=M,p=p,support=c(0,1)) 

dp<-draw.pcf(pcf,pnum=c(N))

# frame
   plot(dp$x,dp$y,type="l",xlab="",ylab="",xaxt="n",cex.axis=1.5)
   points(M,rep(0,mixnum))