similar to: beanplot-Error: sample is too sparse to find TD

Is it possible to force beanplot not to use a log scale? I want to be able to create multiple beanplots of different data on the same specified y-axis. When I try to force a ylim, I get the following... > beanplot(reg5vel,ylim=(c(0,12000))) log="y" selected Warning message: In plot.window(xlim = c(0.5, 7.5), ylim = c(0, 12000), log = "y") : nonfinite axis limits

Hi, I would like to use the beanplot() function from the beanplot package. Unfortunately, I can't find out how to suppress the dashed horizontal line, that shows the overall mean. In the help I've found the argument "overallline", but it only allows for "mean" or "median" . I have tried overallline = F, overallline="n", and

Dear Helpful R Users, I am graphing some data using the beanplot, but I am having trouble getting the output I desire. I have five tanks (A-E) and 2 groups for each tank grp1 or grp2, except tank C where there is only grp1. (I only changed the grouprep to "C grp1" for the example) When I plot them, I would like A B C(only grp1 - half of the bean plot) then D and E (as full beans).

Hi there, I have a set of data that looks like this: As1988<-c(1254.0, 22.0, 4.2, 1081.0, 35.0, 6.0, 1772.0, 192.0, 7.6) The mean of this (as calculated by R) is: 485.9778 The median of this (as calculated by R) is: 35 If I then make a beanplot(As1988), I find that the beanline (average) is now 77.68561 while the beanline (median) is 35.39739 (using the locator function to check the

Dear all, I am trying to create beanplots from a dataset for which boxplot works fine. (MACOS, R 2.8.1 GUI 1.27 Tiger build 32-bit (5301)) I am getting the following error message: Error in shapiro.test(x) : sample size must be between 3 and 5000 I am not even sure why the shapiro.test is being used, but is there any workaround ? Thanks ! Markus [[alternative HTML version deleted]]

Hi, I am writing a function that requires a specific package to be installed. Is there a way of checking if the package is installed and returning a TRUE / FALSE result so my function can return an appropriate error message and exit the function gracefully rather than just bombing out? I'm thinking along the following lines (but want code that works), f_checkpackage <- function() { if

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

Dear R-Users, I have a newbie-question about producing several plots with five variable pairs within in one code (loop). Problem: I define two objects, erveryone has five vectors (columns /variables): dimensionen <- data.frame(meinspss$attr_diff_gesamt, meinspss$finanz_diff_gesamt, meinspss$leist_diff_gesamt, meinspss$soz_diff_gesamt, meinspss$wert_diff_gesamt) gruppe <-

Estimados amigos, mi pregunta es bien directa, a ver si alguien me puede ayudar: Como marcar la media en un grafico boxplot ?? Por lo que veo, los boxplot en R marcan por defecto la mediana dentro de la caja, pero necesito marcar la media y ojala sacar la mediana para que no aparezcan los dos parametros. Por lo que he leido, imagino que se puede hacer con alguna funcion de bajo nivel, pero soy

Dear List, I have a set of distributions recorded at an equal interval of time and I would like to plot them as series of horizontal histograms (with the x-axis representing time, and y-axis representing the bins) since the distribution shifts from unimodal to multimodal in several occasions. What I would like to see is something close to a violinplot, but I do not want a kernel density

Dear list I wish to extract from a population genotypized for 10 SNP a subsample of the same population of size n with similar allele frequencies. Essentially i have a matrix of 200 rows (df) like this Name,Condition,rs1385699_X,rs6625163_X,rs962458_X,Rs4658627_1, sample01,Case,1,1,1,-1 sample02,Control,1,1,1,1 sample06,Control,1,-1,1,0 sample10,Case,1,1,1,0 sample11,Control,1,1,1,1

Good morning ! Today I found a strange, for my poor knowledge of R, behaviour of 'which' on a matrix: HAL9000> str(cluster.matrix) num [1:227, 1:6300] 2 2 2 2 2 2 2 2 2 2 ... HAL9000> class(cluster.matrix) [1] "matrix" HAL9000> ase <- cluster.matrix[1:5,1:5] HAL9000> ase [,1] [,2] [,3] [,4] [,5] [1,] 2 2 2 0 -2 [2,] 2 2 2 0 -2 [3,]

Hi, I'm trying to generate a distance matrix between sample pairs (example below). I'm not very familiar with the loop command which I expect I will need for this. The example below demosntrates what I'd like to get out of the data - essentially, to calculate the proportion of positions where two samples differ. Any help much appreciated! Also, any notes on how the functions work

I have multiple conditions that I would like to plot in a grid. To save space, I don't need the ticks, labels etc on the plots which are not on the outside. I tried split.screen, but if it is clearer to use, it seems pretty rigid since it decides of each screen size apriori while I need more space for the left column and top or bottom row screens than for the other plots. axis(1,tick=FALSE)

This question pertains to setting up a model in the package "dlm" (dynamic linear models, http://cran.r-project.org/web/packages/dlm/index.html I have read both the vignette and?"An R Package for Dynamic Linear Models" (http://www.jstatsoft.org/v36/i12/paper), both of which are very helpful. There is also some discussion at

Given a vector : pvec=(p1,p2,.... p J) with sum(pvec)=1, all the elements are non-negative, that is, they are probabilities a matrix A ( N* J ), with the elements alpha(ij) are 0 or 1 I want to MAXIMIZE THE RESULT RESULT= product( i=1, to N [ sum ( alpha(ij)* pj , j =1,to J ) ] ) thus, I need to get pvec. how should I do ? for example

Hi, Is there a function for raising a matrix to a power? For example if you like to compute A%*%A%*%A, is there an abbreviation similar to A^3? Atte Tenkanen > A=rbind(c(1,1),c(-1,-2)) > A [,1] [,2] [1,] 1 1 [2,] -1 -2 > A^3 [,1] [,2] [1,] 1 1 [2,] -1 -8 But: > A%*%A%*%A [,1] [,2] [1,] 1 2 [2,] -2 -5

Hello, I have a multi-year dataset (see below) with date, a data value and a flag for the data value. I want to find the monthly median for each month in this dataset and then plot it. If anyone has suggestions they would be greatly apperciated. It should be noted that there are some dates with no values and they should be removed. Thanks Emily > print ( str(data$flow$daily) )

--/9DWx/yDrRhgMJTb Content-Type: text/plain; charset=us-ascii Hello everyone, verison: 1.0.1 platform: Gnu/linux mandrake 7.1beta I'm running a user function to step through a frame and plot one column against the next. This crashes the program fairly reliably although not on the first try. Here is a log of the session run through gdb. First I launch R, print out the top of the frame and