[R] which question

Richard Pearson richard.pearson at postgrad.manchester.ac.uk
Fri Jun 6 16:41:49 CEST 2008


I didn't get any attached data, but my suspicion here is that you have somehow got RefSeq IDs in column 8 of human, as well as the gene symbols. Did you read this data in from a text file?

Eleni Christodoulou wrote:
> An example is:
> 
> symbol=human[which(human[,3] %in% genes.sam.names),8]
> 
> The data* human* and *genes.sam.names* are attached. The result of the above
> command is:
>> symbol
>  [1] CCL18                  MARCO                  SYT13
>  [4] FOXC1                  CDH3
>  [7] CA12                   CELSR1                 NM_018440
> [10] MICROTUBULE-ASSOCIATED NM_015529              ESR1
> [13] PHGDH                  GABRP                  LGMN
> [16] MMP9                   BMP7                   KLF5
> [19] RIPK2                  GATA3                  NM_032023
> [22] TRIM2                  CCND1                  MMP12
> [25] LDHB                   AF493978               SOD2
> [28] SOD2                   SOD2                   NME5
> [31] STC2                   RBP1                   ROPN1
> [34] RDH10                  KRTHB1                 SLPI
> [37] BBOX1                  FOXA1                  NM_005669
> [40] MCCC2                  CHI3L1                 GSTM3
> [43] LPIN1                  DSC2                   FADS2
> [46] ELF5                   CYP1B1                 LMO4
> [49] AL035297               NM_152398              AB018342
> [52] PIK3R1                 NFKBIE                 MLZE
> [55] NFIB                   NM_052997              NM_006023
> [58] CPB1                   CXCL13                 CBR3
> [61] NM_017527              FABP7                  DACH
> [64] IFI27                  ACOX2                  CXCL11
> [67] UGP2                   CLDN4                  M12740
> [70] IGKC                   IGKC                   CLECSF12
> [73] AY069977               HOXB2                  SOX11
> [76]                        NM_017422              TLR2
> [79] CKS1B                  BC017946               APOBEC3B
> [82]                        HLA-DRB1               HLA-DQB1
> [85]                        CCL13                  C4orf7
> [88]                        NM_173552
> 21345 Levels:  (2 (32 (55.11 (AIB-1) (ALU (CAK1) (CAP4) (CASPASE ... ZYX
> 
> As you can see, apart from gene symbols, which is the required thing, RefSeq
> ID sare also retrieved...
> 
> Thanks a lot,
> Eleni
> 
> 
> 
> 
> 
> 
> On Fri, Jun 6, 2008 at 1:23 PM, Dieter Menne <dieter.menne at menne-biomed.de>
> wrote:
> 
>> Eleni Christodoulou <elenichri <at> gmail.com> writes:
>>
>>> I was trying to select a column of a data frame using the *which*
>> command. I
>>> was actually selecting the rows of the data frame using *which, *and then
>>> displayed a certain column of it. The command that I was using is:
>>> sequence=*mydata*[*which*(human[,3] %in% genes.sam.names),*9*]
>> ....
>> Please provide a running example. The *mydata* are difficult to read.
>>
>>
>> Dieter
>>
>> ______________________________________________
>> R-help at r-project.org mailing list
>> https://stat.ethz.ch/mailman/listinfo/r-help
>> PLEASE do read the posting guide
>> http://www.R-project.org/posting-guide.html
>> and provide commented, minimal, self-contained, reproducible code.
>>
> 
> 	[[alternative HTML version deleted]]
> 
> ______________________________________________
> R-help at r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
> and provide commented, minimal, self-contained, reproducible code.
> 

-- 
Richard D. Pearson             richard.pearson at postgrad.manchester.ac.uk
School of Computer Science,    http://www.cs.man.ac.uk/~pearsonr
University of Manchester,      Tel: +44 161 275 6178
Oxford Road,                   Mob: +44 7971 221181
Manchester M13 9PL, UK.        Fax: +44 161 275 6204



More information about the R-help mailing list