similar to: workspace question

Hello all, I have matrices V.i of dimension n.i x n.i, where i = 1, ..., J, and the sum of n.i equals N. (and n.i ! = n.j) goal: create one large matrix V, where V has matrices V.i on diagonal. I create each matrix V.i in a for loop (1 to J), so each time I'd like to augment V with the most recently calculated V.i, such that I'll have V after the final iteration of the for loop.

Hello all, X is a list of 20 lists, and each individual list has 65 elements. Y is a list of 65 elements. WANT: subtract Y from each of the 20 lists in X. Here's what I tried and the error messages: > X - rep(Y, 20) Error in X - rep(Y 20) : non-numeric argument to binary operator I tried several methods w/o success. Any suggestions kindly appreciated. Thanks, Dave ps - please copy

All, In the function below I have 24 individuals and 6 calculations per individual. The 6 calculations are collected each time in a 1:24 loop when calculating "delta". I'd like to collect all 144 = 24*6 calculations in one vector ("delta.patient.comb"). The function works as is via indexing, but is there an easier way to collect the measurements via appendinng the 6

All, For data consisting of serial measurements on subjects, one may use the aggregate function to say compute the peak response for each subject for each design condition. Is there a way to alter this or another one-liner to also retain the time at which the peak occurred and thus avoid writing a doing this via a loop? I suppose one could attempt to employ the split function but that's

All, Is there are a wildcard in R for varible names as in unix? For example, rm(results*) to remove all variable or function names that begin w/ "results"? cheers, Dave ps - please respond directly to afshar@miami.edu [[alternative HTML version deleted]]

All, Usually when I save a variable I have ascii = FALSE since I usually load the variable via the load(variable) command and do not need to view the variable outside of R. When I set ascii = TRUE and view the variable outside of R in a text editor, I notice additional characters (starting w/ RDA2 ...) before the first actual value in the vector. Is there are way to save to ascii such that this

All, When saving plots to a pdf file via the pdf function, I would like to be able to automatically expand the graphics device to achieve the same result as when one does this manually (e.g., clicking the green expand button on the upper left of the graph on Mac OS). Consider simple example below: library("lattice") foo.frm = data.frame(Subject = rep(c(1:4), each = 9), Y =

All, When defining the arguments of a function, is it possible to supply a function as an argument? If so, how is this introduced into the function code as well? For example, in the body of the function I have: result = function(x) and I'd like to supply either function.1 or function.2. Please reply directly to afshar@miami.edu Thanks! Dave [[alternative HTML version deleted]]

All, I'm having a problem selecting directly from a vector. I've written ways to do this indirectly, but I'd rather do it directly and didn't see this in the manual. Essentially, I have: > group.label.new [1] 7 9 6 1 10 4 8 3 2 5 > junk [1] 1 2 > group.label.new[junk && 8:10] [1] 7 9 6 1 10 4 8 3 2 5 I'd like to select the elements

All, I'm trying to create an augPred plot in the nlme library, similar to the plot on p.43 of Pinheiro & Bates (Mixed Effects Models in S and S-Plus) for their Pixel data. My data structure is the same as the example but I still get the error msg below. > comp.adj.UKV <- groupedData(adj.UKV ~ Time | Patient_no/Lisinopril, data = comp.adj.UKV.frm, order.groups = F) >

All, I'd like to plot the main relationship of a grouped data object for all levels of a factor in a single panel. The sample code below creates a separate panel for each level of the factor. I realize that this could be done in other ways, but I'd like to do it via plotting the grouped data object. thanks! dave z = rnorm(18, mean=0, sd=1) x = rep(1:6, 3) y =

All, When I take a subset of a factor the reduced factor still maintains all the original levels of the factor when say forming the key in a plot. The data is correct, but the variable still "remembers" the original levels. See below for reproducible code. Does anyone know how to fix this? cheers, dave fact = as.factor(c(rep("A", 3),rep("B", 3), rep("C",

All, I have a repeated measures data frame and was wondering if the covariance matrix can be calculated via some created indexing or built-in R function. Specifically, say there are 3 variables, where potassium concentration is measured 6 times on each patient. Patient number (discrete) Time (1 to 6, discrete) Potassium (continuous variable) I want the covariance/correlation matrix for the

All, Simple question but I don't seem to be able to find the answer in the documentation: When using "plot" within a loop, is there any way to supply the argument to "main" dynamically, i.e., so that the title is Patient k below as the loop cycles through each value of k? plot(x,y, xlim=c(0,250), ylim=c(0,1000), xlab="gamma", ylab="r1",

All, Is there an efficient way to apply say "mean" or "median" to a dataframe according to say all combinations of two variables in the dataframe? Below is a simple example and the outline of a "manual" solution that will work but is not very efficient (could also generalize this to a function). Searched the archives and docs but didn't see anything close to

sent this to the list yesterday but didn't see it listed in the daily summary ... apologies if you receive it twice ... ________________________________ From: Afshartous, David Sent: Tuesday, June 27, 2006 10:02 AM To: 'r-help@stat.math.ethz.ch' Subject: sapply question All: I'm trying to use sapply to break up data within another function. (tapply doens't seem to work

All, Is there one liner way to obtain the max per observation for two vectors? I looked at apply and lapply but it seems that groundwork would have to be done before applying either of those. The code below does it but seems like overkill. Thanks! Dave x = rnorm(10) y = rnorm(10) ind = which(x < y) z = x z[ind] <- y[ind] ## z now contains the max's

All, I've coded a function and it works manually if I copy it line by line into R. However, when I try to "load" (copy and paste) the entire function into R, I get the following error after the listed line of code: + N.j.list = lapply(rej.hyp, length) Error: symbol print-name too long Does anyone you know what this error means? Strangely, when I copy the same line verbatim

All, When specifying colors to xyplot w/ a groups argument, using auto.key no longer maintains the colors properly. I've searched the docs and help but haven't found exactly what I need ... I saw a few examples in the archives involving par.settings but that doesn't seem to do it. I also saw some people using key instead of auto.key, but that didn't seem consistent. Is there a

All, I'm trying to make a basic plot: data points superimposed upon the a line connecting the points w/ a different color. Example below doesn't work as the first xyplot call doesn't remain. Suggestions? David Hour = c(NA,1,2,3,4) y = c(2,2,3,2,1.5) xyplot(y ~ Hour, xlab = list("Hour", font=2, cex=2), ylab= list("U_x * V", font=2, cex=2), type = 'l',