R/DDupLambda.1.R

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

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

  D<-nrow(B)

  M= nrow(Z.a)


  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.1^2 * M / 2 + sum(Y.ia)) * diag(D)  + f1s2To.ia) * M / 2

}