Ergänzen Sie eine DNA-Sequenz

7

Angenommen, ich habe eine DNA-Sequenz. Ich möchte die Ergänzung davon bekommen. Ich habe den folgenden Code verwendet, aber ich verstehe es nicht. Was mache ich falsch?

%Vor%     
Anurag Mishra 04.12.2013, 09:43
quelle

5 Antworten

13

Verwenden Sie chartr , das für diesen Zweck erstellt wurde:

%Vor%

Geben Sie ihm zwei gleich lange Zeichenfolgen und Ihre Zeichenfolge. Auch über das Argument zur Übersetzung vektorisiert:

%Vor%

Hinweis: Ich mache den Ersatz für die String-Repräsentation der DNA und nicht für den Vektor. Um den Vektor zu konvertieren würde ich eine Lookup-Karte als benannte Vektor erstellen und indizieren, dass:

%Vor%     
Spacedman 04.12.2013 10:54
quelle
9

Das Bioconductor Paket Biostrings hat viele nützliche Funktionen für diese Art von Operation. Einmal installieren:

%Vor%

verwenden Sie dann

%Vor%     
Martin Morgan 04.12.2013 17:42
quelle
5
%Vor%

Wenn Sie die komplementären Namen nicht möchten, können Sie sie immer mit unname entfernen.

%Vor%     
42- 04.12.2013 09:58
quelle
5

Es gibt auch ein Paket seqinr

%Vor%

Biostrings hat ein viel kleineres Speicherprofil, aber Seqinr ist auch schön, weil Sie den Fall der Basen (einschließlich gemischt) wählen und sie zu allem ändern können, das Sie wollen, zum Beispiel, wenn Sie eine Mischung aus T und U in der wollen gleiche Reihenfolge. Biostrings zwingt Sie, entweder T oder U zu haben.

    
JeremyS 14.01.2014 04:10
quelle
4

Als Ergänzung können Sie in Groß- und Kleinschreibung chartr() :

verwenden %Vor%

Um die Nukleotidsequenz um einen Schritt weiter und umgekehrt zu ergänzen, können Sie die folgende Funktion verwenden:

%Vor%     
Megatron 14.06.2016 14:31
quelle