similar to: axis for small ranges (PR#390)

I use R Version 0.60 Alpha (September 18, 1997) on a Linux Pentium (Debian 1.3) and on a Sparc-Sun-Solaris 2.5.=20 R> svd(matrix(1:16,4,4)) =09yields on both machines Error: error 4 in dsvdc R> svd(matrix(1:20,4,5)) =09gives a result on the Linux computer $d [1] 0 0 0 NA $u [,1] [,2] [,3] [,4] [1,] 1 0 0 0 [2,] 0 1 0 0 [3,] 0 0 1 0 [4,] 0

I noticed that apply is VERY SLOW when applied to a "large" dimension as for example when computing the row sums of a matrix with thousands of rows. To demonstrate it, I did some benchmarking for different methods of computing the row sums of an nx10 matrix with n =3D 2000, ..., 10000. The first method (M1) I used is the normal apply command: y <- apply(x,1,sum) The second method

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

The following statements yield a seg.fault: R --vanilla --nsize 500K x <- rep(letters,10000) f <- function(x) {z<-paste("\"",x,"\"",sep=""); z} y <- f(x) Segmentation fault If a turn gcinfo on, I get Garbage collection [nr. 1]... 387529 cons cells free (75%) 3807 Kbytes of heap free (62%) Garbage collection [nr. 2]... 273868 cons cells free

Sirs, may be that a fault exists in the file "postscript" under "r documentation" which contains "Notes on R:A Programming..." at http://www.ci.tuwien.ac.at/R/contents.htm| I did not succeed in reading and printing it. Can You help me? Many thanks in advance. F.Petronio petronio at univ.trieste.it

Dear R-Users and Developpers, Currently R does not include functions to compute the skewness and kurtosis. I programmed it myself in the following way, but probably *real* programmers/statisticians can do that better: mykurtosis <- function(x) { m4 <- mean((x-mean(x))^4) kurt <- m4/(sd(x)^4)-3 kurt } myskewness <- function(x) { m3 <- mean((x-mean(x))^3) skew <-

Andreas Weingessel <Andreas.Weingessel@ci.tuwien.ac.at> writes: > When I tried to port the statlib-Pkg rpart to R, I came across the > c-Function prlab which seems to be part of S+, but not of R. > > Does anyone know what prlab does and how difficult it would be to > emulate it in R? Can you show us the context? If I had to guess I would say it was a routine to print

Hello, I want to try the Design and Hmisc library from FE Harrell but, i have trouble with the copyright for the library Hmisc : in the home of the author : http://fharrell.biostat.virginia.edu/s/unix/ in the file Hmisc.README ###COPYRIGHT NOTICE ###You may not port code in the Hmisc library to R. Aie, Why ???? ###You may distribute these functions freely as long as you do so without

I have the following problem. I have some multidimensional data points "x" and a curve "fit" fitted to these points. How can I combine R> pairs(x) and R> pairs(fit,panel=lines) in one plot, so that I can see how good the curve fits the data? A command like R> pairs(x, panel=function(x,y) {points(x,y); lines(lowess(x,y))}) does not work, since I fit the curve in all

Hi Everyone, This example code results in R 'crashing'; that is the R application closes with no warnings or error messages. #----------------------- myD <- read.table(stdin(), header=TRUE, nrows=20) Broth Salt pH Temp N Y Growth 1 310 9.0 2.92 10 90.0 NA 0 2 615 6.0 7.82 30 1.0 2 1 3 217 2.0 7.34 10 7.0 8

Dear friends, when i use glm() to fit my data, i use glm(formula = snail ~ vegtype + mhveg + humidity + elevation + soiltem, *family = Gamma(link = inverse),* data =a,)) It shows: error in eval(expr, envir, enclos) : *gamma distribution don't allow negative value*. But i use result<-glm(formula = snail ~ vegtype + mhveg + humidity + elevation + soiltem, family = poisson, data =a) #this

On Thu, Jan 31, 2013 at 10:29 AM, Weijia Wang <wwang.nyu@gmail.com> wrote: > Hi, > > > > I have a new question about subsetting in R. > > > > Say we have this data frame: > > > > PT_ID Blood_Pressure OBS_TYPE > > 92 1900 90.0 DBP > > 94 1900 90.0 DBP > > 174 2900 140.0 SBP > > 176 2900

The following example from the help page does not work: > f1_ gl(2,3); f1 > f2_ gl(3,2); f2 > f1:f2 # a factor, the ``cross'' f1 x f2 yields Warning in f1:f2 : Numerical expression has 6 elements: only the first used [1] 1 whereas R-0.61.2 yields > f1:f2 # a factor, the ``cross'' f1 x f2 [1] 1:1 1:1 1:2 2:2 2:3 2:3 Note that the new behavior is consistent with

Hi all, I realize this is asking a lot of lattice, but I want to add a second y axis inside a xyplot and have y1 and y2 have different ranges. Given dat below, I can add a second y axis by overlaying a new plot with par(new=T) and label axis 4 with standard graphics. I've seen an example for doing something similar in xyplot even though Deepayan has indicated that lattice isn't the right

Dear all, I can not find/understand the solution to this from the help pages: Say we have the following script: x<-matrix(c(1,1,0,1,0,1,0,1,1),3,3) image(x) How can I change the ranges on the vertical and horizontal axis to a range of my own or at least place a box frame around the image if I choose to use "axes=FALSE"? Apologies for such a bsic question and thanks beforehand

I have a data set like this: > head(FormatedData) ID Target Actual Date Time 1 2030 0 -57.0 12/20/08 17:03:00 2 2030 90 90.0 12/20/08 18:41:00 3 2030 45 43.8 12/21/08 14:36:00 4 2030 0 -23.8 12/21/08 19:30:00 5 2030 90 90.2 12/21/08 21:48:00 6 2030 45 48.6 12/22/08 13:02:00 I wan to convert the format of Date and Time, so I did this: pdate

[This is something like a bug report; maybe somewhat longish & technical ..] As an introduction, just try the following code (it should work both in R and S). I think it screws up the postscript() driver both for S and R, but this is not the issue here. is.R <- function() { ## returns 'TRUE' iff we are using 'R' exists("version") && !is.null(vl

If x is a data.frame containing a column of mode character, this column is changed to a factor after applying fix(x). Example: R> x_data.frame(a=1:3,b=letters[1:3],c=c("hansi","pepi","karli")) R> x$c_as.character(x$c) R> is.character(x$c) [1] TRUE R> is.factor(x$c) [1] FALSE R> fix(x) ## I change "karli" to "sepp" R> x$c [1]

================================================= 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 had a look at this, and worked out most of a function to do it. There are three things I haven't been able to figure out, however. 1) Getting the last N lines of the command history into a vector of strings. I thought something like: last.commands<-history(N) would work, but no go. There doesn't seem to be a file that I can tail, either. 2) Finding continuation lines by