[R] Need some help in R : value more than equals to a row.

Gabor Grothendieck ggrothendieck at gmail.com
Sun Jun 7 19:58:28 CEST 2009


Try this:

t(apply(m, 1, function(x) colMeans(outer(x, x, ">="))))


On Sun, Jun 7, 2009 at 10:09 AM, suparna
mitra<mitra at informatik.uni-tuebingen.de> wrote:
> Hallo,
>  I was trying some code, but couldn't make one step of the code properly.
> Can anybody  please  help me?
>
> I have one matrix like this
>> values
>          [,1]       [,2]       [,3]      [,4]      [,5]
> [1,] 0.7777778 0.36111111 0.22222222 0.1388889 0.0000000
> [2,] 1.0000000 0.00000000 0.53846154 0.0000000 0.5384615
> [3,] 0.5200000 0.48000000 0.64000000 0.0000000 0.8800000
> [4,] 0.8928571 1.00000000 0.00000000 0.8928571 0.1071429
>
>
> And I want to get some matrix like:
>> values.new
>      [,1]  [,2]  [,3]  [,4]  [,5]
> [1,] 0.2  0.4  0.6  0.8  1.0
> [2,] 0.2  1.0  0.6  1.0  0.6
> [3,] 0.6  0.8  0.4  1.0  0.2
> [4,] 0.6  0.2  1.0  0.6  0.8
>
>
> This table should be computed by taking proportion of values in the row that
> are larger or equals to the value being considered with the total no of
> objects in the row.
>
> for example
>> values[1,2]
> [1] 0.3611111
> and
>> values[1,]
> [1] 0.7777778 0.3611111 0.2222222 0.1388889 0.0000000 0.8611111
> So there are two numbers more than equals to values[1,2]
> So
>> values.new[1,2]= 2/length(values[1,])
>
>
>
> With best regard,
> Suparna
> ----------------------------
> Ms Suparna Mitra
> Eberhard-Karls-Universität Tübingen
> Wilhelm-Schickard-Institut
> Algorithmen der Bioinformatik
> Sand 14, 72076 Tuebingen
> Germany
> Tel. ++49-7071-29 70453 (0)
> Fax ++49-7071-29 5148 (0)
> Phone: ++49-176-20361469 (M)
>            ++49-7071-1477169 (R)
> Alternative e-mail: mitra at informatik.uni-tuebingen.de
>
>        [[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