R/is.clock.R

In bomeara/phybase: Basic functions for phylogenetic analysis

`is.clock` <-
function(nodematrix, nspecies,threshold)
{

   clock<-TRUE
   if(dim(nodematrix)[1]==-9)
  	clock<-FALSE
   else{
	for(i in (nspecies+1):(2*nspecies-1)){
		leftson<-nodematrix[i,2]
		rightson<-nodematrix[i,3]
		leftlength<-nodematrix[leftson,4]+node.height(leftson,nodematrix,nspecies)
		rightlength<-nodematrix[rightson,4]+node.height(rightson,nodematrix,nspecies)
		if(abs(leftlength-rightlength)>threshold){
			clock<-FALSE
			break
		}
	}
   }
		
   return(clock)
}