[R] Expanding a VCV Matrix
Maithula Chandrashekhar
m.chandrashekhar1977 at gmail.com
Thu Apr 21 17:08:43 CEST 2011
Thanks Uwe for your insight. However this could not solve my problem.
Actually, I meant to define 'Str' in this way:
> mat <- matrix(c(1,2,0,2,5,0.5,0,0.5,3), 3, 3)
> colnames(mat) <- rownames(mat) <- paste("variable", 1:3, sep="")
> mat
variable1 variable2 variable3
variable1 1 2.0 0.0
variable2 2 5.0 0.5
variable3 0 0.5 3.0
>
> Str<- c(paste("variable", 1:4, sep=""), "variable2"); Str
[1] "variable1" "variable2" "variable3" "variable4" "variable2"
Therefore 'variable2' should appear twice in my modified VCV matrix,
which was not there in your suggestion. Is there any better any on how
can I do that? or it is really imposable to achieve?
Thanks,
2011/4/21 Uwe Ligges <ligges at statistik.tu-dortmund.de>:
>
>
> On 20.04.2011 19:26, Maithula Chandrashekhar wrote:
>>
>> Dear all, I have special task to expand a given VCV matrix, however
>> could not accomplice yet. Let say I have following VCV matrix
>>
>>> mat<- matrix(c(1,2,0,2,5,0.5,0,0.5,3), 3, 3)
>>> colnames(mat)<- rownames(mat)<- paste("variable", 1:3)
>>> mat
>>
>> variable 1 variable 2 variable 3
>> variable 1 1 2.0 0.0
>> variable 2 2 5.0 0.5
>> variable 3 0 0.5 3.0
>>
>> Now, say I have a general string vector like this:
>>
>>> Str<- c(paste("variable", 1:4), "variable2")
>>> Str
>>
>> [1] "variable 1" "variable 2" "variable 3" "variable 4" "variable 2"
>
>
>
> If Str is like the printed one above rather than the code example, you can
> do:
>
> mat[Str[Str %in% rownames(mat)], Str[Str %in% colnames(mat)]]
>
> Uwe Ligges
>
>
>
>
>> Now according to this string, I want my previous VCV matrix also
>> expands. Therefore, as "variable 4" is not there in VCV matrix, so it
>> will be ignored. Therefore final VCV matrix will be of order 4, and
>> will look like:
>>
>> variable 1 variable 2 variable 3 variable 2
>> variable 1 1 2 0 2
>> variable 2 2 5 0.5 5
>> variable 3 0 0.5 3 0.5
>> variable 2 2 5 0.5 5
>>
>> However, I do not think it is just some straightforward expansion,
>> which could be done just by the subsetting mechanism of R. Is there
>> any idea on how can I do it for general case?
>>
>> Thanks for your time,
>>
>> ______________________________________________
>> 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