search for: paulda

--- "Paul, David A" <paulda at BATTELLE.ORG> wrote: > Four questions: > > (1) Is Splus faster (on a PC) under Linux than under > Windows? I am > specifically thinking of intensive loops here, but > would also > solicit > information regarding how Splus stores information > within the Linux &...

R1.7.1/Win2k: Apologies if this posts twice - the first message was not in plain text. I have looked in help.start() and tried typing "crossprod" and "%*%" into the RGui to get an idea for what R is using as internal algorithms for its matrix computations/manipulations... to no avail. Could someone point me in the direction of some documentation? All I get for

I am one of only 5 or 6 people in my organization making the effort to include R/Splus as an analysis tool in everyday work - the rest of my colleagues use SAS exclusively. Today, one of them made the assertion that he believes the numerical algorithms in SAS are superior to those in Splus and R -- ie, optimization routines are faster in SAS, the SAS Institute has teams of excellent numerical

Win2k, R1.6.2. I've been using Splus 6.1 and wanted to try the same regression analysis in R. Using "names( blah.lm )" in R yields [1] "coefficients" "residuals" "effects" "rank" [5] "fitted.values" "assign" "qr" "df.residual" [9] "xlevels"

Undoubtedly a simple question: I've looked at order() and sort() in the help pages for R1.7.1. It doesn't appear that these functions are immediately suited to doing the same thing as PROC SORT DATA = BLAH; BY X Y Z; RUN; in SAS. I have also checked Frank Harrell's Hmisc library. Could someone point me in the right direction so I can sort by the levels of Z within the levels of

Hi All I am really very interested in starting to use R within our company. I particularly like the open source nature of the product. My company is a medical research company which is part of the University of London. We conduct contract virology research for large pharma companies. My question is how do we validate this software? I wonder if anyone else has had the problem and might be able to

...000,000 > > > >Dave > > > > > 20-OP > >________________ > >Date: 04 Jun 2003 08:53:36 -0500 > >From: Marc Schwartz <MSchwartz at medanalytics.com> > >Subject: RE: [R] Rounding problem R vs Excel > >To: "Paul, David ?A" <paulda at BATTELLE.ORG> > >Cc: R-help at stat.math.ethz.ch, "'Duncan Murdoch'" > <dmurdoch at pair.com> > > > >On Wed, 2003-06-04 at 08:09, Paul, David A wrote: > >> I don't have the reference, but a biologist friend of mine > >> once...

I have an older Splus script that generated some regression graphics for me. In it, I used the commands key(10,800,text="<blah1>",border=0) key(10,750,text="<blah2>", border=0, background = 9) text(locator(1),"95% Confidence Limits") locator(n=2,type="l") text(locator(1),"Model for Group 1") locator(n=2,type="l") There is

In a contour plot having numeric values displayed next to the contours, is it possible to modify these numeric values so that they are bolded, enlarged, and/or placed differently? Much thanks in advance, David Paul, Ph.D. Battelle Memorial Institute 614.424.3176

I think you've badly misinterpreted the purpose of the R listserv with this request: https://www.stat.math.ethz.ch/mailman/listinfo/r-help says "The `main' R mailing list, for announcements about the development of R and the availability of new code, questions and answers about problems and solutions using R, enhancements and patches to the source code and documentation of R,

I've been using par() to check the graphics parameters associated with both plot(<fitted linear model>) and plot(<grouped data object>). AFAIK the only differences are in the $cxy, $usr, $xaxp, and $yaxp parameters but the background color for the grouped data plot is grey while the linear model plot has a white background. When I've tried par(bg = "white") prior to

R 1.6.2 for Windows, Win2k: I have fitted a weighted least squares model using the code "wls.out <- gls(y ~ x1 + x2 + x3 + x4 + x5 + x6 - 1, data = foo.frame, weights = varConstPower(form = ~ fitted(.), fixed = list(power = 0.5), const = 1))" The data has 62 rows and the response is zero when the covariates are zero. The purpose of the model was to account for the the fact that

R1.8.0, Win2k: When I run the code lset( list( background = list(col = "white"))) xyplot ( y ~ TIME , data = foo.frame, scales = list(alternating = FALSE), ylim = c(.75,y.max), panel = function(x, y, panel.number, ... ) { panel.superpose(x = foo.frame$TIME[foo.frame$group == 1], y = foo.frame$y[foo.frame$group == 1], subscripts = TRUE, groups =

R1.7.0, Win2k: When I use plot( ) on a groupedData object, if I have 165 subjects I'm supposed to be able to use plot(..., layout = c(5,3,11)) to get all 165 on 11 separate sheets. The graphics window is only displaying the first 10. If I use plot(..., layout = c(5,3,12)) I get the appropriate 11 sheets. Furthermore, yesterday I noticed that when using plot( <linear model object> )

In Splus the code test.lm <- lm(y ~ x, data = test.data) plot(test.lm) generates a graphics window that contains multiple graph sheets that one may choose from via the "page" tabs at the bottom of the window. Is there a way to do this sort of thing in R? As another example, I have some repeated measures data with continuous outcomes and have been working with the nlme library in

For the first error message: > win.metafile(file = "//.../plot1.wmf", + width = 8.5, height = 6.25) > lset( list( background = list(col = "white"))) Error in get(x, envir, mode, inherits) : variable "win.metafile://.../plot1.wmf" was not found > traceback() 4: get(device) 3: trellis.device(device = .Device, new = FALSE) 2: trellis.par.get(item) 1:

Last week (Wed 9/10/2003, "regression questions") I posted a question regarding the use of gnls( ) and its dissimilarity to the syntax that nls( ) will accept. No one replied, so I partly answered my own question by constructing indicator variables for use in gnls( ). The code I used to construct the indicators is at the end of this email. I do have a nagging, unanswered

My thanks to Drs. Armstrong, Bates, Harrell, Liaw, Lumley, Prager, Schwartz, and Mr. Wang for their replies. I have pasted my original message and their replies below. After viewing http://www.itl.nist.gov/div898/strd/ as suggested by Dr. Schwartz, it occurred to me that it might be educational to search for some data repositories on google. I was able to find some,though I'm sure many of

Win2k, R1.7.1: I am currently working with some growth curve data from a biotoxicology experiment. Each of 12 subjects had their blood drawn at 0, 2, 4, 6, 8, and 10 weeks. For the purposes of the project, it would be helpful if I were able to do the following: a. Produce 12 panels, each displaying the *same* data, with the "strip" at the top of a particular panel showing