R/UpdateLambda2.EN.R

Defines functions UpdateLambda2.EN

UpdateLambda2.EN <-
function(survObj, priorPara, ini){

	p		<- survObj$p
	be.ini	<- ini$beta.ini
	sigmaSq	<- ini$sigmaSq 
	
	r2		<- priorPara$r2
	delta2	<- priorPara$delta2
	
	lambda2	<- rgamma(1, shape = r2 + p/2, rate = delta2 + sum(be.ini^2)/(2*sigmaSq))

	return(lambda2)

	}

Try the psbcGroup package in your browser

Any scripts or data that you put into this service are public.

psbcGroup documentation built on May 29, 2024, 4:02 a.m.