similar to: export figure by pdf command

Hi all mailing listers, I wanna get the last several elments of vector. e.g. x <- c(1,2,3,.....,78, 79, 80) How can I implement to assign last three elements to y, y <- c(78, 79, 80) ? In Matlab, It can easily achieve by y=x(end-2:end) Thanks Huaping Wan [[alternative HTML version deleted]]

I would like to plot two data sets (frequency (y-axis) of mean values for 0-1(x=axis)) on a single histogram for comparison. The hist() only allow the overlay of two histograms, and although barplot() allows beside=TRUE, it does not show frequency values (like hist) but rather all of the values. Is there any way that I can use the hist() to plot two data sets similar to the barplot(). Any help or

Hi All, When carrying out hist(samples,breaks=50,probability=TRUE), the column values are considerably greater than 1, which seams very unreasonable. The plot is attached. I think the column value of the hist plot should correspond to x$counts/sum(x$counts) (x=hist(samples,breaks=50,probability=TRUE)). The size of data is a little bit larger, causing failure of uploading. If you need the

I need to repeat a function many times, with differing parameters held constant across iterations. To accomplish this, I would like to create a list (or vector) of parameters, and then insert that list into the function. For example: q<-("l,a,b,s") genericfunction<-function(q){ } ###### The equivalent code would of course be genericfunction<-function(l,a,b,s){ } Any help

Hi, The following script (which I did not develop) is used to calculate and plot a skewed normal curve. The script currently requires the user to input six parameters, rather than reading these directly from a file. I've been spinning wheels here, trying to figure out how to modify the script to automate it. I have four data sets, each in excess of 300 records that I need to process. My

Hello This is Elaine. I am using package lattice to generate boxplots. Using Richard's code, the display was almost perfect except the outlier shape. Based on the following code, the outliers are vertical lines. However, I want the outliers to be empty circles. Please kindly help how to modify the code to change the outlier shapes. Thank you. code package (lattice) dataN <-

Hi Listers, I am producing some graphics that the commands are in a FUNCTION... The problem is that I end up viewing just last graphic and in my FUNCTION there are 4 graphics with the PAR command function... Like those below... How do I view/get the other 3 graphics? Any help? Thanks in advance... histogram<-par(mfrow=c(1,2)) hist(rw_mean_app,main='Bootstrap Method

Hi all. I want to generate a sequence of n plots and save them into a single PDF file, one plot per page. From the R docs and other sources I gather the basic way to do this is save plot 1 into a file then append the 2:n plots to the same file.  This code shows my basic approach, but for some reason only the last plot is saved into the pdf. I've tried different variations (e.g. using onefile

Dear R-help listers, I am new. I just want to get helps on how to plot two histograms overlapped in the same plane coordinate. What I did is very ugly. Could you please help me to improve it? I want to got a plot with semi- transparent overlapping region. And, I want to know how to specify the filled colors of the different histograms. I also prefer other solutions other than ggplot2. Many

Dear listers, A quick question about breaks in hist(). The histogram is highly screwed to the right, say, the range of the vector is [0, 2], but 95% of the value is squeezed in the interval (0.01, 0.2). My question is : how to set the breaks then make the histogram look even? Thanks in advance, Leaf

Hi everyone, I'm working on a script trying to use foreach %dopar% but without success, so I manage to run the code with foreach %do% and looks like this: The code is part of a MCMC model for projects valuation, returning the most important results (VPN, TIR, EVA, etc.) of the simulation. foreach (simx = NsimT, .combine=cbind, .inorder=FALSE, .verbose=TRUE) %do% { MCPVMPA = MCVAMPA[simx]

Hi, I am running a model with count data and one categorical predictor (simple model for me to understand it fully), I did in R a glm like this: glm(Recruitment~Depth, family=poisson). I get the coefficientes and confidence intervals and all is ok. But then I want to do the same model with Bayesian stats, here is my code: model { for (i in 1:232) { Recruitment[i]~dpois(lambda[i])

Dear expeRts, I would like to colorize the backgrounds of a pairs plot according to the respective panel number. Here is what I tried (without success): count <- 0 mypanel <- function(x, y, ...){ count <<- count+1 bg. <- if(count %in% c(1,4,9,12)) "#FDFF65" else NA points(x, y, cex=0.5, bg=bg) } U <- matrix(runif(4*500), ncol=4) pairs(U, panel=mypanel) I

Hi listers, I want to make lattice plots xyplots with the indication of legends inside each panel with only the points and the lines actually ploted inside each given panel according to the group(ing) factor. The code below shows what I have achieved so far and I hope will make clear what I want to have. It seems to me that my solution is a very "dirty hack" and there certainly is

Hi, I need to generate some figure using the Arial font as a requirement for PLoS. Following their guidelines, I have converted the windows font files arial.tff files to .afm using tff2afm (exec file from MikTeX), but when I try to generate a postcript file, the postcript device does not recognize the .afm files. The code I use to test the new font is the following:

> > Leaf Sun wrote: > > The histogram is highly screwed to the right, say, the range > > of the vector is [0, 2], but 95% of the value is squeezed in > > the interval (0.01, 0.2). I guess the histogram is as you wrote. See http://web.maths.unsw.edu.au/~tduong/seminars/intro2kde/ for a short explanation. > -----Original Message----- > From: Berton Gunter

i calculated the density and wanna do something like this separate in 0-19-29-39-49-59-69-79-99 and put in these spaces 8 densities .. 0.something i have the frequency in % and divided already in 20 or 10 to get the density i tried and tried..made breaks vector to separate but couldn't put the other vector with the frequency density onit directly anyone know how to do it?? tks

I'm attempting to refactor an R script that does a lot of plotting, among other things. Ideally I'd like to do something like setup # does pdf(...) for each part of input { plot(process(part)) } cleanup # does dev.off() but have problems: 1 I'm plotting to PDF, so everytime I dev.off() creates a new file, and I want everything in one file (as does my boss :-) 2 I'm

Hi list, When exporting to PDF a graph with a legend, in the final PDF, the text is going beyond the legend box. > dev2bitmap("test.pdf", type="pdfwrite", h=6, w=6) The legend looks OK on the screen. I noticed that the size of the legend box depends on the size of the screen window, which is not the case for other graphical parts (text of the legend, title, axis

Buenos días: Sigo cabezudeando con el tema de los gráficos; a ver si a alguien se le ocurre alguna idea... Supongamos este código: library(plotrix) graphics.off() windows(width=8.27, height=11.69) plot(3:10) #plot(3:10, main="Axis break test") par(xpd=TRUE) text(4, 11,"Axis break test") axis.break() axis.break(2,2.9,style="zigzag") mtext("Texto