similar to: boxplot subsets

Hi, I am looking at the effects of two explanatory variables on chlorophyll. The data are an annual time-series (so are autocorrelated) and the relationships are non-linear. I want to account for autocorrelation in my model. The model I am trying to use is this: Library(mgcv) gam1 <-gam(Chl~s(wintersecchi)+s(SST),family=gaussian, na.action=na.omit, correlation=corAR1(form =~

I''d like to make the following changes (differences are from R1.0.1): boxplot.default() 1c1 < function (x, ..., range = 1.5, width = NULL, varwidth = FALSE, --- > function (x, ..., range = 1.5, width = NULL, varwidth = FALSE, boxwex=0.8, 37c37,38 < bxp(groups, width, varwidth = varwidth, notch = notch, --- > bxp(groups, width, varwidth = varwidth,

Hello, I'm currently using a boxplot to visualize data for three different models. As I have three models, I'm plotting three parallel boxplots for each factor. This works fine - what I need now is a line connecting the medians of each boxplot of each model. I want to do this in order to visualize the trend that one of the models exhibit. Basically, I want to plot a curve for each model

A student asked me if it was possible to draw boxplots where the boxes themselves were grouped: I was able to hack boxplot.formula() to do the right thing, more or less, by incorporating an argument (space) and some code from barplot. with the extended boxplot.formula() below, the following commands "do the right thing" (produce boxes grouped by levels of the second factor): r =

Hi all, I would like to ask a question with regard to R graphics. I am currently trying to add data points in boxplots using points(). In my script, I add up boxplots three times using add=TRUE. I use points() every time I add boxplots, and I see data points. However, I see all data points are clustered at three specific points on x-axis; the locations of the data points are not overlapping with

Thank you for your guys reply for my previous question. But I got one more question about the boxplot. With the code in the R-help: boxplot(len ~ dose, data = ToothGrowth, boxwex = 0.25, at = 1:3 - 0.2, subset = supp == "VC", col = "yellow", main = "Guinea Pigs' Tooth Growth", xlab = "Vitamin C dose mg", ylab =

Hello, I have created two boxplots with following R code. There is one outlier in B group. The outlier is 33. But the all other data are between 0 to 4. How can I skip y-axis around 5 to 25, and expand 0-4 for this case. Also I want keep the outlier in my boxplot. I want my boxplot look like the second one, keep the outlier, and make an interrupt of y-axis from 5 to 25. Thanks, Jianghong

Hello, I am trying to plot two sets of boxplots together. These are estimates of two experiments and?seven?factors. The results of the two experiments I want to plot as boxplots stacked to each other. Therefore I plot first the results of the first experiment; and next with the add option the second set of boxplots. The boxplots are plotted at 'at = 1:7 - 0.15 for the first experiment and

Suppose that I have data on three species for a variable and datasets from two time periods. I want to make a boxplot of the first dataset and then add the second using 'at = ' and 'add = T' as in the example for 'boxplot.' Since the boxes are paired by species, I want to do is have the x labels be centered between the boxes. I'm doing this now with mtext and entering

Dear List, ? I am creating a boxplot with two subsets, very similar to the example by Roger Bivand at ?boxplot (reproduced below).? I am trying to change the labels on the x-axis to have one number to cover both subsets.? I can do this in other plots by using axis=FALSE followed by a separate axis() command.? I have also tried variations in the names= argument but can't get it to work.?

Hi All, I'm using the following code to produce some stacked bar graphs. *setwd("C:\\Users\\Tinus\\Documents\\NMMU\\R\\Seamounts")* *SChla <- read.csv("SM_Chla_data.csv")* * * *#Extract mean values from data file* * * *Coral <- SChla[185:223,] #Reduce SChla to Coral only* *coral <- with(Coral , aggregate(cbind(Pico, Nano, Micro), list(Depth),FUN=mean))*

Dear R users, x-values (EI) = Adw, EG1, LA1, Ad1, LA2, LA3...(14 levels, insect stages) y-valus = antpop within the boxplot function x-values are ordered alphabetically Idea: x-values ranked by list order (insect stage: Egg stage 1 is followed by Larvae 1 and not by Egg stage 2 as it would be in an alphabetically order) Problems with the order(tapply()) function: variable lengths

Hi R-user, I am new with R and I have a problem with the map and mapproj fonctions : with the following code : >test_proj=mapproject(LONG_d01_vec_t1, LAT_d01_vec_t1, projection="lambert", parameters=c(30,60)) >longitude_vec=test_proj$x >latitude_vec=test_proj$y >longitude_mat=matrix(longitude_vec,ncol=e_sn_d01) >latitude_mat=matrix(latitude_vec,ncol=e_sn_d01)

Hello, I would like to evaluate the relationship between flows and phytoplankton abundance (or Chlorophyll a concentrations) using a changepoint analysis.? Specifically, I have two study questions: Study Question 1: Are there certain flow thresholds that result in spikes in phytoplankton abundance? Study Question 2: Are the duration of certain flows important for phytoplankton abundance (e.g.,

Hi, im a newbie with very wobbly coding abilities. Tearing my hair out over getting the boxplot i want... I have a dataset called 'eagle' which consists of year (2011 or 2012), month (jan - dec), roof (TT6, TT13 or BARE) and temp (the continuous variable that i want to plot). So i want boxplots of the three roof treatments in every month organised in chronical order along x axis 2011 -

Good Morning, I am trying to get side by side boxplots of two groups on the same variable. The last item under ?boxplot led me to some useful code. I use "boxwex" to make the boxes narrower, "at" to shift them over and "add" to draw them both on the same graph. Something along the lines of: attach(group1) boxplot(Y~X, col="blue", boxwex=.4,

Sorry for disturbing you with another newbie question! I have a data frame about coastal waters quality parameters: for some parameters (e.g. NH3) I have only 1 observation for each sampling station and each sampling date, while in other cases (chlorophyll) I have 1 obs for each meter-depth for each station and date. How can I select only the max chlorophyll value for each station/date? example

Hi, I would be very grateful for some advice. I have read the help pages for Date, strptime, etc. All examples seem to use some version of day month year as date format. However I have Weekly composite data so ONLY want to input the dates as Week.Year (eg 35.2011). strptime seems to show this is possible using %W for week (UK convention) and %Y for year. My data is in a df called chlorophyll and

Hi, set.seed(45) test1<-data.frame(columnA=rnorm(7,45),columnB=rnorm(7,10)) #used an example probably similar to your actual data apply(test1,2,function(x) sprintf("%.1f",median(x))) #columnA columnB # "44.5"? "10.2" par(mfrow=c(1,2)) lapply(test1,function(x) {b<-

I was poking around in the Gentoo Linux repository yesterday, trying to satisfy my curiosity about how many programming languages there are in "common use" today. I found one called "Eleven" which is quite strangely similar in intent to Rails. The home page is http://eleven.sourceforge.net/. Eleven is much simpler than Rails, and there are some interesting ideas in it. To