R/UpdateLambda.GL.R

Defines functions UpdateLambda.GL

UpdateLambda.GL <-
function(survObj, priorPara, ini){

	p		<- survObj$p
	K		<- priorPara$K
	tauSq	<- ini$tauSq 
	
	r		<- priorPara$r
	delta	<- priorPara$delta
	
	lambda1Sq	<- rgamma(1, shape = (p + K)/2 + r, rate = delta + sum(tauSq)/2)

	return(lambda1Sq)

	}

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.