similar to: plot data with a colour scale - more details!

i'd like to use the C interface to R in a program i'm writing. as a starting point, i'm trying to create a very simple C program that uses R. i've read the R documentation on this, but i'm having trouble figuring out where SEXP is defined and how to use it. i noticed someone else on this list also tried to use the C interface, but they ran into similar problems:

Good Morning. Below I like statement like j<-grep(".r\\b",colnames(mydata),value=TRUE); j with the \\b option which I read long time ago which Ive found useful. Are there more or these options, other than ? grep? Thanks. dstat is just my own descriptive routine. > x ?[1] "age"????????? "sleep"??????? "primary"????? "middle" ?[5]

does anyone know how to use the apply function on a two-variable t-test? i've tried everything. apply(data[1:4],1,t.test) # works, but only a one-variable test apply((data[1:4],data[5:9]),1,t.test) # returns a syntax error i've also tried using the for loop: for (i in 1:1000) t.test(data[i,1:4],data[i,5:9]) doesn't print anything! it processes for awhile, so i know it's doing

Hello Masters, wishing you all a great 2011 I was also going to ask if anyone knows a quick and efficient way to plot a regression plane (z~x*y). I have tried the regr2.plot{HH} function but it is only an educational tool and has poor graphical properties. I also tried to run the following script on a fictitious longitudinal problem, with poor results set.seed(1234)

Hi, Please brace yourselves for some possibly stupid questions. Why are there only 2 codebooks for floors and many more for residue values? What makes the encoding of floor values fundamentally different from encoding residue values? I have a basic understanding of VQ, and looking at the floor decode function, it all makes sense 1) Read the codebook index from the stream 2) Do a

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

?s 02:10 de 02/08/2024, Steven Yen escreveu: > Good Morning. Below I like statement like > > j<-grep(".r\\b",colnames(mydata),value=TRUE); j > > with the \\b option which I read long time ago which Ive found useful. > > Are there more or these options, other than ? grep? Thanks. > > dstat is just my own descriptive routine. > > > x > ?[1]

Javier, Creo que aquí aplica la ley de Linus que dice: "Dado un número suficientemente elevado de ojos, todos los errores se convierten en obvios". La persona que revisa y encuentra un error no necesariamente tiene que ser la misma que la que lo escribe. Una motivación muy importante al compartir un código es la de recibir los beneficios del control de calidad por parte de tus pares.

I found this sentence in the Ogg format specs: http://www.xiph.org/ogg/vorbis/doc/vorbis-spec-res.html "A residue vector may represent spectral lines, spectral magnitude, spectral phase or hybrids as mixed by channel coupling." But where does the spectral phase information come from ? AFAIK MDCT doesn't provide any phase information. And in OGG-encoding, MDCT is taking place a few

Hi I am trying to build a new decoder in Java, I know there are a couple of java decoders already but they don't implement the bit of the spec I am having questions on Question 1) In section 8.6.2 it mentions "However, it is also possible to set [residue_begin] and [residue_end] to specify a range partially or wholly beyond the maximum vector size. Before beginning residue decode,

Dear R People: I have the output from an arima model fit in an object xxx. I want to verify that the ma1 coefficient is there, so I did the following: > xxx$coef ar1 ar2 ma1 intercept 1.3841297 -0.4985667 -0.9999996 -0.1091657 > str(xxx$coef) Named num [1:4] 1.384 -0.499 -1 -0.109 - attr(*, "names")= chr [1:4] "ar1" "ar2"

Dear list, I tried to use the nparcomp to run some post hoc non-parametric comparison and got and error. Error in uniroot(pfct, interval = interval) : f() values at end points not of opposite sign Appreciate any comments. the command line: >nparcomp(Ulceration~Group,data=test,type='Dunnett',control='Non-treated') Jun

Hello, I'm trying to modify the residue vectors in Vorbis audio files at encoding time. For testing, I'm just adding +1 to each residue vector's element, so e.g., if I have an original residue vector of [0,1,2,3] I modify it during the encoding to be [1,2,3,4], and at the decoder's side I should receive the [1,2,3,4] vector. The problem is that, after adding the +1, I still

Hi, I have the following file, and I need to work out the linear regression for each sample. I tried the model(*) and receive the error message (**): > data=split(mydata,rep(1:(nrow(mydata)/6),each=6)) > arrang.linear=lapply(data,lm,formula=KA~PA) Erro em storage.mode(y) <- "double" : invalid to change the storage mode of a factor Al?m disso: Warning message: In

Hi, Does anybody know how codebooks are generated in OggVorbis encoder? We are porting oggorbis encoder on embedded platform for which VQ codebook memory is hugeeee to imagine. How can we reduce that? Can we do VQ with less codebooks and if yes how? If any help available? Parul Embedded Engineer Einfochips Ltd

hello every one, How to function more than one loop in R? I have the following problem to be solved with the a method of three loops, can you help me please? The data is attached with this message. The data is composed of two parts, cleaved (denoted by ?cleaved?) and non cleaved (denoted by ?noncleaved?). ? to access to the ith peptide, you can use X$Peptide[i] ? to access to the ith label,

Dear All, I need to separate one column which is in a format like, "can_region", into two columns like, "can", and "region". Do we have a function in R to do this? Thanks. Regards, Jin [[alternative HTML version deleted]]

Dear All, As i have understood it , the Huffman codebook uses a huffman code present in the bitstream ( lets us say it is decoding a floor 1 Y value using the specifed codebook number ) to walk through the huffman code book tree till it hits a leaf node which is nothing but the entry number. If we consider our floor 1 Y value example , this entry number is directly being used as the Y value.

Dear all, I am a R beginner, and I am looking for a way to do the same thing for all levels of a column in a table. Basically, I have a bunch of protein sequences composed of different amino acid residues, and each residue is represented by an uppercase letter. I want to calculate the ratio of different amino acid residues at each position of the proteins. Here is an example table: Proteins

Hi folks, I'm doing a merge of my current branch onto the mainline (for testing) today. I believe it to be stable. Just a little more vorbisfile testing. After the merge, I have a few more patches to apply, then onto cascading/coupling. New stuff: Floor backend 1 and residue backend 1; both are present, but the mainline modes won't use either yet. Naturally, both are enabled for