Skip to content

Problems dealing with matrices

1 message · halim10-fes

#
Hi Arun,

Thank you very much for your response. Sorry, if I couldn't explain clearly. I 
think, I should restate the problem to get exactly what I want. Here it goes:

I have 2 matrices and 1 vector, namely,

dcmat<-matrix(c(0.13,0.61,0.25,0.00,0.00,0.00,0.52,0.37,0.09,0.00,0.00,0.00,
                0.58,0.30,0.11,0.00,0.00,0.00,0.46,0.22,0.00,0.00,0.00,0.00,
                0.09),nrow=5,ncol=5)

volmat<-matrix(c(100,0,0,0,0),nrow=5,ncol=1)

volinp<-c(1:40)

What I essentially want to do is to multiply 'dcmat' with 'volmat' and dump 
the output in a new matrix 'vol'. But before that, in the first step, I want 
to add volinp[1] with volmat[1,1]. So, the first column of the output matrix 
'vol' matrix will be:

        [,1]
[1,]   13.13
[2,]   61.61
[3,]   25.25
[4,]    0.00
[5,]    0.00

In the 2nd step, I want to add vol[1,1] with volinp[2]. and replace this value 
with volmat[1,1]. The new 'volmat' will look like:

       [,1]
[1,]   15.13
[2,]    0.00
[3,]    0.00
[4,]    0.00
[5,]    0.00

And then multiply 'dcmat' with the new 'volmat', and the 2nd column of output 
matrix 'vol' will look like:

       [,2]
[1,] 1.9669
[2,] 9.2293
[3,] 3.7825
[4,] 0.0000
[5,] 0.0000

Then again, replace the volmat[1,1] with vol[1,2] + volinp[3] and multiply the 
new 'volmat' with 'dcmat'. This addition, multiplication, and dumping will 
continue up to the length of 'volinp' and the final output matrix 'vol' will 
be something like:

      [,1]    [,2]      [,3]  ...length(volinp)
[1,] 13.13  1.9669  0.645697  ...
[2,] 61.61  9.2293  3.029809  ...
[3,] 25.25  3.7825  1.241725  ...
[4,]  0.00  0.0000  0.000000  ...
[5,]  0.00  0.0000  0.000000  ...   

Within my limited capacity, I've tried to come up with a solution but failed. 

I'll appreciate your/others' help with gratefulness.

Regards,

Halim

---------------
Md. Abdul Halim
Assistant Professor
Department of Forestry and Environmental Science
Shahjalal University of Science and Technology,Sylhet-3114,
Bangladesh.
Cell: +8801714078386.
alt. e-mail: xou03 at yahoo.com
 
On Sat, 23 Nov 2013 14:10:12 -0800 (PST), arun wrote
description.
volinp<-