Question: want the piece of the matrix I wanted
0
gravatar for yueli7
15 months ago by
yueli70
China
yueli70 wrote:

Hello,

 

I only want a piece of the matrix, not the whole.

 

I only want the fifth column to the end and the first row to the fourth.

 

That means I only want

 

P

V1 V2 V3

V1 0.3506 0.2943 0.6931

V2 0.7733 0.0383 0.8874

V3 0.4476 0.2438 0.0652

V4 0.1597 0.6114 0.6958

 

 

Thanks in advance for great help!

 

> library("Hmisc")

Loading required package: lattice

Loading required package: survival

Loading required package: Formula

Loading required package: ggplot2

 

Attaching package: ‘Hmisc’

 

The following objects are masked from ‘package:base’:

 

format.pval, units

 

> x<-as.matrix(read.table("tmp01"))

> x

V1 V2 V3 V4

[1,] 1 2 5 4

[2,] 4 4 7 5

[3,] 3 3 2 5

[4,] 3 5 5 3

[5,] 2 4 4 1

> y<-as.matrix(read.table("tmp02"))

> y

V1 V2 V3

[1,] 7 8 1

[2,] 10 1 2

[3,] 7 8 11

[4,] 10 1 3

[5,] 1 3 4

> P<-matrix(0,ncol(x),ncol(y))

> P

[,1] [,2] [,3]

[1,] 0 0 0

[2,] 0 0 0

[3,] 0 0 0

[4,] 0 0 0

> P<-rcorr(x,y)

> P

V1 V2 V3 V4 V1 V2 V3

V1 1.00 0.62 0.27 0.42 0.54 -0.59 0.24

V2 0.62 1.00 0.27 -0.37 0.18 -0.90 -0.09

V3 0.27 0.27 1.00 0.10 0.45 -0.64 -0.85

V4 0.42 -0.37 0.10 1.00 0.73 0.31 0.24

V1 0.54 0.18 0.45 0.73 1.00 -0.23 -0.15

V2 -0.59 -0.90 -0.64 0.31 -0.23 1.00 0.44

V3 0.24 -0.09 -0.85 0.24 -0.15 0.44 1.00

 

n= 5

 

 

P

V1 V2 V3 V4 V1 V2 V3

V1 0.2692 0.6659 0.4822 0.3506 0.2943 0.6931

V2 0.2692 0.6659 0.5436 0.7733 0.0383 0.8874

V3 0.6659 0.6659 0.8745 0.4476 0.2438 0.0652

V4 0.4822 0.5436 0.8745 0.1597 0.6114 0.6958

V1 0.3506 0.7733 0.4476 0.1597 0.7109 0.8040

V2 0.2943 0.0383 0.2438 0.6114 0.7109 0.4595

V3 0.6931 0.8874 0.0652 0.6958 0.8040 0.4595

matrix • 200 views
ADD COMMENTlink modified 15 months ago by Axel Klenk940 • written 15 months ago by yueli70
Answer: want the piece of the matrix I wanted
1
gravatar for Axel Klenk
15 months ago by
Axel Klenk940
Switzerland
Axel Klenk940 wrote:

Dear yueli7,

try

P[1:4, 5:7]

or

P[1:4, 5:ncol(P)]

is that what you want? If yes, it is very basic R subsetting and not related to Bioconductor.

It should be asked on StackOverflow or r-help and you may benefit from reading some introduction to R.

Hope this helps.

ADD COMMENTlink written 15 months ago by Axel Klenk940

Avel Klenk,

 

Thank you so much for your great help!

 

 

 

> library("Hmisc")
Loading required package: lattice
Loading required package: survival
Loading required package: Formula
Loading required package: ggplot2

Attaching package: ‘Hmisc’

The following objects are masked from ‘package:base’:

    format.pval, units

> x <- as.matrix(read.table("tmp01"))
> x
     V1 V2 V3 V4
[1,]  1  2  5  4
[2,]  4  4  7  5
[3,]  3  3  2  5
[4,]  3  5  5  3
[5,]  2  4  4  1
> y <- as.matrix(read.table("tmp02"))
> y
     V1 V2 V3
[1,]  7  8  1
[2,] 10  1  2
[3,]  7  8 11
[4,] 10  1  3
[5,]  1  3  4
> P <- matrix(0, ncol(x), ncol(y))
> P
     [,1] [,2] [,3]
[1,]    0    0    0
[2,]    0    0    0
[3,]    0    0    0
[4,]    0    0    0
> for (i in 1:ncol(x)) {
+  for(j in 1:ncol(y)) {
+ P[i,j] <- rcorr(x[,i],y[,j])$P[1,2]}
+ }
> P
          [,1]       [,2]       [,3]
[1,] 0.3506314 0.29431174 0.69307702
[2,] 0.7732771 0.03833574 0.88741510
[3,] 0.4476207 0.24378172 0.06520707
[4,] 0.1597474 0.61142749 0.69575215
 

 

ADD REPLYlink written 15 months ago by yueli70
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 16.09
Traffic: 159 users visited in the last hour