similar to: Interaction versus combinations

plot.ts has an argument yax.flip, plot.zoo does not. Is there a way to make the yaxis flip in plot.zoo? I tried using a custom panel function: panel.yaxis<-function(...) { npnl<-parent.frame$panel.number if (npnl %% 2 == 0) { axis(side=3) } else { axis(side=2) } } This leads to a blank window. I am stuck with the intricacies of the plotting and axis

Reading a flat text file 138 Mbyte large into R with a combination of scan (to get the header) and read.table. After conversion of text time stamps to POSIXct and conversion of integer codes to factors I convert everything into one data frame and release the old structures containing the data by using rm(). Strangely, the rm() does not appear to reduce the used memory. I checked using

I have to decode raw vectors read from a connection like this: +++++++++++++++++++ fl<-file.choose() fb<-file.info(fl)$size fc<-file(fl) open(fc,open='rb') seek(fc,where=offset) dat<-readBin(fc,'raw',n=(fb-offset),size=1,signed=FALSE) close(fc) ++++++++++++++++++ The decoding performs a combination of a simple bitwise operator and a simple arithmatic operator on

Is there a way to make any of the write functions (.table or .csv) append to the same file? I get this warning message and do not know how to enable the appending. Warning messages: 1: In write.csv(names(Lbys)[c], fo, append = TRUE) : attempt to set 'append' ignored 2: In write.csv(Lbys[[c]], fo, append = TRUE) : attempt to set 'append' ignored Regards, Alex van der Spek

panel.lmline returns intercept and slope of y ~ x subsetted to the combination of conditioning factors given to xyplot in lattice. for instance: xyplot(Xvalues ~ log(Qvalues)|Tfac, data = df7, panel = panel.lmline) I am looking to find the equivalent formulation for lm() proper. If I do this: lmcal <- lm(Xvalues ~ log(Qvalues):Tfac, data = df7) Only one value of the intercept is returned.

I cannot figure out how to tune the minor tic marks on the date axis of a zoo plot. I read hundreds of CSV files from a zip archive transparently. The time/date strings I convert to POSIXct format, order them and then make a zoo object as there may be cases which have unequal time stamping. As follows: ###################### #Transform timestamps and reorder them dat <- transform(dat,

I use by() to generate a summary statistics like so: Lbys <- by(dat[Nidx], dat$LipTest, summary) where Nidx is an index vector with names picking out the columns in the data frame dat. This returns a list of character arrays (see below for str() output) where the columns are named correctly but the rownames are empty strings and the values are strings prepended with the summary

Hello all, I am trying to create a matrix of 1s and -1s without any repetitions for a specified number of columns. e.g. 1s and -1s for 3 columns can be done uniquely in 2^3 ways. -1 -1 -1 -1 -1 1 -1 1 -1 -1 1 1 1 -1 -1 1 -1 1 1 1 -1 1 1 1 and for 4 columns in 2^4 ways and so on. I finally used the function combn([0 1],3) that I found at the following link

Let say, I have a matrix with 8 rows and 6 columns:  >  df1  <- matrix(NA, 8, 4)  > df1       [,1] [,2] [,3] [,4]  [1,]   NA   NA   NA   NA  [2,]   NA   NA   NA   NA  [3,]   NA   NA   NA   NA  [4,]   NA   NA   NA   NA  [5,]   NA   NA   NA   NA  [6,]   NA   NA   NA   NA  [7,]   NA   NA   NA   NA  [8,]   NA   NA   NA   NA  Now I want to get **all possible** ways to fetch 6 cells at a

I have a list (length 750), each element containing a vector of unique strings (unique gene ids), with length up to ~40 (median 15). I want to compile a matrix of all possible triplets and their frequency within gene elements. Using combn and a lot of looping, I am accomplishing this but it is VERY slow. I've tried to figure out a way to vectorize this, using "match" and

Hello I apologise for the length of this entry but please bear with me. In short: I need a way of subsampling communities from all possible communities of n taxa taken 1:n at a time without having to calculate all possible combinations (because this gives me a memory error - using combn() or expand.grid() at least). Does anyone know of a function? Or can you help me edit the combn or

Hello! I have a regular data frame (DATA) with 10 people and 1 column ('variable'). Its cases are people with names ('a', 'b', 'c', 'd', 'e', 'f', etc.). I would like to write a function that would sum up the values on 'variable' of all possible combinations of people, i.e. 1. I would like to write a loop - in such a way that it

Hi, I am trying to use a vector of indices to select some rows from a matrix. But before I can do that I somehow need to convert 'combinations' into a list, since 'mode(combinations)' says it's 'numerical'. Any idea how I can do that? library("combinat") combinations <- t(combn(8,2)) indices <- c(sample(1:length(combinations),10)) # convert ???

Hi, Suppose I have: > summary(manova(plank.man)) Df Pillai approx F num Df den Df Pr(>F) plankton.new[, 1] 1 0.5267 9.8316 6 53 2.849e-07 *** Residuals 58 --- Signif. codes: 0 `***' 0.001 `**' 0.01 `*' 0.05 `.' 0.1 ` ' 1 My understanding is the MANOVA summary returns a list.

Dear all, I need to do a Manova but I have an unbalanced design. I have morphological measurements similar to the iris dataset, but I don't have the same number of measurements for all species. Does anyone know a procedure to do Manova with this kind of input in R? Thank you very much, Naiara. -------------------------------------------- Naiara S. Pinto Ecology, Evolution and Behavior 1

Hello. Me again. Have you tried to encode piano solo? Noticed high bitrate Opus gave? And there's also artefact at 15kHz which wasn't in the original audio. Visible with Spek program. Download FLAC and Opus both files, new link: http://www.filedropper.com/example_3 FLAC full: 1084 kbps; FLAC solo: 465 kbps. with --bitrate 160: Opus full: 158 kbps; Opus solo: 190 kbps. Included also Spek

I cannot understand why a function appears to be available to R (is.loaded('planckwR') returns TRUE) but the call with .C results in an error message: C symbol name 'planckwR' not in DLL for package <path to DLL file here> This is what I do: Loading a homebrewed DLL, compiled with MS VS97. The functions are all void, the C calling sequence is used and all arguments

Dear R users, Can anyone tell me how to get the p value out of the output from summary.manova? I tried all the methods I can think of, but failed. Many thanks Yu-Kang _________________________________________________________________ ¥ß§Y¥Ó½Ð MSN Mobile ªA°È¡G¦b±zªº¤â¾÷¤W¦¬µo MSN Hotmail http://msn.com.tw/msnmobile

Hi, I am trying to generate all unordered combinations of a set of numbers / characters, and I can only find a (very) clumsy way of doing this using expand.grid. For example, all unordered combinations of the numbers 0, 1, 2 are: 0, 0, 0 0, 0, 1 0, 0, 2 0, 1, 1 0, 1, 2 0, 2, 2 1, 1, 1 1, 1, 2 1, 2, 2 2, 2, 2 (I have not included, for example, 1, 0, 0, since it is equivalent to 0, 0, 1). I have

R-2.9.1, Windows7 Dear list, I have a question to you that seems very simple to me, but I just can't figure it out. I have a dataframe called "ratings" which contains the following variables: evalR1, evalR2, evalR3, evalR4, scoreR1, scoreR2, scoreR3, scoreR4, opinionR1, opinionR2, opinionR3, opinionR4. (there are more variables, but this gives an idea of the data structure). What