similar to: acutally a statistics question

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

Dear R Users! I've want to save a lattice "parallel" plot with legend as png: Plotting with legend works well, but when trying to save the plot there occurs a problem with the legend. ################################################### require(lattice) #Some data: data<-matrix(rnorm(100),ncol=5) #Plot works parallel(data)

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

Hello dudes. I'm developing VAR analysis based on suggestions made by Horváth in its paper Canonical Correlation Analysis and Wiener-Granger Causality Tests. That's the reason I'm looking for if there's any R package to develop Wiener - Granger Causality Test. Thanks a lot for your unvaluable help. Regards from Mexico [[alternative HTML version deleted]]

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

I have categories and sub-categories in a table called "categories". Every row has a "parent" field to note if the entry is the sub-category of another (never goes deeper than 1 level). So there might be something like this: id | name | parent 1 Dog 0 2 Pug 1 3 Siamese 5 4 Shih Tzu 1 5 Cat 0 6 Wiener Dog 1 7 Persian

I want to run deconvolution of a time series by an impulse or point-spread function through Wiener filter, regularized filter, Lucy-Richardson method, or any other approaches. I searched the CRAN website and the mailing list archive, but could not find any package for such a deconvolution analysis. Does anybody know an existing R function for deconvolution? TIA, Gang

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

From: tarkhil@over.ru [mailto:tarkhil@over.ru] > On Wed, Aug 20, 2003 at 01:12:51PM +0400, tarkhil@over.ru wrote: > > 10-20 minutes of work. New processes doesn't create anymore, process trying > > to fork looks in top having -20 PRI and "temp" STATE. I was unable to find > > I've got infected computers in my network. When ipnat mapping > table grew to

Hi, all. I'm running R on a beowulf system, with software that can take a list of single-processor jobs and parcel them out over the system. (In case it's of interest, they're dual pentium 866 machines connected by fast ethernet.) The administrators of the system found that I was causing a huge load on the file server, and when they investigated it turned out that each process was

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:

Dear Prof. Lang - I am trying to install the XML library on a 64-bit SUSE linux system (version info below) running 2.4.1. I have gcc version 3.3.3, and libxml2 version 2.6.7. I know this is not current, but I'm on a machine used and administered by others, and updating libxml2 would require updating libc, and things get pretty complicated from there. Trying to install through R (1.9-0),

Hi. I'm still working on figuring out why some C code that interfaced nicely with R-1.1.1 is causing segmentation faults in R-1.2.0. Brian Ripley advised me that R-1.2.0 can be more sensitive to array overruns than was R-1.1.1, so I'm looking for overruns. But in the course of running R with a debugger I've come across something else. Without thinking about it much, I put the

Hi, all. I've got a question that arose when I was switching some row and columns in a matrix. I was expecting the row and column names to follow their rows and columns, and they didn't. I think I see some logic to the behavior, but I'm not sure whether it's intentional or not, so I thought I'd ask the list whether this is a feature or a bug. Here's an example (R-1.3.1,

The following commands illustrate a problem with graphing relatively large data sets using the "line" option to plot: temp <- runif(200000) # or whatever other numbers you like; # same thing happens with a sine wave plot(1:200000, temp) # everything fine plot(1:200000, temp, type="l") # data gets cut off The number of points isn't crucial. The cutoff point

Hi, all. I'm not sure whether this is a problem with ESS or R (or even a feature I don't understand!). Say I've got a function named f.count.these.patterns (as I do -- I'm using actual names in case they're somehow important). I previously had a version called f.count.these.patterns.3, but discarded it because it was worse. So the function is no longer there. If I create a

Christian -- You don't provide enough information (like a call) to answer this. I suspect, though, that you may be subsetting in a way that passes randomForest no data. I'm not aware offhand of an easy way to get this error from randomForest. I tried creating some data superficially similar to yours to see whether something would break if there were only a single value in the variable

It's the "|source" in your formula that tells lattice to separate them. If you drop that, you'll get all points without S and P distinguished at all. If you add a groups argument, you should get them presented with different colors/symbols/etc. depending on your trellis settings (warning: untested code): par.plot(lnvol~lnden, groups =

Hi, all. I'm trying to use some linear regression models in which both the dependent and independent variables are measured with some error. To make things worse, while the errors in the dependent variable are uniform, the errors in the independent (or explanatory, or "x") variables can be heteroskedastic. I've been looking at the book _Measurement Error Models_ by Fuller

Hi, all. I'm trying to cluster 12,500 objects using hclust from package mva. The distance matrix takes up nearly 600 MB. The distance matrix also needs to be copied when being passed to the fortran routine that actually does the clustering (it's modified during the clustering), so that's 1200 MB. I'm actually on a machine with 2.5 GB of memory (and nothing else running), so I