similar to: R listing on SAL

I am embarrassed to have to ask this but can anyone tell me of a Tukey multiple comparisons procedure available for R? I have looked through the search page, through the FAQ, and in the index of V&R (1999), and I still can't find such a thing. I see there is a ptukey function and a qtukey function but that is as far as I got. Do I need to roll my own? -- Douglas Bates

Is there a recommended binary version of perl to use under Windows XP for `Rcmd build'? The R FAQ for Windows just mentions "perl". I think that Kevin Wang's document on R for Windows suggests ActiveState perl but I have not been able to install it. I am running Windows XP Professional on a laptop. The operating system was originally Windows 98 then upgraded to Windows ME

Zhang, Xiantao wrote: >>From 5c70c038c57190144390ae9d30c3d06afba103d4 Mon Sep 17 00:00:00 2001 > From: Xiantao Zhang <xiantao.zhang at intel.com> > Date: Tue, 1 Apr 2008 14:59:30 +0800 > Subject: [PATCH] KVM:IA64 : Add kvm sal/pal virtulization support. > > Some sal/pal calls would be traped to kvm for virtulization > from guest firmware. > Signed-off-by: Xiantao

Hello! I am running openSSH 2.9x on an IRIX 6.5.x platform. This was recently installed using SGI-supplied "freeware" binaries. I find that as time goes on, it takes more attempts to establish an ssh connection from the IRIX platform to another machine, as it fails with "not enough entropy in PRNG." I posted a note asking for assistance, and received a reply suggesting I

Hi, df <- read.table("data.txt", header=T); library(nls); fm <- nls(y ~ a*(x+d)^(-b), df, start=list(a=max(df->y,na.rm=T)/2,b=1,d=0)); I was using the following routine which was giving Singular Gradient, Error in numericDeriv(form[[3]], names(ind), env) : Missing value or an Infinity produced when evaluating the model errors. I also tried the

Would it be possible to include a new command in the ${RHOME}/etc path, which updates the Library index (Rd and html) via "R LIBINDEX". I need it because I use RPM to manage three different R installations (at home, at the institute and in the seminar rooms). I put all libraries into several packages, so it is easy for me to update a single library at all these different places. But

I have installed the binary packages for Debian GNU/Linux release 2.3 (woody) in the U.S. mirror of the CRAN archive. They should propagate to the main CRAN archive within a day and to the other mirrors within two days. These packages have been compiled with gcc-3.0 and g77-3.0. I believe the testing distribution currently provides only a snapshot of gcc-3.0, not the latest released version, so

I have installed the binary packages for Debian GNU/Linux release 2.3 (woody) in the U.S. mirror of the CRAN archive. They should propagate to the main CRAN archive within a day and to the other mirrors within two days. These packages have been compiled with gcc-3.0 and g77-3.0. I believe the testing distribution currently provides only a snapshot of gcc-3.0, not the latest released version, so

I want to create a grid of function values for use in `contour' or `persp'. The function is the log-likelihood for the gamma. The sample is stored as vector of length 20 called `Survival'. A single evaluation of the log-likelihood at, say, scale = 9 and shape = 10 would be obtained by sum(dgamma(Survival, scale = 9, shape = 10, log = TRUE)) (This may work only 0.99.0, I'm not

Windows 2000 (v.5.00.2195), R 1.5.1 I have an lme object, fm0, which I examine with anova() and intervals(). The anova output indicates one of the interaction terms is significant, but the intervals output shows that the single parameter for that term includes 0.0 in its 95% CI. I believe that the anova is a conditional (sequential) test; is intervals marginal or approximate? Which should I

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, Given a LME model (following the notation of Pinheiro and Bates 2000) y_i = X_i*beta + Z_i*b_i + e_i, is it possible to extract the variance-covariance matrix for the estimated beta_i hat and b_i hat from the lme fitted object? The reason for needing this is because I want to have interval prediction on the predicted values (at level = 0:1). The "predict.lme" seems to

Dear list, i am a student of psychology and have to do a multilevelanalysis on some data. About that i have one general and one specific question. This is what i have copied from the help-file on lme: data(bdf) fm <- lme(langPOST ~ IQ.ver.cen + avg.IQ.ver.cen, data = bdf, random = ~ IQ.ver.cen | schoolNR) summary(fm) after summary(fm) i get the following error:

Full_Name: W.B.Kloke Version: 1.8.0 OS: FreeBSD-4.7 Submission from: (NULL) (195.253.22.63) I found a parameter evaluation inconsistency in NLME package. I tried to use gls() inside a function, and I wanted use gls() for different subsets of a data frame: prgls <- function(name){ gls( log10(Y)~(cond-1)+(cond-1):t ,pr,subset=subject==name)} Applying this function with a string as parameter

R-Listers: A very minor -- and maybe silly -- question just for personal enlightenment. In S (either R or S-Plus, AFAIK) when one types or pastes a trellis graphics command into the commands/console window, the graph is automatically produced: e.g., trellis.device(...) xyplot(y~x) If one puts these in a function and calls the function, the same occurs. However, if one sources in these command

I am the maintainer of one of the CRAN mirrors in the U.S., cran.us.r-project.org. It appears that this site has been inaccessible since Friday. This has me puzzled because I can connect to the same machine through ssh and I can see that the http daemon (apache) is running. In fact I can make an http connection from this machine to this machine and from other machines on any of our subnets to

>>>>> "DougB" == Douglas Bates <bates@stat.wisc.edu> writes: DougB> ....... DougB> ....... { time comparisons in testing lme(..) for R } DougB> ....... DougB> This can be expected to run faster when a version of nlm that accepts DougB> gradients and Hessians is available. --- Doug, do I understand properly that it won't

//Referring to the response posted many years ago, copied below, what is the specific criterium used for singularity of the gradient matrix? Is a Singular Value Decomposition used to determine the singular values? Is it the gradient matrix condition number or some other criterion for determining singularity? // //Glenn // / / /> What does the error 'singular gradient' mean

Hello I'm reading up on fitting truncated Weibull distribution to data. There are posts in 2002 that point to this presentation by Prof Bates: http://www.stat.wisc.edu/~bates/JSM2001.pdf but now the file is not there. I can't find it anywhere else, Google doesn't have a cached copy for it. Could someone please give me a copy of this file, if they have it? Thanks and regards,

Hello, I am interested in using nlme to model repeated measurements, but I don't seem to get good CIs. With the code below I tried to generate data sets according to the model given by equations (1.4) and (1.5) on pages 7 and 8 of Pinheiro and Bates 2000 (having chosen values for beta, sigma.b and sigma similar to those estimated in the text). For each data set I used lme() to fit a model,