Missing Parental Data
Opened this issue · 0 comments
dcgerard commented
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)