R/rm_sdt_evaluate_prior_derivative.R

In sirt: Supplementary Item Response Theory Models

## File Name: rm_sdt_evaluate_prior_derivative.R
## File Version: 0.08


rm_sdt_evaluate_prior_derivative <- function(partable, h)
{
    partable <- partable[ partable$est, ]
    m1 <- partable$prior_M
    sd1 <- partable$prior_SD
    y1 <- sirt_dnorm(partable$value+h, mean=m1, sd=sd1, log=TRUE)
    y2 <- sirt_dnorm(partable$value-h, mean=m1, sd=sd1, log=TRUE)
    res <- - ( y1 - y2 ) / (2*h)
    res <- ifelse(is.na(res),0,res)
    return(res)
}