phoenix_tracker: Phoenix Tracker

Y<-event_matrix
diag(Y)<-0
lab<-actors
setwd("~/Documents/git/phoenix_tracker/data/gbme")
gbme(Y=Y, k=2, NS=500,odens=1, ofilename="OUT",fam="poisson", directed=F)

z<-as.matrix(read.table(file="Z", header=F, sep=" "))
#z <- z[-c(1:(length(actors)*200)),]
zz<-array(z,dim=c(dim(Y)[1],500,2))
zzz.mean<-apply(zz,c(1,3),mean)

pdf(height=5,width=5,file="latentnet.pdf")
plot(zz[,,1],zz[,,2],type="n",axes=FALSE,ylab="Latent Dimension 2", xlab="Latent Dimension 1",ylim=c(min(zz[,,2]),max(zz[,,2])),xlim=c(min(zz[,,1]),max(zz[,,1])))
axis(2,las= HORIZONTAL<-1)
axis(1)
points(zz[1,,1],zz[1,,2],col= "red")
points(zz[2,,1],zz[2,,2],col= "green")
points(zz[3,,1],zz[3,,2],col= "blue")
points(zz[4,,1],zz[4,,2],col= "violet")
points(zz[5,,1],zz[5,,2],col= "yellow")
points(zz[6,,1],zz[6,,2],col= "pink")
points(zz[7,,1],zz[7,,2],col= "orange")
dev.off()

pdf(height=5,width=5,file="latentnet.pdf")
plot(zzz.mean,type="n",axes=FALSE,ylab="Latent Dimension 2", xlab="Latent Dimension 1",ylim=c(-0.75,0.75),xlim=c(-0.75,0.75))
axis(2,las= HORIZONTAL<-1)
axis(1)
#points(zzz.mean,col="black")
text(zzz.mean,label=lab,col="black")