similar to: lme

> Date: Fri, 7 Apr 2000 09:55:08 -0400 (EDT) > From: Faheem Mitha <faheem at email.unc.edu> [I have given a more meaningful subject line.] > On Fri, 7 Apr 2000, Prof Brian D Ripley wrote: > > > On Thu, 6 Apr 2000, Faheem Mitha wrote: > > > > > I'm running splus 5 on a solaris platform remotely, and running R on linux > > > on my home machine.

Hello, Is there any way we convert SAS dataset into R dataset? Kindest Regards, Peppy Adi-Purnomo ------ Peppy Adi-Purnomo Energy Market Analyst Energy Link Ltd Dunedin - New Zealand Ph.: +64 3 479 2475 Fax: +64 3 477 8424 Email: s.adi.purnomo at energylink.co.nz www.EnergyLink.co.nz -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read

Hello, I am wondering if someone can direct me to the syntax of the Generalised Additive Model in R? In Splus this is gam(formula, ... inputs would be pretty much the same as glm I suspect..) Any response on that (or what package should I grap) would be appreciated very much. Thanks, Peppy Adi-Purnomo -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help

Hello again, If I want to have a set of add-ons functions (say: nlme and tseries) loaded to R everytime I run R, do I put command: library(nlme) library(tseries) in the rw1021/etc/Rprofile file? thanks, Peppy -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info",

Hello, I am wanting to download the recent R update off the web. May I assume confidently that R package from the website is virus free? I apologise if this is a very silly question, I just need to convince myself if the site is very hygiene. Thank you very much, Peppy -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read

Hi, Suppose I have a by b matrix which a and b are big numbers. How can I convert this matrix into a vector with length of a*b ? By column-wise or row-wise doesnt matter. My R version is 0.64.2. Kindest Regards, Peppy Adi-Purnomo ------- Stephanus (Peppy) Adi-Purnomo Energy Market Analyst - Energy Link Ltd - Dunedin email: s.adi.purnomo at energylink.co.nz http://www.energylink.co.nz

Dear Help desk, What is the difference between glm and lm in R? I am so sorry the help pages dont really help me much since they dont have description how R handles the two. Thanks Peppy -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info", "help", or

hello, I have a very old R version which calculates anova table for a stream of outcomes say: 1,5,4,3,2,4,5,2 which are grouped into 4 blocks, say A,B,C,D, of 2 say namely X,Y What is the command in R 0.64.2 to do this? My very old R command is anova.table but this command is unkown in R 0.64.2 Thanks, Peppy -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help

Hello, I want to update my new R version. Now something has change is that I am now in Windows 2000. Which one should I use whether the R for window 9x or Win NT? Thanks a lot, Peppy -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info", "help", or

Hi, > Say we have a large-sized vector of observation A and a time vector with > the same size. The latter vector basically is the time flag when ith > obervation of A occurs, hence A[i] occurs at T[i]. > > I have no trouble with doing: T[A <= 10 | A >=20] for example > > Now, what I am confused: I tried to pull out when the largest value of A > with: > >

Dear Help desk, I have a model that I believe to be non linear. The relationship is something like: Response = Var1 + Var1^2 + Var1^3 Can I examine this relationship using glm or lm? However, glm or lm gives out pnly 1 coeficient for the model which is on Var1, unless I created some dummy var like (Var2 = Var1^2 and Var3 = Var1^3). Is there any other function handling this relationship?

Dear Help desk, I have a model that I believe to be non linear. The relationship is something like: Response = Var1 + Var1^2 + Var1^3 Can I examine this relationship using glm or lm? However, glm or lm gives out pnly 1 coeficient for the model which is on Var1, unless I created some dummy var like (Var2 = Var1^2 and Var3 = Var1^3). Is there any other function handling this relationship?

Greetings, 1) Is there a nice way of extracting the variance estimates from an lme fit? They don't seem to be part of the lme object. 2) In a series of simulations, I am finding that with ML fitting one of my random effect variances is sometimes being estimated as essentially zero with massive CI instead of the finite value it should have, whilst using REML I get the expected value. I guess

tl;dr anova.lme() claims to provide sums of squares, but it doesn't. And some names are misspelled in ?lme. I can submit all this stuff as a bug report if that's preferred. ?anova.lme says: When only one fitted model object is present, a data frame with the sums of squares, numerator degrees of freedom, denominator degrees of freedom, F-values, and P-values The output of fm1

I am using rw0641. In my continuing quest to understand repeated measures analysis, I again return to lme. I exported the Potthoff and Roy data Orthodont.dat from S-PLUS 4.5 to avoid capture errors and ran the examples in the R help. I imported the data.frame with data <- read.table("Orthodont.dat",header=T) attach(data) and created the objects Orthodont.fit1 <-

When I try to run the example from Variogram with an lme object, I get an error (although summary works): R : Copyright 2005, The R Foundation for Statistical Computing Version 2.2.1 (2005-12-20 r36812) ISBN 3-900051-07-0 ... > fm1 <- lme(weight ~ Time * Diet, BodyWeight, ~ Time | Rat) Error: couldn't find function "lme" > Variogram(fm1, form = ~ Time | Rat, nint =

Dear all, Edmond Ng (http://multilevel.ioe.ac.uk/softrev/reviewsplus.pdf) provides an example to fit the mixed effects meta-analysis in Splus 6.2. The syntax is: lme(fixed=d~wks, data=meta, random=~1|study, weights=varFixed(~Vofd), control=lmeControl(sigma=1)) where d is the effect size, study is the study number, Vofd is the variance of the effect size and meta is the data frame.

Dear R developers, One function in my mkin package [1] returns an object that is originally created by nlme(), but contains some additional information. Its class is c("mmkin.nlme", "nlme", "lme"). Now I would like to use the anova() method for lme objects for comparing such S3 objects. Unfortunately, anova.lme currently does not check for inheritance, but

I'm puzzled by the following problem, which appears when attempting to run an analysis on part of a dataset: If I try: csubset <- dat$Diagnosis==0 cont <- lme(fixed=cform, random = ~1|StudyName, data=dat,subset=csubset,na.action=na.omit) Then I get: Error in eval(expr, envir, enclos) : Object "csubset" not found But if I do

Dear people, A somewhat clueless question follows: I just discovered that the lme function in contrib package nlme for R, while similar to the lme function in Splus, does not use the cluster function option. This difference does not appear to be documented in the V&R `R Complements' file. I have data which is divided into 6 groups The lme model is of the form (simplified from the actual