> check.df<-read.csv("check190.txt",header=T,fill=TRUE,sep="")
> 
> # UAF version
> #H_D: U1
> 
> mixD<-DNAmixtureUAF(list(check.df), k=1,K=character(0),C=list(0.001),database=USCaucasian2,dyes=list(SGMplusDyes))
> startpar<-mixpar(rho=list(2),eta=list(100),xi=list(0.1),phi=list(c(U1=0.9)))
> 
> # unconstrained ML
> 
> mlD<-mixML(mixD,startpar,trace=F)
> 
> mlD$lik/log(10)
[1] -4.304199
> mlD$mle
     rho     eta                          xi   phi.U1
1   69.3   1.375   0.00000000000000000003876        1
> 
> #H_P:suspect
> mixP<-DNAmixtureUAF(list(check.df), k=1,K="suspect",C=list(0.001),database=USCaucasian2,dyes=list(SGMplusDyes))
> startpar<-mixpar(rho=list(2),eta=list(100),xi=list(0.1),phi=list(c(suspect=0.9)))
> 
> mlP<-mixML(mixP,startpar,trace=F)
> 
> mlP$lik/log(10)
[1] -2.679907
> mlP$mle
      rho      eta                          xi   phi.suspect
1   214.3   0.4435   0.00000000000000000009202             1
> 
> #WoE
> (mlP$lik/log(10))-(mlD$lik/log(10))
[1] 1.624292
> 
> # vanilla version
> #H_D: U1
> mixD<-DNAmixture(list(check.df), k=1,K=character(0),C=list(0.001),database=USCaucasian2,dyes=list(SGMplusDyes))
> startpar<-mixpar(rho=list(2),eta=list(100),xi=list(0.1),phi=list(c(U1=0.9)))
> 
> # unconstrained ML
> 
> mlD<-mixML(mixD,startpar,trace=F)
> 
> mlD$lik/log(10)
[1] -4.441308
> mlD$mle
      rho     eta                          xi   phi.U1
1   47.11   2.025   0.00000000000000000009193        1
> 
> #For H_p : suspect
> mixP<-DNAmixture(list(check.df), k=1,K="suspect",C=list(0.001),database=USCaucasian2,dyes=list(SGMplusDyes))
> startpar<-mixpar(rho=list(2),eta=list(100),xi=list(0.1),phi=list(c(suspect=0.9)))
> 
> mlP<-mixML(mixP,startpar,trace=F)
> 
> mlP$lik/log(10)
[1] -2.679907
> mlP$mle
      rho      eta                          xi   phi.suspect
1   214.3   0.4435   0.00000000000000000009202             1
> 
> #WoE
> (mlP$lik/log(10))-(mlD$lik/log(10))
[1] 1.761401