# samples 1-10 only give mean final numbers
D1 <- david[[1]]; D1
# estimate mutation probability
mutestim(D1$mc,mfn=D1$mfn)
# sample 11 has individual final numbers
D11 <- david[[11]]; D11
# with mean final number
mutestim(D11$mc,mfn=mean(D11$fn))
# with true final numbers
mutestim(D11$mc,D11$fn)
Run the code above in your browser using DataLab