Error in Winner's curse correction

[Edwa0506]

The R winner's curse correction uses the below function:

bias <- function(betaTrue, betaObs, se){
z <- betaTrue / se
num <- dnorm(z - c) - dnorm(-z - c)
den <- pnorm(z - c) + pnorm(-z - c)
return(betaObs - betaTrue + se * num / den)
}

The below line is meant to return the bias.
return(betaObs - betaTrue + se * num / den)
But there is a typo. It shoudl really be:
return(betaObs - betaTrue - se * num / den)

[Cloufield]

Hi,
thanks a lot for pointing out the error. Indeed, it should be betaObs - betaTrue - se * num / den for solving.
Just fixed the R codes.