summing matrices to a single matrix in R -


i have following matrices:-

alive <- array(0,c(num_samples,num_time_periods,num_groups)) alive[,1,] <- 100   for(i in 2:num_time_periods){     alive[,i,] <- rbinom(num_samples, alive[,i-1,], exp(-delta[,i,]))} alive    , , 1         [,1] [,2] [,3] [,4] [,5]  [1,]  100   98   94   89   87  [2,]  100   98   96   94   92  [3,]  100   99   95   94   92    , , 2         [,1] [,2] [,3] [,4] [,5]  [1,]  100   98   94   89   87  [2,]  100   98   96   94   92  [3,]  100   99   95   94   92   , , 3         [,1] [,2] [,3] [,4] [,5]  [1,]  100   98   94   89   87  [2,]  100   98   96   94   92  [3,]  100   99   95   94   92

how sum of matrices element give me single matrix?

i have tried write this:-

 totalalive <- array(0,c(num_samples,num_time_periods,num_groups))     for(i in 2:num_groups){  totalalive[,,i] <- sum(alive[,,i])    }

but wrong. want single matrix below:-

sum:-

       [,1] [,2] [,3] [,4] [,5]  [1,]  300   294  ..   ..   ..  [2,]  300   294  ..   ..   ..  [3,]  300   297  ..   ..   ..

you can use apply collase on first 2 dimensions.

so if have 3x5x3 array mm 

mm<-structure(c(0, 0, 0, 4721.565, 4721.565, 4721.565, 4244.95, 4288.055,  4158.742, 3834.17, 3755.25, 3677.222, 3390.485, 3355.014, 3319.538,  0, 0, 0, 4310.424, 4310.424, 4310.424, 3873.18, 3912.528, 3794.482,  3498.017, 3425.959, 3354.792, 3092.831, 3060.493, 3028.149, 0,  0, 0, 3934.859, 3934.859, 3934.859, 3533.586, 3569.504, 3461.75,  3190.963, 3125.172, 3060.271, 2820.944, 2791.468, 2761.986), .dim = c(3l,  5l, 3l))

then

apply(mm, c(1,2), sum)

will return

     [,1]     [,2]     [,3]     [,4]     [,5] [1,]    0 12966.85 11651.72 10523.15 9304.260 [2,]    0 12966.85 11770.09 10306.38 9206.975 [3,]    0 12966.85 11414.97 10092.28 9109.673

Comments

Post a Comment

Popular posts from this blog

commonjs - How to write a typescript definition file for a node module that exports a function? -

consider, toml node module can use: // toml.d.ts declare module toml { export function parse(value:string):any; } declare module "toml" { export = toml; } then: /// <reference path="../../../../../defs/toml/toml.d.ts"/> import toml = require('toml'); toml.parse(...); however, node module export single function, such 'glob' ( https://github.com/isaacs/node-glob ). the node usage of module be: var glob = require("glob") glob("*.js", {}, callback(err, files) { ... }); you'd naively expect this: // glob.d.ts declare function globs(paths:string, options:any, callback:{(err:any, files:string[]):void; ...but because typescripts 'import' semantics bit strange, seems can use 'import .. = require()' statement alias modules . trying call: /// <reference path="../../../../../defs/glob/glob.d.ts"/> import blog = require('glob'); results in: error ts2072: mod...
Read more

openid - Okta: Failed to get authorization code through API call -

i'm integrating okta own idp server using okta's api. i'm implementing authorization code flow following steps below: in own server, use /api/v1/authn endpoint sessiontoken. use sessiontoken obtain authorization calling endpoint: /oauth2/v1/authorize?client_id=" + clientid + "&sessiontoken=" + sessiontoken + "&response_type=code&response_mode=query&scope=openid&redirect_uri=" + redirecturl + "&state=evanyang&nonce=" it's supposed return response status code 302 , location header containing redirect url code value. however, keep getting response status code 200 , without location header, html body saying "you using unsupported browser." , "javascript disabled on browser." according api documentation: http://developer.okta.com/docs/api/resources/oidc.html#authentication-request , sessiontoken parameter sufficient this: an okta one-time sessiontoken. allows api-based ...
Read more

ios - Change Storyboard View using Seague -

i have following code checks connection server, , depending on want change views if no connection found. nsstring *connect = [nsstring stringwithcontentsofurl:[nsurl urlwithstring:@"http://myserver.com"] encoding:nsutf8stringencoding error:nil]; if (connect == null) { [self performseguewithidentifier:@"network_failure" sender:self]; } else { // other code here } my problem not change views. have added storyboard segue between 2 views identifier of network_failure . not have navigation controller or that, problem? thanks help! no not problem. first of all, in objectivec exists nil object oriented null . must use that. so: if(connect == nil) or better: if(!connect) so due fact nil not same of null , probably code does't enter in if scope . you have been able discover breakpoint. anyway: be sure set trigger segue (line) 2 view controller , , name network_failure (that not best name. should instead view1toview2 ) as, m...
Read more