R/neighbor.R

In mdhglm: Multivariate Double Hierarchical Generalized Linear Models

neighbor <-
function() {
nbr<-matrix(0,56,56)

nbr[1,c(5,9,11,19)]<--1

nbr[2,c(7,10)]<--1

nbr[3,c(6,12)]<--1

nbr[4,c(18,20,28)]<--1

nbr[5,c(1,11,12,13,19)]<--1

nbr[6,c(3,8)]<--1

nbr[7,c(2,10,13,16,17)]<--1

nbr[8,c(6)]<--1

nbr[9,c(1,11,17,19,23,29)]<--1

nbr[10,c(2,7,16,22)]<--1

nbr[11,c(1,5,9,12)]<--1

nbr[12,c(3,5,11)]<--1

nbr[13,c(5,7,17,19)]<--1

nbr[14,c(31,32,35)]<--1

nbr[15,c(25,29,50)]<--1

nbr[16,c(7,10,17,21,22,29)]<--1

nbr[17,c(7,9,13,16,19,29)]<--1

nbr[18,c(4,20,28,33,55,56)]<--1

nbr[19,c(1,5,9,13,17)]<--1

nbr[20,c(4,18,55)]<--1

nbr[21,c(16,29,50)]<--1

nbr[22,c(10,16)]<--1

nbr[23,c(9,29,34,36,37,39)]<--1

nbr[24,c(27,30,31,44,47,48,55,56)]<--1

nbr[25,c(15,26,29)]<--1

nbr[26,c(25,29,42,43)]<--1

nbr[27,c(24,31,32,55)]<--1

nbr[28,c(4,18,33,45)]<--1

nbr[29,c(9,15,16,17,21,23,25,26,34,43,50)]<--1

nbr[30,c(24,38,42,44,45,56)]<--1

nbr[31,c(14,24,27,32,35,46,47)]<--1

nbr[32,c(14,27,31,35)]<--1

nbr[33,c(18,28,45,56)]<--1

nbr[34,c(23,29,39,49,42,43,51,52,54)]<--1

nbr[35,c(14,31,32,37,46)]<--1

nbr[36,c(23,37,39,41)]<--1

nbr[37,c(23,35,36,41,46)]<--1

nbr[38,c(30,42,44,49,51,54)]<--1

nbr[39,c(23,34,36,40,41)]<--1

nbr[40,c(34,39,41,49,52)]<--1

nbr[41,c(36,37,39,40,46,49,53)]<--1

nbr[42,c(26,30,34,38,43,51)]<--1

nbr[43,c(26,29,34,42)]<--1

nbr[44,c(24,30,38,48,49)]<--1

nbr[45,c(28,30,33,56)]<--1

nbr[46,c(31,35,37,41,47,53)]<--1

nbr[47,c(24,31,46,48,49,53)]<--1

nbr[48,c(24,44,47,49)]<--1

nbr[49,c(38,40,41,44,47,48,52,53,54)]<--1

nbr[50,c(15,21,29)]<--1

nbr[51,c(34,38,42,54)]<--1

nbr[52,c(34,40,49,54)]<--1

nbr[53,c(41,46,47,49)]<--1

nbr[54,c(34,38,49,51,52)]<--1

nbr[55,c(18,20,24,27,56)]<--1

nbr[56,c(18,24,30,33,45,55)]<--1

nbr<-nbr*-1

return (nbr)
}