similar to: another automation question

I have columns in two tables that need to be compared in order to get a statistic. The function calculating the statistic is somers2. Here are the two columns I must calculate, row by row. name1 name2 1 Pred_pres_a_indpdt[,3,,] population[,23] 2 Pred_pres_b_indpdt[,3,,] population[,24] 3 Pred_pres_c_indpdt[,3,,] population[,25] 4 Pred_pres_d_indpdt[,3,,]

I have a series of columns that need to be evaluated in various tables. I need to apply a function in the following manner somers2(name1,name2). name1 is a vector of x inputs for the function which correspond to a vector of y inputs in name2. y<-rep(c(3,4,5,8),6) z<-rep(c(23,24,25,26,27,28),4) name1<-sprintf("Pred_pres_%s_indpdt[,%s,,]",x,y)

Hi everybody, I'm working on Mac OS X and R-2.1.0 (by the way, sorry for the multiple emails last time, had a technical problem, and thanks for the constructive comments.) I have a .R function which calls a .C function. The R function's name is K_MEANSR.R and the other is K_MEANSC.C. I compile it with 'R CMD SHLIB K_MEANSC.C' in the terminal. I get both .so and .o files.

I am trying to load binary files in the following fashion load("pred/Pred_pres_a_indpdt") load("pred/Pred_pres_b_indpdt") load("pred/Pred_pres_c_indpdt") load("pred/Pred_pres_d_indpdt") load("pred/Pred_pres_e_indpdt") load("pred/Pred_pres_f_indpdt") but I would like to set up a for loop to replace the letters a:f Here is what I have so

Hi, i am working with large data frames with many dependend variables. I want to write some functions that will allow me to quickly select variables from the frame and plot them in various colors depending on factor columns, possibly selecting rows according to factor conditions. In order to do this in a nice function, i need to understand how to work with a column name in the body of a

Dear all, I'm looking for primary help at aggregating table() results and at writing a loop (if useful) My dataset ( http://goo.gl/gEPKW ) is composed of 23k rows, each one representing a point in the space of which we know the land cover over 10 years (column y01 to y10). I need to analyse it with a temporal sliding window of 5 years (y01 to y05, y02 to y06 and so forth) For each period

Dear All, Thank you for the replies to my first thread here: http://r.789695.n4.nabble.com/global-optimisation-with-inequality-constraints-td3799258.html. So far the best result is achieved via a penalised objective function. This was suggested by someone on this list privately. I am still looking into some of the options mentioned in the original thread, but I have been advised that there may

Hi all, I'm trying to use the package read.table within a foreach loop. I'm grabbing 500M rows of data at a time from two different files and then doing an aggregate/tapply like function in read.table after that. I had planned on doing a foreach loop 39 times at once for the 39 files I have, but obviously that won't work until I figure out why the segfault is occurring. The

hi everybody I wonder if there is a simple way, but not simple would be ok too, to customize info/welcome page at session start time? what I'd like to do is to put together simple short howto / dos & don'ts page for users, I'm thinking it would be great if it was possible many thanks [[alternative HTML version deleted]]

I have written a function which returns an "incorrect number of subscriptions error" and am unable to identify why this should happen. Any help would be gratefully received. The problem: I would like to generate combinations/permutations across different selections - for example all the combinations of "3 chosen from 6" combined with "5 chosen from 8" combined with

biplot could benefit from allowing colour palettes, eg. allowing col=list(xcol= xxpal, ycol= yypal) where xxpal and yypal are palettes. For an example, see the robust pricipal components in the principal componetent section of <http://www.statlab.uni-heidelberg.de/data/olive/olive10.pdf> The required change is to replace col[1L] resp. col[2L] by col[1L]] reep. col[[2L]] for example in

I'm finding that readLines() and read.fwf() take nearly two hours to work through a 3.5 GB file, even when reading in large (100 MB) chunks. The unix command wc by contrast processes the same file in three minutes. Is there a faster way to read files in R? Thanks!

Dear all, I'm used to draw bubble plots (gstat) to get a first view of my spatial marked data, but I can't find a way to label legend with marks after having replaced them by a numeric scale. I have a data frame with numeric coords and factors: > names(data) [1] "x" "y" "bloc" "sub" "inoc"

Hi, I do not have much R experience just the basics, so please excuse any obvious questions. I would like to create bubble plot that have Categorical data on the x and y axis and then the diameter if the bubble the value related to x and y. Attached to the email is a pic of what I would like to do. I do hope someone can help me. -- Regards/Groete/Mit freundlichen Gr??en/recuerdos/meilleures

Hi All! I am desperately needing some help figuring out how to calculate LD50 with a GLMM (probit link) or, more importantly, the standard error of the LD50. I conducted a cold temperature experiment and am trying to assess after how long 50% of the insects had died (I had 3 different instars (non significant fixed effect) and several different blocks (I did 4 replicates at a time)=

Dear R users I use Windows XP, R2.5.1 (I have read the posting guide, I have contacted the package maintainer first, it is not homework). In a research project on renal cell carcinoma we want to compute Harrell's c index, with optimism correction, for a multivariate Cox regression and also for some univariate Cox models. For some of these univariate models I have encountered an error

Hello all. I am trying to run PAN, multilevel multiple imputation program, in R to impute missing data in a longitudinal dataset. I could successfully run the multiple imputation when I only imputed one variable. However, when I tried to impute a time-varying covariate as well as a response variable, I received an error message, “Error: subscript out of bounds.” Can anyone tell if my commands

Hi, I am using lmer to analyze habitat selection in wolverines using the following model: (me.fit.of <- lmer(USED~1+STEP+ALT+ALT2+relM+relM:ALT+(1|ID)+(1|ID:TRKPT2),data=vdata, control=list(usePQL=TRUE),family=poisson,method="Laplace")) Here, the habitat selection is calaculated using a so-called discrete choice model where each used location has a certain number of alternatives

Dear all, I am trying to make a series of waffle plot-like figures for my data to visualize the ratios of amino acid residues at each position. For each one of 37 positions, there may be one to four different amino acid residues. So the data consist of the positions, what residues are there, and the ratios of residues. The ratios of residues at a position add up to 100, or close to 100 (more on

Dear alltogether, I tried pan() to impute NAs for longitudinal data. The terminology in the following output follows the pan manpage. No data are attached to this script as this may be too huge. y = 15 responses pred = at first just intercept was tried (later on covariates should follow) subj = 168 different subjects with 4 to 6 observations for each subject at time points t1, t2, ..., t6 #