similar to: Background color(s) for groupedData plot

Sorry if this is obvious, but my S skills aren't great and I haven't been able to find it documented anywhere. I want to write a new function for use with lme objects; the function will simply calculate an ICC (aka "rho") for each level of a mixed-effects model. What I need for this is pretty simple: (c(var1..varn, residual)) / sum(c(var1..varn, residual)) where var1..varn

I have a number of plots that I'm trying to do using the lattice package, unfortunately, they involve Posix values. A small sample of the data is as follows: Time TC.1 <Snipped> 1 2003-04-24 13:29:15 26.61 2 2003-04-24 13:29:30 26.48 3 2003-04-24 13:30:00 26.53 4 2003-04-24 13:30:30 27.85 <Snipped> where str(data) `data.frame': 415 obs. of

Dear r-help, I've got data of the following structure 1979 93.428747 0 1979 87.298608 20 1979 78.506340 40 ... 1979 45.567890 340 1980 60.815289 0 1980 49.630904 20 1980 24.981362 40 ... The first column is year and the last one is the longitude. I need a set of graphs showing the dependence of the middle value on the longitude, for each year, situated one blow the other.

Dear all, Has anyone written functions to compute (and ideally project / plot onto a plane) confidence ellipsoids for q linear model parameters, q > 2 ? Greetings and best wishes, Renaud -- Dr Renaud Lancelot, v?t?rinaire CIRAD, D?partement Elevage et M?decine V?t?rinaire (CIRAD-Emvt) Programme Productions Animales http://www.cirad.fr/presentation/programmes/prod-ani.shtml ISRA-LNERV

To Windows R users: The binary build of R that I uploaded on Tuesday contained errors in the compiled HTML help files and some of the PDF versions of the manuals. In particular, the index was missing from the CHM files, and the table of contents was missing from at least one of the manuals. I've since found the error in my setup, and have rebuilt the distribution. I uploaded it to CRAN

To Windows R users: The binary build of R that I uploaded on Tuesday contained errors in the compiled HTML help files and some of the PDF versions of the manuals. In particular, the index was missing from the CHM files, and the table of contents was missing from at least one of the manuals. I've since found the error in my setup, and have rebuilt the distribution. I uploaded it to CRAN

Hi, Does anyone know how to make the color of the lines all black when plotting groupedData with an outer factor: For example, library(nlme) plot(Dialyzer, outer=~QB, key=F) This generated colored curves in R.2.4.1. How to make all the curves black ? (or how to alter the color (type) of lines for the nlme groupedData plotting function in general?) Thanks Qiong

[I have abbreviated the subject as jitterbug has been having probems with long subjects.] The issue here is that one cannot mix log/non-log axes in the calls to screen(), as the appropriate par() parameter is read-only, but the meaning of yaxp depends on it. But beyond that you can't set x/yaxp for log axes. You should be able to do this: you can in the S original. A simpler version: plot

> version _ platform i386-pc-mingw32 arch i386 os mingw32 system i386, mingw32 status major 1 minor 5.0 year 2002 month 04 day 29 language R OS Win98 I have problems to install some of the packages dowloaded from CRAN (whatever the

Please, forgive my ignorance on statistics, but I have a rather simple question concerning AIC. Small values of AIC mean good fit of the model. How about negative values of AIC? Is a model with AIC=0.5 considered to have better fit than a model with AIC=-500? I couldn't find anything explaining this in the documentation or any elementary statistics textbook that I have at my disposal. Thanks

? stato filtrato un testo allegato il cui set di caratteri non era indicato... Nome: non disponibile URL: <https://stat.ethz.ch/pipermail/r-help/attachments/20101011/7fd9c9d5/attachment.pl>

I've put up a build of the current R-patched (and updated recommended packages) for Windows at http://www.stats.ox.ac.uk/pub/bdr/RWin/Rpre/SetupRpatched.exe This should be able to rename files on Win 9x (unlike the mid-May build) and can build packages accessing Rblas.dll. As ever, this is less tested than an official release, in particular only on Windows XP, but ought to have fewer

On Fri, 3 May 2002 drew.tyre@csiro.au wrote: > Full_Name: Drew Tyre > Version: 1.4.0 > OS: Win 2000 > Submission from: (NULL) (140.253.74.225) > > > When creating a rather complex plot, the resulting figures 'shrink' when either > copied to the clipboard or saved as a metafile using the GUI menu. Also saving > as postscript or pdf causes the same problem, as

Dear all, I am trying to avoid the warnings produced by: > x <- -2:2 > log(x) [1] NaN NaN -Inf 0.0000000 0.6931472 Warning message: production de NaN in: log(x) I thought that using ifelse would be a solution, but it is not the case: > ifelse(test = x < 0, yes = NaN, no = log(x)) [1] NaN NaN -Inf 0.0000000 0.6931472 Warning message: production

Hi, lme() users, Can some one tell me how to do this. I model Orthodont with the same G for random variables, but different R{i}'s for boys and girls, so that I can get sigma1_square_hat for boys and sigma2_square_hat for girls. The model is Y{i}=X{i}beta + Z{i}b + e{i} b ~ iid N(0,G) and e{i} ~ iid N(0,R{i}) i=1,2 orth.lme <- lme(distance ~ Sex * age, data=Orthodont, random=~age|Subject,

Greetings everyone, I am new at R. My questions is rather basic. Looking R manuals looks like there should be a way to read MS excel files into R. Could somebody tell me which library should I use for that? Thanks in advance Victor H. Marín Laboratorio de Modelación Ecológica Depto. de Ciencias Ecológicas, Fac. de Ciencias, Universidad de Chile Casilla 653 Santiago, Chile

Hello R-users, I'm using the boot's library. I have two questions: 1. I execute for example data<-rnorm(100) boot(data,mean,500) and I obtain 500 equal bootstrap samples of data. Why? 2. I need to calculate the bootstrap confidence intervals for a statistic, but when I use the function boot.ci it shows the error message: Error in if

hi, during a for(i in 1:xx]) loop I always newly calculate a vector (e.g. tmp<-c(a,b,c,d,e) ) now i need that vector to be attached at the bootom of a matrix (or data.frame). e.g. m<-matrix() for(i in 1:5]){ #some calculations for a,b,c,d,e a<- ... b<- ... c<- .. .... tmp<-c(a,b,c,d,e) ??? now I need to attach this tmp to the matrix m ??? } i couldn't get rbind

Dear all, 1. How can I compute the intersecting area between 2 polygons ? 2. I have polygons with coordinates in decimal degrees (i.e. 13 deg 30 min = 13.5 decimal degrees). I want to compute their area and get the results in square meters or square kiometers. Can anyone give me a conversion coefficient or a pointer where I can find this information (sorry for this off topic question) ? Thanks

Hi there, I am wanting to create 8 side-by-side histograms which have been rotated 90 degrees clockwise from how they usually sit.. all with the same scales. Is someone able to help me out? Thanks so much, Rachel Cunliffe -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send