similar to: Hmisc or Lattice plot with error bars, Depth (independent variable) on Y axis

I have been using the ECDF function in the Hmisc package to produce cumulative distribution function plots. The problem is that for small datasets the steps "look bad" (not my characterization but from the client). Is there a way to get the same information but smoothed? I have tried the densityplot (lattice), which gives a smoothed line, but this does not give the cumulative density.

I am using R for environmental forensics (determination of the sources and/or groupings in mixtures of organic chemicals in the field). The goal is to determine in there are groups of samples with similar/dissimilar compositions, and to assign samples to a potential source or a mixture of sources based on the composition (unmixing and source allocation). Typically there are 10 to 50 chemicals that

I am looking to create what I would call a "simple variation" on the boxplot. What I would like to do is to be able to plot the upper and lower confidence limits as the "box" and the 10th and 90th percentile as the whiskers. What I have done is write the code to create a dataframe, the columns of which are the mean, sd, 10th percentile, 90th percentile, lower confidence limit

Dear all I have a dataset of how metal concentrations change through time. I have made a plot of date versus metal concentration. However I want to add error bars in the plot. Could you help me? Thanks Maria [[alternative HTML version deleted]]

Greetings, I am running a buch of wilcox tests and need to be able to rapidly export the results into a csv file. I have attached example code as well as my attempts to get what I need. I have tried unlist,cbind,rbind etc but I am obvously missing something simple. FYI I am actually running about 50 WRS tests per dataset, this is just an example. Thanks 10^6 Mike AKCCR <-

I have written a few different summary functions. I want to calculate the statistics by groups and I am having trouble getting the output as a dataframe. I have attached one example with a small dataset that calculates summary stats and percentiles, I have others that calculate upper confidence limits etc. I would like the output to be converted to a dataframe with one of the columns as the

I am debugging some code and found a function that returns and error most of the time. I have issolated the problem to when it raises an argument to the 1/3 power but for the life of me I can not figure out why it is not working. i have gone through the FAQs and have found nothing (not to say I might have missed something). I am highly embarrased with my inability find the problem (my face is

I am using lattice to make figures as pdfs: trellis.device(device = "pdf",file = "Figure6.pdf",color = FALSE) I need to specify some blank space on the left-hand margins (the pages will be bound so we need about 0.5 inch)). I have tried a number of solutions but none seems to work (e.g. par.set). Can this be done when initiating the plotting device? Or is the some other way

I was preparing an e-mail for the help list and ran across a quandary. When asking for help it is useful to include the code/data so others can run your code and test it. I was running code on a data frame and wanted to include a small version of the data frame. The data frame was based on experimental data. What is the best way to do this? I didn't want to send an attachment so a wrote code

I have searched the archives and instructions high and low but have not found what I need. I have a dataframe named Data with columns: AdjResND0 - the numeric result Parameter - the classification (chemical name) Sort - the order I want the chemical names to appear, with leading zeros so all are the same length Client.Name - the sample name I am trying to generate a series of 1 page barcharts,

Hello! In a grouped Dotplot, is there any way to set the color of error bars to be the same as the corresponding symbols? Example data: require(lattice) require(Hmisc) data(barley) Dotplot(variety~Cbind(yield, yield+2, yield-2)|year, groups=site, data=barley) I experimented with changing trellis settings of plot.line (as mentioned in the Hmisc documentation) as well as col.line settings in

Dear R-helpers, (1) I can't figure out how to tell Dotplot that I want the colors of the CI bars to be the same as the colors of the dots. For example: t2a.ci3 <- data.frame(est = c(7, 20, 75), lower = c(-9, 0.5, 42), upper = c(22, 39, 109)) mypal <- c('skyblue3', 'mistyrose3') condNames <- c('noise', 'horn', 'marimba') Dotplot(1:2 ~

Estimados, tengo un grafico al que agregue una linea con abline(), este es el codigo ... library(Hmisc) # agregar una linea horizontal con el valor del angulo de contacto del control trellis.device(color = FALSE) xYplot(Cbind(ca.med[,3],ca.med[,5],ca.med[,6])~ca.med[,2], groups=Sol, data=ca.med, xlab=list(label="Concentration (mM)", cex=1.3), ylab=list(label="Contact Angle

Hello-- I have a data for small population who took 1 drug at 3 different doses. I have the actual drug concentrations as well as predicted concentrations by my model. This is what I'm looking for: - Time vs Concentration by ID (individual plots), with each subject occupying 1 plot -- there is to be 9 plots per page (3x3) - Observed drug concentration is made up of points, and predicted drug

A data frame has one factor, one date, and one numeric column. When I plot these using the default pch of the open circle (first attachment), xyplot(TDS ~ sampdate | she.s, data = sheep.cast, main = 'TDS in Sheep Creek', ylab = 'Concentration (mg/L)', xlab = 'Time') I see the higher concentration points toward the right on several panels. But, when I change to using a

Dear R users, today I've got the following problem. Here you are a dataframe as example. There are some SAMPLES for which a CONCentration was recorded through TIME. The time during which the concentration was recorded is not always the same, 10 points for Sample A, 7 points for Sample B and 11 for sample C Also the initial concentration was not the same for the three samples. I would like

I am looking for procedure that allow one to fit multiple distributions to a variable. For example, based on analysis of the data we suppose that the data can be represented by 3-5 normal distributions added together. I would like to be able to determine the mean, sd, and weight associated with each distribution and examine the improvement of the fit when 3,4 or 5 normal distributions are used as

I am trying to get the same text printed on each page of a multi-page series of bar charts. The text need to appear in the upper left-hand corner of the page, outside of the plot area. A watermark might be the closest analogy to what I am after This is what I have so far: PData <- na.omit(subset(TData,Matrix == "Product")) barchart(AdjResND0 ~ reorder(Compound, Sort) | Label ,

Estimada comunidad, estoy haciendo un grafico con la funcion xYplot (paquete Hmisc) para agregar barras de error a los puntos del grafico. Pero pasa algo que no me gusta: aparece el nombre de las series en el grafico. Alguien sabe como puedo eliminarlas ? no quiero usar ggplot porque ya tengo hechos todos los graficos con lattice. Aqui el codigo que uso y los datos (estan separados por tab, no

r-help lattice adepts: I have a question which is somewhat geospatial, so I posted to r-sig-geo rather than here: https://stat.ethz.ch/pipermail/r-sig-geo/2012-November/016757.html > summary: How to overlay a geographical map on each panel in a lattice > (or Trellis), e.g., of levelplot's? Note I am not inquiring about > creating choropleth maps[,] which Sarkar 2008 covers quite