Displaying 1 result from an estimated 1 matches for "wtmat".
Did you mean:
tmat
2012 Apr 26
0
Problem with constrOptim when hitting boundary
...NA,nBelowCutoff=nBelowCutoff,nZero=nZero,nMut=length(m)))
}
###Make 'weight matrix' of P(X=X[i]|X>0, P=p(k,s)) of size nStages x nCopies for each x --
###an array that is nMuts x nCopies of the prob of X given a certain allele
probX<-function(p){
dbinom(x=x,size=m,prob=p)
}
wtMat<-sapply(possAlleles,probX) #nMuts x nAlleles
if(is.null(dim(wtMat))){
if(length(x)==1) wtMat<-matrix(wtMat,nrow=1)
else stop("Programming error -- single row wtMat")
}
###Need to convert so in terms of pi[-0]
transMat<-rbind(rep(-1,K),diag(K))
Aid<-A%*%transMat
numG...