search for: amgen

...iple columns: Value Value Value GeneID Method1 Method2 Method3 6 123 456 987 7 234 NA 432 8 864 190 34 Is there an elegant way to do this type of reformating in R? Thanks. Mark -- Mark Dalphin email: mdalphin at amgen.com Mail Stop: 29-2-A phone: +1-805-447-4951 (work) One Amgen Center Drive +1-805-375-0680 (home) Thousand Oaks, CA 91320 fax: +1-805-499-9955 (work)

...power = .8913665 coef(model3$modelStruct$varStruct)["const"]; #const = -0.4229219 coef.varFunc(model3$modelStruct$varStruct); #R can't seem to find this function, Splus can Any advice on what I am doing wrong would be appreciated. Hugh Rand Senior Scientist Amgen rand at amgen.com

...#39;, 'Marron' and 'scale-space' didn't return anything that I recognized as useful here, but I am new to the idea of confidence in a lowess line so I may not be using the appropriate vocabulary. Thank-you, Mark ---------------------- Mark Dalphin Dept Comp Biol, M/S AW2/D3262 Amgen, Inc. 1201 Amgen Court W Seattle, WA 98119 Phone: +1-206-265-7951

...t;- pcycle(pars$outcex, pars$cex, par("cex")) Where pars$cex is matching pars$cex.axis. This would also effect other parameter specifications. This is a behavior I did not expect; I'm not sure if this is a design feature or something else. --Matt Matt Austin Statistician Amgen, Inc 800 9AMGEN9 x77431 805-447-7431 "The purpose of computing is insight, not numbers" --Richard Hamming [[alternative HTML version deleted]]

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

...ation, despite being under "Essential Programs" section is described in a manner suggesting that documentation is optional. http://cran.r-project.org/doc/manuals/R-admin.html#Essential-and-useful-othe r-programs-under-Unix ---------------------- Mark Dalphin Dept Comp Biol, M/S AW2/D3262 Amgen, Inc. 1201 Amgen Court W Seattle, WA 98119 Phone: +1-206-265-7951 [[alternative HTML version deleted]]

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

...racle date format remaining static; I prefer to code more 'defensively'. I assume this is a bug in DBI/ROracle processing of functions in the selected variables, but I don't know. Any suggestions are welcome. Thanks, Mark -- Mark Dalphin email: mdalphin at amgen.com Mail Stop: 29-2-A phone: +1-805-447-4951 (work) One Amgen Center Drive +1-805-375-0680 (home) Thousand Oaks, CA 91320 fax: +1-805-499-9955 (work)

Hello there, I am doing analysis on survival data. How do I pick out a probability of survival at a chosen landmark time(for example, 3 years, 4 years) from the result of "survfit"? or any other functions? As I know, the result of 'survfit' only have the probabilities for all event or cencor time. Thank you very much Lisa Wang Princess Margaret Hospital Toronto, Ca tel:

...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 Drive M/S 24-2-C Thousand Oaks CA 93021 (805) 447 - 7431 "Today has the fatigue of a Friday and the desperation of a Monday" -- S. Pearce 2005 -----Original Message----- From: r-help-bounces at stat.math.ethz.ch [mailto:r-help-bounces at stat.math.ethz.ch]On Behalf O...

Hello How do I use function 'MERGE" to combine the FILE A and FILE B below to make FILE C? Thank you FILE A 140 151 167 30.1 11.4 40 FILE B 140 167 5.7 30.3 FILE C 140 151 167 30.1 11.4 40 5.7 NA 30.3

...look like # Do the splom. Notice that the lower panel should have reported that # it was seeing 'censored' data, but it didn't... splom(lung[, c("sex", "age", "OS")], lower.panel = panel.lowerTest) ________________ John Gosink, PhD Senior Scientist, Amgen, Inc. 1201 Amgen Court West Seattle, WA 98119-3105 (w) 206 265-8217

...$modelStruct$varStruct)["const"]; #const = -0.4229219 > coef.varFunc(model3$modelStruct$varStruct); #R can't seem to find this > function, Splus can > > > Any advice on what I am doing wrong would be appreciated. > > Hugh Rand > Senior Scientist > Amgen > rand at amgen.com > > ______________________________________________ > R-help at stat.math.ethz.ch mailing list > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html

...n Princeton, NJ. This year's conference focuses on the techniques and methodologies pivotal to the increased demand for statistics in business organizations. The conference program will include presentations by S-PLUS experts from market-leading companies such as Novartis, Procter & Gamble, Amgen, and MCI. Of particular interest to subscribers to r-help is the panel discussion on future plans for compatibility between S-PLUS and R. The panel discussion is on the topic "Ensuring Compatibility and Portability Between the Two Prime Dialects of the S Language". Martin DeBono from...

Hi, when I generated a survfit() object, I can get number of patients at risk at various time points by using summary(): fit<-survfit(Surv(time,status)~class,data=mtdata) summary(fit) class=1 time n.risk n.event survival std.err lower 95% CI upper 95% CI 9.9 78 1 0.987 0.0127 0.963 1 41.5 77 1 0.974 0.0179 0.940 1 54.0 76

Hello: Let me rephrase my question to attract interest in the problem I'm having. When I appply anova() to two equations estimated using glmmPQL, I get a complaint, > anova(fm1, fm2) Error in anova.lme(fm1, fm2) : Objects must inherit from classes "gls", "gnls" "lm","lmList", "lme","nlme","nlsList", or "nls"

...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 Statistician, Amgen Inc. Example of problem derived from example(layout) code. ##Start windows() x <- pmin(3, pmax(-3, rbeta(5000,1,5))) y <- pmin(3, pmax(-3, c(rnorm(2500, .25, 1), rnorm(2500))) ) trtp.f <- factor(rep(c("Placebo", "Not Placebo"), each=2500), levels=c(&q...

I'm just getting started with R, so feel free to point me to the appropriate documentation if this is already answered somewhere (though I've been unable to find it myself). This does seem like a rather basic question. I want to fold a table into a matrix. The table is formatted like so: Column_Index Value 1 486 2 688 3 447 4 555 5

Hi, Just a quick query - if I'm creating a function to produce a number of histograms per page of output (one per column from a matrix), how can I pass the column name of the matrix into the title (or indeed to form part of the x-axis label)? TIA, Laura Laura Quinn Institute of Atmospheric Science School of Earth and Environment University of Leeds Leeds LS2 9JT tel: +44 113 343 1596 fax: