similar to: Please begin a new subject if your post isn't a response to a previous one

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

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

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

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

You can get help on the whole package by > help(package="multtest") which is likely pretty close to what you want. There's the index etc for the package on the web as well. You can also just look in the package's installation directory. If it's loaded you can do an ls(2) say if it loaded in position 2, to get all objects in the package. Reid -----Original Message-----

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

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

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

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

It''s gotta be Rascal. There''s no fixing Victoria - that''s a re-write from the ground up. And Jake a recipe for disaster; Class Jake < Bat include Pig attr_accessor :one_eye . . end That pig ''ll never fly. ---------- Forwarded message ---------- From: mongrel-users-request at rubyforge.org < mongrel-users-request at rubyforge.org> Date:

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

help.search("Kalman") says to look at help(KalmanLike, package=ts). Andy -----Original Message----- From: Mohamed A. Kerasha [mailto:mohamed at engr.uconn.edu] Sent: Tuesday, November 19, 2002 9:27 AM To: r-help at stat.math.ethz.ch Subject: [R] Kalman Filter Hi all, Does any one know if there is Kalman Filter code or library in R. Thanks, Mohamed.

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

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