similar to: t() converts data.frames ind matrix

hi, I have a matrix (representing original data) that looks e.g. like this (consider it beeing x,y,z coords): 441 447 0 265 407 0 374 223 0 288 574 0 669 309 0 591 195 0 595 475 0 424 351 0 I get a second matrix (subject data) that is similiar to the above matrix but it is scaled, translated and rotated (and of course a little inprecise). Also I have an

hi, i am working on a little R-project with a couple od other guys.we use CVS, but everyone keeps the R-source files in different locations in his home-directory. of course this causes trouble when sourcing R-files. i thought a UNIX environment variable could be the solution, but R doesn't seem to know about the environment variables. e.g. >> source("$PROJECT/xxx.R")

hi, is there a way to source all R-files that reside in a given directory with a single R-call ? greetinx jan -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info", "help", or "[un]subscribe" (in the "body", not the subject !) To:

hi, I am looking for a possibility to open several graphics in independent windows (is that the graphics device ?). up to now, i was only able to generate a graphic, then save it (e.g. as eps) and then generate a new plot that always overwrites the old graphic. I am aware of the fact that I can generate multiple plots in a single graphics window, but this is not what i am looking for ..

hi, another anova question !! how does aov() treat missing values ? since it might depend on my experimental design: i have a within subject design with two factors ... i use aov() like this: aov(performance~type*block+Error(subj/(type+block)),data=anovaData) now it might happen that i have a couple of NaN data points in my result dataFrame. how can i know/control how aov() is treating those

Hi, simple question I guess: the following line works well: aveBehav=c(apply(sdata, 2, mean)) However, I would like to pass an argument to the function mean, namely na.rm=TRUE Does anyone knows how to do this? Thanks in advance, Jan

hi, i am using the barplot function, but i have problems with the plot region when adjusting the bar width manually. if you only want to plot 2 bars the default barwidth is to wide (aesthetically speaking). by adjusting width and xlim one can actually produce narrow bars, but then the plotted graph is not in the center of the plot region, which means that for example the title of the plot

hi, knowing there are a lot of statistics guru's in the R-mailing list, i throw in a statistics question. i have data that i do not know how to statistically test: subjects are repeatedly asked to make a decision (e.g. left-right -> coded as 0 or 1). i have 20 subjects, each subject made 8 decisions. i now want to analyse whether my experimental manipulation induced a systematic bias

hi, during a for(i in 1:xx]) loop I always newly calculate a vector (e.g. tmp<-c(a,b,c,d,e) ) now i need that vector to be attached at the bootom of a matrix (or data.frame). e.g. m<-matrix() for(i in 1:5]){ #some calculations for a,b,c,d,e a<- ... b<- ... c<- .. .... tmp<-c(a,b,c,d,e) ??? now I need to attach this tmp to the matrix m ??? } i couldn't get rbind

hi, when performing a wilcox.test or a wilcox.exact i get results that looks like this: wilcox.exact(x, mu=.5) Exact Wilcoxon signed rank test data: x V = 207, p-value = 0.0006905 alternative hypothesis: true mu is not equal to 0.5 the way i understand the wilcox.test (or wilcox.exact) the V-value represents the summed up ranks of either the positive or negative differences,

hi, i want to plot percentage values using barplot. everything works fine so far. giving a chance level of 50% i want to plot only the area above 50 percent. i can of course cut using ylim=c(50,100) but this does not work properly since i also see some of the area below 50%. thanks for suggestions, jan -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help

hi, i am looking for a way to correct for ties in the wilcoxon signed rank test -> e.g. wilcox.test(x,mu=.5) one way i have heard of is to look up the p value in a table that has been produced by Buck (1975). obviously i need to know the T-value to do that -> how do i get the T-value from the wilcox.test() function. is there any other (already implemented) way to correct for ties in

hi, suppose you have a for-loop like this: for(i in1:x) {} ... now you want to generate a variable in every cycle (since you do not know the size of x in beforehand you have to do this dynamically). the variable should e.g. look like this: variableName1<-c() (if x==1) variableName2<-c() (ifx==2) variableName3<-c() (if x==3) .. .. i tried this(which obviously didn't work):

-- Dr. Jan Malte Wiener Max-Planck-Institute for Biological Cybernetics Spemannstr. 38, 72076 Tuebingen, Germany tel.: +49 7071 601 631 email: jan.wiener at tuebingen.mpg.de url: www.kyb.tuebingen.mpg.de/~malte -------------- next part -------------- An embedded message was scrubbed... From: Jan Wiener <jan.wiener at tuebingen.mpg.de> Subject: Re: [R] reading "special" text files

In Julian Faraway's text on pgs 117-119, he gives a very nice, pretty simple description of how a glm can be thought of as linear model with non constant variance. I just didn't understand one of his statements on the top of 118. To quote : "We can use a similar idea to fit a GLM. Roughly speaking, we want to regress g(y) on X with weights inversely proportional to var(g(y).

hi, i am trying to generate nice barplots with std-errors. do i really have to generate the std-errors myself by the segments() command ? thanks for help, jan -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info", "help", or "[un]subscribe" (in the

Hi, I am looking for ways to donwsample one-dimensional vectors. For example, x=sample(1:5, 115, replace=TRUE) How do I downsample this vector to 100 entries? Are there any R functions or packages that provide such functionality. I did find the zoo package and the aggregate() function, but these appear to be rather specific for time-series. Thanks in advance, Jan

Hi guys I have implemented xapian on a website, and it currently has about 2M items in its index. Its all been working quite nicely so far, until I tried removing some old items from the index (removing items when the index was smaller was no problems at all). When I try to remove them now (using writabledatabase_delete_document() via php), it halfway freezes up the machine, and the apache

Hello frieds, do I need to do 'net rpc join -U Administrator' (to selfjoin the localhost to the domain) if I have samba3.0.4 running in ROLE_DOMAIN_PDC and with ldapsam-Backend and want to use winbindd for ntlm_auth squid authentification? -- Best regards, Malte mailto:malte.woelky@gmx.de _________________ Malte Woelky -=[SkyNet]=-

Dear R-users, I was wondering for the following: Let 'x' be a matrix and 'ind' and indicator matrix, i.e., x <- array(1:20, dim = c(4, 5)) ind <- array(c(1:3, 3:1), dim = c(3, 2)) I'd like to get (as a vector) the elements of 'x' which are not indexed by 'ind'. Since negative indices are not allowed in index matrices I thought of using something like: