search for: infel

Dear list, Please find below an example of odd behaviour of the by function. It occurs both under GNU/Linux R 2.6.2 and Windows R 2.7.0alpha. Respective sessionInfo()'s are given below. I hope I do not overlook anything. testFactor <- factor(sample(LETTERS[1:6], size = 42, replace = TRUE)) testMatrix <- matrix(rnorm(42 * 6), nrow = 42) testDf <- as.data.frame(testMatrix) by(data

...nteractive than in "R CMD check", at least under both Fedora 13 and Windows 7 with Windows 7 sessionInfo() copied below: In interactive, the result is c('A', 'b', 'C'); with R CMD check, it is c('A', 'C', 'b'). This produced the infelicity of a bug in "R CMD check" that I could not replicate with interactive R because a *.Rd file contained the equivalent example of stopifnot(all.equal(sort(c('A', 'b', 'C')), c('A', 'b', 'C'))): It worked just fine interactively but f...

...rmat(as.character(c(1:2, NA)), na.encode=FALSE), na.print ="") Everything about this is documented (if you look carefully enough), but IMO it violates the principle of least surprise https://en.wikipedia.org/wiki/Principle_of_least_astonishment , so I would call it at least an 'infelicity' (sensu Bill Venables) Is there any chance that this design decision could be revisited? cheers Ben Bolker --- Consider dd <- data.frame(f = factor(1:2), c = as.character(1:2), n = as.numeric(1:2), i = 1:2) dd[3,] <- rep(NA, 4) print(dd, na.print = "")...

...lt;- function(x) stop("not this one") UseMethod("mycoef") } so I used ENCLOS(env). This adds some overhead, hopefully no more than searching for methods. BTW, I noticed that R_LookupMethod uses findVar, that is looks for any object not for functions: that must be another infelicity. -- Brian D. Ripley, ripley at stats.ox.ac.uk Professor of Applied Statistics, http://www.stats.ox.ac.uk/~ripley/ University of Oxford, Tel: +44 1865 272861 (self) 1 South Parks Road, +44 1865 272866 (PA) Oxford OX1 3TG, UK Fa...

utils::citation() offers the option auto=(NULL|TRUE) to obtain the citations for a package from the package DESCRIPTION file or a CITATION file, as described below: |auto| a logical indicating whether the default citation auto-generated from the package 'DESCRIPTION' metadata should be used or not, or |NULL| (default), indicating that a 'CITATION' file is used if it exists, or

In theory, dlnorm(x, ...) == dnorm(log(x), ...)/x, x>0. Unfortunately, when sd=0, dlnorm and plnorm return NaN, while dnorm returns (if(x != mean)0 else Inf) and pnorm returns (if(x<mean)0 else 1). [A numerical optimization, maxLik{maxLik}, reported the NaNs for me.] help('dlnorm') says, "dlnorm is calculated from the definition (in ?Details?). [pqr]lnorm are based on the

Dear list, "+" (and "-") being defined for difftime class, I expected mean() to return something sensible. This is only half-true : > mean(c(1:5, 5:1),na.rm=TRUE) [1] 3 > mean(as.difftime(c(1:5, 5:1),unit="mins"),na.rm=TRUE) Time difference of 3 mins Fine so far. However : > mean(c(1:5, NA,5:1),na.rm=TRUE) [1] 3 > mean(as.difftime(c(1:5,

Hello, The main reason for resetting the tagbuf in `print.default()` and other entry points to the print routine is that it is currently not reset on exit. Creating a context to reset it on exit to its last value might work. This should be done in the entry points rather than in print-value-rec though, since callers of the latter might write to the tagbuf. Another solution to this problem is

In using cor(data.frame), it is annoying that you have to explicitly filter out non-numeric columns, and when you don't, the error message is misleading: > cor(iris) Error in cor(iris) : missing observations in cov/cor In addition: Warning message: In cor(iris) : NAs introduced by coercion It would be nicer if stats:::cor() did the equivalent *itself* of the following for a data.frame:

In a paper I'm writing using Sweave, I make use of lattice graphics, but don't want to explicitly show (or explain) in the article text the print() wrapper I need in code chunks for the graphs to appear. I can solve this by including each chunk twice, with different options, as in <<ortho-xyplot1-code, keep.source=TRUE, eval=FALSE>>= library(nlme) library(lattice)

Currently S3 objects nested in generic vectors cause the tag buffer to be reset.? This feels sub-optimal for those objects that don't have a print method: > list(a=list(b='hello')) $a $a$b???????????????? ### <<<< notice "$a$b" [1] "hello" > list(a=structure(list(b='hello'), class='world')) $a $b?????????????????? ###

...uot;Read error" *when the chunk name had accented characters in it* (my system works in UTF8). Plain ASCII chunk names work OK. I do not know ODF specification well enough to be sure that this is an odfWeave problem (some chunk name mangling might be necessary) rather than an OpenOffice "infelicity" (non-implementation of a specification-requested feature) : checking with another ODF-compliant word processor might be useful, but I don't have any on hand. In any case, this problem is not obvious after re-reading the odfWeave documentation, and does not seem to have been already...

Hi all, I'd like to discuss a infelicity/possible bug with gsub. Take the following function: f <- function(x) { gsub("\u{A0}", " ", gsub(" ", "\u{A0}", x)) } As you might expect, in utf-8 locales it is idempotent: Sys.setlocale("LC_ALL", "UTF-8") f("x y")...

Hi In the following, the graph I see on the screen and the .png output coincide. However, in the .pdf file, the fonts seem to be scaled fairly larger, resulting in the label for the top legend disappearing. Is this an infelicity or bug, or is there something I've missed? More generally, how do I control the size of fonts used in legends and axis labels? library(car) library(ggplot2) qp <-qplot (education , income , shape=type , size=women , colour=prestige , xlab="Education" , ylab="Income&q...

...t; > traceback() >2: log(a) >1: f("A") Oops, I get the same behaviour as you now. I think I must have edited the transcript before and changed the order of the lines: now I get the "no traceback" message only if the try(f("A")) comes first. And notice this infelicity: > f("A") Error in log(x) : Non-numeric argument to mathematical function > traceback() 2: log(a) 1: f("A") > try(f("B")) Error in log(x) : Non-numeric argument to mathematical function > traceback() 2: log(a) 1: f("A") # Th...

I've discovered an infelicity (I guess) in qbeta(): it's not a bug, since there's a clear warning about lack of convergence of the numerical algorithm ("full precision may not have been achieved"). I can work around this, but I'm curious why it happens and whether there's a better workaround -- i...

Hi, I have sparse (tridiagonal) matrices, and I use the Matrix package for handling them. For certain computations, I need to either set the first row to zero, or double the last row. I find that neither operation preserves sparsity, eg > m <- Diagonal(5) > m 5 x 5 diagonal matrix of class "ddiMatrix" [,1] [,2] [,3] [,4] [,5] [1,] 1 . . . . [2,] . 1

...uot;pkgType") [1] "both" I think this is somehow related to the fact that install.packages() will now tell you if a more current source package exists and ask if you want to install from source. But this would require RTools to be installed to work. IMHO, I consider this a major infelicity or design flaw, that will cause problems for naive and even experienced users. install.packages("foo") should always work unless package "foo" cannot cannot be found in getOption("repos") -Michael

Hi useRs, A recent coding infelicity along these lines yielded a corrupt data frame. foo <- matrix(1:12, nrow = 3) bar <- data.frame(foo) bar$NewCol <- foo[foo[, 1] == 4, 4] bar lapply(bar, length) > foo <- matrix(1:12, nrow = 3) > bar <- data.frame(foo) > bar$NewCol <- foo[foo[, 1] == 4, 4] > b...

Full_Name: Anne CORI Version: 2.6.1 OS: windows Submission from: (NULL) (134.157.220.19) let A be a mtrix with n lines and p columns, p>36 ; matplot(A) will only draw the plot of the first 36 columns of A. This is due to a bug in the definition of pch in the function : pch <- c(paste(c(1:9, 0)), letters)[1:k] should be replaced by pch <- c(paste(c(1:9, 0)), letters) I hope this