similar to: Repeated Measures

HI There, I was wondering if there was a procedure in R for carrying out greenhhouse geiser type repeated measures analysis and ante-dependence testing? I've checked out the nicely written manuals, but it seems to lack an entry in the index (maybe because this technique has not been implemented yet). THanks, Matt Redding Soil Scientist Intensive Livestock Environmental Management-Research

Hi All R guru's, A couple of simple questions which I cannot find answers to in "an introduction to R" or in the reference manuals. I am running R on Win95, and have plotted up quite a nice box and whisker plot: char<-read.table("n:/reddinm/grifpr~1/survey/char.txt") attach(char) boxplot(Al,Fe,OC,pH, range=0,xaxis=FALSE, ylab="Proportion of Unirigated

> Hi All R guru's, > > A couple of simple questions which I cannot find answers to in "an > introduction to R" or in the reference manuals. > > I am running R on Win95, and have plotted up quite a nice box and whisker > plot: > > char<-read.table("n:/reddinm/grifpr~1/survey/char.txt") > attach(char) > boxplot(Al,Fe,OC,pH,

HI All, I am trying to analyse the results of a heap of samples, to determine how many samples are required to meet a defined confidence interval. I checked out the R manual for bootstrapping and re-sampling methods, but did not find anything. Are there any contributed procedures that deal with this? thanks, Matt Redding Soil Scientist Intensive Livestock Environmental Management-Research

Hi All, I am preparing a few plots for publication. Unfortunately the journal requires figures to have all text in a sans serif fonts. I realise that there is support for the use of Hershey vector fonts in the function TEXT, but how do I get the text on axes to be sans serif? For example, the labels and axis values produced from this line of code: plot(SolPC, SorbPC, ylab="Sorbed P

As an alternative to the effects package, try predict() with type="terms" JM On 7 Oct 2005, at 8:00 PM, Peter Dunn wrote: > From: r-help-bounces at stat.math.ethz.ch > [mailto:r-help-bounces at stat.math.ethz.ch] On Behalf Of Peter Dunn > Sent: Wednesday, October 05, 2005 9:06 PM > To: R-help mailing list > Subject: [R] R/S-Plus equivalent to Genstat

HI All R Gurus, Please reply directly to this email if you have a suggestion, as I am not currently an R-help member. I would like to complete a non-linear regression with groups. Having gone down the path of completing such things under GENSTAT, I would prefer to do it with R. My example data set is... TimeTreat ExtrVol2 SorbedT ET1 39.5 384.5 ET1

Does anyone know if there is a package or other method of reading Genstat files directly into R. Genstat isn't listed in the foreign package. Many thanks, Graham [[alternative HTML version deleted]]

Hi all I'm doing some things with a colleague comparing different sorts of models. My colleague has fitted a number of glms in Genstat (which I have never used), while the glm I have been using is only available for R. He has a spreadsheet of fitted means from each of his models obtained from using the Genstat "predict" function. For example, suppose we fit the model of the type

I would like to announce that my DataLoad utility which can be found at: http://www.vsn-intl.com/genstat/downloads/datald.htm has been updated to support R data.frames. DataLoad reads a large variety of data formats (listed on the web page) and can convert these to ASCII or XDR data.frames. Only the Windows and Linux versions are up to date due to the fact that I no longer has access to a Unix

Hi all - R2.0.1 on OSX;MASS library;nlme library I am trying to emulate the solution to a problem set that has normally been run in Genstat, using R. The problem that I am having at the moment is with the following glmm question (using glmmPQL from the MASS library): "We have two different forest habitats (first rotation thicket, and high forest) which we want to survey for the

Hi, I know that is not possible make a stepwise procedure using REML in R, I can use ML for this. For nested design it may be very dangerous due the difference in variance structure, mainly in a splitplot design. ML make significative variables that REML dont make. I read an article that is made a stepwise procedure using GENSTAT. from article: "Terms were dropped from a model in a

Hello, Has anyone written R functions for applying self-controlled case series methods (http://statistics.open.ac.uk/sccs/). In fact only thing needed is to modify glm function to allow absorption of effect. Eg. in Poisson model individual effect is used as factor, but it is considered as nuisance term where parameter estimates are not needed. Could anyone point how absorbing individual

I would like to do a logistic multiple regression on a binary variate, with what Genstat calls an offset, ie a regression variable with the coefficient constrained=1. I'm far from home and have Splus but not R here, but any information from either source would be very welcome. I'd hate to have to go back to using Genstat . . . Thanks, Ted. --------------- Dr E.A. (Ted) Catchpole

Dear friends. I hope you will bear over with a lengthy question. I have a problem of 76 patients on hemodialysis measured repeatedly once a month for nine months for methylmalonic acid and vitamin B12. The vitamin had been administered previously to these patients but was curtailed in the first part and then started again. There was a very clear response in that as vitamin B12 went down

Hello everybody, I have a question. I would like to get a correlation between constitutive and induced plant defence which I messured on 30 plant species. So I have table with Species, Induced defence (ID), and constitutive defence (CD). Since Induced and constitutive defence are not independant (so called spurious correlation) I should do a randomisation test. I have a syntax of my

Hi R helpers, I'm trying to fit a rather complex model to some simulated data using lme and am not getting the correct results. It seems there might be some identifiability issues that could possibly be dealt with by specifying starting parameters - but I can't see how to do this. I'm comparing results from R to those got when using GenStat... The raw data are available on the

Hi All, I've been trawling through the documentation and listserv archives on this topic -- but as yet have not found a solution. I'm sure this is pretty simple with R, but I cannot work out how without resorting to ugly nested loops. As far as I can tell, grep, match, and %in% are not the correct tools. Question: given these vectors -- patrn <- c(1,2,3,4) exmpl <-

Hi R Gurus, I have a question that some of you may have dealt with before. I would redirect output to the clipboard, rather than to the screen. Is this possible? Is a function already written for it, or is it just an option in print? Any help would be appreciated. Also, thanks to all those who have helped develop this great software, or who make the effort here to support the user-base.

I am changing the title because this is really about the history of anova, and about strata in analysis of variance. As this kind of question has been arising very frequently, an extended comment may be in order. The ideas, and the sums of squares breakdowns, go back to Fisher; see in particular his "Design of Experiments", first published in 1935. This book is still a good