R/DupLambda.0.R

In selenashuowang/jlsm: Joint Latent Space Model for Social Networks with Multivariate Attributes

DupLambda.0<-function( alpha.1, B,  Y.ia, Z.i,Z.a, Lambda.0, Lambda.1 ,p.lambda.0){

  D<-nrow(B)
  N<-nrow(Y.ia)



  rbZ.a=Z.a

  for(each in 2:nrow(Z.i)){rbZ.a=rbind(Z.a,rbZ.a)}

  rbZ.i=Z.i[rep(1:nrow(Z.i), times = rep(nrow(Z.a),nrow(Z.i))), ]



  dnj.ia <- rbZ.a - rbZ.i

  A.ia <-  1 + 1 / sqrt(det(B)) * exp(- alpha.1) * exp(apply(dnj.ia, 1, function(x) x %*% B %*% x))


  f1s2To.ia <- 2 * B %*% (t( dnj.ia  / A.ia) %*% dnj.ia) %*% B -   sum(1 / A.ia) * B


  solve((1 / p.lambda.0^2 * N / 2 +  sum(Y.ia)) * diag(D) +   f1s2To.ia) * N / 2



}