Missing Parental Data

[dcgerard]

I should allow parental data to be encoded as missing with NA rather than only accepting 0 counts.

Generate data

library(updog)
nsamp <- 191
genovec <- rgeno(n = nsamp, ploidy = 4, model = "f1", p1geno = 2, p2geno = 2)
sizevec <- rep(6, nsamp)
refvec <- rflexdog(sizevec = sizevec, 
                   geno = genovec, 
                   ploidy = 4, 
                   seq = 0.001, 
                   bias = 1, 
                   od = 0.001)

Errors if you give it NA

fout <- flexdog(refvec = refvec, 
                sizevec = sizevec, 
                ploidy = 4, 
                model = "f1",
                p1ref = NA, 
                p1size = NA)

OK if you have 0 ref and 0 size.

fout <- flexdog(refvec = refvec, 
                sizevec = sizevec, 
                ploidy = 4, 
                model = "f1",
                p1ref = 0, 
                p1size = 0)
plot(fout)