# NOT RUN {
# Generate a toy z matrix
z <- abs(matrix(rnorm(50), ncol=2))
z <- z/rowSums(z)
# Return the median DBS width
dbs(z, summ="median")$msw
# For real sequence data
data(mvad)
# }
# NOT RUN {
mod <- MEDseq_fit(seqdef(mvad[,17:86]), G=11, modtype="UUN", weights=mvad$weight)
dbs(mod$z, weights=mvad$weight)
# }
Run the code above in your browser using DataLab