search for: indv

Hi R sages, Here is my latest problem. Consider the following toy example: x <- read.table(textConnection("y1 y2 y3 x1 x2 indv.1 bagels donuts bagels 4 6 indv.2 donuts donuts donuts 5 1 indv.3 donuts donuts donuts 1 10 indv.4 donuts donuts donuts 10 9 indv.5 bagels donuts bagels 0 2 indv.6 bagels donuts bagels 2 9 indv.7 bagels donuts bagels 8 5 indv.8 bagels donuts bagels 4 1 indv.9 donuts donuts donuts 3 3 indv.10 bagels...

Hello again R users, I have a devilishly hard problem, which should be very simple. I hope someone out there will have the answer to this on the tip of their tongue. Please consider the following toy example: x <- read.table(textConnection("y x1 x2 indv.1 bagels 4 6 indv.2 donuts 5 1 indv.3 donuts 1 10 indv.4 donuts 10 9 indv.5 bagels 0 2 indv.6 bagels 2 9 indv.7 bagels 8 5 indv.8 bagels 4 1 indv.9 donuts 3 3 indv.10 bagels 5 9 indv.11 bagels 9 10 indv.12 bagels 3 1 indv.13 donuts 7 10 indv.14 bagels 2 10 indv.15 bagels 9 6"), header = TRUE)...

Hello all, I am trying to write a program in R in which I call a function multiple times within a loop. The problem is that sometimes the function breaks down while calling another function, and produces an error message that breaks my loop and the program stops. I would like to keep the loop running when this function breaks down, and just move on to the next iteration in the loop. Is there

Hi, I was hoping someone could help me with the following problem. Consider this toy example. For the input dataset there are four individuals (rows "indv.1" through "indv.4"), measured for two different variables (columns "var.1" and "var.2") at two different levels of a factor (column "factor.level"). I want to calculate a matrix that has the average values for each trait (rows "var1.avg" and &...

Hello again, How do I split a data frame into smaller, completely separate data frames (rather than separate data frames comprising a single "list")? Consider the following data, and my coding attempt: x <- read.table(textConnection("id type number indv.1 bagel 6 indv.2 bagel 1 indv.3 donuts 10 indv.4 donuts 9"), header = TRUE) closeAllConnections() x.split <- split(x, x$type) This is where I'm stuck. Now I have one "list" comprised of different data frames, but what I want is separate data frames. Ideally, I'd like t...

Hi I have the followed code for which I am writing a loop pass. int main() { int i = 0; for (i=0; i<20; i++) { printf ("hello world %d\n", i); } return 0; } In the function runOnLoop, I have the following instruction PHINode *indv = NULL; indv = L->getCanonicalInductionVariable(); However, when I check indv is always set to NULL. Since the code has a canonical induction variable, I was expecting the endv to start pointing to the induction var phi node. Please let me what I'm missing here Thanks Malveeka -----------...

...ave the followed code for which I am writing a loop pass. > > int main() { > int i = 0; > for (i=0; i<20; i++) { > printf ("hello world %d\n", i); > } > return 0; > } > > In the function runOnLoop, I have the following instruction > PHINode *indv = NULL; > indv = L->getCanonicalInductionVariable(); > > However, when I check indv is always set to NULL. > Since the code has a canonical induction variable, I was expecting the endv > to start pointing to the induction var phi node. > > Please let me what I'm missing...

Hi all, I have a very simple problem that I cannot seem to find the answer to. Consider the following toy dataset: x <- read.table(textConnection("V1 apples bananas cherries indv.1 7 8 4 3 indv.2 7 7 4 9"), header = TRUE) How would I change the column name of ONLY the first column, not the others? Surely I should not have to re-specify the names of ALL the columns -- e.g., colnames(x) <- c("pears", "apples", "bananas", "cherries&q...

Hi, I turn to you, the R Sages, once again for help. You've never let me down! (1) Please make the following toy files: x <- read.table(textConnection("var.1 var.2 var.3 var.1000 indv.1 1 5 9 7 indv.210000 2 9 3 8"), header = TRUE) y <- read.table(textConnection("var.3 var.1000"), header = TRUE) write.csv(x, file = "x.csv") write.csv(y, file = "y.csv") (2) Pretend you are starting with the files "x.csv" and "y.csv." Th...

Hello R listserve, I would appreciate someone's help with this problem. Consider the following toy dataset: x <- read.table(textConnection("worldclim.1 worldclim.2 cru.1 cru.2 indv.1 7 8 32 658 indv.2 7 7 39 422"), header = TRUE) How could I create a subset of the data based on the column prefix? For instance, let's say I wanted to subset only the columns with the "cru" prefix. Keep in mind that my real data set is very large and likely to change through t...

Is there a way to combine two columns within a data frame? Example data: id snp AL1 AL2 1500 30 A B 1510 30 A A 1520 30 A B This is what I would like: indv snp AL1AL2 1500 30 AB 1510 30 AA 1520 30 AB Any help is greatly appreciated. Alysta

I just had a quick play with gzip compression on a filesystem and the result was the machine grinding to a halt while copying some large (.wav) files to it from another filesystem in the same pool. The system became very unresponsive, taking several seconds to echo keystrokes. The box is a maxed out AMD QuadFX, so it should have plenty of grunt for this. Comments? Ian

...int main() { > >>   int i = 0; > >>   for (i=0; i<20; i++) { > >>     printf ("hello world %d\n", i); > >>   } > >>   return 0; > >> } > >> In the function runOnLoop, I have the following instruction > >> PHINode *indv = NULL; > >> indv = L->getCanonicalInductionVariable(); > >> However, when I check indv is always set to NULL. > >> Since the code has a canonical induction variable, I was expecting the > >> endv to start pointing to the induction var phi node. > >>...

...wish to allow for different repeated measures covariance parameter estimates for different groups (men and women), each covariance matrix having the same structure. In proc mixed this would be done by specifying group= in the REPEATED statement. Is this simple to do in R? (I've tried form=~time|indv/sex for example but this doesn't seem to do the job). 2) I've read that other correlation structures can be specified. Does anyone have any examples of how toeplitz or (first-order) ante-dependence structures can be specified? Many thanks, Gareth

...me privately] WSt> I do not want to call this a bug, but it bugged me. Please try the WSt> following: WSt> tform <- sqrt(RADAI) ~ sqrt(RADAI.e) + TAGE.ej + SPITAL + ARZT + DAS28 + WSt> SJC + TJC + DGA + HAQ + PGA + PAIN + YEAR.SYM + YEAR.DIA + WSt> ALTER + IndV + m.NSAIDs + m.Stereoids + DM.Cyclosp + DM.Penic WSt> WSt> as.character(tform) --> length-3 vector "~" "<LHS>" "<RHS>". However, the end of <RHS> is TRUNCATED! and that truncation of the third part is a bug in my opinion. The fact that...

...0 62 0 38 7 0 0 19 124 30 3889 15 962 51948 230 44 0 62 0 38 We ran lockstat to see what was causing the high SMTX and found this: root at dcivolume01:~# lockstat -D 10 sleep 2 Adaptive mutex spin: 706601 events in 2.060 seconds (342956 events/sec) Count indv cuml rcnt nsec Lock Caller ------------------------------------------------------------------------ ------- 402581 57% 57% 0.00 13105 0xffffff0d492ad678 kidmap_cache_lookup_uidbys id+0x51 265969 38% 95% 0.00 10882 0xffffff0d492...