similar to: dyn.load the same code more than once?

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

Hello, all. I've got what is probably a simple question about comparison of models using anova, specifically about the situations in which it's valid. I understand, I think, what's going on when the models are strictly nested (as most are in the demo(lm) examples). My question involves what happens when the models aren't strictly nested. In my particular case, I'm doing

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, all. Does anyone have an equivalent to Splus "peaks", which finds local maxima (with locality defined by a parameter "span")? I thought I'd check whether anyone has done it already before trying to put something together myself. Thanks, Matt Wiener -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read

Hi, all. The following does not behave as I think it should, and as it seems to me it has in the past (although I can't check this easily). I know it happens in both R-0.64.2 and R-0.65.0 on an old Power Computing running Linux-PPC 1999, and in R-0.64.2 on an SGI running Irix 6.5. Try the following: t1 <- data.frame(matrix(rnorm(16), nc=4)) > t1 X1 X2 X3 X4 1 -0.7206945

Hi, all. The help file for tapply says that if simplify is true, and the result of the calculation is always a scalar, then tapply will return a vector. Nonetheless: > t1 <- tapply(runif(10), rep(1:5, 2), mean) > is.vector(t1) [1] FALSE > is.array(t1) [1] TRUE > I have found this in version 0.65.1 on an SGI running Irix 6.5, and on a Mac running Linux-PPC. I've also

Hi, all. I've got some C code that I've used successfully, and fairly extensively, in R-1.1.1, which is suddenly causing segmentation faults in R-1.2.0. The original code used calloc (just plain calloc) with a free at the end, but I've tried replacing the calloc with R_alloc (based on the Writing R extensions documentation, which I'm not sure I've understood properly), and

Thanks to John Fox, Brian Ripley, and Peter Dalgaard for responding. The short answer (as in Peter Dalgaard's reply, already posted to the list) is that the models I'm concerned with can in fact be compared using ancova. The key fact is that while the parameters may not be nested, the subspaces I'm examining are. An additional note from Prof. Ripley on AIC and BIC (which I quote in

Hi, all. In trying to compile R v.0.90 on an alpha running Red Hat Linux 6.0, I'm getting the following error: make[3]: Entering directory `/home/mcw/alpha-R/R-0.90.0/src/appl' gcc -I. -I../include -I../../src/include -DHAVE_CONFIG_H -mieee -g -O2 -c cpoly.c -o cpoly.o /tmp/ccPjZNwd.s: Assembler messages: /tmp/ccPjZNwd.s:312: Error: unknown opcode `sqrttsu' /tmp/ccPjZNwd.s:2631:

================================================= This is an automated summary of the status of the R-bugs repository. Note that this may be neither complete nor perfectly correct at any given instance: Not all bugs are reported, and some reported bugs may have been fixed, but the repository not yet updated. Some bug fixes are difficult to verify because they pertain to specific hardware or

Hi, all. Using the patched version of R-0.99.0, I cannot reproduce the following example from the strsplit documentation: (Example) unlist(strsplit("a.b.c", ".")) ## [1] "" "" "" "" "" ## Note that `split' is a regexp! ## If you really want to split on `.', use unlist(strsplit("a.b.c", "\."))

Two more details: Someone else who works here, and hasn't used 0.65.0, thinks he's seen the error before. He restarted R, and it went away. Second, I think the relevant code is in SEXP do_function, in eval.c in the main directory. line 695: WrongArgCount("lambda") Beyond that I'm afraid I'm getting in over my head. Matt

Hi. I'm having a problem with disappearing comments when I transfer files (entire .RData files) from one machine to another, and I'm not sure whether it's an R problem or an ESS problem. Since ESS seems more likely, I'm trying this list first. (However, it seems to me this problem only showed up when I upgraded to R-0.65.0 or R-0.65.1.) Functions originally shows up fine -- if

Full_Name: Matthew Wiener Version: 1.1.1 OS: linux Submission from: (NULL) (156.40.248.102) I've come across a glitch in hist. I can reproduce it on linux for Intel, linux for PPC, and Irix 6.5. t1 <- c(41, 42, 42, 43, 43, 43, 44, 44, 45, 46) hist(t1, breaks = 10) (OK) hist(t1/50, breaks = 10) Error in hist.defauilt(t1/50, breaks = 10): some 'x' not counted; maybe

================================================= This is an automated summary of the status of the R-bugs repository. Note that this may be neither complete nor perfectly correct at any given instance: Not all bugs are reported, and some reported bugs may have been fixed, but the repository not yet updated. Some bug fixes are difficult to verify because they pertain to specific hardware or

Hi all, I just noticed that my two nics that are bonded via mode 0 are showing the same mac addy in ifconfig. Is this normal? In syslog I keep getting; kernel: pbond0: received packet with own address as source address Some guidance is appreciated before I dive in and troubleshoot. If I don't list the mac addy in my ifcfg-th# files, only one of the 2 nic comes up during boot.

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

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)

Matthew Wiener <mcw@ln.nimh.nih.gov> writes: > t1 <- data.frame(matrix(rnorm(16), nc=4)) > t1$X1 <- 1 > t1$X2 <- 2 > print(t1) > Error: dim<- length of dims do not match the length of object Well, it is prototype-compatible. Splus 5.3 does likewise. A way out is t1<-data.frame(unclass(t1)) However, we do seem to have a bug in the area: > t1 <-