similar to: R-alpha: several bugs ...

Two things. TASKS has TASK: Naming with Numeric Values and "unlist" STATUS: Open FROM: <hornik@ci.tuwien.ac.at> R> l <- list("11" = 1:5) R> l $11 [1] 1 2 3 4 5 R> unlist(l) 111 112 113 114 115 1 2 3 4 5 [ Bug or feature ? ] This seems to be a feature (at least, SPLUS does the same), so please close that task. *** However, please add

Attached is an updated version of the CRAN announcement, also to be sent out as soon as 0.50 is out. Again, please have a look if possible. -k ************************************************************************ This is the first announcement of the Comprehensive R Archive Network (CRAN) CRAN is a collection of sites which carry identical material, consisting of the

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

Thanks for fixing formula() and update(). Here are some open problems. * effects: ** There is no help available for effects() ** The lm() help page has effects(lm.obj), but from the code it is clear that the one-arg version cannot work. Shouldn't it simply return z$effects? ** The two-arg version only works for factors (?) ... otherwise, z$model.frame[factors & pattern != 0] gives

[I do think this discussion belongs to R-devel rather than anywhere else .. MM] >>>>> "Kurt" == Kurt Hornik <Kurt.Hornik@ci.tuwien.ac.at> writes: >>>>> Peter Dalgaard BSA writes: >> Kurt Hornik <Kurt.Hornik@ci.tuwien.ac.at> writes: KH>>> While trying to write documentation for data.class(), I came

Dear all, I know there have been various questions posted over the years about loops but I'm afraid that I'm still stuck. I am using Windows XP and R 2.9.2. I am generating some data using the multivariate normal distribution (within the 'mnormt' package). [The numerical values of sanad and covmat are not important.] > datamat <-

Dear R-devel list members, Ben Fairbank draw it to my attention that factanal() (in the stats package) doesn't report factor correlations for oblique rotations. Looking at the source, I see that factanal also doesn't save the factor-transformation (rotation) matrix from which these correlations can be computed. I've modified the source, attached below, so that the transformation

Dear List, First off, my deepest gratitude to the Sweave developers: this tool has improved my quality greatly. A question in my work I use \Sexpr{} statements scalar values and the xtable package for all manner of tables. What I'd like to do is to use a vector inline, rather than a whole separate table. Something like: %%%%%%%%%%%%%%%% begin code % Latex junk % Sweave block:

Hi everyone, I'm fairly new to R, and I don't have a background in statistics, so please bear with me. ;-) I'm dealing with 2^k factorial designs, and I was just wondering if there's any way to analyze more than two factors of a gage R&R study in R. For example, Minitab has an "expanded gage R&R" function that lets you include up to eight additional factors

R version: 1.7.1 OS: Red Hat Linux 7.2 When "covmat" is supplied in princomp(), the output value "center" is all NA's, even though the input matrix was indeed centered. I haven't read anything about this in the help file for princomp(). See code below for an example: pc2$center is all NA's. Jerome Asselin x <- rnorm(6) y <- rnorm(6) X <- cbind(x,y)

In summary.glm I'm trying to get a better feel for the z output. The following lines can be found in the function 1 if (p > 0) { 2 p1 <- 1:p 3 Qr <- object$qr 4 coef.p <- object$coefficients[Qr$pivot[p1]] 5 covmat.unscaled <- chol2inv(Qr$qr[p1, p1, drop = FALSE]) 6 dimnames(covmat.unscaled) <- list(names(coef.p), names(coef.p))

Hi, I have a set of covariance matrices but not the original data. I want to carry out some exploratory factor analysis. So, I am trying to construct a covariance matrix list as the input for factanal. I can construct a list which includes the cov, the centers, and the n.obs. But it doesn't work. I get an error that says "Error in sqrt(diag(cv)) : Non-numeric argument to mathematical

--=-hiYzUeWcRJ/+kx41aPIZ Content-Type: text/plain; charset="UTF-8" Content-Transfer-Encoding: 8bit Hi, on March 10 I filed a wishlist bug report asking for the inclusion of some changes to factanal() and the associated print method. The changes were originally proposed by John Fox in 2005; they make print.factanal() display factor correlations if factanal() is called with rotation =

Dear friends, I'm using R to produce the following Factor Analysis: > matriz.cor<-hetcor(matrix(as.factor(data), ncol=variables, byrow=T))$correlations > factanal(x=data, factors=2, covmat=matriz.cor, scores='regression') Then the screen output shows the following message: Error en factanal(x = data, factors = 2, covmat = matrix, : requested scores without

Hi, Try: dat1<- t(dat[,-1]) ?colnames(dat1) <- dat[,1] covmat <- cov(dat1) A.K. I'm running into this error and I'm not sure how fix it Error: is.numeric(x) || is.logical(x) is not TRUE This is my data frame geneExpression_Lab10.txt This is the homework instructions in case you want to see Lab10.pdf This is my code dat <-read.delim("~/Dropbox/Homework/R

Dear all, I've sent this question 2 days ago and got response from Sarah. Thanks for that. But unfortunately, it did not really solve our problem. The main issue is that we want to use our own (manipulated) covariance matrix in the calculation of the mahalanobis distance. Does anyone know how to vectorize the below code instead of using a loop (which slows it down)? I'd really appreciate

Dear all, We have a data frame x with n people as rows and k variables as columns. Now, for each person (i.e., each row) we want to calculate a distance between him/her and EACH other person in x. In other words, we want to create a n x n matrix with distances (with zeros in the diagonal). However, we do not want to calculate Euclidian distances. We want to calculate Mahalanobis distances, which

On 4 Jan 2003 at 12:51, Wolfgang Lindner wrote: > Please excuse me, if the following questions are *too* off-topic, but I found it > interesting. In inspecting your code I came across an R feature, I could not > find in the online manuals: > > Q1. Looking at the left-handside in your function def: > > "cov.cor" <- function ( covmat ) { >

I'm attempting the heruclean task of teaching myself Introductory Statistics and R at the same time. I'm working through Peter Dalgaard's Introductory Statistics with R, but don't understand why the answer to one of the exercises works. I'm hoping someone will have the patience to explain the answer to me, both in the statistics and R areas. Exercise 6.1 says: The zelazo data

The function tapply is not working in the Windows version of R=20 (Version 0.50 Beta (Sept 29, 1997)) In tapply <- function (x, INDEX, FUN=3DNULL, simplify=3DTRUE, ...)=20 ... The part: if (simplify && all(unlist(lapply(ans, length)) =3D=3D 1)) { ans <- unlist(ans, recursive =3D FALSE) names(ans)<-namelist[[1]] return(ans) } should be replaced by if (simplify