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_maur...@yahoo.fr :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.