Hank, Attached is the dataframe that can be supplied for the argument comp.CAND.frm. The argument for comp.LIS.frm can be deleted since it currently isn't used. The other arguments can be set as: delta.0.Y.0 = 50 gamma = 40. cheers, dave
-----Original Message----- From: Martin Henry H. Stevens [mailto:[EMAIL PROTECTED] Sent: Wednesday, June 14, 2006 1:23 PM To: Afshartous, David Cc: [email protected] Subject: Re: [R] appending Hi David, It would be helpful if you supply a little data, upon which this would operate. Hank On Jun 14, 2006, at 12:31 PM, Afshartous, David wrote: > All, > > In the function below I have 24 individuals and 6 calculations per > individual. > The 6 calculations are collected each time in a 1:24 loop when > calculating "delta". > > I'd like to collect all 144 = 24*6 calculations in one vector > ("delta.patient.comb"). > The function works as is via indexing, but is there an easier way to > collect the measurements via appendinng the 6 measurements each time > to the current set? I couldn't find anything in Intro to R on > appending. > > cheers, > Dave > ps - please respond directly to [EMAIL PROTECTED] > > > > creatine.function.new = function(delta.0.Y.0, gamma, comp.LIS.frm, > comp.CAND.frm) { ## function to calcuate the delta.i, i.e. the percent > ## leftover ## gamma = rate of Cr going into bucket, e.g., mg/hr ## > delta.0.Y.0 = product of delta.0 and Y.0 at baseline ## > Y.1 = delta.0.Y.0 + gamma > delta = numeric(6) > delta.patient = numeric(24) > delta.patient.comb = numeric(144) > ## > for (k in 1:24) { ## each patient > patient.k.CAND = which(comp.CAND.frm$Patient_no == k) > Ucr.CAND.patient.k = comp.CAND.frm$Ucr[patient.k.CAND] > C = Ucr.CAND.patient.k ## 6 observed creatanine values for each > patient > delta[1] = (Y.1 - C[1])/Y.1 > Y.i = Y.1 > delta.i = delta[1] > for (i in 1:5) { ## six measurments per > patient > Y.i.plus.1 = delta.i * Y.i + gamma > delta.i.plus.1 = (Y.i.plus.1 - > C[i+1])/Y.i.plus.1 > delta[i+1] = delta.i.plus.1 > delta.i = delta[i+1] > Y.i = Y.i.plus.1 > } > delta.patient[k] = list(delta) > delta.patient.comb[(6*(k-1)+1):(6*(k-1)+ 6)] = delta > } > list(delta.patient, delta.patient.comb) } > > ______________________________________________ > [email protected] mailing list > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide! http://www.R-project.org/posting- > guide.html Dr. M. Hank H. Stevens, Assistant Professor 338 Pearson Hall Botany Department Miami University Oxford, OH 45056 Office: (513) 529-4206 Lab: (513) 529-4262 FAX: (513) 529-4243 http://www.cas.muohio.edu/~stevenmh/ http://www.muohio.edu/ecology/ http://www.muohio.edu/botany/ "E Pluribus Unum"
______________________________________________ [email protected] mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
