similar to: R-alpha: R 0.50.a3 family.lm()

A few problems in R-0.50-a3 (which were also in R-0.49): 1) y <- c(10,11,12,13,14,NA,NA) n <- length(y) missed <- (1:n)[is.na(y)] notmissed <- (1:n)[!is.na(y)] blocks <- cut(missed,breaks=c(0,notmissed,n+1)) a <- function(v) { q <- range(v) c(q[1]-1,q[2]+1) } brackets <- tapply(missed,blocks,a) This codes gives the following in S: > brackets $"0+ thru 1":

We are having problems building R on Red Hat Linux 6.0/Sparc. R compiles with no problem but then it fails "make tests". Specifically, the line eigen(cbind(-1,c(1:2,0),0:2)) # complex values causes a segmentation fault. Previous versions of Red Hat Linux did not give this problem. The eigen() function is one of the rare cases in the base library where Fortran code is called via the

Not sure if I already reported this ... it seems that as of 0.50.a2, the adj graphics parameter is completely ignored. Try e.g. plot(1:10) title(main = "This is a title", adj = 0) title(main = "This is a title", adj = 1) No difference ... -k =-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=- r-devel mailing list -- Read

This is a bug both in R and S-plus (in R it's "worse"), I think. Look at this : ma <- cbind(1,1:8); logi <- rep(c(T,F),4) ma[logi,] ##> [,1] [,2] ##> [1,] 1 1 ##> [2,] 1 3 ##> [3,] 1 5 ##> [4,] 1 7 ##-- now the same with [ , drop = FALSE] ## R (-0.50-a3) : ma[logi, drop = F] ##>> Error: invalid subscript type ##

Not really important, but ... `make tests' fails if a user's .Rprofile asks for libraries which are not found by the binary in the compilation tree (yes I know, we're here perhaps the only ones who use private libraries ...). Perhaps we should have options --no-site-file --no-init-file which inhibit the loading of a global and a user's init file? -k

Full_Name: Ricardo Luiz de Andrade Abrantes Version: 2.1.1 OS: Debian Linux, kernel 2.6.8 Submission from: (NULL) (201.6.83.153) The problem can be well explained with the following example: Suppose I made a program in fortran, and a C interfacece to it. Now I want to use this C interface in R to call my fortran program. Then I modified my C file to deal with SEXPs and compile it as a shared

It is very wierd... Can some of you confirm the following behavior ? It is a new bug (feature ?) which was not yet in 0.49 ... noquote <- function(obj) { ## constructor for a useful "minor" class if(!inherits(obj,"noquote")) class(obj) <- c(class(obj),"noquote") obj } "[.noquote" <- function (x, subs) structure(unclass(x)[subs], class =

hi all i have a question that about the eigen analysis found in R and in eviews. i used the same data set in the two packages and found different answers. which is incorrect? the data is: aa ( a correlation matrix) 1 0.9801 0.9801 0.9801 0.9801 0.9801 1 0.9801 0.9801 0.9801 0.9801 0.9801 1 0.9801 0.9801 0.9801 0.9801 0.9801 1 0.9801 0.9801 0.9801 0.9801 0.9801 1 now > svd(aa) $d [1] 4.9204

.. R is now configured for sparc-sun-solaris2.6 Source directory: . Installation directory: /unsup/R-0.64.2 C compiler: gcc -g -O2 FORTRAN compiler: g77 -O2 Gnome support: no $ gcc --version egcs-2.91.66 $ g77 --version GNU Fortran 0.5.24-19981002 "make check" later fails. The failure is in the test of the eigen() function. The

Batch: Putting q(save=F) at the end of my file does not work in my context because I can no longer source the file without quitting. I have that quit statement in my .First so that I always quit that way interactively. The problem is that it is ignored in batch. eigen: The crash occurs on my 586 running Red Hat Linux 2.0.27 but not on my son's 486 running SLackware Linix 2.0.29. We both

Hi, with small matrices eigen works as expected: > eigen(cbind(c(1,4),c(4,7)), only.values = TRUE) $values [1] 9 -1 $vectors NULL > eigen(cbind(c(1,4),c(4,7))) $values [1] 9 -1 $vectors [,1] [,2] [1,] 0.4472136 -0.8944272 [2,] 0.8944272 0.4472136 > eigen(cbind(c(1,-1),c(1,-1))) $values [1] -3.25177e-17+1.570092e-16i -3.25177e-17-1.570092e-16i $vectors

Hi, All: What tools exist for diagnosing singular gradient problems with 'nls'? Consider the following toy example: DF1 <- data.frame(y=1:9, one=rep(1,9)) nlsToyProblem <- nls(y~(a+2*b)*one, DF1, start=list(a=1, b=1), control=nls.control(warnOnly=TRUE)) Error in nlsModel(formula, mf, start, wts) : singular gradient matrix at initial

Dear All i have a little puzzle about eigenvector in the R. As we know that the eigenvector can be displayed on several form. For example A=matrix(c(1,2,4,3),2,2) if we want to get the eigenvalue and eigenvector, the code followed eigen(A) $values [1] 5 -1 $vectors [,1] [,2] [1,] -0.7071068 -0.8944272 [2,] -0.7071068 0.4472136 however, we also can calculate the vector matrix

Hi list, It appears there is a problem with the OpenMPI I/O library on CentOS 6.2 & 6.3 (package openmpi-1.5.4-1.el6.x86_64). When I compile the attached program it ends up in the error path since MPI_File_open returns 16. The corresponding (unhelpful) message is: MPI_ERR_OTHER: known error not in list I couldn't find any pointers on the net and the same program works with OpenMPI

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

Hi, I'm trying to compile two functions into one function. the first funtion is called 'fs' which is self-made function, another function is from the built-in 'integration' function that is copy-paste-edited. If built separatey, these functions work well. However that is not the case if combines together, where certainly I made mistake somewhere when constructing the code.

This message is in MIME format --_=XFMail.1.1.p0.Linux:970825095458:252=_ Content-Type: text/plain; charset=us-ascii Here is a patch which fixes some problems with time series functions. Some examples of what goes wrong... > x <- ts(rnorm(100),start=1,deltat=2) > start(x) Error in ts[1] : object is not subsettable > end(x) Error in ts[2] : object is not subsettable > y <-

I've only been using R on and off for 9 months and started using the glasso package for sparse covariance estimation. I know the concept is to shrink some of the elements of the covariance matrix to zero. However, say I have a dataset that I know has some underlying "baseline" covariance/correlation (say, a value of 0.3), how can I change or incorporate that into to the

Hello everyone, i just updated from 4.2.1 to 4.2.2. If i try to fire up my win2k8 domU with a pci device attached, the dom0 machine hardcrashes. my system log (idrac) shows the following: CPU 2 has an internal error (IERR). A bus fatal error was detected on a component at bus 0 device 0 function 0. CPU 1 machine check detected. and plenty of other entries. The machine hardresets then. If i

Problems in R but not in S: --------------------------- 1) 'unlist' seems to have several other problems than the ones reported up to now. For instance, 'unlist' can be used on almost any object in S without much trouble. Eg.: S> unlist(c(2)) [1] 2 S> --- R> unlist(c(2)) Segmentation fault (core dumped) This occurs in R-0.49 and in R-0.50-a1. 2) Problem with the