similar to: error in bargraph.CI {sciplot}

Hello, When using bargraph.CI in package sciplot can the bars for each group be different colors? How do I select the color for each group? When I use this instead of the default (SD vs SE): bargraph.CI(x.factor = dose, response = len, data = ToothGrowth, ci.fun= function(x) c(mean(x)-sd(x), mean(x) + sd(x)) ) Am I getting 95% CI bars? Thank you kindly, Michael Just

Hi there, I am a beginner. I would like to change the error bars in the bargraph.CI function from the default (se) to (sd). The help file says ci.fun= function(x) c(fun(x)-se(x), fun(x)+se(x)) Is there a simple way of telling the function what (x) precisely is - I already define in in the of the bargraph.CI function and assume that is should be able to use that information. cheers, Herwig --

Hi there, I wonder if anyone can help me. I'm trying to use bargraph.CI in the Sciplot package when there is a missing combination of the factor levels. Unfortunately the standard errors on the plot do not appear to be correct. Consider an analysis consisting of two factors A and B. When all factor level combinations are present all appears fine: library(sciplot) #all data

I'm trying to create a bargraph of means with standard error bars using the function bargraph.CI (in the sciplot package). Like this: bargraph.CI(x.factor, response,data,xlab, ylab, par(family="serif"),font=11) However, an error message comes up when I try to use the par funtion. Does the character font/style need to be changed in a different way for this kind of plot? Or am I

hi List and Manuel, I have encounter the following problem with the function "lineplot.CI".? I'm running R 2.10.1, sciplot 1.0-7 on Win XP.? It seems like it's a scoping issue, but I couldn't figure it out. Thanks! ...Tao > lineplot.CI(x.factor = dose, response = len, data = ToothGrowth)??? ## fine > lineplot.CI(x.factor = dose, response = len, data = ToothGrowth,

Is it possible to add standard error bars to the means on interaction plots? Thanks Murray -- Murray Pung Statistician, Datapharm Australia Pty Ltd 0404 273 283 [[alternative HTML version deleted]]

Hi all, I I need to manipulate the x axis in a stripchart. I will use one of the data sets included in R to explain what I need to do. attach(ToothGrowth) stripchart(len[supp=='VC']~dose[supp=='VC'], vertical=TRUE, group.names=c('A','A','A')) stripchart(len[supp=='OJ']~dose[supp=='OJ'], add=TRUE, vertical=TRUE, at=c(1:3)+.1,

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 =

I'm trying to plot a scatterplot of two variables using pch to plot different characters based on a third factor. Here is my example > data("ToothGrowth") > attach(ToothGrowth) > levels(supp) [1] "OJ" "VC" > plot(len ~ dose,pch=as.numeric(supp)) > legend(locator(1),pch=as.numeric(supp),legend=levels(supp)) The command as.numeric(supp) returns 2 2

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.?

Hello: I'm having quite a bit of difficulty adjusting the x-labels using bargraphCI. I've tried using text and srt=45 to rotate the labels or mtext for 2 lines to break up the labels. However, using either method, I cannot line up the labels with the midpoints of the bars (they line up with some sort of tick mark that is off the midpoint of the bars). Any suggestions would be greatly

Dear list, How can I add the mean and standard deviation to each of the boxplots using the example provided in the boxplot function? 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",

Dear all How does one resolve the variance between the actual availability of data sets in the default implementation with those mentioned in the documentation? I am unable to attach some of the datasets, even though help() is available for the same datasets. For example, "ToothGrowth" is a dataset that is widely mentioned in the documentation as being part of the base install.

Dear all I have written a function to perform a very simple but useful task which I do regularly. It is designed to show how many values are missing from each variable in a data.frame. In its current form it works but is slow because I have used several loops to achieve this simple task. Can anyone see a more efficient way to get the same results? Or is there existing function which does

Hello, I'm trying to sort my bargraph.CI plot (function like barplot in the SCIPLOT package) by the mean of the response variable. Does somebody have a trick for it. Thank you. Romain Mayor, PHD student.

Hi, sorry for bothering your guys. I will trying to make some nice graph using boxplot. when I check the help file of boxplot, there is a sample code as: boxplot(len ~ dose, data = ToothGrowth, add = TRUE, boxwex = 0.25, at = 1:3 + 0.2, subset = supp == "OJ", col = "orange") legend(2, 9, c("Ascorbic acid", "Orange juice"),

Hello, I have a couple more questions about boxplots. In the books I've read on statistics (I'm not a statistics expert but just a poor engineer trying to help his wife with the analysis of the data she gathered for her medieval history PhD), I seen two kinds of box plots: one using quartiles to determine the height of the rectangles and whiskers which is what the standard R boxplot does,

Hi, I would like to have lineplot.CI and barplot.CI to actually plot confidence intervals , instead of standard error. I understand I have to use the ci.fun option, but I'm not quite sure how. Like this : > qt(0.975,df=n-1)*s/sqrt(n) but how can I apply it to visualize the length of the student's T confidence intervals rather than the stdandard error of the plotted means

I'm using xyplot to generate a trellis plot with each panel containing a scatterplot and a best fit line. Is it possible to write the slope of the best fit line in each panel? -- Alan S Barnett <asb at mail.nih.gov> NIMH/CBDB

Dear all, I want to plot 8 chemical variables in the same plot. I have used the interaction.plot command: > interaction.plot(speclong$Date,speclong$time2,speclong$outcome,xaxt="n",type="l",pch=20,xlab="", log="y", ylab="Concentration (ug/L)",trace.label="",col=rainbow(8)) However, what it gives me is dotted lines and the colours are