similar to: Solution to "Can anyone help me reproduce this SAS Mixed output??"

I asked this before and I am going to try again in more applied terms. I am trying to use R to extract variance components for a two-factor random effects model with both factors crossed. It would also be nice to generate some confidence intervals as well. For example, a data set using SAS Proc Mixed is below followed by the four variance component estimates and the respective confidence

Hello, I work with a mixed model with 4 predictor variables Time, Size, Charge, Density and Size, Charge, Density are factors, all with two levels. Hence I want to put their interactions with Time into the model. But, I have two data sets (Replication 1 and 2) and I want that Replication is random effect. Here is my code: knots <- default.knots(Time) z <- outer(Time, knots, "-")

Following lme model runs fine in general under R.2.1.1 but only for 9 out of my 11 response variables under R.2.2.0. model for one of my response variables: lme(Yresp~F1fix,random=list(const=pdBlocked(list(~F2mix-1,~Ass:F1fix-1,~F3mix-1,~F1fix:F3mix-1,~F2mix:F3mix-1),pdClass="pdIdent"))) Yresp is my response variable, F1fix is a fixed effect factor whereas F2mix and F3mix are random

Does anyone have any R code or suggestions on how to generate a population that follows a Neyman-Scott process? More specifically, I'd like to randomly generate both the parents and the children from Poisson distributions and the children locations from a bivariate normal onto an N by N grid. If the locations spill over outside the grid, then I'd like to reflect those locations back into

I am trying to use lme to fit a mixed effects model to get the same results as when using the following SAS code: proc mixed; class refseqid probeid probeno end; model expression=end logpgc / ddfm=satterth; random probeno probeid / subject=refseqid type=cs; lsmeans end / diff cl; run; There are 3 genes (refseqid) which is the large grouping factor, with 2 probeids nested within each refseqid,

All - I'm having an issue with trying to plot a model derived from nls() onto a simple plot.? I have included a sample data set and the code that I've been using. ?? year month day?????? date location mileage? cost gallon????? cpg ? mpg????????? x 2009???? 1?? 4?? 1/4/2009????? BZN? 124585 19.39? 14.37 1.349339 10.71677 2009-01-04 2009???? 1? 15? 1/15/2009????? BZN? 124888? 23.2? 16.12

I am trying to use lme to fit a mixed effects model to get the same results as when using the following SAS code: proc mixed; class refseqid probeid probeno end; model expression=end logpgc / ddfm=satterth; random probeno probeid / subject=refseqid type=cs; lsmeans end / diff cl; run; There are 3 genes (refseqid) which is the large grouping factor, with 2 probeids nested within each refseqid,

Dear Rhelp folks, I think I read this in recent Rhelp, but I can't find it in the archive. I need to know how to get lg2c in order to build the nlme package. I'm trying to install nlme_3.1-2.tar.gz on a Linux machine (6.1 Redhat) under R Version 0.99.0 (February 7, 2000). Thanks, Jim Jim Robison-Cox ____________ Department of Math Sciences | |

Just after R-0.63.0 came out I reported a problem to the list. I was unable to get ./configure to complete it's job, because it couldn't find a working fortran. It stopped after trying to test for underscore in subroutine names. This was especially perplexing, since I could still configure the older version (I had 62.2) with no error messages. I tried suggestions of B.D. Ripley and

Dear R-helpers, I just upgraded from version 62.1 to 62.2 and now am faced with the same problem I've run into before when upgrading: Is there an easy way to take the libraries I have been using in the older version and pull then into the new version. My organization looks like this on Solaris 2.5.2: /opt -- /R -- /R-0.62.1 --/library -- subdirectories for lots of

To R-devel: If anyone is testing the aov() function and it's relatives which I posted to the list, thanks for the effort, but I would ask you to hold off now. Brian Ripley is developing better and more complete versions, so you should save your effort for looking over his functions. His preliminary version is located at: http://www.stats.ox.ac.uk/pub/R/aov.tar.gz Prof. Ripley says

When I run configure for the new release, R-0.63, I get an error: loading cache ./config.cache checking for a BSD compatible install... aux/install-sh -c checking whether ln -s works... yes checking for ranlib... ranlib checking for bison... bison -y checking for ar... ar checking for ratfor... no checking for latex... /psoft/TeX/bin/latex checking for dvips... /psoft/TeX/bin/dvips checking for

When I run configure for the new release, R-0.63, I get an error: loading cache ./config.cache checking for a BSD compatible install... aux/install-sh -c checking whether ln -s works... yes checking for ranlib... ranlib checking for bison... bison -y checking for ar... ar checking for ratfor... no checking for latex... /psoft/TeX/bin/latex checking for dvips... /psoft/TeX/bin/dvips checking for

This message is in MIME format. The first part should be readable text, while the remaining parts are likely unreadable without MIME-aware tools. Send mail to mime@docserver.cac.washington.edu for more info. --0-488534781-901753860=:6003 Content-Type: TEXT/PLAIN; CHARSET=US-ASCII Content-ID: <Pine.NEB.4.00.9807291716501.6003@mathfs.math.montana.edu> Attached is R code for an aov

Dear R-help-ers (Currently running Version 0.65.1 Release (October 07, 1999) on Solaris 2.6) My problem has to do with updates of the R language. When I install a new version, I would like to just change one directive, or a soft link to the current doc/html/ directory, to update a bookmark I have pointing to the current R doumentation directory. I've tried building a soft link:

R-crew, I'm trying to install the R-Snapshot which I downloaded a few hours ago and I'm having a problem I can't trace. My system is Sun SPARC, SunOS 5.3 : 1996 >From running ./configure:*************************************** loading cache ./config.cache checking for a BSD compatible install... /opt/local/R/R-snapshot/etc/install-sh -c checking whether ln -s works... yes

Dear R-help folks, I am seeing unexpected behaviour from the function mean with option na.rm =TRUE (which is removing a whole column of a data frame or matrix. example: testcase <- data.frame( x = 1:3, y = rep(NA,3)) mean(testcase[,1], na.rm=TRUE) [1] 2 mean(testcase[,2], na.rm = TRUE) [1] NaN OK, so far that seems sensible. Now I'd like to compute both means at once:

R developers, Last week I asked several questions which I think I can now answer myself in large part. > > Q 1) Is aov supposed to avoid qr decomposition and matrix inversion? > In S it is claimed to be faster than lm() for large datasets. Is > that due to avoidance of qr()? > (My main goal is to get the Error strata working.) My use of matrix

I am trying calculate a probability using numerical integration. The first program I ran spit out an answer in a very short time. The program is below: ## START PROGRAM trial <- function(input) { pmvnorm(lower = c(0,0), upper = c(2, 2), mean = input, sigma = matrix(c(.1, 0, 0, .1), nrow = 2, ncol = 2, byrow = FALSE)) } require(mvtnorm) require(adapt) bottomB <- -5*sqrt(.1) topB <-

R developers, 2 things: a bug in model.frame and a question about setup of generic functions. I don't understand the following behavior for generic functions: Suppose I'm working with the cats data in the MASS library and I want to create a formula object to model Hwt on Sex: # This works: > formula(Hwt ~ Sex) Hwt ~ Sex # But the following does not: > formula(Hwt ~ Sex,