similar to: logarithmically scaled y-axis in vioplot

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

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

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 some questions for the wireframe function of the lattice package. My dataset's "x-data" are sampled logarithmically and as such I would like to have a semilogarithmic 3D plot when plotting a time series. Does anyone know how to change the example in http://addictedtor.free.fr/graphiques/sources/source_40.R such that x or y would be logarithmically scaled? Also, I would

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

dbox format needs at least two tools: 1. Expunging only marks the message having refcount=0. The actual disk space is freed only after calling a tool pointing to the user's dbox dir. 2. Ability to force index resync/rebuild in case it doesn't figure out itself that there's something wrong. I'm starting to think about building a more generic mailutil tool, maybe something similar

Dear all, I have a problem with the command 'getURL' from the RCurl package, which I have been using to obtain a ftp directory list from the MOD16 (ET, DSI) products, and then to download them. (part of the script by Tomislav Hengl, spatial-analyst). Instead of the list of files (from ftp), I am getting the complete html code. Anyone knows why this might happen? This are the steps i

Dear R-Help, I've never had any trouble importing data into R until I had to import an .sdd file for a class. The file can be found here: http://www.math.umt.edu/steele/Math%20549/Farms.sdd. It begins with the line "## Dump S Version 4 Dump ##". I first attempted read.S which issued the message "not an S object". I then checked the Import/Export manual which seemed to

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

Hi, does anyone know how to make pairs of boxplots? E.g. 5 pairs of boxes. I figured a not very elegant way, by inserting a dummy NA field to the list where the data is given. So, where I have the dummy I got more distance between boxes, as a result it looks like groups. But I am having problems with putting the names under the groups, and the thick marks. Any better idea? Thanks, Kati ___

Dear All, I need to calculate the variance components in a mixed effect model (one fixed and one random effect) with REML (maximizing the proportion of the likelihood that does not depend on the fixed effects). In S+ there is the varcomp function, but I would like to do it in R. Is there a way to do that? Thanks! Katalin ___ Katalin Csillery Division of Biological Sciences University of

Hi I'm currently prowling for some IMAP tools to help me maintain our mailboxes. I'd prefer "pure" IMAP, i.e. no messing with file system directly, but would settle for more crude ways if need be. Specifically I am looking for: - Expunge tool: Outlook does not seem to auto-expunge messages in any way shape or form. So I'd like to run a daily cron to expunge all

Hi, I have a coding problem. I have data pairs, and I want to calculate a parameter for the data pairs as a whole data set (so not using the information that thely are pairs) I am looking for bias in that certain parameter so I want to bootstrap these data pairs(!), and than calculate the parameter. I run into a coding problem. How should I define theta in order to resample the pairs, but than

Hello, I am interested in fitting a generalized nonlinear regression (gnlr) model with negative binomial errors. I have found Jim Lindsay's package that will do gnlr, but I have having trouble with the particular model I am interested in fitting. It is a threshhold model, where below a certain value of one of the parameters being fitted, the model changes. Here is a sample: Cones:

Dear All, I am trying to use the bca.ci function on a boot.out object which consists a few NaN's and I want to ignore those NaN's, and get a ci only for the "normal" values. boot.out$t has R number of values for 3000 different statistics, so when I use boot.ci(boot.out, index=i) and i happens to be a column in boot.out$t with some NaN's in there I get an error message. I