similar to: Bug in read.table() (PR#1477)

R 1.7.1 on Windows XP The "normal Q-Q plot" produced by plot.lm() mislabels points when the model is fitted using na.action=na.exclude. Example: x <- 1:50 y <- x + rnorm(50) y[c(5,10,15)] <- NA # insert some NA's y[40] <- 50 # add an outlier plot(lm(y ~ x, na.action=na.omit)) # outlier correctly labeled in all # four plots

I am using R 1.7.1 on Windows XP Pro. The 'overwrite=FALSE' argument appears to be ignored in file.copy(): > file.exists(c("file1.txt", "file2.txt")) [1] TRUE TRUE > file.copy(from="file2.txt", to="file1.txt", overwrite=FALSE) [1] TRUE > and sure enough, file1.txt has been overwritten. Rich Raubertas Biometrics Research, RY33-300 Merck

Dear list, Utils::modifyList() drops NULL components in its second argument, instead of adding them to the first argument. Compare: > modifyList(x=list(A=1), val=list(B=2, C=3)) $A [1] 1 $B [1] 2 $C [1] 3 > modifyList(x=list(A=1), val=list(B=NULL, C=3)) $A [1] 1 $C [1] 3 To me this seems inconsistent with the documentation ("Elements in 'val' which are missing from

Hi, It is my understanding that the eigenvectors of a circulant matrix are given as follows: 1,omega,omega^2,....,omega^{p-1} where the matrix has dimension given by p x p and omega is one of p complex roots of unity. (See Bellman for an excellent discussion on this). The matrix created by the attached row and obtained using the following commands indicates no imaginary parts for the

R 1.6.1 on Windows NT4: The boxplot() function appears to draw its own tick marks and axis values even when called with add=TRUE. As a toy example, try x <- rnorm(100) f <- factor(rep(1:4, each=25)) plot(c(0,4), c(-3,3), type="n", xaxt="n", yaxt="n") boxplot(x ~ f, add=TRUE) My expectation is that a high-level plotting function will not mess with the axes

https://bugzilla.samba.org/show_bug.cgi?id=4855 Summary: skipped local filters cause segfault while deleting (- rFR --delete) Product: rsync Version: 3.0.0 Platform: x86 OS/Version: Linux Status: NEW Severity: normal Priority: P3 Component: core AssignedTo: wayned@samba.org

Hi, I am having some problems mentioned below, the box is in production live environment with traffic around 30 - 100 calls. I am running T/E410P in a Dual P4 xeon with HT disabled. I am using zaptel 0.9.0 and asterisk stable 1 release. There is no gui, just mysql, perl (small script) and asterisk. System runs very smoothly if the calls are around 40-50 and comes one by one , however sometimes

R 3.5.1 on Windows 7 The documentation for 'var' says: "These functions return 'NA' when there is only one observation (whereas S-PLUS has been returning 'NaN'), and fail if 'x' has length zero." The function 'sd' (based on 'var') has similar documentation. However, I get: var(numeric(0)) [1] NA rather than an error. Personally I

Hi, Is it possible to add labels to variables in R (so as to have a better description of what the variables represent)? Thanks, Neil

Dear R-helpers, I have found a slightly annoying problem when trying to plot lines on graphs. I first created my data using tapply, thus:- y1=as.vector(fit1$coef$random$id) x1=tapply(o1,id,median,na.rm=T) x2=tapply(o2,id,median,na.rm=T) #then I plot the data, thus:- plot(x1[x2==0],y[x2==0]) #if I now try to fit the linear regression, R 'hangs up'

I am currently the only use-R at my company, but they are considering buy a more powerful server and letting multiple people use it. They asked me if R supports client-server setups. I know S+ has a server version that does that. I didn't find anything about that on CRAN, but hopefully someone can correct me. I did see some stuff about R web servers

I would like to call C routines from Fortran under linux as suggested in section 5.6 of the "Writing R extensions" documentation. I'm familiar with Fortran but not with C. I understand the example provided in Fortran: subroutine testit() double precision normrnd, x call rndstart() x = normrnd() call dblepr("X was", 5, x, 1) call rndend() end but I don't understand

Hi, Does anyone know if R has the functionality to calculate a simple moving average. I cant seem to find it in the help menu. thanks, Wayne Dr Wayne R. Jones Statistician / Research Analyst KSS Group plc St James''s Buildings 79 Oxford Street Manchester M1 6SS Tel: +44(0)161 609 4084 Mob: +44(0)7810 523 713 KSS Ltd A division of Knowledge Support Systems Group plc Seventh

Assume we have a function like: foo <- function(x, y) how is it possible to define a binary indexing operator, denoted by $, so that x$y functions the same as foo(x, y)

Dear all, For some time now, R is becomming more and more popular in more and more countries. France is for sure one of them, but french people being french one of the obstacle they might tackle is the lack of documentation and support in their native language. To offer this support in french an IRC channel (#Rfr on irc.freenode.net) was created some months ago, beside the official english

Hi guys im having a problem getting R to numerically integrate for some function, say f(bhat)*optimize(G(bhat)), over bhat. Where id like to integrate this over some finite range, so that here as we integrate over bhat optimize would return a different optimum. For instance consider this simple example for which I cannot get R to return the desired result: f <- function(bhat) exp(bhat) g

Duncan Murdoch wrote: > Huntsinger, Reid wrote: > >>Sorry, I was looking at 2.0.1 when I meant to be looking at 2.1.0. The line >>numbers for the latter are 161-179 and line 164 is the one with what I think >>is a typo. >> >>Reid Huntsinger >> >>-----Original Message----- >>From: r-devel-bounces@stat.math.ethz.ch

Dear helpers in this forum, This is a clarified version of my previous questions in this forum. I really need your generous help on this issue. > Suppose I have the following data set: > > id x y > 023 1 2 > 023 2 5 > 023 4 6 > 023 5 7 > 412 2 5 > 412 3 4 > 412 4 6 > 412 7 9 > 220 5 7 > 220 4 8 > 220 9 8 > ...... > Now I want to compute the

Hi, I am trying to run a very computationally expensive procedure in R-2.0.0. and the process always gets killed after approx 8 minutes. This procedure calls some of my own C++ code - in case it was this code causing a memory leak I unload and then reload the .so file every time, however I still get the same problem. The procedure is run 16000 times and always calls the lm() function. My

Could anyone help me with the following issue. Using the GSL library in R I define the following code: ######### library(gsl); S<-function(n) { r<-0:n; ans<-sum(gsl_sf_choose(n,r)*(-1)^r*2^(2*r)*gamma_inc(6-2*r,2)) ans } ######### >SS(10) yields 34.91868 >SS(40) yields 5.340422 >SS(60) yields 180.3162 Doing the same computations in maple I get