[R] Problem with rowMeans()
Erik Iverson
iverson at biostat.wisc.edu
Fri Jun 13 01:34:51 CEST 2008
ss wrote:
> Hi Wacek,
>
> Yes, data is data frame not a matrix.
>
>> is.numeric(data[3])
> [1] FALSE
what is class(data[3])
>
> But I looked at the column 3 and it looks okay though. There are few NAs and
> I did find
> anything strange.
>
> Any suggestions?
>
> Thanks,
> Allen
>
>
>
> On Thu, Jun 12, 2008 at 7:01 PM, Wacek Kusnierczyk <
> Waclaw.Marcin.Kusnierczyk at idi.ntnu.no> wrote:
>
>> ss wrote:
>>> Thank you very much, Wacek! It works very well.
>>> But there is a minor problem. I did the following:
>>>
>>>> data <-
>>> read.table('E-TABM-1-processed-data-1342561271_log2_with_symbols.txt',
>>> +row.names = NULL ,header=TRUE, fill=TRUE)
>> looks like you have a data frame, not a matrix
>>
>>
>>>> dim(data)
>>> [1] 23963 85
>>>> data[1:4,1:4]
>>> Probe_ID Gene_Symbol M16012391010920 M16012391010525
>>> 1 A_23_P105862 13CDNA73 -1.6 0.16
>>> 2 A_23_P76435 15E1.2 0.18 0.59
>>> 3 A_24_P402115 15E1.2 1.63 -0.62
>>> 4 A_32_P227764 15E1.2 -0.76 -0.42
>>>> data1<-data[sapply(data, is.numeric)]
>>>> dim(data1)
>>> [1] 23963 82
>>>> data1[1:4,1:4]
>>> M16012391010525 M16012391010843 M16012391010531 M16012391010921
>>> 1 0.16 -0.23 -1.40 0.90
>>> 2 0.59 0.28 -0.30 0.08
>>> 3 -0.62 -0.62 -0.22 -0.18
>>> 4 -0.42 0.01 0.28 -0.79
>>>
>>> You will notice that, after using 'data[sapply(data, is.numeric)]' and
>>> getting
>>> data1, the first sample in data, called 'M16012391010920', was missed
>>> in data1.
>>>
>>> Any further suggestions?
>>>
>> surely there must be an entry in column 3 that makes it non-numeric.
>> what does is.numeric(data[3]) say? (NAs should not make a column
>> non-numeric, unless there are only NAs there, which is not the case
>> here.) check your data for non-numeric entries in column 3, there can
>> be a typo.
>>
>> vQ
>>
>
> [[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.
More information about the R-help
mailing list