Thanks you all!!!!

 
Maurice TEKO

 
Biostatisticien,Doctorant.

Université de Liège
Département des Sciences et Gestion de l'environnement
185 avenue de Longwy
6700 Arlon (Belgique)
Mail :teko_maurice@yahoo.fr
        :anoumou.tekoahatefou@ulg.ac.be
 http://www.ulg.ac.be
  



________________________________
 De : Rui Barradas <ruipbarradas@sapo.pt>
À : Jim Lemon <jim@bitwrit.com.au> 

Envoyé le : Mercredi 5 décembre 2012 16h26
Objet : Re: [R] Help for a function

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@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.
	[[alternative HTML version deleted]]


