Hi Ken, Your command works fine and your comment was right as well, here is the final command (inspired from yours) directly on DATA1.x and DATA1.y foo3<-list() for (i in 1:1000) foo3[[i]]<-nls( Qe~LgmFla(a,b,Ce), data=list(Ce=DATA1.x[,i],Qe=DATA1.y[,i]), start=c(a=0.5,b=0.1) ) Thanks a lot :) Regards/Cordialement
Benoit Boulinguiez _____ De : Kenn Konstabel [mailto:[EMAIL PROTECTED] Envoyé : lundi 27 octobre 2008 16:59 À : Benoit Boulinguiez Cc : r-help@r-project.org Objet : Re: [R] create list of data frames Hi, > To realize the data frame I've tried this > > for (i in 1:1000) > { > foo<-list(c(foo[],data.frame( Ce=DATA1.x[,i],Qe=DATA1.y[,i]))) > } I think the following would do it: foo <- list() for(i in 1:1000) foo[[i]] <- data.frame(Ce = DATA1.x[,i], Qe=DATA1.y[,i]) But then again, do you really need a list of data frames? You can run your analyses directly on DATA1.xyz... Best regards, Kenn On Mon, Oct 27, 2008 at 5:44 PM, Benoit Boulinguiez <[EMAIL PROTECTED]> wrote: Hi all, I need to realize nonlinear regression on a thousand data sets. I guess the lapply function would help me on that thus I'd like to create a list of data frames, each data frame containing the data as follows: Ce Qe 1 1.849147 0.1958672 2 10.054250 0.5771036 3 18.077246 0.7718514 4 27.576468 0.8079606 5 35.146862 0.8500489 6 43.245078 0.8366673 7 51.745760 0.8879672 8 59.767086 0.8860329 The data are in arrays DATA1.x which contains Ce and DATA1.y which contains Qe as follows: ...... [,988] [,989] [,990] [,991] [,992] [,993] [,994] [1,] 2.28584 1.464201 2.252296 2.050368 1.820814 1.460908 2.489860 [2,] 10.69722 10.060908 10.394847 10.305939 9.313467 10.600574 9.684784 [3,] 19.02983 19.025505 18.452735 18.764500 18.474322 18.906141 18.410527 [4,] 27.33536 27.083139 26.741224 26.943614 26.706248 26.888886 27.268291 [5,] 34.86732 34.636166 35.185439 35.407979 34.849281 34.195723 34.669640 [6,] 42.67302 43.202754 43.179338 43.530424 43.861326 43.888251 42.663388 [7,] 51.51916 52.331856 50.661131 51.272102 51.938743 51.866432 50.892845 [8,] 60.73692 60.468978 60.582377 59.610627 60.857182 59.293317 59.654958 [,995] [,996] [,997] [,998] [,999] [,1000] [1,] 2.303937 2.959057 2.471204 2.089243 2.673633 1.849147 [2,] 10.687650 10.647979 9.720253 10.546309 10.144749 10.054250 [3,] 18.787858 19.015047 18.086981 18.854964 18.577910 18.077246 [4,] 27.248525 26.615157 27.876889 27.137564 26.607100 27.576468 [5,] 35.461289 34.686118 35.122716 34.873995 34.953602 35.146862 [6,] 43.915466 42.534826 43.264155 43.105342 42.877609 43.245078 [7,] 51.738493 52.289541 51.456639 52.464800 51.406434 51.745760 [8,] 59.549702 59.817789 59.471866 60.149036 60.411239 59.767086 To realize the data frame I've tried this for (i in 1:1000) { foo<-list(c(foo[],data.frame(Ce=DATA1.x[,i],Qe=DATA1.y[,i]))) } but the program craches when I want to see foo. Regards/Cordialement ------------- Benoit Boulinguiez Ph.D Ecole de Chimie de Rennes (ENSCR) Bureau 1.20 Equipe CIP UMR CNRS 6226 "Sciences Chimiques de Rennes" Avenue du Général Leclerc 35708 Rennes CEDEX 7 Tel 33 (0)2 23 23 80 83 Fax 33 (0)2 23 23 81 20 <http://www.ensc-rennes.fr/> http://www.ensc-rennes.fr/ [[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. [[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.
