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 [EMAIL PROTECTED] 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: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] Namens Ingrid Tohver Verzonden: donderdag 29 mei 2008 3:23 Aan: [email protected] 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 ______________________________________________ [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 and provide commented, minimal, self-contained, reproducible code. ______________________________________________ [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 and provide commented, minimal, self-contained, reproducible code.
