Re: [R] cbind results to original data frame

From: ONKELINX, Thierry <Thierry.ONKELINX_at_inbo.be>
Date: Thu, 29 May 2008 09:46:12 +0200

Have you tried to merge the two datasets? See ?merge for the details.

HTH, Thierry




ir. Thierry Onkelinx
Instituut voor natuur- en bosonderzoek / Research Institute for Nature and Forest
Cel biometrie, methodologie en kwaliteitszorg / Section biometrics, methodology and quality assurance
Gaverstraat 4
9500 Geraardsbergen
Belgium
tel. + 32 54/436 185
Thierry.Onkelinx_at_inbo.be
www.inbo.be

To call in the statistician after the experiment is done may be no more than asking him to perform a post-mortem examination: he may be able to say what the experiment died of.
~ Sir Ronald Aylmer Fisher

The plural of anecdote is not data.
~ Roger Brinner

The combination of some data and an aching desire for an answer does not ensure that a reasonable answer can be extracted from a given body of data.
~ John Tukey

-----Oorspronkelijk bericht-----
Van: r-help-bounces_at_r-project.org [mailto:r-help-bounces_at_r-project.org] Namens Ingrid Tohver
Verzonden: donderdag 29 mei 2008 3:23
Aan: r-help_at_r-project.org
Onderwerp: [R] cbind results to original data frame

I have the following script, which gives me a prediction for each of my observations (> 49k), and subsets my data frame by site (n = 183):

ts <- (by(dem16,dem16['Site'],function(dat)

		try(predict(nls(Tw ~ mu + ((alpha - mu)/
  	                   (1 + exp(gamma*(B - Mean_air)))),
  	         data = dat,
  	         start = list(mu = 0.8, alpha = 21.8, gamma = 0.22, B =

12.8),
		 algorithm = "port",
		 lower = list(0,0,0,0))))))

What I would like to do is apply the resulting "ts" list to my original data frame so I can perform a goodness of fit test. However, the resulting "ts" is a list only as long as the number of sites (183), so I can't force it to cbind to my original data frame with 49,000 observations. How can I format ts to cbind with dem16 (the original)?

Thank you,
Ingrid



R-help_at_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.

R-help_at_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. Received on Thu 29 May 2008 - 08:07:04 GMT

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.2.0, at Thu 29 May 2008 - 08:30:42 GMT.

Mailing list information is available at https://stat.ethz.ch/mailman/listinfo/r-help. Please read the posting guide before posting to the list.

list of date sections of archive