# \donttest{
# Subset example data to one isotype and sample as a demo
data(ExampleDb, package="alakazam")
db <- subset(ExampleDb, c_call == "IGHA" & sample_id == "-1h")[1:80,]
# Create model using only silent mutations and ignore multiple mutations
model <- createTargetingModel(db, model="s", sequenceColumn="sequence_alignment",
germlineColumn="germline_alignment_d_mask",
vCallColumn="v_call", multipleMutation="ignore")
# View top 5 mutability estimates
head(sort(model@mutability, decreasing=TRUE), 5)
# View number of silent mutations used for estimating mutability
model@numMutS
# }
Run the code above in your browser using DataCamp Workspace