similar to: printing the name of the arguments passed to a function

On all current R versions, afaik: When using save(....., ascii = TRUE), the resulting file is not always readable by load(.). The following is executable R code exhibiting two (slightly different) examples of which I have seen even more [~/R/MISC/saveload-bug.R]: ### BUG : save(..., ascii = TRUE) produces load() - unreadable files : tstf1 <- function(n = 500) { ## Purpose: ##

Dear all How can I obtain the name of the argument passed in a function? Here is a simplistic example of what I would like to obtain: myfunction= function(name) { print(paste("The parameter name was",unknownFunction(name)) } myfunction(myobject) [1] "The parameter name was myobject" Thanks Francisco

I am trying to install R.2.1.1 under Linux FC4 using RPM. During the installation I get the following message: The following package could not be found on your system. Installation cannot continue until it is installed. The package in question of libblas.so.3. Does anyone know where I can get the package? Can I use RPM to intall the package? Thanks, John John Sorkin M.D., Ph.D. Chief,

Is there anyway I can convert a vectors objectname to a string to be used in fx: Monkey<-c(0,0,0,1,1,1) Wax<-c(1,0,1,0,1,0) f<-function(x,y){ table(x,y) } f(Monkey,Wax) so that the printout is not y x 0 1 0 1 2 1 2 1 but Wax Monkey 0 1 0 1 2 1 2 1 -- View this message in context:

R 2.1.1 Win 2k Would someone suggest a method (or methods) that can be used to determine ntile cutpoints of a vector, i.e. to determine values that can be used to divide a vector into thirds such as 0-33 centile, 34-66 centile, 67-100 centile. If for example I had a vector: 1,2,3,4,5,6,7,8,9 and wanted to divide the vector into thirds I would have cut-points of 3, and 6. Thanks, John John

I am using read.csv to read a CSV file (produced by saving an Excel file as a CSV file). The columns containing dates are being read as factors. Because of this, I can not compute follow-up time, i.e. Followup<-postDate-preDate. I would appreciate any suggestion that would help me read the dates as dates and thus allow me to calculate follow-up time. Thanks John John Sorkin M.D., Ph.D. Chief,

I want to place four plots on a page, and I would like to have all four plots share a common title. I have tried the following code, but the title is centered over the fourth graph and not centered across all four plots. Does anyone have any suggestions? R 2.1.1 windows xp oldpar<-par(mfcol =c(1,4),ask=TRUE) plot(p,varp) plot(p,SEp) plot(p,CVp) plot(p,ppval) title(paste("P and 95%CI

Can anyone suggest a good non-parametric test, and an R implementation of the test, that allows for two or more factors? Wilcoxon signed rank allows for only one. Thanks, John John Sorkin M.D., Ph.D. Chief, Biostatistics and Informatics Baltimore VA Medical Center GRECC and University of Maryland School of Medicine Claude Pepper OAIC University of Maryland School of Medicine Division of

I am having trouble selecting rows of a dataframe that will be included in a regression. I am trying to select those rows for which the variable Meno equals PRE. I have used the code below: difffitPre<-lm(data[,"diff"]~data[,"Age"]+data[,"Race"],data=data[data[,"Meno"]=="PRE",]) summary(difffitPre) The output from the summary indicates that

Is there an R function that can be used to calculate a sample size for a linear regression: i.e. Given a linear regression: y=f(x,z) i.e. fit1<-lm(y~x+z) is there a function which can be passed coef(z), SE(z), and other parameters to determine a sample size based on z, SEz and other information from the regression? Thanks, John John Sorkin M.D., Ph.D. Chief, Biostatistics and Informatics

I think I am doing something wrong when I try to bootstrap R square obtained from lm. My code is included below. No matter how many times I run the simulation, I always get exactly the same result, the bias and std.error are always zero. I would think that these values should be non-zero. I would appreciate any suggestions as to what I am doing wrong, or perhaps what I fail to understand. R 2.1.0

R 2.7.2 Windows XP I am using apply to compute a series of Student's t-test from two matrices, sample1 and sample2. boo<-apply(sample1,1,t.test,sample2) I want to pick of the p-values from the tests, but can't seem to get it to work. I have tried several methods to get the values including: boo<-apply(sample1,1,t.test$t.test,sample2) boo<-apply(sample1,1,t.test,sample2)$t.test

I am trying to model data in which subjects are followed through time to determine if they fall, or do not fall. Some of the subjects fall once, some fall several times. Follow-up time varies from subject to subject. I know how to model time to the first fall (e.g. Cox Proportional Hazards, Kaplan-Meir analyses, etc.) but I am not sure how I can model the data if I include the data for those

Is there any way to load a package using R code rather than the Load Package command in the Package menu? I would like to load the boot package. R.2.1.0 Patched Win 2k. Thanks, John John Sorkin M.D., Ph.D. Chief, Biostatistics and Informatics Baltimore VA Medical Center GRECC and University of Maryland School of Medicine Claude Pepper OAIC University of Maryland School of Medicine Division of

I am having trouble with apply. Could someone suggest changes to the code below that will eliminate the for loops? R 2.1.1 patched Windows 2k Thanks, John function () { NaCon<-array(dim=c(10,10)) for (i in 1:10) {NaCon[i,]<-rnorm(10,10.77,0.02)} MeanNaCon<-vector(mode="numeric",length=10) for (j in 1:10) {MeanNaCon[j]<-mean(NaCon[j,1:10])} print(MeanNaCon)

I am trying to fit lme using R 2.1.1 under Windows 2k. I am getting the following message noted below. Any suggestions that would help me correct my error would be greatly appreciated. Thanks, John > fit1lme<-lme(Velocity~time,data=gate) Error in getGroups.data.frame(dataMix, groups) : Invalid formula for groups John Sorkin M.D., Ph.D. Chief, Biostatistics and Informatics

R 2.1.1 on Windows 2K I hope one and all will allow both an R question and a general stats question: (1) Is there any function that will perform non-parametric tests such as a sign test or a signed rank test? I know I could program both but I would prefer not re-inventing the wheel. (2) When performing a sign test, the usual practice is to drop zero values. Is there any theoretical reason

R 2.2.0 windows XP How can I perform a regression analyses using a vector of means, a variance-covariance matrix? I looked at the help screen for lm and did not see any option for using the afore mentioned structures as input to lm. Thanks, John John Sorkin M.D., Ph.D. Chief, Biostatistics and Informatics Baltimore VA Medical Center GRECC, University of Maryland School of Medicine Claude D.

R 2.3.1 Windows XP Question: How can I obtain the values of BIAS and STD. ERROR from a bootstrap. Background: I am running a bootstrap: result2<-boot(1:400,regSEvssample,R=5000) and obtain the following results: Bootstrap Statistics : original bias std. error t1* 1.876602 -0.0001368616 0.1630380 I would like to get the values of ORIGINAL and BIAS. I can get the value of

Given a square matrix of variables X, is there any way to compute a multivariable correlation among all the variables in the array? It is possible to calculate the correlation of all pairs of variables in the array, but I want to know the correlation of all the variables taken together, i.e. for the matrix X=x1 x2 x3 x4 x5 x6 x7 x8 x9 I don't want the pair-wise correlations