similar to: a very preliminary lme for R

Sorry, by mistake I sent this to Professor Bates instead of r-help. Han -------- Original Message -------- Subject: Re: [R] degrees of freedom for t-tests in lme Date: Thu, 25 Apr 2002 09:16:16 -0700 From: Han-Lin Lai <Han-Lin.Lai at noaa.gov> To: Douglas Bates <bates at stat.wisc.edu> References: <3CC6E87F.5400277D at noaa.gov> <6rg01lottu.fsf at franz.stat.wisc.edu>

Dear R-users, I am a newbie to this site and a relative new-comer to S/R, so please tread lightly, for you tread... There have been several posting relating to problems with augPred() from the nlme library. Here is a "fix" for one of these problems which may lie at the root of others. In my case the problem with augPred() lay in gsummary(), which augPred() uses, causing it to fail.

Full_Name: Volker Franz Version: 1.5.1 OS: Debian-Linux Submission from: (NULL) (134.176.77.64) Hi, working with the panel.superpose function, I found out that this function does not pass the subscripts and groups arguments to panel.groups functions. In my view, this seems an unnecessary restriction, because the subscripts-mechanism which allows to access the original data should also work if

Could someone give me a hand with the format of the gsummary function? Basically, I have a large set of xyz coordinates generated by LiDAR data (>37 million points) and I am trying to derive various summary statistics on the z-coordinates by a grid cell. I wrote a function to do this by creating factors from the x- and y- coordinates and then using gsummary. However, I want the function to

Full_Name: Ben Saville Version: 2.1 OS: Windows XP Submission from: (NULL) (152.2.94.145) I'm using the gsummary function to calculate a sum of V1 (column one) from my data 'mytest' by group (V2,or column 2). If V1 (the variable of interest) is all the same value (in this case all 2's), I do not get back the correct summation. If there is at least one difference in V1 (all

Dear R list members, I've turned up a strange discrepancy between results obtained from the lme function in the nlme package in R and results obtained with lme in S-PLUS. I'm using version 3.1-24 of nlme in R 1.4.1 under Windows 2000, and both S-PLUS 2000 and 6.0, again under Windows 2000. I've noticed discrepancies in a couple of instances. Here's one, using data from Bryk

--Hipsgkxbeg Content-Type: text/plain; charset=us-ascii Content-Description: message body text Content-Transfer-Encoding: 7bit This one was sent privately to me. --Hipsgkxbeg Content-Type: message/rfc822 Content-Description: forwarded message Content-Transfer-Encoding: 7bit Received: from tuvok.kom.tuwien.ac.at (tuvok.kom.tuwien.ac.at [192.35.241.66]) by fangorn.ci.tuwien.ac.at (8.9.3/8.8.5)

Martin Maechler <maechler@stat.math.ethz.ch> writes: > This caused 3 postings (1 to R-help, 2 to R-devel) to be aborted. > For some reason, it also seems majordomo did not save the postings in > a 'dead.letter' or something. > So could you please re-submit them? > These are the sendmail entries: > Sep 17 21:14:21 .. sendmail[9364]: VAA09363:

On Wed, Jul 8, 2009 at 12:16 PM, Eric Christopher<echristo at apple.com> wrote: >> For reference, here is the corresponding patch against GDB: >> http://web.mit.edu/rnk/www/jit-patch.diff >> > > This appears to be missing some hunks... Ah, crap, there's a gdb/jit.[ch] but I don't know how to get cvs diff to include them. 'cvs add' wants write access

torsten@www:/raid/home/torsten > ping rsync.r-project.org PING franz.stat.wisc.edu (128.105.174.95): 56 data bytes --- franz.stat.wisc.edu ping statistics --- 12 packets transmitted, 0 packets received, 100% packet loss Torsten -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-devel mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send

Hi all! I am working on a nested lme model with one fixed effect ("treatment", which 3 levels) and two random effects for "Individuals" (four of them) within "treatment" and "replicate -2 levels-" within "individual" within "treatment". For doing so, I´ve been trying to create a factor for Individual%in%Treatment, say IT by

This may have been announced before but there are now both Debian GNU/Linux and RedHat Linux packages for R-0.61 on the CRAN archive sites. The master site for CRAN is http://www.ci.tuwien.ac.at/R/ That site provide a list of available mirror sites. Within the CRAN archive sites the pre-compiled packages for Linux are located in the bin/i386-linux directory. Under RedHat you can install

Dear Douglas Bates. I just downloaded the compiled version (I'm a poor Windows devil, not yet having found the time to move to a more advanced platform...) from NT- the files are dated 30.5-1999 so they are not old - and the problem persisted....wonder what I did wrong ? R : Copyright 1999, The R Development Core Team Version 0.64.0 Patched (May 3, 1999) R is free software and comes with

Dear listers, I am trying to fit a model using nlsList() using alternately a SSfol() selfstart function or its developped equivalent formulae. This preliminary trial works well mydata<-groupedData(Conc~Tps|Organ,data=mydata) mymod1<-nls(Conc~SSfol(Dose,Tps,lKe,lKa,lCl),data=mydata) as well as a developped form: mymod2<-nls(Conc~Dose * exp(lKe+lKa-lCl) *

Right now, GDB has no way to be told about JITed code. I'm working on adding such an interface, and the LLVM JIT would be the first client. The interface is evolving, and I'll admit that right now it's downright terrible. Here's how it works right now: - The JIT generates the machine code and DWARF call frame info (.eh_frame/.debug_frame) for a function into memory. - The JIT

Hello, I want to use lme to fit two (or more) models, and then compare the fits on each individual. I know how to write my own code to do this (for each individual, plot the raw data, followed by lines() to plot each fitted curve) but I would like to use plot(augPred(... as it produces a nice trellis plot. I thought I could do this with par(new=T) but it does not seem to work.

(please keep sebastian cc'd) Doug - we are managing to trigger, both on R 1.1.1 and R 1.2.0 (proper versions of NLME), the following: > warnings(model1) Warning message: Singular precision matrix in level -1, block 1 Error in cat(list(...), file, sep, fill, labels, append) : argument 2 not yet handled by cat > suggesting either an "inappropriate" use of cat(), or a

I am trying to install a library for my linear models class. (The name "st849" comes from the course identification, Statistics 849.) # R INSTALL st849 Installing package `st849' ... funs data help >>> Building help index for package `st849' ... done. >>> Building help pages for package `st849' NONE:0: m4: ERROR: EOF in argument list RWC

If I used a groupedData object, if I do fit<-lme(blah) then plot(augPred(fit)) produces a nice trellis plot of the data along with the fitted lines However I find that the lines and the data points are in the same colour (light blue against a medium grey background). Is there a way to make the lines in a different colour (e.g. black)? It would also be nice if the line were plotted after the

With the assistance of Saikat DebRoy we have produced an nlme package (linear and nonlinear mixed-effects models) for R (see http://www.r-project.org/ for details about R). The nlme package for R requires R-0.90.0, which was released earlier today. Also released today was NLME 3.2 for S-PLUS. We will be working on incorporating the changes from NLME 3.1 to 3.2 into the R version but we wanted