Take the following DNAStringSet foo as follows

x <- DNAStringSet('AGAT')
y <- DNAStringSet('GGGTTCC')
z <- DNAStringSet('GGGCC')

foo <- c(x,y,z)
names(foo) <- c('x','y','z')

I'm interested in creating a new DNAStringSet of just x and y.  How do I go
about doing this?

At first, I thought you could do something like

c(foo['x'], foo['y'])

but this generates an error message, "cannot subset a DNAStringSet by
names."

Thanks,
Andrew

	[[alternative HTML version deleted]]

