Displaying 1 result from an estimated 1 matches for "union_ch".
2010 Jul 07
2
how to define method for "+" function in a new class
...# Check both objects are correct
checkMolecule(x)
checkMolecule(y)
# Extract chains information
ch1 <- getMoleculeChains(x)
ch2 <- getMoleculeChains(y)
union_ch <- unique(c(ch1,ch2))
not_used <- .ALPHABET
for (i in seq(along =3D union_ch))
{
idx <- which(not_used !=3D union_ch[i])
if (length(idx) !=3D 0) not_used <- not_used[idx]...