similar to: possible spam alert

Most mornings, I have over a screenful of messages mostly from R-help and it's very useful to have them threaded. However, the usefulness of threading is lost when posters reply to a message and then change the subject instead of creating a new message. People who don't have a mail client that can display email in threads are probably unaware that this sort of thing can happen in ones

All, I am using a sum of squared differences in the objective function of an optimization problem I am doing and I have managed to speed it up using the outer function versus the nested for loops, but my suspicion is that the calculation could be done even quicker. Please see the code below for a simple example. If anyone can point out a faster way I would appreciate it greatly. Thanks, Roger

I have 3 years of experience with R and have an interest in becoming a better programmer so that I might someday be able to contribute packages. Other than R, my only experience was taking Lisp from Daniel Friedman in the 1970's. I would like to learn either C or C++ for several reasons: To gain a better concept of object oriented programming so that I can begin to use S4 methods in R. To

Dear all, I have two vectors and connect the points by line segments: x <- c(1990,1994,1995,1997) y <- c(30,40,80,20) plot(x,y,type="l") I want to get the values of y's on these lines at missing x's, i.e., at 1991,1992,1993,and 1996. Is there a simple way to do this? Best regards, Remigijus mailto:remigijus.lapinskas at maf.vu.lt

Hello R-help, My vector aa has some NA's: > aa_c(NA,1) > aa=="NA" [1] NA FALSE > aa==NA [1] NA NA How can I detect NA's? In other words, I would like my answer to be TRUE FALSE. Best regards, Remigijus mailto:remigijus.lapinskas at maf.vu.lt -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help

Dear All, There is an example for the priceIts function (the its package) which does not work for me as expected. > ?priceIts > x1 <- priceIts(instrument = c("^ftse"), start = "1998-01-01", + quote = "Close") Error in validObject(.Object) : invalid class "its" object: Missing values in dates > x2 <-

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<-

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

Dear all, I have a function MYFUN which depends on 3 positive parameters TETA[1], TETA[2], and TETA[3]; x belongs to [0,1]. I integrate the function over [0,0.1], [0.1,0.2] and [0.2,0.3] and want to choose the three parameters so that these three integrals are as close to, resp., 2300, 4600 and 5800 as possible. As I have three equations with three unknowns, I expect the exact fit, i.e., the SS

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]]

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

Dear all, I got this error message > library(dprep) > mardia(Savg) Error in cov(data) : 'x' is empty But with the same data, I got > library(mvnormtest) > mshapiro.test(Savg) Shapiro-Wilk normality test data: Z W = 0.9411, p-value = 0.6739 What does the error message "Error in cov(data) : 'x' is empty" mean? Thanks a lot! Jiao

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

Hi, I am having some trouble trying to plot a histogram in more than one colour. What I want to do is, plot two vectors in the same histogram, but with different colours, for instance: > x <- rnorm(1000,20,4); > y <- rnorm(1000,10,2); Then I'd like to have x and y ploted on the same hist (I can do that already doing w <- c(x,y) then hist(w)) but the bars

I am running ess/emacs on Linux and have an R instance running on a remote Unix server. I would like to be able to direct input from my ess buffer to R (hope I am using the right lingo, I am new to emacs). The ess manual contains a section describing how to do just that, but a prerequisite is to install a lisp file ssh.el , which can be found at

Dear List, There has just been this question (as well as the answer to it) in the S list: > I would like to put some text in plots... but in french. How can I > write accented letters? > I suppose that there is a possibility similar as the one which > works in LaTeX, isn't it? What about R (for Windows and in East European languages)? Specifically, how can I give a title like

Is it possible to fit a lagged regression, "y[t]=b0+b1*x[t-1]+e", using the function "lag"? If so, how? If not, of what use is the function "lag"? I get the same answer from y~x as y~lag(x), whether using lm or arima. I found it using y~c(NA, x[-length(x)])). Consider the following: > set.seed(1) > x <- rep(c(rep(0, 4), 9), len=9) > y <-

Hi to all on the samba lista. I'm new on Samba and on the mailing list so forgive me for any misunderstading that I may commit. I have a Samba daemon running on a OpenBSD machine with only a directory shared. Anyone can read and write on that directory. However, on the same local network there are various Windows machines and the users of those machines are users of a Windows Domain. Is

Dear All, R 1.9.1, Windows When copying and pasting a few lines from the 'predict.Arima' help, I get an error message: > data(lh) > predict(arima(lh, order = c(3,0,0)), n.ahead = 12) Error in eval(expr, envir, enclos) : Object "xreg" not found On the other hand, the following is OK: > data(lh) > predict(arima0(lh, order = c(3,0,0)), n.ahead = 12) $pred Time

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