search for: indiciations

Hello All, Searched around, haven't found a decent solution. I'd like to translate a vector of numbers to a matrix (or to a list of vectors) such that the vector values would serve as indicies of the 1's in an otherwise-zero-filled matrix (or list of vectors). For example: > a = c(1,3,3,4) # perform operation [,1] [,2] [,3] [,4] [1,] 1 0 0 0 [2,] 0 0 1

Hi -- I have found a bug in an HPC code where llvm is vectorizing a loop on Skylake that has an obvious recurrence. I derived a small test case based on the original benchmark below: /*****************************************************************/ static void __attribute__ ((always_inline)) one( const int *restrict in, const int *const end, const unsigned shift, int *const restrict index,

Hello, I would like to plot a series of math symbols on a graph with ordered indicies : expression(delta[1]) expression(delta[2]) etc. Is there a mean to write this in a for loop ? for(i in 1:6) text(i,i,expression(delta[i])) does not work ('i' is not evaluated). Thanks in advance, Y. Noel

The real weird thing is I wrote "weird" and it comes out "wired"... The background is that I did a lot of sparse matrices calculations. Therefore need some indices tweaks, e.g subtract a block from a big matrix. Then I just create an indicies matrix and use it directly as a vector since a matrix is just a vector with dimensions attributed in R. I know I can convert them to a

-----BEGIN PGP SIGNED MESSAGE----- Hash: SHA1 I have two matrices: > m1 <- matrix(1,4,4) > m1 [,1] [,2] [,3] [,4] [1,] 1 1 1 1 [2,] 1 1 1 1 [3,] 1 1 1 1 [4,] 1 1 1 1 > m2 <- matrix(0,3,3) > diag(m2) <- 1 > m2 [,1] [,2] [,3] [1,] 1 0 0 [2,] 0 1 0 [3,] 0 0 1 I want to get indicies from m2

Dear R, I found a very wired behavior for a 2-by-2 matrix, see this example > A <- matrix(1:4, 2) > idx4A <- matrix(1:4, 2) > A[idx4A] Error in A[idx4A] : subscript out of bounds But other matrices are fine, > B <- matrix(1:9, 3) > idx4B <- matrix(1:9, 3) > B[idx4B] [1] 1 2 3 4 5 6 7 8 9 I can reproduce this for both 32bit windows and 64bit linux with R 2.12.1

Hello, I''m using R2HTML library to make a HTML page output (both data frames and graphics): HTMLStart(outdir=paste(getwd(),"/prove html", sep = ""),filename="index", echo = F, HTMLframe = T, withprompt = "HTML> ", CSSFile = "R2HTML.CSS", Title = "Indici di attivit?") ...

I sent this to the list last week, and haven't seen it pop up. Either I deleted it when it did appear, or possibly it was destroyed as spam...? If it did appear and I somehow missed it, appologies. In a nutshell, can the function FUN supplied to by() deduce what level of factor by() was on when FUN was called? I've been digging through the functions, and can't see where the

I would like to run a for loop with an index going from 0 to 499 but the following seems to miss out the first value: C <- 499 for (i in 0:C) The alternative is: C <- 500 for (i in 1:C) { #Then every time I use i, I replace it with i-1 } Is this a good way to do it or is tere a better way? Thank you, ThomasThis message and any attachment are intended solely for the addressee and may

A table without names displays like a vector: > unname(table(2:3)) [1] 1 1 1 and preserves the table class (as with unname in general): > dput(unname(table(2:3))) structure(c(1L, 1L), .Dim = 2L, class = "table") Does that make sense? R is not consistent in its treatment of such unname'd tables: In plot, they are considered erroneous input: >

I am new to R and Im some trouble with the following question... Generate 100 standard normal N(0,1) samples of size 100, X1(k),...,X100(k) where k=1,...,100 (The k is and indicie in brackets) Calculate the sample mean for each sample. For each sample mean Xbark the 0.95-confidence interval for the mean mew=0 is given by... Ik= ( Xbark plus or minus 1.96/10) Find the number of intervals such

Hello: I have a vector with 120,000 reals between 0.00000 and 0.9999 They are not sorted but the vector index is the time-order of my measurements, and therefore cannot be lost. How do I use R to find the starting and ending index of ANY and ALL the "series" or "sequences" in that vector where ever there are 5 or more members in a row between 0.021 and 0.029 ? For

Hi all, Using the FLAC++ level 2 metadata interface the following function doesn't do what you would expect: cuesheet->get_track(i).get_num_indices(); (where cuesheet is a FLAC::Metadata::CueSheet *) It returns a bool instead of an int, why is this? To get the mumber of indicies in a track I need to use: cuesheet->get_track(i).get_track()->num_indices; ie. retrieve the C struct

I'm a newbie, so please forgive me: I have a data frame that includes two long columns of x and y coordinate values thusly: ... ... ... x y ... ... ... ... ... 115 127 ... ... ... ... ... 120 129 ... ... ... ... ... 0 0 ... ... ... ... ... 135 162 .... ... ect... i need to ask R tell me how many x-y pairs contain coords that are

Hello, I am trying to plot SADIE red-blue plots of cluster indicies using filled.contour. I want a plot which only has three bins for the data: <-1.5, -1.5 - 1.5, >1.5, but I am having trouble getting there. example X1 X2 X3 X4 X5 1 -5 -4.5 1.0 4.5 6 2 -3 -2.0 1.2 -1.0 3 3 0 0.0 0.0 -0.5 -1 4 -2 -3.0 1.0 1.5 3 5 -6 -2.0 0.5 3.0 2 example<-as.matrix(example)

I am using the PAM algorithm in the CLUSTER library. When I allow PAM to seed the medoids using the default __build__ algorithm things work well: > pam(stats.table, metric="euclidean", stand=TRUE, k=5) But I have some clusters from a Hierarchical analysis that I would like to use as seeds for the PAM algorithm. I can't figure what the mediod argument wants. When I put in the

Jean-Marc, What exactly does the speex code do to calculate the granulepos? With vorbis it's the 'count of decodable samples including this packet'. So if you had a packet that allowed you to decode 1024 samples since the beginning of playback, the granulepos is 1024, not for example 1023 which would be the *index* of the last sample assuming C-style array indicies starting at

[If I knew who to report this to privately, I would. Sorry to embarrass anyone who's just trying to contribute to the R-project.] There seems to be some oddities with the RSiteSearch web page. When I enter 'RSiteSearch("console")' I'm taken to http://search.r-project.org/cgi-bin/namazu.cgi?query=console&max=20&resu

Hi, I posted around a week ago and I will try this time with a more detailed description. I am converting R code to C code to boost speed as currently it takes a few days to run. Now I will try describe this with as detail as possible. I have a list of size 47. Each element of the list is itself a list of size 7, and each of these lists are arbitrarily long lists of integers. I.e. If my main

Hello, I'm most grateful for your time to read this. I have a uber size 30GB file of 6 million records and 3000 (mostly categorical data) columns in csv format. I want to bootstrap subsamples for multinomial regression, but it's proving difficult even with my 64GB RAM in my machine and twice that swap file , the process becomes super slow and halts. I'm thinking about generating