similar to: How to calculate GLM least square means?

Hi, I have some data (lots of year,distance pairs), which I can straightforwardly boxplot: dists <- read.table("movedists.dat") with(dists,plot(as.factor(V1),V2)) If I want to plot these data as violin plots using vioplot, this approach doesn't work: > with(dists,vioplot(as.factor(V1),V2)) Error in Summary.factor(c(1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, : min not

I am trying to make horizontal violin plots. I have tried both vioplot and simple.violinplot, but both of them seem to not be willing to take the horizontal option. Is this correct, or am I just bungling it somehow? For instance, for vioplot (from the example shown, with the horizontal modification): > vioplot(bimodal,uniform,normal, horizontal=TRUE) Error in median(data) : need numeric data

Hi I am using violin plots (type of boxplots) and  I am trying to increase the font size in the plots. It looks like that the violin plots do not work as "normal" plots as the cex parameters are ignored. You can have a loot at the code below require('vioplot') data1<-rnorm(100) data2<-rnorm(10) data3<-rnorm(1000) vioplot(data1,data2,data3,cex.lab=2,cex.axis=2)

Hi, I have a list of vectors and am trying to coerce them into something that vioplot will take as groups of data to be plotted independently. Can someone nudge me into the right direction? Thanks, Joh

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",

I am creating boxplots from a dataframe and would like to add to the standard output a marker representing the value from a particular row in the dataframe. .....And I apologize if the solution is as trivial as it seems it should be. Thanks for any assistance. Keith

Dear Rgurus, To my knowledge the best way to visualize the distribution of a discrete variable X is plot(table(X)) The problem which I have is the following. I have to discrete variables X and Y which distribution I would like to compare. To overlay the distribution of Y with lines(table(Y)) gives not satisfying results. This is the same in case of using density or histogram. Hence, I am

Dear all I am trying to make a graphic with the "vioplot" package. I use the following code: library(vioplot) x1 <- GSMrxDL x2 <- WIFI x3 <- UMT vioplot(x1, x2, x3, ylim=c(0, 10), names=c("GSMrxDL", "WIFI", "UMT"), col="gold") title("NIS Strahlung", xlab="Sender", ylab="V/m") Now I want to scale the y-axis

We are using GAM in mgcv (Wood), relatively new users, and wonder if anyone can advise us on a problem we are encountering as we analyze many short time series datasets. For each dataset, we have four models, each with intercept, predictor x (trend), z (treatment), and int (interaction between x and z). Our models are Model 1: gama1.1 <- gam(y~x+z+int, family=quasipoisson) ##no smooths Model

I have been using the package vioplot to make boxplot/densityplots. Now I am looking for a variation of this, and I was wondering if someone could give me any tips on how to do what I want. I have nine groups of values, each group containing two sets of numbers which show the difference between two values, such as this: Start dev Stop dev -1 10 5 -2 0

Hi,? I have made violin plot with both ggplot2 and vioplot package with same data.? The results and code are as follows.? ################### ?Loading data? data.melt <- dget("https://gubox.box.com/shared/static/rirth0eym114afwyjxwe128sjzipzdym.txt") ################### ?Violin plot by library(vioplot) library(vioplot) data.use_11<- data.melt$value[data.melt$ident == 0]

I have read and reread the boxplot and the boxplot stats page, and I still cannot understand how and what boxplot shows. I realize that this might be due to me not knowing enough statistics, but anyway... First, how does boxplot determine the size of the box? And is the line inside the box the mean or the median (or something completely different?) And how does it determine how long out the

I'm attempting to produce something like a violin plot to display how y changes with x for members of different groups (My specific case is how floral area changes over time for several species of plants). I've looked at panel.violin (in lattice), which makes nice violin plots, but is really set up to work on a single variable - the area trace represents the frequency of each value of x

Hi First post and a relative R newbie.... I am using the vioplot library to produce some violin plots. I have an input CSV with columns off irregular length that contain NAs. I want to strip the NAs out and produce a multiple violin plot automatically labelled using the headers. At the moment I do this Code:? ds1 = read.csv("http://www.lecturematerials.co.uk/data/spelling.csv")

Using Trellis, am successfully setting up a number of panels (25) in which I have two box and violin plots. I would like to colour - one plot as RED and the other as BLUE (in each panel). I can do that with the box plots, but the violin density areas just take on one colour. My basic call is as follows: bwplot(rmsd ~ file | code, data=spread_data.filtered, panel = function(...,

Hello, I am trying to graph a violin plot like the example in panel.violin(lattice) with the singer data - see below Does anyone know how to change the color of the boxplot on the graph from blue to green? I tried a number of changes, but to no avail. Thanks for your help, John bwplot(voice.part ~ height, singer, panel = function(..., box.ratio) {

Is it possible to create weighted boxplots or violin plots in lattice? It seems that you can specify weights for panel.histogram() and panel.densityplot(), but not for panel.bwplot or panel.violin(). Please let me know if I've missed something in the package documentation. Thanks! -- Raphael D. Mazor Freshwater Biologist Southern California Coastal Water Research Project 3535 Harbor Blv

Hi, Can anyone tell me why I am not getting the correct intervals for fixed effect terms for the following generalized linear mixed model from HPDinterval: > sessionInfo() R version 2.4.1 (2006-12-18) i386-pc-mingw32 locale: LC_COLLATE=English_United States.1252;LC_CTYPE=English_United States.1252;LC_MONETARY=English_United States.1252;LC_NUMERIC=C;LC_TIME=English_United States.1252

I'm trying to sort a DATAFRAME by a column "ID" that contains alphanumeric data. Specifically,"ID" contains integers all preceeded by the character "g" as in: g1, g6, g3, g19, g100, g2, g39 I am using the following code: DATAFRAME=DATAFRAME[order(DATAFRAME1$ID),] and was hoping it would sort the dataframe by ID in the following manner g1, g2, g3, g6, g19,

Dear R users, I am having problems trying to fit quasipoisson and negative binomials glm. My data set contains abundance (counts) of a species under different management regimens. First, I tried to fit a poisson glm: > summary(model.p<-glm(abund~mgmtcat,poisson)) Call: glm(formula = abund ~ mgmtcat, family = poisson) . . . (Dispersion parameter