nullestimation function

nullestimation calculates the estimation under the null hypothesis

nullestimation(A, d)

Arguments

• A: List of symmetric A matrices
• d: Dimension of the latent space

Returns

Returns a list of the following

Examples

#simulate data
U <- matrix(0, nrow=20, ncol=3)
U[,1] <- 1/sqrt(20)
U[,2] <- rep(c(1,-1), 10)/sqrt(20)
U[,3] <- rep(c(1,1,-1,-1), 5)/sqrt(20)

L<-list(diag(c(11,6,2)),diag(c(15,4,1)))
A <- list()
for(i in 1:2){
  P <- U%*%L[[i]]%*%t(U)
  A[[i]] <-apply(P,c(1,2),function(x){rbinom(1,1,x)})
  A[[i]][lower.tri(A[[i]])]<-t(A[[i]])[lower.tri(A[[i]])]
}

#fit model
nullestimation(A,3)

See Also

multiRDPG

Author(s)

Agnes Martine Nielsen (agni@dtu.dk)