[R] Help for a function

Rui Barradas ruipbarradas at sapo.pt
Wed Dec 5 16:26:13 CET 2012


Hello,

Also, t1 and min(xt) do not vary inside the loop so if it enters the 
loop it never exits.
And

     res1[j] <-(a*h)
     res2 <-sum( res1[j])

is equivalent to

     res2 <- a*h

so the inner-most loop is not needed at all.

Hope this helps,

Rui Barradas
Em 05-12-2012 04:20, Jim Lemon escreveu:
> On 12/05/2012 01:01 AM, anoumou wrote:
>> Hello all,
>> I need a help.
>> I am modeling a disease and a create a R function like that:
>> ...
>> But i do not get the results,i try by all means but i d'ont 
>> understant the
>> problem.
>
> Hi anoumou,
> Your function provides almost no indication of what two of its five 
> arguments are supposed to be. If, with a certain degree of optimistic 
> inference, we suppose "x" to be a data frame organized as shown below 
> the function. "t", "i" and "CONTAGIEUX" must be the appropriately 
> named columns of the data frame. This leaves two columns, "Symptomes" 
> and "Incubation". Say we flip a coin to decide which of these to 
> assign to "r", and with all of our degrees of freedom gone, we assume 
> that the other is "h". We are forced to the conclusion that "a" is a 
> nuisance argument, added to throw us off the scent.
>
> Peering within the function, we notice that the date format is wrong, 
> braces are unmatched and that our data frame is alphabetically ordered 
> by name of country to no purpose whatsoever. The best I can do here is 
> to make the function potentially able to do something if you can work 
> out what to do with it.
>
> Lambda<-function (x,date1,r,h,a) {
>  ndate1 <- as.Date(date1, "%Y-%m-%d")
>  t1 <- as.numeric(ndate1)
>  x[order(x$i),]
>  xt <-x[,"t"]
>  xi <-x[,"i"]
>  CONTAGIEUX <-x[,"CONTAGIEUX"]
>  while ( t1 < min(xt) ){
>   for (i in 1:length(xi) ){
>    for (j in 1:CONTAGIEUX[length(CONTAGIEUX)]){
>     res1[j] <-(a*h)
>     res2 <-sum( res1[j])
>    }
>   }
>  lambda[i] <- r*res2
>  }
>  x<-data.frame(x,lambda)
>  x
> }
>
> Jim
>
> ______________________________________________
> R-help at 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.




More information about the R-help mailing list