similar to: lattice question: how to change the dot on boxplot to line

Hello, I am drawing a boxplot using bwplot in lattice. I would like to learn how to draw a thicker median line instead of the default setting. Thank you. Code bwplot(........, pch=rep("|",2)) Elaine [[alternative HTML version deleted]]

Dear friends, I have two questions regarding the use of lattice. First some code: ## begin code z <- cbind(rep(c("BIC", "ICL", "s_v", "Q_v", "sig-q", "s_lsk", "s_lML", "s_mlsk", "s_mlML", "s_la8", "s_haar"), each = 250), rep(c(5, 10, 20, 30, 50), each = 50)) z

Dear all, I would like to have the same color for the all boxplots from the same panel, but my code below shows the two colors alternating. Thanks! set.seed(42) D1 <- rnorm(200) D2 <- factor(sample(letters[1:2],200,TRUE)) D3 <- factor(sample(letters[3:5],200,TRUE)) DF <- data.frame(x=D1,a=D2,b=D3) print(bwplot(b~x|a,data=DF,col=c("black","black"),

Dear list, Does anyone have a suggestion (or better still) code for sampling from the uniform distribution over the convex hull of a set of points? Many thanks and best wishes, Ranjan

On Tue May20'25 05:56:59AM, Rui Barradas wrote: > From: Rui Barradas <ruipbarradas at sapo.pt> > Date: Tue, 20 May 2025 05:56:59 +0100 > To: Ranjan Maitra <mlmaitra at gmx.com>, r-help at r-project.org > Subject: Re: [R] how to read a PSB file in R? > > ?s 03:17 de 20/05/2025, Ranjan Maitra via R-help escreveu: > > I have come across this file (warning,

Dear all, I have a question with regard to the use of hclust. I would like to be able to specify my own distance matrix instead of asking R to compute the distance matrix for me. It is computationally easier for me this way. My question is: How can I get hclust to accept this? Thanks, Ranjan -- *************************************************************************** Ranjan

Thanks! Yes, however, this seems a bit wasteful. Just wondering if there are other, more efficient options possible. Best wishes, Ranjan On Thu, 29 Mar 2018 22:20:19 -0400 Boris Steipe <boris.steipe at utoronto.ca> wrote: > If one is equal to the reverse of another, keep only one of the pair. > > B. > > > > > On Mar 29, 2018, at 9:48 PM, Ranjan Maitra

Dear All, is any way I can change the colours of the box and umbrella of a bwplot without having to go to the length of: box.rectangle<-trellis.par.get("box.rectangle") box.rectangle$col<-"black" trellis.par.set("box.rectangle", box.rectangle) etc... but straight from the call: bwplot(y ~ x | z, mydata) ? Regards, Federico Calboli

I don't know if this is more efficient than enumerating with distinct directions and weeding... it seems kind of heavyweight to me: ####### library(gtools) directionless_circular_permutations <- function( n ) { v <- seq.int( n-1 ) ix <- combinations( n-1, 2 ) jx <- permutations( n-3, n-3 ) x <- lapply( seq.int( nrow( ix ) ) , function( i ) {

Jeff, I wanted to let you know that your function is faster than generating the directional circular permutations and weeding. Here is the time for n = 10. I compared with just doing the permutations, there is no point in proceeding further with the weeding since it is slower at the start itself. system.time(directionless_circular_permutations(10)) user system elapsed 1.576 0.000

Dear friends, I would like to get all possible arrangements of n objects listed 1:n on a circle. Now this is easy to do in R. Keep the last spot fixed at n and fill in the rest using permuations(n-1, n-1) from the gtools package. However, what if clockwise or counterclockwise arrangements are the same? I know that half of the above (n - 1)! arrangements are redundant. Is there an easy way to

Dear friends, Is there a way to read data files written in lisp into R? Here is the file: https://archive.ics.uci.edu/ml/machine-learning-databases/university/university.data I would like to read it into R. Any suggestions? Thanks very much in advance for pointers on this and best wishes, Ranjan -- Important Notice: This mailbox is ignored: e-mails are set to be deleted on receipt. Please

Thanks! I am trying to use it in R. (Actually, I try to give my students experiences with different kinds of files and I was wondering if there were tools available for such kinds of files. I don't know Lisp so I do not actually know what the lines towards the bottom of the file mean.( Many thanks for your response! Best wishes, Ranjan On Wed, 17 Jan 2018 20:59:48 -0800 David Winsemius

Dear friends, Sorry for this somewhat generically titled posting but I had a question with using contrasts in a manova context. So here is my question: Suppose I am interested in doing inference on \beta in the case of the model given by: Y = X %*% \beta + e where Y is a n x p matrix of observations, X is a n x m design matrix, \beta is m x p matrix of parameters, and e is a

New function below is a bit faster due to more efficent memory handling. for-loop FTW! directionless_circular_permutations2 <- function( n ) { n1 <- n - 1L v <- seq.int( n1 ) ix <- combinations( n1, 2L ) jx <- permutations( n-3L, n-3L ) jxrows <- nrow( jx ) jxoffsets <- seq.int( jxrows ) result <- matrix( n, nrow = factorial( n1 )/2L, ncol = n ) k

Dear list, I have been using the Rmath library for quite a while: in the current instance, I am calling dnt (non-central t density function) repeatedly for several million. When the argument is small, I get the warning message: full precision was not achieved in 'pnt' which is nothing unexpected. (The density calls pnt, if you look at the function dnt.) However, to have this happen a

Dear colleagues, I have an image which I can display in the greyscale using image. On this image, for some pixels, which I know, I want to display their activity based on a third measure. One way to do that would be to color these differently, and use an opacity measure to display the third measure. An example of what I am trying to do is at:

?s 03:17 de 20/05/2025, Ranjan Maitra via R-help escreveu: > I have come across this file (warning, massive, 4.3 GB) https://esahubble.org/media/archives/images/original/heic1502a.psb and it appears to be a filetype I was not aware of previously. > > Is it possible to read the file in R using any tool? It is an image and I am looking for the RGB of the file. > > Many thanks and

Dear colleagues, I had a question wrt the car package. How do I evaluate whether a simpler multivariate regression model is adequate? For instance, I do the following: ami <- read.table(file = "http://www.public.iastate.edu/~maitra/stat501/datasets/amitriptyline.dat", col.names=c("TCAD", "drug", "gender", "antidepressant","PR",

Possibly because: panel.hist is not an existing R function -- you have to first create it so pairs() can use it. ?pairs shows you how in the Help examples, i.e. panel.hist <- function(x, ...) { usr <- par("usr") par(usr = c(usr[1:2], 0, 1.5) ) h <- hist(x, plot = FALSE) breaks <- h$breaks; nB <- length(breaks) y <- h$counts; y <- y/max(y)