search for: maustin

Hi, Is it possible to add labels to variables in R (so as to have a better description of what the variables represent)? Thanks, Neil

I have generated some some survival times and censoring indicators. Thus I have an ordered pair for each observation. How do I sample these ordered paris? I only know how to sample from a vector? I would appreciate any help I could get. Thanks Matt

...mailto:r-help-bounces at stat.math.ethz.ch]On Behalf Of Kunal Shetty Sent: Friday, October 08, 2004 5:16 AM To: Austin, Matt Cc: R-help Subject: RE: [R] Read.Table Reading a Text file thanks austin, it worked..it was exactly what I was looking for regards Kunal "Austin, Matt" <maustin at amgen.com> wrote: > x.1$V1 > or > x.1[,1] > or > x.1['V1'] > > and you shouldn't need to call print.default() directly, just call > print(). > > --Matt > -----Original Message----- > From: r-help-bounces at stat.math.ethz.ch > [mailto:r-hel...

I have installed the binary packages for Debian GNU/Linux release 2.3 (woody) in the U.S. mirror of the CRAN archive. They should propagate to the main CRAN archive within a day and to the other mirrors within two days. These packages have been compiled with gcc-3.0 and g77-3.0. I believe the testing distribution currently provides only a snapshot of gcc-3.0, not the latest released version, so

I have installed the binary packages for Debian GNU/Linux release 2.3 (woody) in the U.S. mirror of the CRAN archive. They should propagate to the main CRAN archive within a day and to the other mirrors within two days. These packages have been compiled with gcc-3.0 and g77-3.0. I believe the testing distribution currently provides only a snapshot of gcc-3.0, not the latest released version, so

Dear R users and Helpers I am beginner with using R and interested in carrying out certain task for my statistical research. I am reading data for a text file, which could contain data in following pattern x y 8 10 11 14 16 16 18 15 6 20 4 4 20 18 As per the example I have two columns and 7 rows of data in each. However is real life data situation I may not know how many columns are present

Greetings, I have a meta-analysis problem in which I have fixed effects regression coefficients (and estimated standard errors) from identical models fit to different data sets. I would like to use these results to create pooled estimated regression coefficients and estimated standard errors for these pooled coefficients. In particular, I would like to estimate the model \beta_{i} = \mu +

Dear R users; Ive got two questions concerning nlme library 3.1-65 (running on R 2.2.0 / Win XP Pro). The first one is related to augPred function. Ive been working with a nonlinear mixed model with no problems so far. However, when the parameters of the model are specified in terms of some other covariates, say treatment (i.e. phi1~trt1+trt2, etc) the augPred function give me the following

...ot;calibration" set from the original dataframe in order to get my "validation" set.....Any hint will be greatly appreciated. >TT > > > >------------------------------ > >Message: 20 >Date: Sat, 14 Aug 2004 17:41:01 -0700 >From: "Austin, Matt" <maustin at amgen.com> >Subject: RE: [R] calibration/validation sets >To: "'Peyuco Porras Porras .'" <levin001 at 123mail.cl>, > R-help at stat.math.ethz.ch >Message-ID: > <E7D5AB4811D20B489622AABA9C53859101F1111D at teal-exch.amgen.com> >Content-Type: text...

Lattice graphics may be the answer depending on your exact problem. Here is an example of four traditional plots without space: par(mfrow=c(2,2), omi=c(.5, .5, .5, .5)) par(mar=c(0, 2, 2, 0)) plot(rnorm(10), rnorm(10), axes=FALSE) box(); axis(2); axis(3) par(mar=c(0, 0, 2, 2)) plot(rnorm(10), rnorm(10), axes=FALSE) box(); axis(3); axis(4) par(mar=c(2, 2, 0, 0)) plot(rnorm(10),

If you used all 28 animals to find the model out of a group of candidate models, I would have my reservations about this 'validation'. Any confidence intervals you get from the final model are bound to be overly optimistic because you haven't accounted for the degrees of freedom chewed up during the model fitting/finding process. Matt Austin Statistician Amgen One Amgen Center

The following works, you need to include x=TRUE in the call to coxph. Passing the time and status variables as additional arguments is a matter of personal preference. f.coxph.zph<-function(x, timeVar, statusVar) { cox.fit <- coxph(Surv(timeVar, statusVar) ~ x, na.action = na.exclude, method = "breslow", x=TRUE) fit.zph<-cox.zph(cox.fit) fit.zph$table[,3] } time.cox <-

Try factor(vec2) in your bwplot() call. -----Original Message----- From: r-help-bounces at stat.math.ethz.ch [mailto:r-help-bounces at stat.math.ethz.ch]On Behalf Of Ernesto Jardim Sent: Friday, July 09, 2004 9:41 AM To: Mailing List R Subject: [R] Problem with bwplot Hi, I'm ploting some box-and-whisker plots with bwplot but I'm not getting any box-and-whiskers ... just dots.

When plotting using the cex.axis argument to boxplot(), the size of the plotting symbols beyond the whiskers of the boxplot are being changes. Example: par(mfrow=c(2,1)) boxplot(c(rnorm(10), 10), horizontal=TRUE, main="Test", las=2, cex.axis=2) boxplot(c(rnorm(10), 10), horizontal=TRUE, main="Test", las=2, cex.axis=1) This is from the following line in bxp()

I don't get the same result, do you have a package loaded that would change the default behavior (such as Hmisc)? > list(grp.1) [[1]] [1] 1 2 Levels: 1 2 > list(grp.1[mask]) [[1]] [1] 1 Levels: 1 2 > library(Hmisc) <<snip>> > list(grp.1[mask]) [[1]] [1] 1 Levels: 1 --Matt > version _ platform i386-pc-mingw32 arch i386 os

A correction. You either need to open the plotting device prior to the simulation that includes your plotting commands and close it after the simulation or have the name change dynamically in your simulation so that the runs go in separate files. An example of the first method would be postscript(file="/where/to/put/file/filename.ps") ## your simulation commands dev.off() for the

I read in a thread in r-help today that the windows device in 2.6 supports transparency, so I tried an example and had some issues. The density plots should be filled with transparent color in the following example (similar to the points), however the color is "fully" transparent. This works in the Cairo device, but not in the windows device. Thanks, --Matt Matt Austin

How about: x <- data.frame(matrix(rnorm(1550),c(50,31))) model <- step(lm(x[,1] ~ as.matrix(x[,2:31]))) --Matt -----Original Message----- From: r-help-bounces at stat.math.ethz.ch [mailto:r-help-bounces at stat.math.ethz.ch]On Behalf Of rongguiwong Sent: Monday, September 20, 2004 20:52 PM To: r-help at stat.math.ethz.ch Subject: [R] how to take this experiment with R? This message uses

Could you just use lines(newX, myPred, col=2) -----Original Message----- From: r-help-bounces at stat.math.ethz.ch [mailto:r-help-bounces at stat.math.ethz.ch]On Behalf Of Paul Johnson Sent: Thursday, September 23, 2004 10:3 AM To: r help Subject: followup: Re: [R] Issue with predict() for glm models I have a follow up question that fits with this thread. Can you force an overlaid plot