similar to: FW: [R] The Quality & Accuracy of R

Hi All, We have all had to face skeptical colleagues asking if software made by volunteers could match the quality and accuracy of commercially written software. Thanks to the prompting of a recent R-help thread, I read, "R: Regulatory Compliance and Validation Issues, A Guidance Document for the Use of R in Regulated Clinical Trial Environments (http://www.r-project.org/doc/R-FDA.pdf).

Hi Folks, SAS Institute is adding official support for R: http://support.sas.com/rnd/app/studio/Rinterface2.html Cheers, Bob ========================================================= Bob Muenchen (pronounced Min'-chen), Manager, Research Computing Support U of TN Office of Information Technology Stokely Management Center, Suite 200 916 Volunteer Blvd., Knoxville, TN

Hi All, I'm using the Hmisc summarize function and used list instead of llist to provide the by variables. It generated an error message. Is this a bug, or do I misunderstand how Hmisc works with lists? The program below demonstrates the error message. Thanks, Bob x<-1:8 group <- c(1,1,1,1,2,2,2,2) gender<- c(1,2,1,2,1,2,1,2) mydata<-data.frame(x,group,gender)

Dear HelpeRs, I'm confused about the role of quotes around package names on the library and detach functions. Books on R use both approaches: library(Hmisc) describe(mydata) detach(package:Hmisc) and library("Hmisc") describe(mydata) detach("package:Hmisc") The help file for detach says "quoted or unquoted" and the help file for library says about the

Hi All, I searched around to find the number of R packages currently available, but didn't find anything, so I choose all repositories & told it to install. The list contained about 2,856 (correcting roughly for those installed). But the list includes repetitions such as 19 names that begin with "bvbovine". Selecting only CRAN and CRAN(extras) I get 1,344. Is there an easier

Hi All, I'm learning JGR 1.4-15 with R 2.4.1 in Windows XP (all patches applied). JGR looks great but I'm having trouble getting the data editor to save my results. I don't see anything in R-help about it. Here are the steps I followed: 1. I chose "Tools>Object Browser" & double-clicked on a data frame, "mydata". 2. A spreadsheet editor popped up and

Hi All, If I need to stack two data frames, I can use rbind, but it requires that all variables exist in both sets. I can make that happen, but other stat packages would figure out where the differences were, add the missing variables to each, set their values to missing and stack them. Is there a more automatic way to do that in R? Below is an example program. Thanks, Bob # Top data frame has

Hi All, I've been fiddling around with various ways to estimate the popularity of R, SAS, SPSS, Stata, JMP, Minitab, Statistica, Systat, BMDP, S-PLUS, R-PLUS and Revolution R. It's not an easy task. You can see what I've come up with so far at http://r4stats.com/popularity . I'm sure people will have plenty of ideas on how to improve this, so please let me know what you think.

Hi All, I'm using the subset function to select a list of variables, some of which are contiguous in the data frame, and others of which are not. It works fine when I use the form: subset(mydata,select=c(x1,x3:x5,x7) ) In reality, my list is far more complex. So I would like to store it in a variable to substitute in for c(x1,x3:x5,x7) but cannot get it to work. That use of the c function

Dear R-Helpers, Why does R show character missing values in vectors as NA and when stored in a data frame as <NA>? I've searched but did not find an explanation. Thanks, Bob > gender <- c("f","f","f",NA,"m","m","m","m") > gender [1] "f" "f" "f" NA "m" "m"

Hi All, I'm perplexed by the way the unclass function displays a factor whose labels have been swapped with the relevel function. I realize it won't affect any results and that the relevel did nothing useful in this particular case. I'm just doing it to learn ways to manipulate factors. The display of unclass leaves me feeling that the relevel had failed. I've checked three books

Dear R-Helpers, I suspect I'm about to ask a FAQ, but I haven't been able to find an answer in the FAQ, AItR or an R Site Search. When I look at the methods of summary (below) it says, "Non-visible functions are asterisked". I looked at the help file for summary.princomp, which did not comment on it being non-visible. I ran its help file example, which printed visible output. I

Hi All, When I cut & paste help file examples into a script window, about half the time it pastes as a single long line. The steps I follow are: 1. Open a help file e.g. ?data.frame. 2. Select the examples at the bottom. 3. Choose File: Copy. 4. Return to the console. 5. Choose File: New script. 6. Choose File: Paste or do CTRL-V. The examples frequently paste as a single long line. I

Hi All, I'm fiddling with an program to read a text file containing periods that SAS uses for missing values. I know that if I had the original SAS data set instead of a text file, R would handle this conversion for me. Data frames do not allow missing values in their indices but vectors do. Why is that? A search of the error message points out the problem and solution but not why they

Greetings, I am pleased to announce the availability of the document, "R for SAS and SPSS Users", at http://oit.utk.edu/scc/RforSAS&SPSSusers.doc . It presents an introductory view of R for people who already know SAS and/or SPSS. Included are 27 programs written in all three languages (i.e. 81 total) so that people can see how R works compared to the other two, task by task. I

Hi All, I'm working on a class example that demonstrates one way to deal with factors and their labels. I create a function called myLabeler and apply it with lapply. It works on the whole data frame when I subscript it as in lapply( myQFvars[ ,myQFnames ], myLabeler ) but does not work if I leave the [] subscripts off. I would appreciate it if anyone could tell me why. The program below

Hi All, I'm trying to understand the difference between do.call and lapply for applying a function to a list. Below is one of the variations of programs (by Marc Schwartz) discussed here recently to select the first and last n observations per group. I've looked in several books, the R FAQ and searched the archives, but I can't find enough to figure out why lapply doesn't do what

Dear R-Helpers, I'm fiddling with my .Rprofile in Windows XP & R 2.7.0 Beta. I prefer to manually save my workspace but automatically save my command history via the .Rprofile. That is working fine once I found that "utils::" was required before the loadhistory & savehistory functions. What I would like to do is add a separator line with a date between the histories of each

Hi All, A few weeks ago I suggested that it would be nice to be able to submit lines from the help files for execution. You can cut and paste them into the console, or enter example(function) to run them all. However, I often find myself wanting to run just a line or two, or even parts of a line to see an intermediate result. It turns out that this is one of the many nice features of the JGR

Hi All, I'm stumped on something that must be trivial. I created a correlation matrix on 4 variables (6 correlations) using Hmisc's rcorr function. I wanted to correct the P-value matrix for the number of tests done, so I ran it through the p.adjust function. That function adjusted for the 12 p-values it saw, rather than 6. I added the argument n=6 to p.adjust but it requires that n be