similar to: Writing a function, want a string argument to define the name of the excel sheet to be called

Hi all, My problem is as follows: I want to run a loop which calculates two values and stores them in vectors r and rv, respectively. They're calculated from some vector x with length a multiple of 7. x <- c(1:2058) I need to difference the values but it would be incorrect to difference it all in x, it has to be broken up first. I've tried the following: r <- c(1:294)*0 rv

I'm using the TSA package (along with all prerequisites) to do some GARCH work and for some reason, something which used to work for me has decided to up and stop. The code is as follows, after loading the package: " gs <- garch.sim(alpha=c(1.9,0.1), beta=c(0.700001, -0.0800003, -0.016),rnd = rnorm, n = 400, ntrans=500) gs1 <- garch.sim(alpha=c(1.9,0.1), beta=c(0.7, -0.08,

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

Has anyone plotted or is it possible to plot CIE *xy* chromaticity diagram http://en.wikipedia.org/wiki/File:CIE1931xy_blank.svg I need this plot in color. ishida [[alternative HTML version deleted]]

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

grDevices::convertColor has arguments 'from' and 'to' which can take on value 'XYZ'. Can someone confirm that 'XYZ' is the same as the CIE chromaticity coordinates that are also sometimes refered to as 'xyY' in the literature? Or are these the CIE tristimulus values? It looks to me like the first case is true, but I would appreciate hearing from one of

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

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

Still playing with Lattice... I want to use panel.text(x, y etc) but with x and y in plot coordinates (0,1), not user coordinates. I think if I had this problem with traditional graphics, I could use grconvertX to make the change. I did come across convertX {grid} but this doesn't seem to be what I need. Is there a function like grconvertX in Lattice, or is there a flag or some other

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

Hello R Folks... Using the example below, I¹d like two of the panels to be plotted with type = ³p² but the third to be done with type = ³h². I can¹t use type = c(³p², ³p², ³h²) because this syntax applies all given types to every panel. I don¹t think I can use groups and distribute.type because these are intended for different styles of plotting within a single panel. As you can see, I tried

Hello R Masters and the Rest of Us: The first of these works fine, the 2nd is accepted but too literal (the "%->%" is shown in the plot label and in the wrong position). The 3rd throws and error due to "unexpected SPECIAL". Would someone recommend a way to format this? I want the two phrases connected by a right arrow. TIA, these things always elude me. Bryan

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 =

On this feed, which I think is the place we should monitor upcoming changes: http://developer.r-project.org/blosxom.cgi/R-devel/NEWS What is the significance of the green and pink highlighting? Thanks, Bryan

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;

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