similar to: [OFF] stepwise using REML???

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

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

Hi, I have this: lme(y~x1+x2,random=~1|x1/x2) How to make this random effect using lmer? I try this: lmer(y~x1+x2+(1|x1/x2) But it dont work. Any idea? Thanks Ronaldo -- System halted! -- |> // | \\ [***********************************] | ( ?? ?? ) [Ronaldo Reis J??nior ] |> V [UFV/DBA-Entomologia ] | / \ [36570-000 Vi??osa -

Hello all, I am trying to figure out an optimal linear model by using stepwise regression which requires partial f-test, I did some Googling on the Internet and realised that someone seemed to ask the question before: Jim Milks <jrclmilks at joimail.com> writes: > Dear all: > > I have a regression model that has collinearity problems (between > three regressor variables). I

Hi R Gurus, I am more or less a beginner with R. I have searched your FAQ, manuals, and the archives of this list for an answer to this question, unsuccessfully. I commonly do analysis of chemical characteristics of soils, by depth intervals (e.g. 0-100 mm, 100-200, 200-300, etc), comparing two soil profiles that have received some treatment (effluent). In Genstat I can handle this problem

Hi, I make this model using lme m.lme <- lme(Glycogen~Treatment,random=~1|rTrt/Liver) How to make this using lmer? I try > m.lmer <- lmer(Glycogen~Treatment+(1|rTrt/Liver)) Erro em lmer(Glycogen ~ Treatment + (1 | rTrt/Liver)) : entry 0 in matrix[0,0] has row 2147483647 and column 2147483647 Al??m disso: Mensagem de aviso: / not meaningful for factors in: Ops.factor(rTrt, Liver)

Hi! How to calculate the correct SE of mean in a nested or spliplot anova? Nested example: --------------------- m <- aov(Glycogen~Treatment+Error(Treatment/Rat/Liver)) > m Call: aov(formula = Glycogen ~ Treatment + Error(Treatment/Rat/Liver)) Grand Mean: 142.2222 Stratum 1: Treatment Terms: Treatment Sum of Squares 1557.556 Deg. of Freedom 2 Estimated

Hi all, We've just had a paper accepted for publication in Nature. We used R for 95% of our analyses (one of my co-authors sneaked in some GenStat when I wasn't looking.). The preprint is available from the Nature web site, in the open peer-review trial section. I searched Nature for previous references to "R Development Core Team", and I received no hits. So I tentatively

Hi all I have another GLIM question. I have been using R as well as Genstat (version 6) in order to fit GLIM models to the data (displayed below). The same models are fitted but the answers supplied by the two packages are not the same. Why? Can anyone help? A discription of the data and the type of model/s fitted can be found below. Regards Allan The

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 guys I've just had a thought about something we do with our search on ReportBuyer.com: we cater for both American and British English in our searches and we have had plans for a while now to implement something that allows users to find 'colour' and 'color', 'tap' and 'faucet' by doing some clever programming. Looking at the Xapian docs, though, it appears

Is there a dataload utility for linux. The link in genstat is down but I managed to find the utility at: http://gurukul.ucc.american.edu/econ/gaussres/UTILITYS/DATALOAD.HTM but this is a dos/windows version. Thank you Jean

<URL: http://finzi.psych.upenn.edu/R/library/datasets/html/morley.html > "The classical data of Michaelson and Morley on the speed of light" Can you provide more information about the data? How were they obtained, etc.? I do not have the book "Genstat Primer" and the nearest location where it is available is University of York which is rather far from my location.

Hello everybody, I am wondering if it is possible to create a gui to plot a time series that is very big, it's an EEG signal of 20mins. What I would like to do is plot the first 5mins, then have a button on the gui that plots the next 5mins when pushed. Is it possible? Thanks in advance ! Gael.

I've been asked about "Canonical Variance Analysis" (CVA). I don't see any reference to it searching the R site. Does it go by other names? Genstat describes it thus: Canonical variates analysis operates on a within-group sums of squares and products matrix, calculated from a set of variates and factor that specifies the grouping of units. It finds linear combinations of the

Hello r-help, I'm trying to fit birds counts over years using glm. In fact I'm trying to reproduce an analysis already perform with genstat (attach document). I have done (with Estate and year as factors): Model1 <- glm(Females~Estate+Year+offset = log(area)), family = quasipoisson(link = log), na.action = "na.exclude") After I have calculated the prediction using: Pred1

Good day, I'm ran aov () and lme() on a split-plot using a 64-bit machine. For aov() I don't see the values for ErrorA, F-value and p-value in the output. For lme(), output is different from results from a 32-bit. Please see codes used and corresponding output. Is my code wrong and/or not sufficient or is this a compatibility issue? ************** model1<-aov(Y~Main*Sub +

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