RE: [R] panel ordering in nlme and augPred plots
Hi Petr try plot(...,as.table=T,) Christian -- Mensaje Original -- From: Petr Pikal [EMAIL PROTECTED] To: r-help@stat.math.ethz.ch Date: Mon, 25 Apr 2005 16:40:40 +0200 Subject: [R] panel ordering in nlme and augPred plots Dear all I am trying nlme together with Pinheiro/Bates book. I constructed grouped data object with suitable plotting layout (according to some common factor, panels from bottom to top are in increasing order). When I do nlme(... some stuff...) I get fitted object which I can plot with plot(augPred(fit.nlme6, level=0:1)) but it results in completely different ordering. Is there any way how I can plot panels in some defined order e.g. ord.f-order(my.1.fac, my.2.fac) plot(augPred(fit.nlme6, level=0:1), ord.f) The only thing I found out is that if I order grouped.data object gr.dat-gr.dat[ord.f,] and do nlme fit, then the ordering in augPred plot is OK. BTW is there a way how to specify in plot row1 - 3 panels row2 - 3 panels row3 - 3 panels row4 - 2 panels row5 - 4 panels Thank you Best regards Petr Pikal [EMAIL PROTECTED] __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
RE: [R] mixed effects model:how to include initial conditions
Christoph, If you take a look at journal articles related to this topic published elsewhere, you will find that the most common analysis is: model-lme(growth~initialsize+block+treatment+) maybe it would be important to check if your model really needs the interactions you pointed out. I would suggest to try simpler models firstwhich generally work OK in analysis of growth of plants, trees, etc CM -- Mensaje Original -- Date: Fri, 21 Jan 2005 12:30:07 +0100 From: Christoph Scherber [EMAIL PROTECTED] To: r-help@stat.math.ethz.ch Subject: [R] mixed effects model:how to include initial conditions Dear R users, I am analyzing a dataset on growth of plants in response to several factors. I am using a mixed-effects model of the following structure: model-lme(growth~block*treatment*factor1*factor2, random=~1|plot/treatment/initialsize) I have measured the initial size of the plants (in 2003) and thought it might be sensible to include this (random) variation into the random effects term of the model. Is that correct? Or should initialsize rather be included as a covariate into the fixed effects term, as in: alternative-lme(growth~block*initialsize*treatment*factor1*factor2, random=~1|plot/treatment) I would very much appreciate any suggestions on how to analyze these data correctly. Best regards Chris. __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
RE: [R] nls question
Mike nlsList from nlme library can fit nonlinear models for dataset grouped by some specification, e.g. by specie in your case Regards Christian -- Mensaje Original -- From: Mike Saunders [EMAIL PROTECTED] To: R Help [EMAIL PROTECTED] Date: Thu, 16 Dec 2004 10:40:00 -0500 Subject: [R] nls question Just a quick question. Is there a way to easily specify factor levels in a function definition within nls? For example, I am trying to fit a 3 parameter, nonlinear Weibull function to tree height and I would like to have results by species (or down the road a bit, by plot). I am hoping there is someway to do it easily, similar to the gas example in the ANCOVA chapter in Modern Applied Statistics with S, and not running it with a by function or making a very complex function statement with all factor levels combinations. Can someone point me in the right direction? Thanks, Mike Mike Saunders Research Assistant Forest Ecosystem Research Program Department of Forest Ecosystem Sciences University of Maine Orono, ME 04469 207-581-2763 (O) 207-581-4257 (F) [[alternative HTML version deleted]] __ [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 __ [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
[R] plot.lm basic question
Hi all, Does anyone know how can I change the color of the 'captions' and 'point labels' in the plots obtained from plot.lm? Thanks for any hint Christian __ [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
[R] proxy
Hi all; I?m working with the latest version of R under Win2000. My internet connection (using a proxy) requires my username and password. In this situation I?m not able to update the packages from CRAN (unless I download and install them from a local drive). I tried looking at information on the list servers but I couldn?t figure out how to change 'environment variables' or something like that. I also tried renaming the files internet2.dll to internet.dll and so on. Any 'new' idea how to solve this situation. Thanks in advance Christian Mora __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
[R] Selecting subsamples
Hi all, I?m working with a dataset with 9 columns and 2000 rows. Each row represents an individual and one of the columns represents the volume of that individual (measured in cubic meters). I?d like to select a sample from this dataset (without considering any probability of the rows) in which the sum of the volume of the individuals in that sample = 100 cubic m. I?ll appreciate any suggestion Thanks CM __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help
