This looks like homework. If it is, you should really tell us along with what
your teacher's policy is on getting help over the internet is (and note that
many teachers monitor this list and can see if you are getting help).
You have done the first part yourself, much better than some who have tried to
get us to do the whole thing for them, so a possible hint: the new problem
really has 3 groups, never sick, currently sick, and healed. Just modify your
current code to allow for people to move from the currently sick to the healed
group.
--
Gregory (Greg) L. Snow Ph.D.
Statistical Data Center
Intermountain Healthcare
greg.s...@imail.org
801.408.8111
-Original Message-
From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-
project.org] On Behalf Of Piotr Arendarski
Sent: Tuesday, April 13, 2010 4:27 PM
To: r-help@r-project.org
Subject: [R] Simulation problem.
Hi,
I have problem with simulating.
This is my task...
Suppose that there are N persons some of whom are sick with influenza.
The
following assumptions are made:
* when a sick person meets a healthy one, the chance is á that the
latter
will be infected
* all encounters are between two persons
Write a function which simulates this model for various values of
N (say, 10 000) and á (say, between 0.001 and 0.1). Monitor the
history of this process, assuming that one individual is infected at
the beginning.
The code is:
*
simulation - function(number, prob){
cumulative.time - 0
current.time - 0
number.sick - 1
while(number.sicknumber){
current.time - current.time + 1
meetings - rhyper(nn=1, m=number.sick, n=number-number.sick, k=2)
if(meetings==1){
one.sick - rbinom(n=1, size=1, prob)
if(one.sick==1){
cumulative.time - c(cumulative.time, current.time)
number.sick - number.sick + 1
}}}
cumulative.time
}
number - 1000
prob - .05
model - simulate(number, prob)*
But than add the assumption that *each infected person has a 0.01
chance
of recovering at each time unit* Do you have idea how to modify the
code
?
Piotr Arendarski
[[alternative HTML version deleted]]
__
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.
