lzip=function(parm,xsam){
p=parm[1]
lam=parm[2]
pdf=(1-p)*dpois(xsam,lam)+p*(xsam==0)
lpdf=-log(pdf)
lzip=sum(lpdf)
lzip
}
gzip=function(parm,xsam){
p=parm[1]
lam=parm[2]
pdf=(1-p)*dpois(xsam,lam)+p*(xsam==0)
gp=((xsam==0)-dpois(xsam,lam))/pdf
gp=-sum(gp)
gl=(dpois(xsam,lam)*(1-p)*(xsam/lam-1.0))/pdf
gl=-sum(gl)
gzip=c(gp,gl)
gzip
}