Rui , I discover all this when running and rerunning my function one and thousand times and per line. Now it's perfect, and work so well. Thanks.
________________________________ À : Rui Barradas <ruipbarra...@sapo.pt>; Jim Lemon <j...@bitwrit.com.au> Cc : "r-help@r-project.org" <r-help@r-project.org> Envoyé le : Mercredi 5 décembre 2012 16h45 Objet : Re: [R] Help for a function 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) :anoumou.tekoahate...@ulg.ac.be http://www.ulg.ac.be ________________________________ De : Rui Barradas <ruipbarra...@sapo.pt> À : Jim Lemon <j...@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]]
______________________________________________ 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.
