similar to: Terminology and canonical statistical user literature

I am looking for a function like my.blockwisesum(vector, n) that computes sums of disjoint subsequences of length n from vector and can work with vector lengths that are not a multiple of n. It should give me for instance my.blockwisesum(1:10, 3) == c(6, 15, 24, 10) Is there a builtin function that can do this? One could do it by coercing the vector into a matrix of width n, and then use

Hello all, trying to prepare my first package for submission to CRAN I am stumbling over the "checking examples" step of R CMD check. I have some examples that produce more than one plot. I currently separate those plot calls by readline("Press <Return> for a plot including a density plot") or some such to have R wait before producing the next plot. This works OK

Hello everybody, I am just trying to put together my first own CRAN-able package. I have completed the programming and my code works OK. But now when I perform R CMD check agsemisc_1.0-2.tar.gz and R gets to the step * checking examples ... I receive the following error message: > print(bwplot(Species~Sepal.Length, data=iris, panel=panel.bwstrip)) Error in rep(T, length(x)) : T used

(was: AW: [Rd] R CMD check: checking examples: how to (not) pause execution) Paul, Brian, > > I think the more usual way to solve this is to have the user set > > par(ask=TRUE) is they want to be prompted. thank you. Looking back, 'par' would have been a likely candidate for looking. One reason I did NOT look there is that I need the functionality for print.trellis in

Hello everybody, I have a bunch of strings like this: "IBM POWER4+ 1.9GHz" "IBM RS64-III 500MHz" "IBM RS64-IV 600 MHz" "IBM RS64 IV 750MHz" "Intel Itanium 2 Processor 6M 1.5GHz" "Intel Itanium2 1 Ghz" "Intel Itanium2 1.5GHz" "Intel

Dear Lattice bit-meddlers, while within a panel function for xyplot, how can I find out the values of (effectively) xlim and ylim -- no matter whether they have been set explicitly or chosen by Lattice itself? I have just tried for an hour to find out, with no success whatsoever. I looked in Grid for something that would return such data, but got lost. I looked at tracebacks of panel calls,

Full_Name: Lutz Prechelt Version: 2.2.1 OS: WinXP SP2 Submission from: (NULL) (130.133.8.114) This command plot(0, 0, xlim=c(3, 5), ylim=c(0, 10, 17)) should complain about 'ylim' (because it has three elements). However, it does in fact (in my german error message at least) complain about 'xlim' instead: Fehler in plot.window(xlim, ylim, log, asp, ...) : ung?ltiger

Full_Name: Lutz Prechelt Version: 2.4.1 OS: Windows XP Submission from: (NULL) (160.45.111.67) I stack a number of data.frames using rbind. Each of these dataframes has a column 'authorname', which is a factor and a column author = unclass(authorname) as piecewise pseudonyms. When using rbind to stack these dataframes, R warns about invalid factor levels and inserts all NAs in the author

I'm very new to R and utterly blown away by not only the language but the unbelievable set of packages and the documentation and the documentation standards and... I was an early APL user and never lost my love for it and in R I find most of the essential things I loved about APL except for one thing. At this early stage of my learning I can't yet determine if there is a way to

NIST maintains a repository of Statistical Reference Datasets at http://www.itl.nist.gov/div898/strd/. I have been working through the datasets to compare R's results to their references with the hope that if all works well, this could become a validation package. All the linear regression datasets give results with some degree of accuracy except one. The NIST model includes 11 parameters,

Dear R-devels, I am looking for a package (or some other infrastructure) for the following situation: I am doing a larger data evaluation. I have several dozen data files and do not want to keep the persistent data in the R workspace (for robustness reasons). I have several dozen R files, some for reading and preprocessing data files, others for doing plots or analyses. I will make frequent

Dear Users! I encountered with some problem in data reading while I challenged R (and me too) in a validation point of view. In this issue, I tried to utilize some reference datasets ( http://www.itl.nist.gov/div898/strd/index.html). And the result departed a bit from my expectations. This dataset dedicated to challenge cancellation and accumulation errors (case SmLs07), that's why this

Hi all, I am a freshman of R,but I am interested in it! Those days,I am learning pages on NIST,with url http://www.itl.nist.gov/div898/handbook/eda/section3/probplot.htm, I am meeting a problem about probability plot and I don't know how to plot a data set with R. Could somebody tell me the answer,and a example is the best! I will look forward to your answer. Thank you very much.

I have a large computation that performs various bootstrap resampling loops over pairs of samples using lapply(). After a few seconds it stops like this > do.summarize.effects(compute.power=T) ... Error: Unimplemented feature in rep > .Traceback [[1]] [1] "rep(no, length = length(ans))" [[2]] [1] "ifelse(t.howto == \"n\", boot.result, l.meanboot.test(onew.x,

Hola a todos: Estoy intentado realizar un modelo senoidal de unos datos de radiación solar con el fin de afrontar el relleno de la serie y aplicar la prueba de Thom para verificar su homogeneidad [0]. De momento me encuentro con los siguientes problemas: 1- ¿Existe la prueba de Thom en R? ¿O debo crearme mi propia función? 2- Para la realización del modelo senoidal estoy siguiendo los pasos

-----BEGIN PGP SIGNED MESSAGE----- Hash: SHA1 Hello all - As I progress in R I am trying to automate functions I would have normally farmed out to Excel, SPSS or Statistica. Single factor anova is one of them. For example, a dataset from NIST StRD (http://www.itl.nist.gov/div898/strd/anova/AtmWtAg.html) has two groups: 1 2 107.8681568 107.8681079 107.8681465 107.8681344

Hello. I'm a new user of R, and I have a question regarding the use of aov and lm-functions. I'm doing a fractional factorial experiment at our production site, and I need to familiarize myself with the analysis before I conduct the experiment. I've been working my way through the examples provided at http://www.itl.nist.gov/div898/handbook/pri/section4/pri472.htm

Hello, I am trying to test the precision of R on datasets from The Statistical Reference Datasets Project http://www.itl.nist.gov/div898/strd/index.html and I don't manage to understand how R is storing its results. For example, I calculate a mean on the michelso dataset (100 values) and find: > m=mean(michel) > m V1 299.8524 > print(m,digits=15) V1 299.8524

The function wtd.var(x,w) in Hmisc calculates the weighted variance of x where w are the weights. It appears to me that wtd.var(x,w) = var(x) if all of the weights are equal, but this does not appear to be the case. Can someone point out to me where I am going wrong here? Thanks. Tom La Bone [[alternative HTML version deleted]]

Hi, > [mailto:r-help-bounces at stat.math.ethz.ch] On Behalf Of Corey Powell > > Do you know of any references that verify the accuracy of R > for basic statistical calculations and tests. The results of > these studies should indicate that R results are the same as > the results of other statistical packages to a certain number > of decimal places on some benchmark