similar to: least summed square distance + fit

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, 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, 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, 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 just recognized that the t() function converts a data frames into a matrix (R 1.5). is that a bug or a feature ? 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

Hi! I'm trying to figure out how to use optim... I get some really strange results, so I guess I got something wrong. I defined the following function which should be minimized: errorFunction <- function(localShifts,globalShift,fileName,experimentalPI,lambda) { lambda <- 1/sqrt(147) # error <- abs(errHuber(localShifts,globalShift, #

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, 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, we are just about to evaluate R as a standard statistics package for our institute. We looked around the help - manual, the FAQ, etc. and did not find anything for the topics of random factors in ANOVA, repeated measures in ANOVA, or within subjects designs in R. Could anyone point us to any information concerning - univariate approach to repeated measures anova - multivariate

-- 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

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, 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, 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

I have a collection fo mp3s that I wish to stream. They came from a windows environment and when ripped they used spaces in the files names. When I prepare a file list with the spaces on the linux box it works fine with mpg123 and muse buts spits out error messages from ices and continues to cycle through the playlist until I break out of ices. <p>Error opening /music/d/m4/ANNOINTED -

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

how can i sort a matrix by sorting the first colloum in ascending order and have the values of the other colloums in the relative position to the sorted values. thanks for help -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info", "help", or

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, I am trying to create 3 mds plots: one with axis 1 vs axis 2, one with axis 2 vs axis 3, and one with axis 1 vs axis 3. When inputting my code, I only end up with one diagram and when inputting mdsg.mds$dims, the program returns 2 for 2 dimensions. How can I create the other two plots? Any help would be greatly appreciated, Calla Carbone The table I use is number of each taxa on each

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, I have been trying to create a function to generate a Procrustes plot, generated from package "vegan" in lattice. standard vegan code as follows library(vegan) pro=protest(P1, P8, permutations=4999,choices=1:4) plot(pro) Now, here is the code for the function that I have failed to get to work properly. panel.procrustes=function(x,y) }Pro=protest(x,y,permutations=4999,choices=1:4)