Question: Reduce saturation by deleting the 3rd position in each codon of DNA seqs. R ape package
0
gravatar for dago
3.3 years ago by
dago2.5k
Germany
dago2.5k wrote:

I am performing a phylogenentic analysis using specific marker genes.

I selected the aa sequences, aligned them, used as template for aligning the respective nucleotide sequences. 

I created a saturation plot with the "ape" package in R, and I can see that same seq are saturated.

In order to have a more robust results I would like, now, to delete the third position in each codon.

Is there a way/program to do this?

I guess this could be done with

seq

test3[seq(1, length(test3), by = 3)] # take only the first nucleotide in the codon

but I cannot find a way of taking only the first one and the second one.

I am also open to any suggestion for different way of doing this task

 

dna alignment R phylogeny • 1.6k views
ADD COMMENTlink modified 3.3 years ago by David W4.7k • written 3.3 years ago by dago2.5k
3
gravatar for David W
3.3 years ago by
David W4.7k
New Zealand
David W4.7k wrote:

In R you can use negative indices to exclude elements. 

 

Here's an exmaple using one of ape's datasets (not this is an alignment (matrix), not a single sequence (vector) as in your example

 

ADD COMMENTlink modified 3.3 years ago • written 3.3 years ago by David W4.7k

Ok, really good thanks!
 

ADD REPLYlink written 3.3 years ago by dago2.5k

I might not understand the command AlignmentToFilter[, -seq(1, ncol(AlignmentToFilter), by=3) ] .. but it is apparently skipping the first codon position in my hand while command AlignmentToFilter[, -seq(3, ncol(AlignmentToFilter), by=3) ] worked for me and omitted the 3rd codon position.

ADD REPLYlink written 9 months ago by al-ash30
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1630 users visited in the last hour