EBSeq/R/Likefun.R

   1 Likefun <-
   2 function(ParamPool, InputPool)
   3 {
   4
   5 NoneZeroLength=InputPool[[5]]
   6 AlphaIn=ParamPool[1]
   7 BetaIn=ParamPool[2:(1+NoneZeroLength)]
   8 PIn=ParamPool[2+NoneZeroLength]
   9 ZIn=InputPool[[4]]
  10 Input=InputPool[[3]]
  11 Input1=matrix(InputPool[[1]],nrow=nrow(Input))
  12 Input2=matrix(InputPool[[2]],nrow=nrow(Input))
  13 RIn=InputPool[[6]]
  14 RInSP1=matrix(InputPool[[7]],nrow=nrow(Input))
  15 RInSP2=matrix(InputPool[[8]],nrow=nrow(Input))
  16 NumIn=InputPool[[9]]
  17 ##Function here
  18 #LikelihoodFunction<- function(NoneZeroLength){
  19         F0=f0(Input, AlphaIn, BetaIn, RIn, NumIn, log=T)
  20         F1=f1(Input1, Input2, AlphaIn, BetaIn, RInSP1,RInSP2, NumIn, log=T)
  21                 F0[F0==Inf]=min(!is.na(F0[F0!=Inf]))
  22                 F1[F1==Inf]=min(!is.na(F1[F1!=Inf]))
  23
  24         -sum((1-ZIn)*F0+ (1-ZIn)* log(1-PIn) + ZIn*F1 + ZIn*log(PIn))
  25 }
  26