similar to: histogram

Hello, I'm doing some graphics for a paper and a need customize such with filled region above the density curve. My attempts I get something very near what I need, but I don't solve the problem of use semitransparent filled. Below a minimal reproducible code. Someone has any idea? require(lattice) # toy data... dt <- expand.grid(A=1:2, B=1:3, y=1:50) dt$y <- rnorm(nrow(dt), dt$B,

I am trying to plot multiple densityplots on each panel, and using panel.groups to do some additional plotting (not included in the example) as in this example. library(lattice) thedata <- data.frame(x=rnorm(1200),class=rep(1:3,each=200,times=2), group=rep(1:2,each=100,times=6)) densityplot(~x | class,groups=group, xlab='x', panel.groups =

Hello Everyone,   Learning about R and its wonderful array of functions. If it's not obvious, I usually try to find out what a function stands for. I think this helps me remember better.   One function that has me stumped is "sink." Can anyone tell me if this stands for something?   Thanks,   Paul         __________________________________________________ [[alternative HTML

Hi I implemented the age-structure model in Gove et al (2002) in R, which is a nonlinear statistical model. However running the model in R was very slow. So Dave Fournier suggested to use the AD Model Builder Software package and helped me implement the model there. ADMB was incredibly fast in running the model: While running the model in R took 5-10 minutes, depending on the

I am not new to R, but I am new to .grib files. I am downloading some climate data and I would like to analyze it in R. R has a nice netcdf package, but I don?t see any package available to deal specifically with grib files. I see a few posts from other people using grib files in R. However, I was unclear if they used grib files in a different software program and then imported the data

Hi R users, What is the command for adding (or how do we add) country boundaries on r-spatial plots. Thanks ZABLONE OWITI GRADUATE STUDENT Nanjing University of Information, Science and Technology College of International Education Add: 219 Ning Liu Rd, Nanjing, Jiangsu, 21004, P.R. China Website: www.nuist.edu.cn ==================================================== [[alternative

... Both readLines() and scan() produce a number_of_lines x 1 vector; trying paste(s, collapse = NULL) leaves it unaffected. How can I concatenate vector elements (lines) into a single string? Thank you. -- View this message in context: http://r.789695.n4.nabble.com/How-to-read-contents-of-a-text-file-into-a-single-string-tp2069303p2069303.html Sent from the R help mailing list archive at

Hello, I have downloaded several times your program R-2.11.1 32.exe, after the download when I click on the exe file I always receive the following error: "The setup files are corrupted. Please obtain a new copy of the program". Could you help me to download English or Russian version for Windows (intel865). You may send a link. Best regards, -- Gulbanu [[alternative HTML

I am trying to create a heap of boxplots, by looping though a series of factors and variables in a large data.frame suing paste to constrcut the facto and response names from the colnames I thought I could do this using get() however it is not working what am I doing wrong? thanks Nevil Amos > sp.codes=levels(data.all$CODE_LETTERS) > > for(spp in sp.codes) { > > >

In converting some older code: ## builds an OWIN to generate a set of points on a grid within the OWIN object: bdry <- read.shape( "../../data/boundary.shp" ) outline.verts <- bdry$Shapes[[2]]$verts[1:bdry$Shapes[[2]]$Pstart[2],] outline.order <- nrow( outline.verts ) : 2 bdry.poly <- vector( 2, mode="list" ) bdry.poly[[1]] <- list(

Dear All I am using the histogram() function to plot two subsets from my data on the same output display, i.e. there are two histogram plots My code is: histogram( ~ Age | Date, layout = c(1, 2), xlab = "Age (years)", strip = FALSE, strip.left = TRUE, col = "black", border = "white", cex.axis = 1.1, family = "serif", cex.lab =

Hello all, I am using R to perform certain calculations on huge amounts of data. In short I need a function that does the mode function, ie returns the most common element. I looked at the mode function in R but it seems to return the type of the data element you give it. Does such a method exist? I have tried googling this to no avail as all the results lead me back to the mode function I do

hi I've used R for many years on windows machines, but have now acquired an Asus eee 1000 linux machine. In order to get the best out of the machine, I used the 'pimpmyeee.sh' script, to get the full KDE desktop. The version of Linux is Xandros, which I believe is a close relative of Debian, but sadly I have only a nodding acquaintance with Linux at present. Naturally I want to

Dear R experts, I have plot the PDF of t distribution with df = 74. curve(dt(x,df=74),from=-4, to=4) how can I shade the area under curve (for example, col="red") from t=+- 1.996? Thank you. Regards, CH -- CH Chan Research Assistant - KWH http://www.macgrass.com

What is the recommended way of requiring a certain version when loading a package (or, indeed, from R itself)? Perl has the require module version use module version require version use version constructs which is kind of what I am looking for (especially 'use' which is evaluated at compile time), but R seems to have lost the version= argument to require(). The best I have been able

Hello Everyone,   My company purchased S-Plus before it was bought out by Tibco. My understanding is that we own version 7.0 outright.   So far, I've been learning in R but thought I might also try working in S-Plus. My understanding is that S-Plus has some useful extra features. Another potential benefit would be the ability to purchase technical support, which I thought might help me to

Hey, I have a data frame x which consists of say 10 vectors. I essentially want to find out the best fit exponential smoothing for each of the vectors. The problem while I'm getting results when i say > lapply(x,ets) I am getting an error when I say >> myprint function(x) { for(i in 1:length(x)) { ets(x[i],model="AZZ",opt.crit=c("amse")) } } The error message is

Hi, I am encountering a strange error in POSIXlt... anyone got a clue? > as.POSIXlt("1982-01-01") Error in as.POSIXlt.character("1982-01-01") : character string is not in a standard unambiguous format > as.POSIXlt("1992-01-01") [1] "1992-01-01" > as.POSIXlt("1972-01-01") [1] "1972-01-01" > as.POSIXlt("1962-01-01")

Hi, Can we run EGARCH in R. If yes, I would be grateful if someone could tell me the R codes for running EGARCH model. Thanks. [[alternative HTML version deleted]]

Dear All, Is it possible to replace function with my own? I want to apply pca clustering, but to use some strange correlation function. I'm asking about replacing, say, mean() with new content of mean() and use standard other functions, which might use mean() as part. karsar