comID<-unique(dat[,1])

output<-data.frame(K.obs=0,K.null.mean=0,
		K.null.lwr=0,
		K.null.upr=0,
		n.rand=0)
		
n.species<-NULL		

for (i in 1:length(comID)) {
	
	tmp<-dat[dat[,1]==comID[i],]
	out<-Krand(phy,tmp,n.rand=1000,sub=TRUE,fig=FALSE)
	n.species[i]<-nrow(tmp)
	output[i,]<-out
	}
	
output<-cbind(comID,n.species,output)