I’m trying to follow the 454 SOP protocol with my data and run unifrac.weighted with a group=myfile.design. But I get back an error that the sequences it’s looking for aren’t in my group file.
I’ve tried running it a few ways.
unifrac.weighted(tree=LaSelva.final.pick.phylip.tre, group=nest.design, subsample=900, processors=2)
unifrac.weighted(tree=LaSelva.final.pick.phylip.tre, name=LaSelva.final.pick.names, group=nest.design, subsample=900, processors=2)
etc
Nothing I’ve done has worked. Am I missing something in how I should be generating my .tre file to use with a design file? Is there another file I should be adding into the command?
Thanks!