[R] re sults from "do.call" function

[kathie]

Dear R users...

I made this by help of one of R users.
_________________________________________________________________

X=matrix(seq(1,4), 2 , 2)
B=matrix(c(0.6,1.0,2.5,1.5) , 2 , 2)

func <- function(i,y0,j) { y0*exp(X[i,]%*%B[,j]) }
list1 <- expand.grid( i=c(1,2) , y0=c(1,2) , j=c(1,2) )
results <- do.call( func , list1 )
_________________________________________________________________

> results
              [,1]          [,2]         [,3]         [,4]          [,5]         
[,6]          [,7]          [,8]
[1,]  36.59823  36.59823  36.59823  36.59823   1096.633   1096.633  
1096.633   1096.633
[2,] 181.27224 181.27224 181.27224 181.27224  59874.142  59874.142 
59874.142  59874.142
[3,]  73.19647  73.19647  73.19647  73.19647   2193.266   2193.266  
2193.266   2193.266
[4,] 362.54448 362.54448 362.54448 362.54448 119748.283 119748.283
119748.283 119748.283
[5,]  36.59823  36.59823  36.59823  36.59823   1096.633   1096.633  
1096.633   1096.633
[6,] 181.27224 181.27224 181.27224 181.27224  59874.142  59874.142 
59874.142  59874.142
[7,]  73.19647  73.19647  73.19647  73.19647   2193.266   2193.266  
2193.266   2193.266
[8,] 362.54448 362.54448 362.54448 362.54448 119748.283 119748.283
119748.283 119748.283
_________________________________________________________________


However, I expected only this below as you imagine.
              [,1] 
[1,]  36.59823
[2,] 181.27224
[3,]  73.19647
[4,] 362.54448
 
                [,5]
[1,]    1096.633
[2,]   59874.142
[3,]    2193.266
[4,]  119748.283

That is, from "list1",

y0*exp(X[i,]%*%B[,j])
------------------------------------
1*exp(X[1,]%*%B[,1]) = 36.59823
1*exp(X[2,]%*%B[,1]) = 181.27224
2*exp(X[1,]%*%B[,1]) = 73.19647
2*exp(X[2,]%*%B[,1]) = 362.54448

1*exp(X[1,]%*%B[,2]) = 1096.633
1*exp(X[2,]%*%B[,2]) = 59874.142
2*exp(X[1,]%*%B[,2]) = 2193.266
2*exp(X[2,]%*%B[,2]) = 119748.283
------------------------------------

How could I find only 8 values?  Any suggestion will be greatly appreciated.

Kathryn Lord
-- 
View this message in context: http://www.nabble.com/results-from-%22do.call%22-function-tp20532136p20532136.html
Sent from the R help mailing list archive at Nabble.com.