search for: arnima

The write.csv() function is currently implemented as function (..., col.names=NA, sep=",", qmethod="double") { write.table(..., col.names=NA, sep=",", qmethod="double") } Surely, it should be function (..., col.names=NA, sep=",", qmethod="double") { write.table(..., col.names=col.names, sep=sep,

Dear r-bugs, I'd like to be able to draw boxplots with solid or dotted whisker lines. The lty argument is currently ignored by boxplot(), and indeed hardwired in bxp(). It's not very difficult to customize these functions to generate various types of boxplots, but generally it might be an idea to allow the user to pass a few more graphical arguments to boxplot(). The approach might as

Dear r-bugs, When I create EPS files, they sometimes appear rotated in my LaTeX PDF document and sometimes they don't. Two examples: ## x1.eps is not rotated in LaTeX x <- seq(-1, 1, length=100) postscript("c:/x1.eps", height=3, width=4, horizontal=FALSE, onefile=FALSE, paper="special") plot(x, dnorm(x), type="l") dev.off() ## x2.eps is not

The audit trail indicates that PR#3663 is not reproducible and/or fixed, but I believe windows(record=T) is still not working as documented: windows(record=T) for(i in 1:5) plot(0, 0, cex=4, pch=as.character(i)) After cycling through the plots with PgUp and PgDn, the user finds out that plot 5 is lost. Indeed, .SavedPlots reports that only 4 plots were saved. I'm running the example in

...se, this is a not-so-great solution to a marginal problem. In the best case, this might be worth a short note in the Journal of Statistical Software. Thanks for your time, Arni P.S. In case the mailing list doesn't handle attachments, I've placed the same files on http://www.hafro.is/~arnima/ for your convenience. -------------- next part -------------- \name{round_multinom} \alias{round_multinom} \encoding{UTF-8} \title{ Round Multinomial Proportions (or Allocate Seats from Election Results) } \description{ This function can round multinomial proportions to a given number of d...

Full_Name: Arni Magnusson Version: 1.6.2 OS: Windows XP Submission from: (NULL) (210.48.49.68) Hi there. I'm experiencing unexpected behaviour from as.POSIXct: > as.POSIXct("1969-12-24") [1] "1969-12-23 23:00:00 New Zealand Standard Time" > as.POSIXlt("1969-12-24") [1] "1969-12-24" > as.POSIXlt("1969-12-24")+1 [1] "1969-12-23

Dear r-bugs, The 'see also' link in help(area,html=T) is broken, see line 91 in library/MASS/html/area.html: "../../integrate/html/integrate.html" should be "../../base/html/integrate.html" Regards, Arni OS: Windows XP R: 1.7.1 MASS: 7.1-8

...nction to compute the summary statistics # ### subset is a logical vector specifying which part of the # ### data to summarize, or TRUE to include all data # ### # ### Author: Arni Magnusson <arnima@u.washington.edu>, inspired by # ### an R-help message from Christophe Pallier # ### # ### Returns: Data frame containing summary statistics # ###...

I'd like to suggest changes to three help pages, regarding the use of pch=NA to suppress plotting symbols. See below. Arni R 2.0.1 on WinXP === help(bxp) === The argument outpch=" " needs to be replaced with outpch=NA in two places. I actually wrote this part of the documentation myself at one point, but have now realized that pch=NA and pch=" " are not the same:

Unlike R 2.1.1, version 2.2.0 generates warnings when an 'oma' argument as passed to pairs(): A <- rnorm(100) B <- rnorm(100) pairs(cbind(A,B)) # no warning pairs(cbind(A,B), oma=c(6,8,10,12)) # warnings in R 2.2.0 I think pairs() should draw the plot quietly, without warnings. Can't see anything in the documentation indicating that an

In Windows XP, the matrix() function crashes the program when 'dimnames' is an empty list: matrix(1:4, nrow=2, dimnames=list()) # R has encountered a problem and needs to close ... This bug is specific to WinXP, as Linux64 handles this situation more gracefully: matrix(1:4, nrow=2, dimnames=list()) Error in matrix(1:4, nrow = 2, dimnames = list()) : invalid type

The weighted.mean() function replaces NA values with 0.0 when the user specifies na.rm=TRUE: x <- c(101, 102, NA) mean(x, na.rm=TRUE) # 101.5, correct weighted.mean(x, na.rm=TRUE) # 67.66667, wrong weighted.mean(x, w=c(1,1,1), na.rm=TRUE) # 67.66667, wrong weighted.mean(x, w=c(1,1,1)/3, na.rm=TRUE) # 67.66667, wrong The weights are

Dear r-bugs, I'm having problems loading the 'grid' package when it has been detached earlier in the same session: > library(grid) > detach("package:grid") > library(grid) Error in .Call("L_initGrid", PACKAGE = "grid") : .Call function name not in load table Error in library(grid) : .First.lib failed This bug was mentioned on r-help

One descriptive statistic that is conspicuously missing from core R is the statistical mode - the most frequent value in a discrete distribution. I would like to propose adding the attached 'statmode' (or a similar function) to the 'stats' package. Currently, it can be quite cumbersome to calculate the mode of a distribution in R, both for experts and beginners. The lack of a

The model frame shows the response and predictors in a data frame with nicely labelled columns: fm <- lm(wt~qsec+log(hp)+sqrt(disp), data=mtcars) model.frame(fm) # ok When the left hand side consists of more than one response, those response variables still look good, inside a matrix: fm <- lm(cbind(qsec,hp,disp)~wt, data=mtcars) model.frame(fm)[[1]] # ok A problem arises when

...stringhash{$s}) .. I wonder if it might be more efficient using regular expression instead, if ($s =~ /^string[1-3]$/) { ... it will depends on how 'regular' the strings are, of course. Michael --__--__-- Message: 9 Date: Fri, 21 Feb 2003 21:06:33 -0800 (PST) From: Arni Magnusson <arnima@u.washington.edu> To: ripley@stats.ox.ac.uk cc: r-devel@stat.math.ethz.ch, <R-bugs@biostat.ku.dk> Subject: Re: [Rd] POSIX problem in New Zealand (PR#2570) On Fri, 21 Feb 2003 ripley@stats.ox.ac.uk wrote: > What exactly is the problem? Those appear to be the same time, in > diffe...

I have been reading about autocorrelation in linear models over the last couple of days, and I have to say the more I read, the more confused I get. Beyond confusion lies enlightenment, so I'm tempted to ask R-Help for guidance. Most authors are mainly worried about autocorrelation in the residuals, but some authors are also worried about autocorrelation within Y and within X vectors

...part of the R base package. The tasks are neither specific nor marginal, so rather than creating one more 'misc' package, I would be happy if the R Development Core Team would adopt these functions and hammer them into shape. The functions are available from http://students.washington.edu/arnima/s, and a short demo session below explains their behaviour. Sorry about the length of this message, but it should be quick read. > env() Environment Objects Kb 1 .GlobalEnv 19 982 2 package:lme4 180 992 3 file:f:/gnu/h...

Creating a plot with reverse Y axis is easy enough with traditional graphics: > x <- 1:3 > y <- 1:3 > plot(y~x, ylim=c(3,1)) +-----------+ 1 | o | | | 2 | o | | | 3 | o | +-----------+ 1 2 3 But xyplot doesn't grasp my idea: > xyplot(y~x, ylim=c(3,1)) +-----------+ 3 | o | | | 2 | o |

I'd like to iterate my earlier request (#7737) to change the documentation for bxp(). The argument outpch=" " needs to be replaced with outpch=NA in two places. I actually wrote this part of the documentation myself at one point, but have now realized that pch=NA and pch=" " are not the same: x <- split(rlnorm(26e4), letters) ## NA generates small file