similar to: bagplot()

Has anyone ported the Bagplot function by Rousseeuw, Ruts, and Tukey from S to R? The S function comprises a script and a FORTRAN function. I assume porting is relatively uncomplicated, but since I have not done any porting before I would not want to invest the effort if a port is readily available.

Hi, I am using the bagplot function within the aplpack package and often get this error: Error in chull(pg[, 1], pg[, 2]) : NA/NaN/Inf in foreign function call (arg 2) Here's an example of data for which this happens: bagplot(c(-4.66, -2.62, -3.65 ,-3.07 ,-4.91, -4.56 ,-3.79, -3.10, -3.01), c(2.14, -7.18, -0.15 , 2.67, -2.49, -0.89 , 0.96 , 0.59, -4.59)) Can anyone help? --

Hello, I would like to simulate randomly from a multivariate normal distribution using a correlation matrix, rho. I do not have sigma. I have searched the help archive and the R documentation as well as doing a standard google search. What I have seen is that one can either use rmvnorm in the package: mvtnorm or mvrnorm in the package: MASS. I believe I read somewhere that the latter was

When I type library() to see what is installed the following list in RED comes up. Packages in library '/home/jason/R/i686-pc-linux-gnu-library/2.13': abind Combine multi-dimensional arrays aplpack Another Plot PACKage: stem.leaf, bagplot, faces, spin3R, and some slider functions biglm bounded memory linear and

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

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

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

I have two vectors x and y, which I would like to plot against each other. I am also displaying other data in this plot. However, I have about 1 million points to plot, and just plotting them x againt y is not very informative. What I'd like to do is to do sort of a continuous box plot. My x values goes from -1 to 1 and my y values from 0 to 1, so I?d like to plot the median and quantiles,

Dear R People: Does anyone have any code for Fractals, chaos, or anything like that, please? This is strictly for demo purposes...decorative only. This is R version 1.5.1 for Windows. Thank you in advance! Have a great weekend! Sincerely, Erin Hodgess Associate Professor Department of Computer and Mathematical Sciences University of Houston - Downtown 1 Main Street Houston, TX 77002

I am recieving a connection unexpectedly closed error while trying to sync a directory on two machines. Rsync is in my path on both machines. I am using ssh and I can connect to the server fine through ssh. I am using a vpn but have no problems connecting to any of the servers with it. This is the full error message I get: Remote Machine: Connection refused rsync: connection unexpectedly

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

I can take the results of a simulation with one random variable and generate an empirical interval that contains 95% of the observations, e.g., x <- rnorm(10000) quantile(x,probs=c(0.025,0.975)) Is there an R function that can take the results from two random variables and generate an empirical ellipse that contains 95% of the observations, e.g., x <- rnorm(10000) y <- rnorm(10000) ?

I have a simple problem: I need to load a ERSI shapefile of US states and check whether or not a set of points are within the boundary of these states. I have the shapefile, I have the coordinates but I'm having a great deal of difficulty bringing the two together. The problem is the various GIS packages for R do not play well with each other. sp, shapefiles, maptools, etc all use different

Hello, I am new to network routing and I need help configuring a linux box with two ethernet cards. In this case it''s a Linux RH 7.3 box, in a cabinet that already has a couple of Windows servers. The Windows server routing is below as an example. The Linux box has an out-of-band interface at 10.130.36.38 and a public eth at 62.50.8.84. I had to add a route for the private interface so

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 have a quick question: It appears that in rgdal v0.7-12 (R version 2.15.1, OSX 10.6.8) writeOGR will not write a shapefile the the current directory. Is this correct? An earlier version of rgdal must have allowed this because I have a older script that used to work, but doesn't now. So, as an example, here is what I get today: > shape = readOGR('.',