similar to: Interpolation

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

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

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

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

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

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

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

Hi, I have the install guide and the admin guide....nothing in either that I can see from the contents pages tells me how to create the first "server" ( I assume that is what I have to do?) Is there another doc Im missing? or a good URL for a howto on a redhat based machine? Also from what I can see the "free" version is cli only? and there is no virtual (vmware) appliance?

>From the <samba>/var/log/<client>.log file: doing parameter domain controller = AZPHXNT01 Unknown parameter encountered: "domain controller" Ignoring unknown parameter "domain controller" >From the smb.conf file: # Domain Master specifies Samba to be the Domain Master Browser. This # allows Samba to collate browse lists between subnets. Don't use

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

Hi, Iwrote 'pdfc' in analogy to 'lpr': pdfc <- function (file = "Rplot.pdf") { current.device <- dev.cur() dev.off(dev.copy(device = pdf, file = file, ...)) dev.set(current.device) print(paste(file, "generated.")) } lpr <- function (object, file = "Rplotlpr.ps", ...) { if (missing(object)) { current.device <-

Has anyone achieved a current successful build of R 1.8 or 1.9 or earlier for HPUX Tru 64 OSF 5.1? There is no binary version of R for HPUX Tru 64 OSF 5.1. The R admin manual mentions that the native make fails on "Alpha/OSF (aka Tru 64)" and gnu make must be used instead. There are problems reported with building R versions around 1.4 from 2002 on R-devel, some of which were fixed with

Hi, I'm looking to create a few basic world maps in R (for Windows). A while back someone asked for the same thing and was referred to the (*.tar) files at: ftp://ftp.mcs.vuw.ac.nz/pub/statistics/map/ which I downloaded. However, I am unable to install them, and clearly I'm not doing something properly (I'm your basic moron when it comes to installation, so please understand).

Dear all, I am sorry because my question is perhaps trivial and is certainly a detail. It's related to the output of some clustering package as mclust and e1071. The output is as follows : [1] 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 [38] 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 [75] 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1

Here is the answer to an old question of mine, regarding maps in R Hope this help. claudia ---------- Forwarded message ---------- Date: Mon, 25 Jun 2001 13:48:05 -0400 From: Kieran Healy <kjhealy at Princeton.EDU> To: Claudia Tebaldi <tebaldi at ucar.edu> Subject: Re: [R] maps in R Hi Claudia -- you can get the maps() library from here:

A student asked me about simulation he is doing. This simulation presented an interesting R programming puzzle that we did solve but I am not happy with the solution. I feel there should be a better way of doing this. Briefly, he is simulating births and deaths in a linear array of cells. Births and deaths are required to alternate. For a birth, one of the cells is selected at random to be the

Hello, I recall reading a thread months ago on this mailing list about handling very datasets under R, but I can't seem to find it. This has become particularly important recently, because I've been playing with a dataset containing information about every fatal car accident in the U.S. since 1975; in total, the relevant files are about 120 megs. I'd like to load all of these into R

Does anyone know what the units are for projected coordinates obtained using mapproj's mapproject function with an Albers projection? Thanks for any and all help! Buck Stockhausen *************************************************** * Dr. William T. Stockhausen * *************************************************** * Resource Ecology and Fisheries Management * *