[R] Help for a function

anoumou teko_maurice at yahoo.fr
Tue Dec 4 15:01:51 CET 2012


Hello all,
I need a help.
I am modeling a disease and a create a R function like that:

Lambda<-function (x,date1,r,h,a){
  ndate1 <- as.Date(date1, "%d/%m/%Y")
  t1 <- as.numeric(ndate1)
  x[order(x$i),]
  t <-x[,"t"]
  i <-x[,"i"]
  CONTAGIEUX <-x[,"CONTAGIEUX"]
  while ( t1 < min(t) ){
  for (i in 1:length(i) ){
    {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
}

on such data :

 	DATE	i	Symptomes	t	Incubation	CONTAGIEUX
1	2009-04-29	 Canada 	13	14363	13	13
2	2009-05-01	 Israel 	2	14365	2	2
3	2009-05-09	argentina	1	14373	1	1
5	2009-05-09	australia	1	14373	1	1
6	2009-05-10	australia	1	14374	2	2
7	2009-04-29	Austria 	1	14363	1	1
8	2009-04-30	Austria 	1	14364	2	2
9	2009-05-01	Austria 	1	14365	2	3
10	2009-05-02	Austria 	1	14366	2	4
11	2009-05-03	Austria 	1	14367	2	5
17	2009-05-09	Austria 	1	14373	2	7
18	2009-05-10	Austria 	1	14374	2	7
19	2009-05-08	brasil	4	14372	4	4
20	2009-05-09	brazil	6	14373	6	6
21	2009-05-10	brazil	6	14374	12	12
22	2009-05-02	canada	51	14366	51	51
23	2009-05-03	canada	85	14367	136	136
24	2009-05-04	canada	101	14368	186	237
31	2009-04-27	Canada	6	14361	6	6
32	2009-04-28	Canada 	6	14362	6	6
33	2009-04-30	Canada 	19	14364	25	25
34	2009-05-01	Canada 	34	14365	53	59
35	2009-05-01	China,HongKong, SAR	1	14365	1	1
36	2009-05-02	China,HongKong, SAR	1	14366	2	2
37	2009-05-03	China,HongKong, SAR	1	14367	2	3
38	2009-05-04	China,HongKong, SAR	1	14368	2	4
44	2009-05-10	China,HongKong, SAR	1	14374	2	7
45	2009-05-04	colombia	1	14368	1	1
46	2009-05-05	colombia	1	14369	2	2
47	2009-05-06	colombia

But i do not get the results,i try by all means but i d'ont understant the
problem.
Thanks for your help.



--
View this message in context: http://r.789695.n4.nabble.com/Help-for-a-function-tp4652054.html
Sent from the R help mailing list archive at Nabble.com.



More information about the R-help mailing list