R/genYY_AR.R

In SAcmu/LLSM: Package to Fit Longitudinal Latent Space Model

#'@title Function to simulate longitudinal network data under LLSM-AR model
#'
#' @description
#' \code{genYY_AR} generate longitudinal networks under LLSM-AR model
#' @details
#' \code{genYY_AR} generate longitudinal networks under LLSM-AR model
#'
#' @export


genYY_AR <-
function(Phi,Tau,Beta,TT,dd,nn)
{
	ZZ = genZZ_AR(TT=TT,Tau=Tau,Phi=Phi,dd=dd,nn=nn)
	distMat = genDD(ZZ=ZZ)
	Prob = getProb(distMat = distMat, Beta=Beta)
	YY = list()
	for(tt in 1:TT ){
		YY[[tt]] =  array(0,dim=c(nn,nn))
        ll = 1
            for(cc in 1:nn){
                for(rr in cc:nn){
                    if(rr == cc){YY[[tt]][rr,cc] = 0}else{
                        YY[[tt]][rr,cc] = rbinom(1,1,Prob[[tt]][ll])
                        YY[[tt]][cc,rr] = YY[[tt]][rr,cc]
                        ll = ll + 1
                   }
               }
        }
	}
	return(list(Phi=Phi,ZZ=ZZ,distMat=distMat,Prob=Prob,YY=YY))
}