Alroy_Forbes<-function(occs)
{
	output<-matrix(nrow=ncol(occs),ncol=ncol(occs))
	rownames(output)<-colnames(occs)
	colnames(output)<-colnames(occs)
	for(i in 1:nrow(output))
	{
		for (j in 1:ncol(output))
		{
			species.A<-rownames(occs)[which(occs[,rownames(output)[i]]!=0)]
			species.B<-rownames(occs)[which(occs[,colnames(output)[j]]!=0)]

			n<-length(unique(c(species.A,species.B)))

			a<-length(intersect(species.A,species.B))
			b<-length(which(species.A%in%species.B==F))
			c<-length(which(species.B%in%species.A==F))
			
			output[i,j]<-1-((a*(n+sqrt(n)))/((a*(n+sqrt(n)))+(3/2*b*c)))
			
		}
	}
	output<-as.dist(output)
	return(output)
}