It is quite clear that you are confused, because your function is full of
variables that are not defined in your email. Unfortunately, that also makes
helping you quite difficult. Also, your (accidental?) use of HTML format means
that your code is garbled on our end (see below). You have to learn how to make
your email program send plain text to avoid this problem. (This point and
others are in the Posting Guide mentioned at the bottom if this and every email
on the list.
Usually your functions will be easier to understand and run faster if they
return a vector, and at the point where you call the function you assign that
return value to a new column in your data table. Normally we just put an
expression (often as simple as a variable name) alone on the last line of the
function in order to return it.
plusXY <- function( x, y ) {
result <- x + y
result
}
DF <- data.frame( a=1:5, b=5:9 )
str( DF )
DF$c <- plusXY( DF$a, DF$b )
str( DF )
The above steps are a reproducible example. Without you giving us a
reproducible example of the kind of code AND data you have it is very difficult
to understand your problem at the receiving end of an email. The link [1] below
gives advice on making such examples.
I also recommend that you (more carefully) read the Introduction to R document
that came with the software. In particular, study indexing or the various
members of the lapply family of functions, since rbind is quite an inefficient
way to add data to a data frame. Patrick Burns' "The R Inferno" [2] may provide
some useful guidance as well.
[1]
http://stackoverflow.com/questions/5963269/how-to-make-a-great-r-reproducible-example
[2] http://www.burns-stat.com/
---------------------------------------------------------------------------
Jeff Newmiller The ..... ..... Go Live...
DCN:<jdnew...@dcn.davis.ca.us> Basics: ##.#. ##.#. Live Go...
Live: OO#.. Dead: OO#.. Playing
Research Engineer (Solar/Batteries O.O#. #.O#. with
/Software/Embedded Controllers) .OO#. .OO#. rocks...1k
---------------------------------------------------------------------------
Sent from my phone. Please excuse my brevity.
On May 18, 2015 4:32:20 AM PDT, "Ariadna García Sáenz" <qui...@msn.com> wrote:
>Good afternoon,
>I'm working in a disease spread model, and I have to use a function to
>put information on a table.
>My question is that I don't know how to get the output of this function
>in a table. Thanks in advance!
>Here I paste my script:
>tbmodel <- function(time,S,E,I,De,Di,beta,alpha,phi,rho){ ts <-
>data.frame(time=sim.time,S=susceptibles,E=exposed,I=infectious,De=De,Di=Di)
>next.time <- ts
>while(next.time$time<limit.time&next.time$S>0&next.time$E>=0&next.time$I>=0)
># Conditions under which within-herd spread continues { next.time
><-gillespie(next.time$time,next.time$S,next.time$E,
>next.time$I,next.time$De,next.time$Di,beta,alpha,phi,rho) ts <-
>rbind(ts,next.time)
>if(ts$time[nrow(ts)]>limit.time|ts$E[nrow(ts)]<0)ts <- ts[-nrow(ts),]
># Elimination of last column (time> limit.time) }
>
>
> [[alternative HTML version deleted]]
>
>______________________________________________
>R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
>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.
______________________________________________
R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
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.
