matrix - Connecting matrices in R -


i show how many endosymbionts present in surviving hosts indicated final timestep in hosts matrix. corresponding columns in endos matrix @ final timestep added represent respective surviving hosts. if has thoughts on how show please let me know! here code far: 

set.seed(123) time  <-10  hosts <-matrix(na,nrow=(time+1),ncol=10)  endos <-matrix(na,nrow=(time+1),ncol=10) hosts[1,] <- seq(1,10) endos[1,] <- rep(10,10)  for(tt in 1:time) {    indexhost <- sample(1:10,size=10,replace=true,prob=null)   hosts[tt+1,] = hosts[tt,indexhost]    for(indexhost in 1:10) {     endos[tt+1,indexhost]<-rbinom(1,(2*endos[tt,indexhost]),0.5)    }  }

to find host survivors, can do:

survivors        = unique(hosts[time+1, ])

which remove duplicates in survivor vector. then, can index number of endos in survivors. i'm assuming dynamics endos independent of survivorship, can do:

endosinsurvivors = endos[time+1, survivors]

to calculate sum of endos id of host @ last time step, can use

tapply(endos[time+1, ], index=hosts[time+1, ], fun=sum) # using set.seed(1234)   2  7  79 38

Comments

Post a Comment

Popular posts from this blog

php - SPIP: From Tag directly to an article -

i work spip , plugin called tagsphere works fine tags in spip. have tag groupe put 1 article tag. when click on tag in tagsphere send me before menu can see articles linked tag. want send me directly article without see tag menu before. here code model of thagsphere: <div id="tagsphere-#env{id_article}"> <ul> <boucle_mot(mots){id_groupe ?}> <li>[<a href="#url_mot">(#titre)</a>]</li> </boucle_mot> </ul> </div> anybody idea? if want skip tag page, have link directly article. <b_article> block prevents display tag if no article linked. <div id="tagsphere-#env{id_article}"> <ul> <boucle_mot(mots){id_groupe ?}> <b_article> <li> <a<boucle_article(articles){...
Read more

jquery - isAjaxRequest always return false -

i got error message 404,'model has not been saved' following code my action controller public function actionextend($id) { $model=$this->loadmodel($id); if(yii::app()->request->isajaxrequest) { $model->todate = date('y-m-d h:i:s a',(strtotime('next month',strtotime(date('y-01-d'))))); if($model->save()) { yii::app()->end(); } else { throw new chttpexception(404,'model has not been saved'); } } } the jquery call $("#buttonextend<?php echo $this->post_row;?>").click(function(e){ $.ajax({ type: "post", url: "<?php echo yii::app()->createurl('post/extend', array('id' => $data['id'])); ?>", }); }); why not ajax call ? thank help i found response, missed } , function ...
Read more

ruby on rails - In a controller spec, how to find a specific tag in the generated view? -

i know integration tests preferred need ran in controller test, i'm testing gem injecting html code in view, xhr can't run in feature spec (if can, please explain me how :) ) so rspec controller tests can assert selector present (with capybara) : response.body.should have_selector('#foobar') has_selector? call all method capybara find selector. what want last child of body , assert id in particular. afaik it's not possible have_selector. what : all('body:first-child').first.id.should == '#foobar' however, capybara dsl, defined (more or less): def all(*args) page.all(*args) end and page empty unless use visit it's integrations specs. how can use capybara all method inside rspec controller test ? i can't test right after googling seems trick def page capybara::node::simple.new(response.body) end source
Read more