similar to: GLME

Dear All, is any way I can change the colours of the box and umbrella of a bwplot without having to go to the length of: box.rectangle<-trellis.par.get("box.rectangle") box.rectangle$col<-"black" trellis.par.set("box.rectangle", box.rectangle) etc... but straight from the call: bwplot(y ~ x | z, mydata) ? Regards, Federico Calboli

Dear All, using R is noticed that I would benefit by going back to the books and "refresh" my linear algebra and calculus. Has anybody suggestions on which books are better as reference on the two these two topics? I was thinking to get Lang's "A First Course in Calculus", but I I have no idea what could be a good one for linear algebra. Something detailed and clear

Dear All, I would like to ask you how to accomplish a little tricky data manipulation. I have a large dataset, looking something like: temp line cage number 18 18 1 6678.63 18 18 1 7774.458 18 18 1 7845.902 18 18 1 9483.578 18 18 1 8983.555 18 18 1 9181.052 18 18 1 9458.696 18 18 1 8138.616 18 18 1 7981.994 18 18 1 7556.491 18 18 1 7672.137 18 18 1 6607.776 18 18 1 8383.65 18 18 1 7129.852 18 18

Dear All, I would like to ask you to help me with my memeory. I remember using some function that would list all the possible methods I could apply to an object. Say, if I had an object of class=lme, it would tell me that that I could do stuff like qqnorm(myobjct), or VarCorr(myobject). In general, a very complete list. I though this list of all possible methods would pop out by typing

Dear R community, I need to plot the results of some simulations I did using QTL Cartographer. I am plotting LOD scores over three chromosomes. The three plot have to be one next to the other. The procedure I am using is: par(mfrow=c(1,3)) plot(x$x, x$y, ylim=c(0,35), type="l", col="blue", las=1, xaxs="i", yaxs="i", xlab="X Chromosome",

Dear All, I am analysing some data for a colleague (not my data, gotta be published so I cannot divulge). My response variable is the number of matings observed per day for some fruitlies. My factors are: Day: the observations were taken on 9 days Regime: 3 selection regimes Line: 3 replicates per selection regime. I have 81 observations in total The lines are coded A to I, so I do not need

Hi All, I would like to ask a question on fixed and random effecti in lme. I am fiddlying around Mick Crawley dataset "rats" : http://www.bio.ic.ac.uk/research/mjcraw/statcomp/data/ The advantage is that most work is already done in Crawley's book (page 361 onwards) so I can check what I am doing. I am tryg to reproduce the nested analysis on page 368:

Dear All, I am writing to ask a clarification on what R, and in particular lme, is doing. I have an experiment where fly wing area was measured in 4 selection lines, measured at 18 and 25 degrees. I am using a lme model because I have three replicated per line (coded 1:12 so I need not use getGroups to creat an orederd factor). The lines are called: "18"; "25";

Hi All, I woul like to ask you a couple of questions on chisq.test. First, I have 40 flies, 14 males and 26 females and I want to test for an a priori hypothesis that the sex ratio is 1:1 sex<-c(14,26) pr<-c(1,1)/2 chisq.test(se, p=pr, correct=TRUE) Chi-squared test for given probabilities data: sex X-squared = 3.6, df = 1, p-value = 0.05778 If my calculations are correct, this is

Dear All, recently, while setting me on the straight and narrow about linear contrasts for a linear mixed effect model, Prof Ripley pointed out that my interpertation of the call anova(mymodel.lme) was not correct, because I was meant to add type = "marginal", as in anova(mymodel.lme, type = "marginal") I tried to look deeper in the issue, asking people, checking on the

Hello! could you give me some advice where i can finde out/recieve/download a GLME package, written by Jose Pinheiro. i couldn'ti find it in the package lists (probably becouse it is still in a beta version). thank you for your replay

Thanks to those who responded to my inquiry about generalized linear mixed models on R and S-plus. Before I summarize the software, I note that there are several ways of doing statistical inference for generalized linear mixed models: (1)Standard maximum likelihood estimation, computationally intensive due to intractable likelihood function (2) Penalized quasi likelihood or similar

Dear All, I suspect this is kind of dumb, but when I was under the thrall of the dark lord (read, using a W2K box), all my work in R files came out as foo.RData. I moved on to GNU/Linux, and all the old .RData files keep on working as they used. No problems in loading and stuff. But I use R from the terminal. Assuming I decide to switch to emacs, do I need to save my work as foo.r? what about my

I've noticed quite a few messages relating to large data sets bedeviling R users, and having just had to program my way through one that actually caused a "Bus error" when I tried to read it in, I'd like to ask two questions. 1) Are there any facilities for aggregation of data in R? ( I admit that this will not do much for the large data set problem immediately) 2) Is there any

Dear All, when I try to: install.packages("Design"); install.packages("Hmisc") I get the following error messages: * Installing *source* package 'Design' ... ** libs g77 -mieee-fp -O2 -fomit-frame-pointer -pipe -march=i586 -mcpu=pentiumpro -O2 -fomit-frame-pointer -pipe -march=i586 -mcpu=pentiumpro -c lrmfit.f -o lrmfit.o make: g77: Command not found make: ***

Dear List, While trying to define a customised CITATION file for a package, following R-exts, I realised that if I use only one 'citEntry' I got both a text description and a BibTex entry for the package, as for the 'nlme' package: -------------------------- citation('nlme') To cite package 'nlme' in publications use: Jose Pinheiro, Douglas Bates, Saikat

Function in the base library of R 1.8.0 seems to use a search path with base coming before the local environment. I think this is intentional and is related to base being a namespace. > log <- function(x, base) 200 > log(1) [1] 200 > log10 function (x) log(x, 10) <environment: namespace:base> > log10(1) [1] 0 > This is most problematic when you are creating a

Dear All, I would like to ask a couple of questions on a LME model. I tested 4 selection lines at 4 food concentrations against a standard competitor stock. I had 3 replicate cages per selection line. In each cage I have 10 vials. I counted the number of wild type flies and competitor stock emerging in each vial. My main question is: is there any difference between selection lines? I did fit

If .../<pkg>/man/ contains a filename beginning with [, such as [.datetimes, Rdconvlib.pl cannot convert that help file into latex at INSTALL time. [Using --no-latex is a workaround.] The error given is: # R INSTALL ssCatalogues : help >>> Building/Updating help pages for package `ssCatalogues' Formats: text html latex example /^[.datetimes$/: unmatched [] in regexp at

Forgive me if this isn't relevant, but: are you stripping off the output you need (e.g., parameter estimates) and throwing away what you don't (e.g., residuals)? Reid Huntsinger -----Original Message----- From: Saikat DebRoy [mailto:saikat at stat.wisc.edu] Sent: Wednesday, March 21, 2001 12:15 PM To: Han Lai Cc: R-help Subject: Re: [R] Is there any help for "Cannot allocate vector