search for: cdrh

I would like to use R for submissions to FDA/CDRH (the medical device company I work for currently uses only SAS). Previous postings to the list regarding R and 21 CFR 11 compliance have been very helpful. However, reluctance to using open source software for statistical analyses and reporting remains high here at my company. Has anyone used R...

I checked the help and the mailing list archives, but I can find no mention of a routine that calculates higher moments like skewness and kurtosis. Of course, these are easy enough to write myself, but I was thinking that they MUST be in here. Am I wrong? Thanks. -Frank

...-- 7 values [1] 2.002566e-308 2.222849e-308 2.225074e-308 2.225074e-308 2.225074e-308 2.227299e-308 2.447581e-308 Does anyone really care about being correct to 1 unit of machine precision? If you do, you have a bad algorithm. -----Original Message----- From: Samuelson, Frank* [mailto:FWS4 at cdrh.fda.gov] Sent: Thursday, August 19, 2004 12:11 PM To: 'r-help at stat.math.ethz.ch ' Subject: [R] precision problems in testing with Intel compilers I compiled the 1.9.1 src.rpm with the standard gnu tools and it works. I tried compiling the 1.9.1 src.rpm with the Intel 8 C and FORTRAN c...

It appears that the maximum number of connections available in R is about 48. Can anyone tell me how to bump this number up? I've been perusing the source, but any info would speed things up. Is there a reason that it was set to such a low number? Thanks for any help. -Frank

Does anybody know if R is FDA or ICH (or EMEA...) compliant? AFAIK S-Plus is but that means nothing... -- Trebate bolji pristup internetu? Nazovite IskonInternet na 0800 1000 ili pogledajte http://www.iskon.biz/individualni/usluge/dialup/

Dear All, discussing with a statistician of a pharmaceutical company I received this answer about the statistical package that I have planned to use: As R is not a validated software package, we would like to ask if it would rather be possible for you to use SAS, SPSS or another approved statistical software system. Could someone suggest me a 'polite' answer? TIA Giovanni -- dr.

I'm running R-2.0.0 (yesterday's snapshot)in its own directory, and everything works great, except: > .Random.seed Error: Object ".Random.seed" not found Does 2.0.0 not use .Random.seed for saving, etc, like it says in the help page? Thanks for any help. -Frank

Is there a way to tell how much memory the computer running R has? -Frank -----Original Message----- From: David Smith [mailto:dsmith at insightful.com] Sent: Monday, February 09, 2004 1:32 PM To: Ross Boylan Cc: r-help Subject: RE: [R] R does in memory analysis only? Ross Boylan writes: > R works only on problems that fit into (real or virtual) memory. > ... does S-Plus have the same

Hi again, thanks a lot for the quick answer. I just forgot the comma, always these stupid mistakes... Anyways, as I said before, I have two data.frames containing about 1000 rows and I would like to avoid looping through all of them... In each data.frame are coordinates (x,y,z), so every row is giving the information on one single point. I would like to calculate the distance from each point

If anyone has a few extra CPU cycles to spare, I'd appreciate it if you could verify a problem that I have encountered. Run the code below and tell me if it crashes your R before completion. library(lme4) data(bdf) dump<-sapply( 1:50000, function(i) { fm <- lme(langPOST ~ IQ.ver.cen + avg.IQ.ver.cen, data = bdf, random = ~ IQ.ver.cen | schoolNR); cat("

Hi all there Can some one clarify me on this issue, features wise which is better R or SAS, leaving the commerical aspect associated with it. I suppose there are few people who have worked on both R and SAS and wish they would be able to help me in deciding on this. THank you for the help --------------------------------- [[alternative HTML version deleted]]

Is there a routine in the boot package to get the jackknife-after- bootstrap influence values? That is, the influence values of a jackknife of the bootstrap estimates? I can see how one would go about it from the jack.after.boot code, but that routine only makes pretty pictures. It wouldn't be hard to write, but I find it hard to believe this isn't part of the package already. Thanks

How does one loop over factors? Perhaps this is a newbie question. I tried: > b [1] caseX caseY caseZ Levels: caseX caseY caseZ > length(b) [1] 3 > > for (i in b) { + print (b == i) ; + print (i); + } [1] FALSE FALSE FALSE [1] 1 [1] FALSE FALSE FALSE [1] 2 [1] FALSE FALSE FALSE [1] 3 > But that strangely doesn't work. I must protest the implications of the above. i , as

Anybody know where I can get the Pinheiro/Bates book? I can't find a bookstore w/ stock and the publisher says they don't know when they'll have it again. Thanks. -Frank

1. How do I get a grid in the background of my xyplots, like Figure 1 of http://www.ci.tuwien.ac.at/Conferences/DSC-2001/Proceedings/Murrell.pdf 2. In my scatter xyplot I have multiple columns and rows of panels. I would like to have different y scales for panels in different rows (panels in same rows have same y scales). In such a scenario, only the outside axes are necessary,

I originally posted this to R-help, but here is probably a better place. I have a rather complex routine which works great in R. Then I bootstrap it. After running a few minutes it typically segfaults, though it happens very randomly. This happens w/ 1.9.0 and 1.9.1 on Fedora1 and a Scyld (RH7.3) system. Traces always point back to routines in src/main/memory.c, though often not the same

I compiled the 1.9.1 src.rpm with the standard gnu tools and it works. I tried compiling the 1.9.1 src.rpm with the Intel 8 C and FORTRAN compilers and it bombs out during the testing phase: comparing 'd-p-q-r-tests.Rout' to './d-p-q-r-tests.Rout.save' ...267c267 < df = 0.5[1] "Mean relative difference: 5.001647e-10" --- > df = 0.5[1] TRUE