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: xo...@yahoo.com On Sat, 23 Nov 2013 14:10:12 -0800 (PST), arun wrote > Hi, > Could you show your expected output? It is a bit unclear from the description. > > On Saturday, November 23, 2013 2:00 PM, halim10-fes <halim10- > f...@sust.edu> wrote: Dear R-friends, > > Hope you doing well. I've been trying to deal with the following > problem for the couple of days but couldn't come up with a solution. > It would be great if any of you could give some insight into it. > > I have three matrices like: > > dcvol<-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) volinp<- > matrix(c(100,0,0,0,0),nrow=5,ncol=1) > > scvol<-matrix(c(1:40),nrow=5,ncol=8) > > What I essentially want to do is to add each value in scvol[1,] with > the volinp[1,1] and then multiply each new volinp with dcvol and > finally put the outputs in a new matrix. > > Thanks in advance. > > 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: xo...@yahoo.com > > -- > This message has been scanned for viruses and > dangerous content by MailScanner, and is > believed to be clean. > > ______________________________________________ > R-help@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. -- This message has been scanned for viruses and dangerous content by MailScanner, and is believed to be clean. ______________________________________________ R-help@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.