source("http://www.stat.washington.edu/hoff/Code/GBME/gbme.r")


###
Y<-as.matrix(read.table(
  "http://www.stat.washington.edu/hoff/Code/GBME/Class/Y.dat"))

n<-dim(Y)[1]

sex<<-c(rep(1,13),rep(0,14))               #1-13 are boys, 14-27 are girls
Xss<-array(outer(sex,sex)+outer(1-sex,1-sex),dim=c(n,n,1))#indicator of samesex 
###

gbme(Y,Xss,fam="binomial",N=matrix(1,27,27),k=2)

#for binomial responses, the "number of trials" per pair is 
#given in the matrix N. For binary data, this is just a matrix 
#of all 1s.