similar to: Introductory Resources

dum is a simple data frame transferred to Splus using the dump() command in Splus and the source() in R. All fields are numeric. There are no missing data. The data frame looks like it is should: > apply(dum,2,mode) yrcl sland s02 s234 "numeric" "numeric" "numeric" "numeric" > apply(dum,2,is.vector) yrcl sland s02 s234

Hi, Does anyone know how I can retrieve a function name, for example If I have a function f as follows: f <- function( myfunc ) { print( name_of(myfunc) ); } I want to know what I should have as "name_of" such that I could call this with : f( median ) and it would print "median" or f( my_function_name ) and it would print "m_function_name". So far all I

Does R have a bootstrap command like SPlus, or do we have to form our own code if we want to do bootstrapping in R? Jason Parcon Western Michigan University e-mail addresses: s0parcon at wmich.edu parcon at stat.wmich.edu -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read

I'm trying to create a function which will allow me to subset a data set based on values of various specified variables. I also want to then apply some other function(s) (e.g., summary). This is what I've tried so far.... > test.fx <- function(dta, expvar, expval) { + newdta <- subset(dta, eval(expvar)>expval) + summary(newdta$eval(expvar)) + } > >

OK, I am forwarding this to the r-help. In the original message I simply hit reply and I thought I send to ALL r-help list. (apparent not) also I recieved another piece of mail.....about the RNG in Splus 4.0 or later indicating it was modified and getting worse. I only have Splus 3.4 for HP. So I cannot try it. This makes the idea of using the pre-generated RN on CDROM more appealing.. Mai

I'm updating some (very) old code, and one particular option of its plot method depends on a once-was-true trick polygon(x, y, border=TRUE, col=0) polygon(x, y, border=TRUE, density=0) would draw the polygon AND erase whatever was underneath it back to background color. Is there a reliable way to do this with the current R (standard graphics)? Terry Therneau PS For the

RNG in R and Splus 3.4 Prof. Ripley asked the details of the example. We were doing parametric bootstrap, so it is similar to simulation. Anyway here is the details. We start with a sample of 19 positive numbers. We know the sample is from truncated exp(0.3)...only the truncation point, theta, is unknown. In other words, the sample can be generated from something like x1 <- rexp(100,

RNG in R and Splus 3.4 Prof. Ripley asked the details of the example. We were doing parametric bootstrap, so it is similar to simulation. Anyway here is the details. We start with a sample of 19 positive numbers. We know the sample is from truncated exp(0.3)...only the truncation point, theta, is unknown. In other words, the sample can be generated from something like x1 <- rexp(100,

I have found that specifying different "sizes" in the sample command has a funny effect on the random sampling. The code below is a condensed version of a function I wrote to simulate a bootstrap method. For simplicity, I eliminated the internal bootstrap loop, but kept a statement to draw one bootstrap sample, because this is where the problem occurs. The output (mean(y)^2) should be

I have been unsuccessful in extracting the dispersion parameter in SPLUS 6.1 using summary or summary.glm(modelobj$dispersion) from a glm object in which the family was set to quasi. This is the syntax given in the manual. I want to write a script to bootstrap the estimate of the dispersion parameter, but cannot seem to access that value. Any suggestions? Thanks, Ed ----------

Hi, In SPlus there is a function called "twoway" for median polishing gridded data. Is there an equivalent function in R? I have been searching for it in R help without much success. Your help is much appreciated. Cheers Sorn [[alternative HTML version deleted]]

Hye every one, My question will certainly seem stupid as I am quite a beginner in R. I would like to trace a curve which fits these two vectors: x<-c( 2,3,4,5,6,7,8,10 ) y<-c( 20, 12, 8, 6, 5, 4.5, 4, 3.8) It seems to follow a non linear model. Could anyone help me because I could'nt find the answer I am looking for in the FAQs. In advance thank you for your time. G. Lefebvre

Hello, how can I get a variable driven summary of one column of my data.frame? Usually I would do > summary(data$columnname) to get a summary of column named "columnname" of my data.frame named "data". In my case the columnname is not static but can be set dynamically. So I save the chosen columname in something like variable <- "columnname" but how can I

Hello, I learnt there's an advanced R/Splus course in Boston this january. Anyone got the announcement? please kindly forward it to me. Best, Eugene

Dear R-Help I am running apply on a data.frame containing factors and numeric columns. It appears to convert are columns into as.character? Does it convert data.frame into matrix? Is this expected? I wish it to recognise numerical columns and round numbers. Can I use another function instead of apply, or should I use a for loop in the case? > summary(xmat) A B

Dear R Help List, I have a vector of n and a vector of n-1 and I want to use boot() to bootstrap the ratio of their respective medians. I want to eventually use boot.ci() to generate confidence intervals. Because the vectors are not equal in length, I tried a few things, but have yet to be successful. Try 1: > x <- runif(20) > > y <- c(runif(19), NA) > > median(x) [1]

Hi Listers, I am producing some graphics that the commands are in a FUNCTION... The problem is that I end up viewing just last graphic and in my FUNCTION there are 4 graphics with the PAR command function... Like those below... How do I view/get the other 3 graphics? Any help? Thanks in advance... histogram<-par(mfrow=c(1,2)) hist(rw_mean_app,main='Bootstrap Method

(I'm up too late so this might come through garbled...) I've just been doing some bootstrapping on data frames and I discovered that S-plus 6.0r1 was a *lot* faster than R 1.3.1 at the task. Splus was completing 100 bootstrap iterations in about 4 seconds while R was taking about 15 seconds. However, doing bootstrapping on equivalent *matrices* R was slightly faster, 1.5 seconds verses

I apologize for this post. I am new to R (two days) and I have tried and tried to calculated confidence intervals for medians. Can someone help me? Here is my data: institution1 0.21 0.16 0.32 0.69 1.15 0.9 0.87 0.87 0.73 The first four observations compose group 1 and observations 5 through 9 compose group 2. I would like to create a bootstrapped 90% confidence interval on the difference of

Dear R-experts, I am trying to bootstrap (and average) the median squared error evaluation metric for a robust regression. I can't get it. What is going wrong ? Here is the reproducible example. ############################# install.packages( "quantreg" ) library(quantreg) crp <-c(12,14,13,24,25,34,45,56,25,34,47,44,35,24,53,44,55,46,36,67) bmi