similar to: nls() and lines()

I have been working through the book "Applied longitudinal data analysis: modeling change and event occurrence" by Judith D. Singer and John B. Willett. I have been working examples using SAS and also using it as an opportunity for learning to use R for statistical analysis. I ran into some difficulties in chapter 8 which deals with using structural equation modeling. I have tried to

Hello all, In lme4 if you want to model two non-nested random effects you code it like this: mod1 <- lmer(y~x + (1|randomvar1) + (1|randomvar2)) How would you go about to model something similar in nlme? In my database I have two variables for which I have repeated measures, lets call them "individual" and "year". But none of the "individuals" were measured in

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:

Dear All, I am working with linear mixed-effects models using the lme4 package in R. I created a model with the lmer function including some main effects, a two-way interaction and a random effect. Now I am searching how I could incorporate an interaction between the random effect and one of the fixed effects. I tried to express the interaction in:

Hello all, I want to extract the survival function estimate from a model fitted by survreg(). Using predict.survreg(..., type="quantile", p=seq(0,1,0.001)), gives the quantiles, which I managed to turn around into a survival function estimate (Prob{T > t} as function of t). Is there a more straightforward way of doing this? I have had difficulties using pweibull() with the

A non-text attachment was scrubbed... Name: not available Type: text Size: 997 bytes Desc: not available Url : https://stat.ethz.ch/pipermail/r-help/attachments/19980618/ee08ba8d/attachment.pl

Hello- Does anyone have a recommendation on how to call R from perl? I'm using the IPC::Open2 module, and running R with the --slave and --quiet options. The problem is that I can't predict how many lines of output I should try to read for each command-- if any! The ultimate goal is to use perl to provide a form-driven web interface, but have R do the underlying calculations. Thanks

Using diag() to extract the diagonal of a matrix loses all but the first dimname (R 0.62.3). The problem seems to be in [ ]: > > x <- matrix(1:9,3,3) > dimnames(x) <- list(c("a", "b", "c"), c("a", "b", "c")) > x a b c a 1 4 7 b 2 5 8 c 3 6 9 > diag(x) a NA NA 1 5 9 > x[c(1,5,9)] a NA NA 1 5 9 > Paul

To those who might be interested -- following is the solution to my previous post regarding reproducing output from SAS Proc Mixed for a two-factor crossed random effects ANOVA model. I am graciously endebted to the kind replys from two statisticians for this solution whose names I will refrain from mentioning for the sake of privacy. I hope this helps someone?! -- Phil Turk > hw7 <-

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

Content Inspecion SMTPMAIL could not deliver the e-mail below because of unreachable host Please check the recipients e-mail address before you try again. Received: from Unknown (61.113.174.142) by ETRUST-SMTP (10.3.0.22) From: r-bugs@r-project.org To: aikins@pixie.udw.ac.za Subject: Mail Delivery System (aikins@pixie.udw.ac.za) Date: Thu, 1 Apr 2004 00:56:02 -0800 MIME-Version: 1.0 Content-Type:

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

Hi, there, S+ has a function leaps which can be used for criterion based model selection. Is there an equivalent function in R? thanks in advance. Kenny Ye Assistant Professor Department of Applied Math and Statistics SUNY at Stony Brook Stony Brook, New York 11794-3600 (516)632 9344 -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read

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

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