similar to: R/S-Plus equivalent to Genstat "predict": predictions over "averages" of covariates

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

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 am converting an S-Plus function into R. The S-Plus code uses some of the glm families, and family objects. The family objects in S-Plus and R have many different features, for example: In R: > names(Gamma()) [1] "family" "link" "linkfun" "linkinv" "variance" [6] "dev.resids" "aic"

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 all I'm getting a NaN returned on using dffits, as explained below. To me, there seems no obvious (or non-obvious reason for that matter) reason why a NaN appears. Before I start digging further, can anyone see why dffits might be failing? Is there a problem with the data? Consider: # Load data dep <-

Hi all I'm using Sweave and LaTeX, and love how they interact. But here's a different interaction I'm not sure how to achieve; I hope someone can help. I use a simple example, of course, to demonstrate. Suppose in my LaTeX document I have this: Here is a linear equation: \begin{equation} y = -1 + 3 x \label{EQ:example} \end{equation} This appears in the final product

Hi all I am doing some stats work for a group of biologists who require windows metafiles (*.wmf) for their publications. To create these, I appear to have two choices: 1. Restart my machine in Windows and use savePlot 2. Keep my machine in linux, save as another format, then convert. I'd rather stay in linux; but how do I get wmf files? I looked at using ImageMagick's convert,

Hi all After posting what follows, Duncan Murdoch suggested perhaps a bug in formatC, or an error on documentation. Any comments? In particular, bug, error or not, any ideas about how I can consistently get two significant figures to print? P. ---------- Original Message ---------- Hi all I am not a C programmer, but I am trying to understand formatC to get consistent printing of

Hi all I'm using Sweave and find it a treat. But one question: I use Sweave to create my pictures which are automatically included into LaTeX. For example, in the file test.Snw, I may have: % LaTeX stuff \begin{figure} <<fig=true,width=5,height=5>>= x1 <- seq(1,5, length=10) x2 <- sin(x1) plot(x1,x2) @ \caption{Plot} \end{figure} % More LaTeX This produces the LaTeX

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 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

Hi all I am trying to package a library in R 1.6.1 (Windoze XP). I have read the document "Writing R extensions" and think I have done things correctly (though apparently not). I have searched the mail archives for help to no avail. I also posted a few days ago and thanks to Peter Dalgaard I have managed to focus my search for errors. I am trying to create a library in R. I run the

Hello, 1. After I generate a 100x50 matrix by x3<-matrix(0,100,50);for (i in 1:100) {x1<-rpois(50, mu[i]);x2<-x1; x2[runif(50)<.01]<-0; x3[i,]<-x2}, 2. I want to calculate means and sample variances of each row and create a new matrix 100x2; 3. I then want to repeat above procedure 500 times so that eventually I will have 500 100x2 matrices. Would someone mind helping me to

Hi, Hope this does not sound too ignorant . I am trying to detrend and transform variables to achieve normality and stationarity (for time series use, namely spectral analysis). I am using the boxcox transformations. As my dataset contains zeros, I found I need to add a constant to it in order to run "boxcox". I have ran tests adding several types of constants, from .0001

Hi there, I'm using the tweedie distribution package and I cant figure out how to run a model using an identity link. I know I can use a log link by having link.power=0 and I think identity would be link.power=1, but I'm not sure. Furthermore when I try running it with link.power=1 it requires starting values which I cant manage to give appropriately so I'm not sure if its

Hi all, I am trying to package a library in R 1.6.1 (Windoze XP). I have read the document "Writing R extensions" and think I have done things correctly (though apparently not). I have searched the mail archives for help to no avail. When I try to attach the library using, eg > library( libname, lib.loc=path.to.library) I get this message: Warning message: NAs introduced by

Hi all Suppose I do the following: set.seed(1000) library(xtable) x <- runif( 10 ) y <- 1 + 2*x + rnorm( length(x) ) test.lm <- lm( y ~ x ) summary( test.lm ) xtable ( summary( test.lm ) ) The final xtable output follows: % latex table generated in R 1.8.1 by xtable 1.2-2 package % Tue Jun 22 09:56:36 2004 \begin{table}[ht] \begin{center} \begin{tabular}{rrrrr} \hline &

Hi all Two questions: 1. I note that help for the postscript device claims "The postscript produced by R is EPS (Encapsulated PostScript) compatible...". It does not say it is EPS *compliant*. Indeed, the EPS produced by R includes a \showpage command which I believe (not that I'm an expert!) non-standard EPS. My question: Is there any reason why it is EPS compatible, but

Hi all In two days, I am giving a small, informal workshop in our Department about using R. Initially, it was just for statistician, but surprisingly (to me anyhow) many mathematicians are also coming who have a MATLAB background. They are coming at the Workshop from a teaching perspective. They are considering using R to avoid licensing issues with MATLAB. One thing they were hoping me to