Skip to content
Permalink
Browse files

small fix in fitMk

  • Loading branch information...
liamrevell committed Sep 27, 2017
1 parent 3f4ea6d commit 554a3abafd708325fd827f92e87afeca36ab18a5
Showing with 1 addition and 1 deletion.
  1. +1 −1 R/fitMk.R
@@ -84,7 +84,7 @@ fitMk<-function(tree,x,model="SYM",fixedQ=NULL,...){
el<-pw$edge.length[ii]
v<-vector(length=length(desc),mode="list")
for(j in 1:length(v))
v[[j]]<-matexpo(Q*el[j])%*%liks[desc[j],]
v[[j]]<-expm(Q*el[j])%*%liks[desc[j],]
vv<-if(anc==root) Reduce('*',v)[,1]*pi else Reduce('*',v)[,1]
comp[anc]<-sum(vv)
liks[anc,]<-vv/comp[anc]

0 comments on commit 554a3ab

Please sign in to comment.
You can’t perform that action at this time.