similar to: Unique :key not maintained after add_indexes?

Hi Terry, On 03/16/2016 08:03 AM, Therneau, Terry M., Ph.D. wrote: > Is the phrase "index <- match(x, sort(unique(x)))" reliable, in the > sense that it will never return NA? This is assuming that match() and unique() will never disagree on equality between 2 floating point values. I believe they share some code internally (same hashing routine?), so maybe it's reliable.

Dear List, Let's say I have a list whose components are 2 matrices (as exemplified in the "mylist" object below). I'd like to create a list with components being 4 matrices based on an logical index vector. is there a way to simplify what I'm doing to obtain the results in "mylist2"? I'd like something that would work on an arbitrary number of elements in

Hi, I am a Ph.D. student from Québec, Canada. I’m a beginner with R and Bioconductor. Until now the only experience I have is in analyzing microarray data using affy and limma packages. Now I am trying to analyze Rat Gene 10 st arrays and I would like to run RMA analysis and Smyth moderated t test on those arrays. Since no cdf official package is available for those arrays, after reading many

Hi everyone, I have 2 environments (2 different R sessions) as described below: Session 1: Name of the environment: "CrlmmInfo" Objects in the environment: index1: logical index - length 238304 index2: logical index - length 238304 priors: list of 4 - (matrix 6x6, 2 vectors of length 6, vector of length 2) - all num params: list of 4: centers [238304 x 3 x

Hi everyone, I have 2 environments (2 different R sessions) as described below: Session 1: Name of the environment: "CrlmmInfo" Objects in the environment: index1: logical index - length 238304 index2: logical index - length 238304 priors: list of 4 - (matrix 6x6, 2 vectors of length 6, vector of length 2) - all num params: list of 4: centers [238304 x 3 x

Hi, I'm trying to plot bar graphs with stacked columns marked with askterisks for certain columns, which follow certain criteria. I've gotten the stacked bar plots with Data set A (please refer to code below). However, I haven't figured how to put asterisks on columns based on some extra information (such as "used/unused") as in Data set B i.e. asterisks at the

Hey, Thanks for your great inputs. While "index = apply(mat == 0, MARGIN = 1, any)" gives you an idea of the rows containing zero(s), "index<-data[!apply(data==0,MARGIN=1,any),]" does the actual job of removing the rows with zeros. Kind regards Ogbos On 9 March 2010 15:12, Paul Hiemstra <p.hiemstra@geo.uu.nl> wrote: > Dimitris Rizopoulos wrote: > >>

Hi All, I have a whole lot of *.raw files in my working folder and I am doing the same analysis on each of those and want to save all the results in a single file. I am making some mistake here and can't figure out how to solve it. Say, the *.raw files are ABCD.raw, EFGH.raw, IJKL.raw ... The files are of this format ID PHI?? aa1? aa2? aa3 .... 1??? 1???? 1.3?? 2.0?? 1.0 2??? 0????

Hi, I'm using "by" to summarize by multiple groups, and want to extract the returned into a pretty dataframe. I'm trying to find a simple way to name the rows of the data frame. I'd like it to be something like index1.val1.index2.val2 where the index1 and index2 are the names of the indices and the val1 & val2 are names of possible values of the index. (the

Buenas, He creado dos andamios para dos tablas a la que accedo por index1 e index2 respectivamente. En cada index puedo hacer las tareas de añadir, borrar, ver lineas... pero me gustaria saber como se hace una busqueda mas compleja...es decir si tengo una tabla Productos y otra tabla descripciones y quiero mostrar en los productos con sus descripciones... he intentado en el index1 acceder a la

I do it because: 1. Different retrieved records use different variations of the same template, and may want more/fewer/different queries and present them differently. Plus templates can display whatever they want and not have to rely on the controller to assign it data. e.g. index1.rhtml can retrieve and display categories, news, calendar; index2.rhtml can retrieve and display categories, news,

Dear listers: As a part of a intermediate process, I need to use and modify a very big matrix (some 30000x30000) inside the body of a function. If the matrix is defined in the function, R shows a error message "Cannot allocate vector of size 6.7 Gb". But if I define the matrix before the function is called (as soon as the dimension can be calculated) R can allocate it. The problem is

Is the phrase "index <- match(x, sort(unique(x)))" reliable, in the sense that it will never return NA? Context: Calculation of survival curves involves the concept of unique death times. I've had reported cases in the past where survfit failed, and it was due to the fact that two "differ by machine precision" values would sometimes match and sometimes not,

This is a proposal to enhance TableGen's ability to analyze and manipulate DAGs. Hopefully this will allows more complex DAGs to be built in TableGen. 1. Add a new value suffix. value(index) The value must be a DAG. The index specifies the operator or an operand, whose value is produced. The index can be 0 produce the operator 1...n produce operand by

rly there would be a heavy penalty to launching a shell so you would want t= o carefully evaluate the frequency this is executed on your system.<br /> <div class=3D"container"> <div class=3D"line number1 index0 alt2"><code class=3D"sql plain"><br /> DELIMITER @@</code></div> <div class=3D"line number2 index1

Hi all, I''ve recently started working with Ferret and I''m getting what seems to be slow searches. I have about 10000 documents in the index, with several fields per document, with some fields having an array of several values that are indexed. I am using a RAMDirectory to store the index for searching. When doing testing, I find that searches are reasonable at around .2 to

Dear Srinidhi, You are trying to fit 1 random intercept and 2 random slopes per individual, while you have at most 3 observations per individual. You simply don't have enough data to fit the random slopes. Reduce the random part to (1|ID). Best regards, Thierry ir. Thierry Onkelinx Statisticus / Statistician Vlaamse Overheid / Government of Flanders INSTITUUT VOOR NATUUR- EN BOSONDERZOEK

Dear Srinidhi, You are trying to fit 1 random intercept and 2 random slopes per individual, while you have at most 3 observations per individual. You simply don't have enough data to fit the random slopes. Reduce the random part to (1|ID). Best regards, Thierry ir. Thierry Onkelinx Statisticus / Statistician Vlaamse Overheid / Government of Flanders INSTITUUT VOOR NATUUR- EN BOSONDERZOEK

I am running a multilevel growth curve model to examine predictors of social anhedonia (SA) trajectory through ages 12, 15 and 18. SA is a continuous numeric variable. The age variable (Index1) has been coded as 0 for age 12, 1 for age 15 and 2 for age 18. I am currently using a time varying predictor, stress (LSI), which was measured at ages 12, 15 and 18, to examine whether trajectory/variation

I am running a multilevel growth curve model to examine predictors of social anhedonia (SA) trajectory through ages 12, 15 and 18. SA is a continuous numeric variable. The age variable (Index1) has been coded as 0 for age 12, 1 for age 15 and 2 for age 18. I am currently using a time varying predictor, stress (LSI), which was measured at ages 12, 15 and 18, to examine whether trajectory/variation