search for: gestating

Dear R What does NaN mean? I recently did a correlation on a batch of data for some reason it didn't like one column cor(sleep,use="complete.obs") BodyWt BrainWt SlowSleep ParaSleep TotalSleep BodyWt 1.00000000 0.95584875 -0.3936373 -0.07488845 -0.3428373 BrainWt 0.95584875 1.00000000 -0.3867947 -0.07427740 -0.3370815 SlowSleep -0.39363729

Hi, As part of the analysis of a longitudinal study, I have a list, called wg, on R 1.3.1 for windows created as follows - wg <- by(xx, xx$mrn, function(xx) lm(xx$maternal.wt ~ xx$gestation, data=xx)) xx is a data.frame, and so wg contains the results of 293 linear regressions, as lm objects. Now what I want is to extract the 293 intercepts and regression coefficients as a vector, and

We are trying to read a sas export file into R. (Works fine on Mac by ftping the export file from solaris box and then importing to R with read.xport, this gives rational numbers etc.) Trying to do this on an IBM Power p655 running linux, reading the same sas export file as used on MAC, in this case the results are not rational. Many of the numbers are interpreted as INF etc. Is it possible

Hi Stuart, I often use this small function standardize <- function(x) ( x - mean(x, na.rm=T) ) / sqrt(var(x, na.rm=T)) to standardize variables. You should be able to use this to do what you want by splitting the data frame into sections based on the factor level, using standardize() to create a new variable in each section, then paste the data frame back together. Something like: #

Well, as Andrew shouted at me again this weekend it's getting to be about that time.... :-) :-). Yep, we're thinking we should do an official 2.2 Samba release early next week, so if you have anything you can't live with, now would be a good time to tell us.... :-). It's been gestating longer than an elephant, so it's really time to let it out into the world and do its own thing :-). Feedback welcome (nay, *required* ! :-) :-). Cheers, Jeremy Allison, Samba Team. -- -------------------------------------------------------- Buying an operating system without source is li...

Hi, I would be grateful for help with a problem which is irritating me. I am quite sure that I am doing something stupid, but I can't see what it is. I am running R 1.7 on Windows 2000. The graphics device is the PC screen. The graphics from the nlme demonstration in Bates an Pinheiro's manual work just as advertised. The CO2 data and the Orthodont data dsiplay

>From owner-r-help at stat.math.ethz.ch Wed Jan 9 10:08:44 2002 > >Does anyone know how to access/list the value labels stored in a SPSS .sav >file? >At the moment I can read in variables using read.spss() (foreign package), >but have to power up SPSS to check value labels. > >Stuart There may be a solution to this in R, and I would love to hear it, but I was just faced

Hello, I'm trying to create birth centile charts from about 90000 birth weights using the package lmsqreg. As this is not a CRAN package I have problems finding any examples, which explain how to estimate the LMS parameters needed for the equivalent degrees of freedom. (I only have some basic understanding of statistics, I roughly understand the formulas, but then I don't know how to

Brian Ripley pointed out that there is already an R function scale() that does the work of my standardize(). -Greg -----Original Message----- From: Warnes, Gregory R Sent: Friday, January 18, 2002 9:38 AM To: 'Stuart Leask' Cc: 'r-help at stat.math.ethz.ch' Subject: RE: z-scores for different factor levels Hi Stuart, I often use this small function standardize <-

I'm using an R program (which I did not write) to simulate multilevel data (subjects in locations) used in power calculations. It uses lmer to fit a mixed logistic model to the simulated data based on inputs of means, variances, slopes and proportions: ? (fitmodel <- lmer(modelformula,data,family=binomial(link=logit),nAGQ=1)) where modelformula is set up in another part of the program.?

There is currently a bug in the arima function. Namely, for arima models with differencing or seasonal differencing, the innovation variance estimator uses the wrong denominator whenever xreg is non-null. This is the case, for example, when fitting an ARIMA(p,1,q) model with a drift term (common in financial applications). I reported the bug (and a fix) at

The bug repository is like an elephant: It doesn't forget, but the gestation period is long. In the present case, it is clear that something is not right, but someone needs to have sufficient recall and insight to check that your proposed fix is not unfixing a deliberate change. We should get to it eventually. (For some value of "we" not including "me"...) -pd On 20 Apr

After an extremely long gestation, a new version of the coxme package is now available on CRAN. This function fits very general mixed-effects Cox models using an lmer-like syntax. For instance here is a model that fits a random treatment effect, using data from a multi-study collaboration. fit <- coxme(Surv(pg.time, pg.stat) ~ stage + trt + (trt|site) + stata(site), data=

Dear All, We are happy to announce, after a long gestation, the release of the new version of FRAILTYPACK (version 2.2-9.5) which is now available from CRAN. The package fit general frailty models using penalized likelihood estimation, for clustered or recurrent events. For instance : -- ADDITIVE FRAILTY MODELS for proportional hazards models with two correlated random effects (intercept

After an extremely long gestation, a new version of the coxme package is now available on CRAN. This function fits very general mixed-effects Cox models using an lmer-like syntax. For instance here is a model that fits a random treatment effect, using data from a multi-study collaboration. fit <- coxme(Surv(pg.time, pg.stat) ~ stage + trt + (trt|site) + stata(site), data=

Dear All, We are happy to announce, after a long gestation, the release of the new version of FRAILTYPACK (version 2.2-9.5) which is now available from CRAN. The package fit general frailty models using penalized likelihood estimation, for clustered or recurrent events. For instance : -- ADDITIVE FRAILTY MODELS for proportional hazards models with two correlated random effects (intercept

Hi, I am very confused with constructing the wilcox.test in R. I have two populations 'original' and 'test'. I want to know if the 'test' is generally 'lower' than original. I use alpha of 0.05. So do I write the function as wilcox.test(original, test, alternative="l")? or wlcox.test(original, test, alternative = "g")? or wilcox.test(test,

Does R, or will R, be integrating the changes to the Chambers/Lucent S language under their version 5.0? If not already, then when? John Thaden Little Rock, Arkansas, USA -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info", "help", or "[un]subscribe"

I noticed that the 3.2.1 release cycle is about to start. Is there any chance that this fix will make it into the next version of R? This bug is fairly serious: getting the wrong variance estimate leads to the wrong log-likelihood and the wrong AIC, BIC etc, which can and does lead to suboptimal model selection. If it's not fixed, this issue will affect every student taking our time series

Dear All, I'am looking for examples showing that correlation does not imply causality, the targeted audience consists of undergraduate students (their first year at the university but in the BioMathStat track). All practicals are under R. I was able to extract this from R datasets: ### begin data(sunspots) data(lynx) spots <- window(sunspots, freq = 1, start = 1880, end = 1900) lnx <-