%*%beta_A.hat)^2))
#mse_PR
#clacualte g1, g2, and g3

####	gi
g1i=u1i*zi0	# defined in (24)
vTi=vech(zi0%*%t(zi0))

- x_mtx%*%beta_A.hat)^2))
#mse_PR
#clacualte g1, g2, and g3

)), 1/(A.BPE + Di)))
G2i = 2*(1-Bi)^2*(Wi - deli)
return(list(Bi = Bi, fi = fi, Wi = Wi, G2i = G2i))

= "ROW")
g2i <- runCompareComplex(goMatrix, intactBGM)
##Those GO complexes which are redundant:

(R2i)
R[2,2]=sum(R4i)
g1i=y*R1i

g1i=y*R1i
g2i=x*y*R1i+(ei*R1i)^2*(slope_CBLS*varX-cov_CBLS)
g[1,1]=sum(g1i)

(R1i)
R[1,2]=sum(R2i)
R[2,1]=sum(R2i)

g1i=y*R1i
g2i=x*y*R1i+(ei*R1i)^2*(slopes_all[j,1]*varX-CBLS.cov)
g[1,1]=sum(g1i)

## gene-index
g2i <- data.table(gene_id=core$genes$gene_id, idx=1:length(core$genes))
## cytoband-index

## gene-index
g2i <- data.table(gene_id=core$genes$gene_id, idx=1:length(core$genes))
## cytoband-index

R[1,2]=sum(R2i)
R[2,1]=sum(R2i)
R[2,2]=sum(R4i)

*R1i
R[1,1]=sum(R1i)
R[1,2]=sum(R2i)

R[1,2]=sum(R2i)
R[2,1]=sum(R2i)
R[2,2]=sum(R4i)

=vech(zi0%*%t(zi0))
g2i=u2i*vTi-sigma	# defined in (25)
gi=rbind(g1i,g2i)

=vech(zi0%*%t(zi0))
g2i=u2i*vTi-sigma	# defined in (25)
gi=rbind(g1i,g2i)

##     ## //         = g2i   if j is in new scan
##     ## //    and the distribution of g2i is specified as;
##     ## //    g2i  = g1i

####	gi
g1i=u1i*zi0	# defined in (24)
vTi=vech(zi0%*%t(zi0))