Rather than deal with the clustered OTUs for this post, I just uploaded all 1120 sequences, and carried out the analysis. The first time through (today) I forgot to include the outgroup! So that gives us a chance to see how much difference it makes.
temp
dental_1120.fna
dental_1120_rdp.fna
It looks like this:
![](https://blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEiJYiZ1hl08Ifrm4GXvowlA-NINcYxnLGD347UXyt-f33hqpjUJ_nIDY0hac5qKA9XQgkNDke0U5QGW7kPzky-AOF1W4NhytdocthAwApIvEMmYRWSCTHOKPdGfZTihRjidJpzs3IGbn_kQ/s320/Screen+shot+2011-03-10+at+12.31.40+PM.png)
pca_web.txt
Here it is:
![](https://blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEibvsH45WCKushnBHJoF3IIU_2rb3tUvUJ20O5Et5cFXoEljQ7ZrpHPJ3P08nX6TOlbcajtPeN6lV53xri1XPmGt8p52W2EK_W03j053IGyTfLZ8FAE9a53GclpQxian-4RatgKDS60KMIk/s320/pca.web.png)
From the UniFrac FAQ:
My tree was not rooted, but I was able to upload my file and perform an analysis. Are the results valid?
There is no way to tell based on a Newick string alone whether a tree is rooted or not. If an unrooted tree is input, UniFrac will usually assign an arbitrary root and allow you to perform the analysis on that tree. How the tree is rooted can affect the results of both UniFrac tests and the P test. You should redo the analysis with a tree that is rooted with an appropriate outgroup.
It turns out to be easy enough..go back to RDP and browse to find SL7 and then add it to the sequence cart. Repeat the download to
dental_1120+_rdp.fna
. Load the last 5 sequences into clustalx.app and look to make sure that SL7 is really properly aligned.![](https://blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEi7tLQ9Occ7PbV8naSi52vKbZW2-qrxynLg9HQbCNOGZIGn0dkds6qTDqyrHlpJmcqPCSmWZhU1ZAKkcbsue6GYq4L1tQQb9NdfpdFDQkqAmwidgV2VePBZ2kiQHXXT7AKcEr3y410m0t6X/s320/Screen+shot+2011-03-10+at+12.20.43+PM.png)
![](https://blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEjRF7QgpCOzgKU8fcrTi-aZq5vdMh-HxxpgNilv2KIoyvwaUXYqPlxk9VTDwGf9sDuKoynyRrkdxmivkkTWCPKQKlpnF29ncSs5HUAX5qSdyN0p2VcZKVDzbTFp-kJIvqo30CVPYuTOzaMJ/s320/Screen+shot+2011-03-10+at+12.30.24+PM.png)
Repeat the PCA. You can look at the data in a spreadsheet app:
![](https://blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEhJml7rcFZhNXxCOgum2uzqW1uBndmyi9vq1V2gaExPO_IO_n2SUyND_fAQnSA4gS9KtTrp7XUe2v8gdEy9ZA1VAcg5t2pYBE-Jr7tSjgFV3BhKoC8fRenU26wbHvL5AZMGu8r4WW2abgob/s320/Screen+shot+2011-03-10+at+12.11.43+PM.png)
Now I plot it in matplotlib. The first image is what I plotted today for the rooted tree. The second is from the paper. Looks pretty good to me. Also, note some minor differences from the previous graphic where the tree we used was unrooted (and UniFrac rooted it for us however it does when it's not properly rooted).
![](https://blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEhuJtYtgQ4za96Ub1NLOe8bL4_0e6gxpJRdKLNXKJoYfhlTHA4g8IM3gtI_qaau35lSaTPAdtpOAX58CaULGW0rCC6vdNIEtf6W48ycswe2EF7TxpWkrI363-g5kJjLQtA6E5ftZLtPAGLT/s320/pca_web2.png)
![](https://blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEhG__xZ_mwp2TdJHV44Sgn1sP2N3tcrjjE7M44xk0ijapXTK-iHb4RuNKibJfGOeM8rUdCF0gU-mEd8MJ2SrCv-u46MxjY-Kx1n2oBsZKLqtIYax06sG60GAYk19nULH3-b594pnP00FQ5V/s320/Screen+shot+2011-03-10+at+12.09.36+PM.png)
plot_web.py