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

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.