

Dear all: 

I am writing the following small function for a probit likelihood.
As indicated, in order to avoid p=1 or p=0, I defined some precisions.
I feel however, that there might be a better way to do this.
Any help is greatly appreciated.

##########################################################

##set limits to avoid px=0 or px=1
precision1   <- 0.999999
precision0   <- 0.000001

logpost <- function(par, data){
px    <- pnorm(b0 + b1x)
# to avoid px=1 or px=0
px[px >  precision1] <- precision1
px[px <  precision0] <- precision0
loga  <- sum( y*log(px)+(1-y)*log(1-px) )
loga
}

#########################################################



Best,
 

________________________________
Keramat Nourijelyani, PhD
Associate Professorof Biostatistics
Tehran University of Medical Sciences
http://tums.ac.ir/faculties/nourij

	[[alternative HTML version deleted]]

