similar to: plot mean + S.E. over time

Dear all! After hours of trying around, I gave up: I have a 2-dimensional array, and I know how to split it into its rows and how to get the mean for every row using 'sapply'. But what I want is to calculate the mean over the first n rows, and then the second n rows, etc., so that I get a vector like: v == mean1(row 1:5), mean2(row6:10),... (trivial, you might say. I find it rather

Dear all, this is a rather basic question; i am not sure how to structure my data well: I want to extraxt various measures from my raw-data. These measures are of different sizes, so I decided to store them in a list, like: run1 <- list(Dom = (my_vector), mean = (my_single_number)) I can do that in a for loop for 40 runs, ending up with 40 lists: run1, run2, ..., run40. To have all the

Dear all, I would like to use the biological symbols for male and female as plotting symbols in a scatterplot (ideally filled and non-filled). R does not seem to have these symbols using pch= in plot() nor are they implemented via expression() or at least I did not find them. I found that the symbols are e.g. available in the wasysym and the marvosym package in LaTeX. I have coded two very rough

Dear all, I would like to use the biological symbols for male and female as plotting symbols in a scatterplot (ideally filled and non-filled). R does not seem to have these symbols using pch= in plot() nor are they implemented via expression() or at least I did not find them. I found that the symbols are e.g. available in the wasysym and the marvosym package in LaTeX. I have coded two very rough

Hi I get a matrix with the columns named like "X13 X22 X1 X14 ...", i.e. not successively. That has it's good reasons, but now how can I rename the columns to "X1 X2 X3 ..."? Thanks a lot, David

I'm using an SVM as I've seen a paper that reported extremely good results. I'm not having such luck. I'm also interested in ideas for other approaches to the problem that can also be applied to general problems (no assuming that we're looking for spirals). Here is my code: library(mlbench) library(e1071) raw <- mlbench.spirals(194, 2) spiral <-

Hi all I am very astonished that R generates a "syntax error" when I want to split up a line with a backslash, which usually works in any shell script. R itself generates the "+" symbols at the beginning of following lines in a splitted line so I've tried with them as well, but also without success. Searching in the h-help archive did not reveal any answer either.

Hi I am trying to replicate a vector in n rows for comparison purposes with another matrix. foo <- c(1,2,3) bar <- rbind(foo,foo) # does the trick for 2 rows bar <- rbind(rep(foo,2)) # does something else How do I generate a matrix with all rows=foo without writing 'foo' n times as arg? Thanks, David

Hi cor(x,apply(x,1,sum)) gives me the correlations of each column with the sums of each row (correct me if I'm wrong, please). What I need are the correlations of each column with the sums of each row except the entry in the given column. It seems that for any one column i I get it by doing: cor(x[,i],apply(x[,-i],1,sum)) But I struggle to get it for all the columns. I was trying things

I thought the "apply" functions are faster than for loops, but my most recent test shows that apply actually takes a significantly longer than a for loop. Am I missing something? It doesn't matter much if I do column wise calculations rather than row wise ## Example of how apply is SLOWER than for loop: #rm(list=ls()) ## DEFINE VARIABLES mu=0.05 ; sigma=0.20 ; dt=.25 ; T=50 ;

Hi Irucka, I tried it and was able to plot it without any errors.? Here, your code indicates you need two lines. temper[[i]][1] ?temper[[1]][1] # which is the column 1. ? Month 1???? 1 2???? 2 3???? 3 ?temper[[1]][2] #? Data1 #1?? 1.5 #2? 12.3 #3? 11.4 Suppose I use names(temper) instead of seq_along(temper) pdf("irucka.pdf") ?lapply(names(temper),function(i)

Hi all, I recently (today) updated the Matrix package and installed the latticeExtra package, but then when I tried to load the lme4 package, I got the following error message:- > library(lme4) Loading required package: Matrix Loading required package: latticeExtra Error in importIntoEnv(impenv, impnames, ns, impvars) : object(s) '.__C__lmeRep' are not exported by

If there is a helpfile open (f.e ?glm) and it is the top window, then an exception error occurs (closing RGUI) when I hit the open file button. If the helpfile is not the top window (of the RGUI) I am able to open a new script without any error. The RGUI is not closing complete there is a blank screen left which I have to close with the X Button or Taskmanager Windows XP Home - German Version

Dear all, I would like to compare two distance matrices to find the correlation between those and create a XY plot.could you please kindly help me and let me know how I can do that using R.(Mantel test?). Your help is very much appreciated!! Many thanks! Regards, Sima Fakheran ............................................................... Sima Fakheran Esfahani PhD Student Institute of

Hi, I work on DDC, the compiler of a research Haskell dialect, Disciple (disciple.ouroborus.net (http://disciple.ouroborus.net)). We are looking to make use of LLVM's type-based alias analysis metadata to encode non-aliasing information between variables. We have found that the tbaa structure is somewhat limited in its expressivity. In particular we couldn't encode intransitive

Hi all, my uni has recently switched to using a proxy (with password) for accessing Web pages outside the uni. I have since found that "update packages" etc no longer seem to work from R. I have tried to fix this problem in two ways:- 1) renaming internet2.dll to internet.dll in the "Modeules" folder and proceeding as usual. When I tried "update packages" from the

Dear list I am comparing the results of several different experimental setups. With kruskal.test() I can test if there is any difference at all in any of them, if I understand it correctly. But now, when there is a difference, how do I generate a (half-) table of pairwise comparisons, using e.g. wilcox.test(), to find the ones where the difference actually occurs. I guess I don't have to

Full_Name: Roberto Ugoccioni Version: 2.2.0 OS: Windows 2000 Submission from: (NULL) (193.203.232.5) Running Windows 2000 Professional, all patches up to nov 8, 2005. > version _ platform i386-pc-mingw32 arch i386 os mingw32 system i386, mingw32 status major 2 minor 2.0 year 2005

Full_Name: Roberto Ugoccioni Version: 2.2.0 OS: Windows 2000 Submission from: (NULL) (193.203.232.5) Running Windows 2000 Professional, all patches up to nov 8, 2005. > version _ platform i386-pc-mingw32 arch i386 os mingw32 system i386, mingw32 status major 2 minor 2.0 year 2005

Hi there, I am running MCMCpack (MCMCirt1d model) on some files (26 items x about 800 – 1200 individuals). I have a problem. When I am working on “big” files, the R program crashes. More precisely I got the following Microsoft Warning: ******************************************* Microsoft Visual C++ Runtime Library Runtime Error! Program: C:\Program Files\R\R-2.4.1\bin\Rgui.exe