similar to: Is R really an open source S+ ?

I'm just reading Modern Applied Statistics with S, 4th Ed., Venables and Ribley I'm typing in their examples in both R and S+. I need insight in the difference in the class() statement shown in Chap. 2. Example from book: > names(powers.of.pi) <- -2:2 > powers.of.pi -2 -1 0 1 2 0.1013212 0.3183099 1.0000000 3.1415927 9.8696044 >

I need some guidance on what is the best way to automate R I am aware of Rterm. First question: Is it true that R is not currently OLE accessible under Windows 2000? Second question: Is there an R command that echoes all R Console results to a run log system file? My first attempt to write a run log file is the example from the R Data Inport/Export documentation. > zz<-

All Is there an R command to produce error bars around a plotted point. Crawley's book uses the command error.bar() but my version of R rejects it. Must be an S+ command(?). Thanks REX [[alternative HTML version deleted]]

To all I have some simple questions pertaining to weights used in regression. If the variability of the dependent variable (y) is a function of the magnitude of predictor variable (x), can the use of weights give an appropriate answer to the regression parameters and the std errors? Assume that y at x=1 and 6 has a standard deviation of 0.1 and at x=11 it is 0.4 Then according to a web page on

Any idea why I can no longer plot two graphs on the same graphics device using the par(fig=) parameter? A simpler par(mfrow=c(1,2)) does work, showing the two plots side-by-side, but I would like the first to be larger. This simple example fails: x<-c(1,1,NA,2,2,NA,3,3) y<-c(2,4,NA,3,5,NA,1,4) par(fig=c(0,2/3,0,1)) plot(x,y) par(fig=c(2/3,1,0,1)) qqnorm(x) When plotted, the last

To All How does one get a list of functions from a loaded package so that one can then get the appropriate help for each of the functions. Currently my method is based on a lot of trial-and-error. Here's an example of what I mean... >From this forum I learn that an interesting package called "multtest" exists on Bioconductor. I then use R Console's "Packages" --

Is there an eloquent solution to re-assign vector element values? I have a vector which contains chemical data, some of them are "flagged" as non-detected values by their negative values. I can find the statistics on the positive values in vector "v" simply by typing: >v<- c(5,5,-3,-3,7,8,10) > v[(v>0)] [1] 5 5 7 8 10 I can also convert to positive values by

The last two times I have originated message threads on R or Bioconductor I have received the message included below from someone named Patrick Connolly. Both times I was the originator of the message thread and used what I thought was a unique subject line that explained as best I could what my question was. Patrick seems to be implying that I am abusing the R and BioC help newsgroups in this

How do I post text-string dates along the x-axis instead of the excel_date position values? I would like to plot a time-series of chemical values changing with time The first column is the internal date from excel. This should be the x-axis position of a graph of the plotted data. The second column is a text-string of date. This should be what is posted (instead of the excel_date). The third

I'm new to Linux/Ubuntu and need some help getting Legions: Overdrive (http://www.legionsoverdrive.com/) with PlayOnLInux and WINE. Legions: Overdrive is a sci-fi FPS game, Freedom of movement returns in this spiritual successor to Tribes. Travel by jetpack and ski over foreign landscapes, dodging grenades and shooting enemies in deathmatch mode. Sharpen your shooting accuracy or prepared to

R-wizards I have successfully run with S+ a statistical power calculation for non-normal distributions as presented in M. Crawley?s new book.? When I tried the newest version of R on the same code, the $ parse statement doesn't seem to retrieve the appropriate number from a table. Note that some of the cosmetic differences in the two tables have to be dealt with by the parser. Any ideas

Has anyone on this list got sendmail and cyrus-imapd configured so that email addressed in the form username+folder at example.com gets delivered directly to cyrus' user.username.inbox.folder? I have been tinkering with this on and off for about a week and no matter what I have tried, all traffic addressed in the +folder format ends up in user.username.inbox instead. If anyone has gotten

Hi R-experts, I want to use the jackknife function from the bootstrap package onto a linear model. I can't figure out how to do that. The manual says the following: # To jackknife functions of more complex data structures, # write theta so that its argument x # is the set of observation numbers # and simply pass as data to jackknife the vector 1,2,..n. # For example, to jackknife #

I''m trying to run mongrel on port 80 and switch to a different user/group from root. As root, this works: mongrel_rails start -p 3000 -u rails --group rails This: mongrel_rails start -p 80 -u rails --group rails dies with: --------------------------- ** Changing group to rails. ** Changing user to rails. ** Starting Mongrel listening at 0.0.0.0:80

Hi dear all, Can aynone help me about delete-d jackknife usually normal jackknife code for my data is: n <- nrow(data) y <- data$y z <- data$z theta.hat <- mean(y) / mean(z) print (theta.hat) theta.jack <- numeric(n) for (i in 1:n) theta.jack[i] <- mean(y[-i]) / mean(z[-i]) bias <- (n - 1) * (mean(theta.jack) - theta.hat) print(bias) but how i can apply delete-d jackknife

Below is a SIMEX object that was generated with the "simex" function from the "simex" package applied to a logistic regression fit. From this mountain of information I would like to extract all of the values summarized in this line: .. ..$ variance.jackknife: num [1:5, 1:4] 1.684 1.144 0.85 0.624 0.519 ... Can someone suggest how to go about doing this? I can extract the

While doing pls I found the following problem > BHPLS1 <- plsr(GroupingList ~ PCIList, ncomp = 10, data = PLSdata, jackknife = >FALSE, validation = "LOO") when not enabling jackknife the command works fine, but when trying to enable jackknife i get the following error. >BHPLS1 <- plsr(GroupingList ~ PCIList, ncomp = 10, data = PLSdata, jackknife = >TRUE,

Dear R friends I´m interested into apply a Jackknife analysis to in order to quantify the uncertainty of my coefficients estimated by the logistic regression. I´m using a glm(family=’binomial’) because my independent variable is in 0 - 1 format. My dataset has 76000 obs, and I´m using 7 independent variables plus an offset. The idea involves to split the data in let’s say 5 random subsets and

Hi, First, thanks to those who helped me see my gross misunderstanding of randomForest. I worked through a baging tutorial and now understand the "many tree" approach. However, it is not what I want to do! My bagged errors are accpetable but I need to use the actual tree and need a single tree application. I am using rpart for a classification tree but am interested in a more unbaised

Dear all I have used the lda and predict functions to classify a set of objects of unknown origin. I would like to use a jackknife reclassification to assess the degree to which the outcomes deviate from that expected by chance. However, I can't find any function that allows me to do this. Any suggestions of how to generate the jackknife reclassification to assess classification accuracy?