Forum

Nome Utente:
Password:
Riconoscimi automaticamente
 Tutti i Forum
 MolecularLab
 Bioinformatica
 [R] determinare se una sequenza è codificante
 Nuova Discussione  Nuovo Sondaggio Nuovo Sondaggio
 Rispondi Aggiungi ai Preferiti Aggiungi ai Preferiti
Cerca nelle discussioni
I seguenti utenti stanno leggendo questo Forum Qui c'è:

Aggiungi Tag Aggiungi i tag

Quanto è utile/interessante questa discussione:

Autore Discussione  

ki4rett4
Nuovo Arrivato



35 Messaggi

Inserito il - 02 febbraio 2011 : 13:25:32  Mostra Profilo  Visita l'Homepage di ki4rett4 Invia a ki4rett4 un Messaggio Privato  Rispondi Quotando
Buongiorno ragazzi per piacere aiutatemi su quest esercizio in R. C'ho provato, come sempre, ma non mi escono mai!!
L' esercizio è:
Scrivere una funzione, codificante(rna), che restituisce TRUE se all’interno della sequenza di rna il codone di start (AUG) precede quello di stop (UAA, UGA, o UAG), FALSE altrimenti.

Io l ho svolto in questo modo:

codificante<-function(rna){
  for(i in 1:(length(rna)-1)){
    if(rna[i]=="A"&rna[i+1]=="U"&rna[i+2]=="G"&(rna[length(rna)-1]=="A"|rna[length(rna)-1]=="G")&(rna[length(rna)-2]=="A"|rna[length(rna)-2]=="G")&rna[length(rna)-3]=="U"){
condizione<-TRUE
    }else{
      condizione<-FALSE
    }
  }
return(condizione)
}



Grazie a tutti per l' aiuto e buona giornata!!

chick80
Moderatore

DNA

Città: Edinburgh


11491 Messaggi

Inserito il - 02 febbraio 2011 : 15:18:52  Mostra Profilo  Visita l'Homepage di chick80 Invia a chick80 un Messaggio Privato  Rispondi Quotando
Dai un'occhiata all'help per grep e funzioni simili...

In particolare quella che serve a te è la funzione gregexpr.


Sei un nuovo arrivato?
Leggi il regolamento del forum e presentati qui

My photo portfolio (now on G+!)
Torna all'inizio della Pagina

ki4rett4
Nuovo Arrivato



35 Messaggi

Inserito il - 03 febbraio 2011 : 20:30:00  Mostra Profilo  Visita l'Homepage di ki4rett4 Invia a ki4rett4 un Messaggio Privato  Rispondi Quotando
Grazie chick80 ma come ho detto anche a globus il professore non ci ha fatto studiare tutte le funzioni, ragion per cui non posso utilizzarle. L' esercizio dovrei svolgerlo con un ciclo for ma come hai visto non mi esce correttamente!!
Torna all'inizio della Pagina

chick80
Moderatore

DNA

Città: Edinburgh


11491 Messaggi

Inserito il - 03 febbraio 2011 : 20:55:09  Mostra Profilo  Visita l'Homepage di chick80 Invia a chick80 un Messaggio Privato  Rispondi Quotando
L'esercizio domanda di controllare se il codone di start (AUG) precede quello di stop, tu controlli se c'è un codone di start e se l'ultimo codone è un codone di stop.

Se proprio vuoi usare un ciclo for allora passa ciascun elemento dell'array fino a che trovi un codone di start e metti la posizione del codone in una variabile.

Dopodichè continui da quella posizione fino a trovare un codone di stop.

Se lo trovi allora restituisci TRUE altrimenti FALSE.

PS: sappi tuttavia che in una situazione reale la soluzione di questo tipo di problemi non è decisamente quella che ho scritto di sopra, che è di difficile comprensione, lenta ed inefficiente. Due chiamate a gregexpr basterebbero per avere la soluzione. :)
R è un linguaggio molto potente, ma decisamente i cicli for sono da evitare il più possibile.

Sei un nuovo arrivato?
Leggi il regolamento del forum e presentati qui

My photo portfolio (now on G+!)
Torna all'inizio della Pagina

ki4rett4
Nuovo Arrivato



35 Messaggi

Inserito il - 04 febbraio 2011 : 13:12:53  Mostra Profilo  Visita l'Homepage di ki4rett4 Invia a ki4rett4 un Messaggio Privato  Rispondi Quotando
Non ho capito molto bene come si procede chick80....ma proverò a mettere in pratica quello che mi hai detto!! E' intuitivo che sei un esperto di programmazione, sia da come mi esponi la soluzione, sia da come riesci a capire al volo quale sia la soluzione più rapida ed efficiente...Ma io purtroppo è da poco che ho iniziato a programmare, causa esame, per cui non ne capisco molto ma ce la sto mettendo tutta per apprenderne qualcosa e andar al meglio all' esame!! Inoltre, come ti ho detto, non posso utilizzare tutte le funzioni in R, in quanto il prof. vuole che utilizziamo solo ciò di cui lui ci ha parlato...non è colpa mia!
Torna all'inizio della Pagina

chick80
Moderatore

DNA

Città: Edinburgh


11491 Messaggi

Inserito il - 04 febbraio 2011 : 17:28:56  Mostra Profilo  Visita l'Homepage di chick80 Invia a chick80 un Messaggio Privato  Rispondi Quotando
Certo, certo capisco!

Comunque sia prova a scrivere un'altra funzione e poi ne possiamo discutere insieme!

Ricorda, inizia a cercare il codone di start e, una volta trovato passi a cercare lo stop, non cercarli entrambi allo stesso tempo.

Sei un nuovo arrivato?
Leggi il regolamento del forum e presentati qui

My photo portfolio (now on G+!)
Torna all'inizio della Pagina
  Discussione  

Quanto è utile/interessante questa discussione:

 Nuova Discussione  Nuovo Sondaggio Nuovo Sondaggio
 Rispondi Aggiungi ai Preferiti Aggiungi ai Preferiti
Cerca nelle discussioni
Vai a:
MolecularLab.it © 2003-18 MolecularLab.it Torna all'inizio della Pagina