similar to: R equivalent for Splus "peaks"

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

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 lower bound in a barplot (or the left bound if you're making a horizontal barplot) is -0.01 no matter what. This causes problems if you're making a barplot of small values (say, < .001), as most of the plot is taken up with blank space beneath the axis, and the bars are squeezed in at the top. The fix seems to be simple: replace lines 41: xlim <- range(-0.01,

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

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

If there something available in R that has the functionality of the S-PLUS peaks() function? Thanks, Mirka

In 2004 there was this R-Help posting from Jan 2004: http://finzi.psych.upenn.edu/R/Rhelp02a/archive/33097.html R equivalent of Splus peaks() function? The peaks function there has worked well for me on a couple of projects, but some code using "peaks" failed today, which had worked fine in the past. I was looking for a peak in a test case that was a sine curve over one

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:

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

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

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

I am looking for a function in R that is equivalent to the function "peaks" in Splus. This function gives the local maxima of a vector. I do not find it using the help and I 'm not sure it exits in R. Could someone help me ? Thanks, Herve Cardot -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read

Hello all I'm trying to do a simulation that involves identifying the minimum point between two peaks of a (usually) bimodal distribution. I can do this easily if there are only two peaks: CnBdens<-density(Ys/Xs) #probability density function for ratio of Ys to Xs for(p in 1:512) ifelse(CnBdens$y[p]>CnBdens$y[p-1],peak1<-p,break) #identifies first peak in probability

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

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