R/UpdateSigma.GL.R

Defines functions UpdateSigma.GL

UpdateSigma.GL <-
function(survObj, priorPara, ini){
	
	p			<- survObj$p
	be.normSq	<- ini$be.normSq
	tauSq		<- ini$tauSq 

	sh.sig     <- p/2
	rate.sig   <- 1/2*sum(be.normSq/tauSq)

	sigmaSq     <- rigamma(1, a = sh.sig, b = rate.sig)
	
	return(sigmaSq)
	
	}

Try the psbcGroup package in your browser

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

psbcGroup documentation built on June 24, 2021, 9:08 a.m.