similar to: ggplot2 - regression statistics how to display on plot

Happy Friday fellow R users. I need some help - I am trying to make a graph using ggplot 2 of some lead isotope ratios. Normally, the isotope mass number appears as a superscript before the chemical symbol. However, I cannot figure out how to do this ggplot2's axis labels. The closest I have come is using "." In front of the numbers I am superscripting e.g: ggplot(....

Greetings - I am working on some data that contain chemical names with air concentrations, and I am creating a data frame with date/time and each chemical having its own column. However, these are organic chemicals (e.g. 1-butene, 2,3,4-trimethylbenzene etc). The package I am going to be using the data with is openair, and many of the great functions require you to specify a column name which

Hi - I am trying to find and describe missing data in a time series. For instance, in the library openair, there is a data frame called "mydata": library(openair) head(mydata) date ws wd nox no2 o3 pm10 so2 co pm25 1 1998-01-01 00:00:00 0.60 280 285 39 1 29 4.7225 3.3725 NA 2 1998-01-01 01:00:00 2.16 230 NA NA NA 37 NA NA NA 3 1998-01-01 02:00:00

Hello, I was interested in trying to write an R script to calculate a UCL for a lognormal distribution using the Chebyshev Inequality — MVUE Approach (based on EPA’s guidance found in http://www.epa.gov/oswer/riskassessment/pdf/ucl.pdf). This looks like it should be straight forward, but I am need to calculate an MVUE for the population mean and an MVUE for the population variance, which requires

Hi all- So sorry to bother you all with something pretty basic. I am trying to add the lines method output from svysmooth to a svyplot with style="grayhex". However, the line either appears in the wrong place or if I am running in R Studio it causes the system to crash. I know this is something to do with Lattice graphics, but for the life of me I can not figure out how. Dr. Lumley

Hello All, I have plotted a scatter plot in ggplot2 and added a regression line and a regression equation. But the processing is very very slow. One reason might be because I have so many data pairs. Is there any way to speed up this code ? I need to create a multiple layout as well. The code I have used is as follows: setwd("C:/Users/jzd0009/Documents/R software") mydata <-

This is what I have now so far. p=ggplot(data = test, aes(x = YEAR, y = TOTAL, colour = TREATMENT)) + geom_point() + geom_smooth(method = "lm", se=FALSE) + facet_wrap(~COUNTRY) > p +scale_x_continuous(limits=c(1,4)) http://r.789695.n4.nabble.com/file/n3626510/graph.gif I would also like to: 1.) change the headline for the faced wraps ?high? and ?low?. Is there any other way of

Hi there, I would like to make a 'bumps chart' like the ones described e.g. here: http://junkcharts.typepad.com/junk_charts/bumps_chart/ Purpose: I'd like to plot the proportion of people in select countries living for less then one USD pr day in 1994 and 2004 respectively. I have already constructed a barplot - but I think a bumps chart would be better # The barplot and data

hi there.. i got a problem with ggplot2. here my example: library (ggplot2) v1 <- c(1,2,3,3,4) v2 <- c(4,3,1,1,9) v3 <- c(3,5,7,2,9) gender <- c("m","f","m","f","f") d.data <- data.frame (v1, v2, v3, gender) d.data x <- names (d.data[1:3]) y <- mean (d.data[1:3]) pl <- ggplot (data=d.data, aes (x=x,y=y)) pl

How do I add to a trellis plot the best fit line from a robust fit? I can use panel.lm to add a least squares fit, but there is no panel.rlm function. -- Alan S Barnett <asb at mail.nih.gov> NIMH/CBDB

Hi - So when I run the following, I get a strange formatting output with MikTeX, and I am unsure if the behavior is due to R, Hmisc, or MikTeX or both: dfr <- data.frame(x=rnorm(400),y=sample(c('male','female'),400,TRUE)) latex(describe(dfr)) What happens is that the x column is summarized in the output of latex correctly. Then the y summary ends up looking like:

VR, James James T. Durant, MSPH, CIH Statistics and Data Analysis Team (proposed) Science Support Branch (proposed) Division of Community Health Investigations (proposed) US Agency for Toxic Substances and Disease Registry Mailing Address: 4770 Buford Highway, MS F-62 Atlanta, Georgia 30341 Telephone: 770-488-0668 (office) _____________________________________________ From: Lewis,

Putting names on a ggplot Can anyone suggest what I am doing wrong here. I am plotting daily temperatures at Ottawa Ontario for 2008 broken down by months, I seperate them by lines and want to put the names of the months at the top of the chart ( with in the graphing area) Everything is working as I want until I try to add the names of the months. I did something similar a few days ago and I

Hello! I am using ggplot2: ggplot(myplotdata, aes(x=att_levels, y=WTP)) + geom_bar(stat="identity",fill="dark orange",colour="black", alpha = 1,position = "identity") + geom_text(aes(label=WTP),colour="black",size=4,hjust=1.1,position='dodge') + coord_flip() + xlab("") +

Hello I'm having trouble figuring out how to use the output of "segmented()" with a new set of predictor values. Using the example of the help file: ??set.seed(12) xx<-1:100 zz<-runif(100) yy<-2+1.5*pmax(xx-35,0)-1.5*pmax(xx-70,0)+15*pmax(zz-.5,0)+rnorm(100,0,2) dati<-data.frame(x=xx,y=yy,z=zz) out.lm<-lm(y~x,data=dati) o<-## S3

Hi there, I'm trying to find out the command to stop clipping to plot region in ggplot. I have a bar chart (axis flipped) with labels on the bars, but the labels are clipped at the plot region box. I know it's possible to turn this off for base and lattice, but how about ggplot? par(xpd=NA) The ggplot2 manual doesn't seem to mention it, and a search for ggplot and

Hi all-- I am in the process of helping colleagues write up a ms in which we fit zero-inflated Poisson models. I would prefer plotting the rate ratios and 95% CI (as I've found Gelman and others convincing about plotting tables...), but our journals usually like the numbers themselves. Thus, I'm looking at a recent JAMA article in which both numbers and dotplot of RR and 95% CI are

?? Hola. Estoy intentando ordenar el gráfico tal que las columnas se apilen en el siguiente orden: debajo, los EPI; los SemiEPI en el medio y los NoEpi arriba de todo. Estos son mis datos: > tabla.barras Rangos Clase Asistentes 1 De 45d EPI 405 2 De 1a EPI 812 3 De 2a EPI 1639 4 De 3a EPI 2202 5 De 4a EPI 1994 6 De 5a EPI

Hola a tod en s, Este es el código principal que utilizo para diseñar mi zona de estudio (cuenca): png("mapa_zona_estudio.png", units="in", width=9, height=8.5, res=300) # Vamos a pintar los distintos instrumentos de medida sobre el mapa (junto con éste) y añadimos leyenda ggmap(mapa, extend='device', legend="left", base_layer=ggplot(datos, aes(x=lon,

Hello! I am using ggplot2 (see the code below) to plot the data in 'myplotdata'. The first column of 'myplotdata' is called "att.levels" and contains strings; the second column is called "WTP" and contains numeric values. Notice - I use 'coord.flip()' The command aes(x=att_levels, y=WTP), if I understand correctly, sorts things alphabetically based on