Bonjour,
J'ai un fichier texte contenant une liste de rsIDs pour de nombreux snp. Je veux faire une même demande de recherche de nom de gène et de transcrit sur chacune des ligne de ce fichier. Ainsi, je souhaiterai lire mon fichier ligne par ligne et indiquer la ligne en cours (boucle ?) dans ma requête. Comment lire mon fichier ligne par ligne et mettre dans "value" (de la fonction getBM) chaque ligne de mon fichier (ici à la place de Data_vec il s'agir d'avoir chaque ligne du fichier).
Voici le code que j'ai essayé :
mais j'ai toujours ce message :
Code : Sélectionner tout - Visualiser dans une fenêtre à part
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15 Library(biomaRt) Data <- readLines ("rsID.txt") longueur(Données) #Data_vec <- as.vector(Data) Data_vec<-c("Data") require(biomaRt) ensembl <- useMart("ENSEMBL_MART_SNP", dataset = "hsapiens_snp") getBM(attributes=c( "refsnp_id", "ensembl_gene_stable_id", "ensembl_transcript_stable_id"), filters="snp_filter", values="Data_vec", mart=ensembl, uniqueRow=TRUE)
Comment puis-je faire pour lie le fichier ligne par ligne et indiquer à value chaque ligne ligne après l'autre ?
Code : Sélectionner tout - Visualiser dans une fenêtre à part
1
2
3 1] refsnp_id ensembl_gene_stable_id ensembl_stable_id ensembl_transcript_stable_id <0 lignes> (ou 0-longueur de ligne.noms)
Partager