similar to: Command Line Expressions

Suppose I have a nx2 matrix of data, X, the following code generate density estimation for each column and plot them denlist <- apply(X, 2, density) par(mfrow=c(1,2)) lapply(denlist, plot) Does anyone know how to change the main title of each density plot to "var 1", "var 2" by passing optional argument "main"? I've tried lapply(denlist, plot,

---- Original message ---- >Date: Thu, 14 Jul 2005 19:50:41 +0200 >From: Uwe Ligges <ligges at statistik.uni-dortmund.de> >Subject: Re: [R] East Asian language >To: Nan Lin <nlin at math.wustl.edu> >Cc: r-help at stat.math.ethz.ch > >Nan Lin wrote: > >> Dear all, >> >> I just installed R 2.1.1. The installation program automatically

Hello, I maintain a self-made R package under my own home directory "~/.R/library" on our university computer net (SunOS 5.8 system). After updating R to 1.8.0, I found other people cannot access my package any more. They got this error message: Error in library(tb, lib.loc = "/home/fs.hpc/43/zhengp1/.R/library") : This is not a valid package -- no DESCRIPTION exists I set

Hello all, Would anybody tell me how to present spatial-temoral point processes in R, for example, I'd like to plot the spatial points in the sequence of their time domain? Cheers -- Pingping Zheng Department of Mathematics and Statistics Fylde College Lancaster University Lancaster LA1 4YF UK

Dear R people The function below should be decreasing, convex, and tend to zero when x tends to infinity. curve((1-pnorm(x))/dnorm(x),from=0, to=9) >From the plot we see that for x between 8.0 and 8.3 the function is fluctuating. As far as I understand, this is due to the function pnorm() not being sufficiently accurate in the tails. I am using pnorm() in a way that has probably not been

Hi, Suppose I have created something like this in R: foo <- "myfoo" and I want to find out the number of character in foo (in other words, R should return 5 since "myfoo" has 5 charactors. How can I do it? I tried: length(foo) but it returned 1. Cheers, Kevin ------------------------------------------------------------------------------ Ko-Kang Kevin Wang

After "image" and "persp" would it be interesting to als have contour as a generic function? Cheers Paulo Justiniano Ribeiro Jr Dept Maths & Stats - Fylde College Lancaster University Lancaster LA1 4YF - U.K. e-mail: Paulo.Ribeiro@est.ufpr.br http://www.maths.lancs.ac.uk/~ribeiro -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-

Derar R-people I have troubles understanding what the function sort.list() is doing. On the homepage it says that it returns a permutation which rearranges a vector into ascending or descending order (like order() but on a vector instead of a sequence). > sort.list(c(0, 2, 10, 11, 4)) [1] 1 2 5 3 4 which does not make sense to me. In fact I am getting the same (non-sensical) result using

Dear R people We have two packages, where the first package (geoR) is required by the second (geoRglmm). Both packages have functions calling C-code via .C(). We would like to call C-functions included in the first package from within the C-code in the second package. Is this possible? An appropriate header file was included in the src directory of the second package but this alone did not

Dear R-devel During the last half a year I have several times encountered the following problem with optim() when using method= "L-BFGS-B". The function return a value which is clearly not the maximum (seen from printing the value each time the function is called). Some output is shown below. A few things I have observed (as I remember it): a. The problem seems to occur when the

Dear R-users, I wondered if anyone has extended GLMMGibbs to include correlated random effects, and if so, whether they would be willing to let me use their code? Jonathan Myles has no plans to extend glmm in this manner within the foreseeable future. With thanks, Patty -- -------------------------------------------------------------------------------- Assoc Prof Patty Solomon

Hello, everybody: I will like to know if anybody else is working with the sgeostat library and where can I find some other fuctions to fit more variogram models (like a kernel one). Thank you very much for you help Kenneth Cabrera Universidad Nacional de Colombia Sede Medellin -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read

Dear all, I wish use the "est.variogram" function under R. My question concerns the determination of number of lags. To make this calculation, i determined the biggest and the smallest distance between 2 pairs of points. I built the following tolerance : tol=mindistance/2 and i defined 1 lag as : 1lag=mindist+tol. Then i made a buckle the step of stop of which the maximal distance.

Dear all, I'm new on R and this mailing list. We work on spatial rainfall estimation with R and Grass. We have a problem with the krige function from the sgeostat package. We would like to limit the estimated area with a polygon limit. I use a 50 points polygon to describe my work area. The krige function work quiet well without limit. But if I use this option I have the following error

Dear R-developers There is a discrepancy between the variables found by typing names(options()) and the variables described on the help page for the function ``options''. The following variables are in names(options()) but not described on the help page : "error.messages" "repositories" "show.coef.Pvalues" "help.try.all.packages"

Hi, I have problems when calling fft_factor and fft_work directly in my C-code. Martin Maechler already told me that I should not do this -- anyway, the code works fine if the length of the vector is a multiple of 2,3,5. If the length of the vector has different form I get, in some cases, a segmentation fault when work and iwork are freed (I use C's malloc to allocate memory for work and

Hello R-Users, I have I simple question that I could not solve: How to extract the elevation and the slope values from a linear model (lm) separately? coef(model.lm) gives both of them. Thanks, Antonio Olinto -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info",

Hi, Is there a simple command to pause a script until the user hits return? Alternatively, is there a command which prompts the user for some input and assigns the input into a variable? dave +-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+ |David Richmond It works on a | + Dept. of Sociology complex scientific + |Saint Mary's College

I'd like to make some polar plots with R, but I can't seem to find anything designed to help with that. Is there anyway to get a real polar mode, where instead of x and y you have r and theta? How about a way to draw a circle? -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send

Dear All, I am trying to do some simple mixture analyses. For instance, I have a sample of n observations and I suspect they come from two different exponential distributions with parameters rate1 and rate2, respectively. So, I want to estimate rate1, rate2, and the proportions of both kinds of individuals in the sample. I had a look at the packages mda and mclust, but they do not seem to do this