similar to: xyplot {lattice} are different types possible for each panel?

Hi R Folks... Are the results of a fit carried out by panel.lmline readily available for use in a lattice plot? I¹d like to put r^2, m, and b on each panel. I can certainly write something that does this manually and then use it with panel.text, but if it¹s already available, that would be preferable, especially as lattice permits condition and subsetting so readily. Looking at panel.abline, I

Gurus: I keep seeing other people?s code that contain ideas like If (x == 2L) X[-1L] X - 1L I have some idea of what?s going on, but where is the use of concepts like ?2L? documented? Thanks, Bryan ************* Bryan Hanson Acting Chair Professor of Chemistry & Biochemistry DePauw University, Greencastle IN USA

I think I more or less understand what a ?wrapper? is, but I?d like to hear how more experienced R users define it, and especially I'd like to know if there is a formal definition. In my reading, it seems like there are a fairly wide range of meanings, but they are all conceptually similar. I've looked in a couple of the classic R texts, the extensions and developers' manuals, and R

Hello Folks. This must be a silly question with a (not) obvious (to me) answer. Consider this: tmp <- matrix(1:200, nrow = 20) vec <- 300:309 tmp[9,] <- vec # replacing one row works fine p <- c(3, 11, 17) tmp[p,] <- vec # replacing multple rows pastes the values down a column and recycles vec. What I want to do is replace multiple rows simultaneously at once. I suppose I can

Hi all... My students were conflating grepping for a value in a vector to get the index, and then removing it with [-index], for instance like this: set.seed(17) v <- rnorm(20) s <- v[-1.18] They were trying to remove the 12th value in v, which is -1.18 or so. But the result is, as documented in ?Extract, to coerce 1.18 to the next lowest integer, and remove the 1st value of v, not the

I?m trying to build a simple formula interface to work with a function using ggplot2. The following scheme ?works? up until the plot(p) request, at which point there are complaints about xlim?s and a blank graphics window. Looking at str(p) I do see the limits are NULL, plus layer 1 claims to have an empty data frame (but df is reproduced correctly). I'm sure I'm missing something really

Hello Again R Folks: I?m trying to clean up some code. Suppose I have an object like this: > str(test) List of 2 $ G:List of 2 ..$ cls:'data.frame': 101 obs. of 2 variables: .. ..$ V1: num [1:101] -0.0019 -0.0019 -0.00189 -0.00188 -0.00186 ... .. ..$ V2: num [1:101] 0.000206 0.000247 0.000288 0.000329 0.000371 ... ..$ rob:'data.frame': 101 obs. of 2

One for the ggplot2 gurus... I have a function which makes a plot just fine if the response vector (res in the example; fac1 is a factor) has no NA in it. It plots the data, then makes a little annotation at the bottom with the data counts using: p <- p + geom_text(aes(x = fac1, y = min(res) - 0.1 * diff(range(res)), label = paste("n = ", ..count.. , sep =

Hi Folks... No code to troubleshoot here. I need some suggestions about the right terminology to use in further searching, and any suggestions about R pkgs that might be appropriate. I am in the planning stages of a project in which IR, NMR and other spectra (I'm a chemist) would be collected on various samples, and individual samples would be followed over time. The spectra will be feature

Hello Wise Ones... I need a clever way around a problem with findInterval. Consider: vec1 <- 1:10 vec2 <- seq(1, 10, by = 0.1) x1 <- c(2:3) a1 <- findInterval(x1, vec1); a1 # example 1 a2 <- findInterval(x1, vec2); a2 # example 2 In the problem I'm working on, vec* may be either integer or numeric, like vec1 and vec2. I need to remove one or more sections of this vector;

Hello i need some precision about nipals in the chemometrics package in R . When i use nipals in chemometrics i obtain T and P matrix. I really don't understand what to do with these two matrix to obtain the scores for every the component (like in spss fo example) Comp1 Comp2 Comp3 quest1 0,8434 0,54333 0,3466 quest2 0,665 0,7655 0,433 Thank you very

My question is this: is there a way I can make one of the arguments of the function be a string the user can enter, and then have that be the excel filename? ie, foo <- function(x,y,NAME){ #make a matrix with x rows and y cols M <- matrix(nrow=x,ncol=y) #write the matrix write.table(M, file = "result.csv",append=TRUE, sep = ",") } I've had a look but I

Hello Kindred R Spirits... I'm trying to get a file (csv) from Dropbox using their shareable link concept. They issue a short URL that goes to a web page where you see a button that says "Download File". They don't really give you the URL of the file itself, just this page. Is there a way to coax R into getting such a file? I don't even really want the file per

Hello Everyone. I?m on a Mac, using rgl. Thanks to all the developers and maintainers on it! If I drag the mouse, I can rotate the view, if I hold down the option key or the ctrl key and drag, I get scaling (though at times I seem to alter the perspective?). Adding the shift key doesn't seem to modify anything. So, what I can't find is a key combination that translates the object.

Hello Again... I?m making a faceted plot of a response on two categorical variables using ggplot2 and having troubles with the coloring. Here is a sample that produces the desired plot: compareCats <- function(data, res, fac1, fac2, colors) { require(ggplot2) p <- ggplot(data, aes(fac1, res)) + facet_grid(. ~ fac2) jit <- position_jitter(width = 0.1) p <- p +

Hello R Folks... Not a technical question, but I need some advice and perspective. I?ve got a set of functions I?m planning to put together into a package. The main hunk of data that gets used by different functions is currently an S3 list. I?ve been reading about S4 objects, and I see the (numerous) advantages of them. I have seen the recommendation that all new packages be done with S4.

Hi All... I?m trying to build a small demo using gWidgets which permits interactive scaling and selection among different things to plot. I can get the widgets for scaling to work just fine. I am using gcheckboxgroup to make the (possibly multiple) selections. However, I can?t seem to figure out how to properly define the gcheckboxgroup; I can draw the widget properly, I think my handler would

Hello R Folks... I have a list with the following structure: > str(df) List of 3 $ y : num [1:4, 1:1242] -0.005379 0.029874 -0.023274 0.000655 -0.004537 .. $ x : num [1:1242] 501 503 505 507 509 ... $ names: Factor w/ 4 levels "PC Loading 1",..: 1 2 3 4 I want to plot each row of df$y against df$x, and have each plot in it¹s own panel according to the levels of df$names.

Hello Again Gurus and Lurkers: I?m trying to build a very user-friendly function which does aov without having the user type in a formula (which would be tedious in this case). The idea is to take the response from a PCA score matrix, and the factors from a list. A simple example is the function given below, along with test data and a sample call to the function. I'm certainly having

Hello All. I must be missing the really obvious here: mm <- function(d, i) median(d[i]) b1 <- boot(gravity$g, mm, R = 1000) b1 b2 <- boot(gravity$g, mm, R = 1000, strata = gravity$series) b2 Both b1 and b2 seem to have done (almost) the same thing, but it looks like the strata argument in b2 has been ignored. However, str(b1) vs str(b2) does show that the strata have been noted