similar to: mva :: prcomp

Dear R Team, Some minor typos in help pages for package:mva Thank you. Rashid Nassar 1. help(kmeans) Details: [k-means? not sure about this] v The data given by `x' is clustered by the k-Means algorithm. When this terminates, all cluster centres are at to the mean of their

I am using R mainly for multiple linear regression and principal components analyses and I am quite happy with it. I think it is a worderful work you have done. But I am having a problem when using Varimax for rotating loadings obtaines from a princomp (I use the cor = TRUE option): the variance across the components in the variables does not mantain. Let me explain it with an example in which

Hi All Not sure if this a bioconductor question or general R mailing list so apologies if this has gone to the wrong one................. When plotting dendrograms created by hclust you can "identify" clusters by clicking on the graphics and returning a list of what is contained in each cluster. However I'd like to be able to "zoom in" on specific clusters and plot

Hello, I am using the following code to plot an MVA plot. library(affy) library(Biobase) library(limma) library(gcrma) pd<-read.phenoData("Clk.targets.2.txt",header=TRUE, row.names=1,as.is=TRUE,sep="\t") Data <- ReadAffy(filenames=pData(pd)$FileName,phenoData=pd) Print(Data) eset <- gcrma(Data) write.exprs(eset,

Per the discussion about the problems with prcomp() when n << p, which boils down to a problem with svd() when n << p, here is a patch to prcomp() which substitutes La.svd() instead of svd(). -Greg (This is really a feature enhancement, but submitted to R-bugs to make sure it doesn't get lost. ) *** R-1.6.0/src/library/mva/R/prcomp.R Mon Aug 13 17:41:50 2001 ---

Dear List, although the library() command tells me that the pcakage "mva" is installed on my machine, I cannot use its functions or get help() about them. And, strange enough, I never installed the package manually. Has it become a part of R-base? I can't find the package among the package sources on CRAN, so I can't (re)install it manually. I'm using R 1.4.0. Can anyone

Full_Name: Laurent Gautier Version: 1.3.0-patched OS: IRIX 6.5 Submission from: (NULL) (130.225.67.199) Hello, The function La.eigen, called by cmdscale in the package mva behaves an unexplicable way (for me). The following lines show what happened. I tried the very same on linux, and it worked fine. >a <- matrix(c(1,2,3,2),3,3) >a [,1] [,2] [,3] [1,] 1 2 3 [2,]

Hello! I have a question about the factanal function. This function returns at the end test statistics like this: Test of the hypothesis that 4 factors are sufficient. The chi square statistic is 4.63 on 2 degrees of freedom. The p-value is 0.0988 Is it possible to get the chi square statistic and the p-value as variables, not the text on the screen? An object of class "factanal"

Vince Carey found this (thank you!). Since the fix to the problem is not entirely obvious, I post this to R-devel as RFC: help(dist) says: >> Missing values are allowed, and are excluded from all computations >> involving the rows within which they occur. If some columns are >> excluded in calculating a Euclidean, Manhattan or Canberra >> distance, the sum is

Hi All: Package TSERIES is stated to depend on MVA. However, there is no MVA package to be found under the list of package sources. Best wishes, ANDREW tseries: Package for time series analysis Package for time series analysis with emphasis on non-linear and non-stationary modelling Version: 0.7-6 Depends: ts, mva, quadprog Date: 2001-08-27 Author: Compiled by Adrian

Inicio del mensaje reenviado: > De: Arnau Mir <arnau.mir@uib.es> > Fecha: 14 de noviembre de 2011 13:24:31 GMT+01:00 > Para: Martin Maechler <maechler@stat.math.ethz.ch> > Asunto: Re: [R] How to compute eigenvectors and eigenvalues? > > Sorry, but I can't explain very well. > > > The matrix 4*mp is: > > 4*mp > [,1] [,2] [,3] > [1,]

Full_Name: Pierre Legendre Version: 2.1.1 OS: Mac OSX 10.4.3 Submission from: (NULL) (132.204.120.81) I am reporting the mis-behaviour of the function 'eigen' in 'base', for the following input matrix: A <- matrix(c(2,3,4,-1,3,1,1,-2,0),3,3) eigen(A) I obtain the following results, which are incorrect for eigenvalues and eigenvectors 2 and 3 (incorrect imaginary portions):

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

Hello. Consider the following matrix: mp <- matrix(c(0,1/4,1/4,3/4,0,1/4,1/4,3/4,1/2),3,3,byrow=T) > mp [,1] [,2] [,3] [1,] 0.00 0.25 0.25 [2,] 0.75 0.00 0.25 [3,] 0.25 0.75 0.50 The eigenvectors of the previous matrix are 1, 0.25 and 0.25 and it is not a diagonalizable matrix. When you try to find the eigenvalues and eigenvectors with R, R responses: > eigen(mp) $values [1]

I was surprised by this (in R 2.0.1): > a <- ordered(-1:1) > a [1] -1 0 1 Levels: -1 < 0 < 1 > model.matrix(~ a) (Intercept) a.L a.Q 1 1 -7.071068e-01 0.4082483 2 1 -9.073800e-17 -0.8164966 3 1 7.071068e-01 0.4082483 attr(,"assign") [1] 0 1 1 attr(,"contrasts") attr(,"contrasts")$a [1]

How can I get the result of, e.g., poly(1:3. degree=2) to give me the unnormalized integer coefficients usually used to explain orthogonal polynomial contrasts, e.g, -1 1 0 -2 1 1 As I understand things, the columns of x^{1:degree} are first centered and then are normalized by 1/sqrt(col sum of squares), but I can't see how to relate this to what is returned by poly(). >

For ordered factor the natural contrast coding would be to parametrize by the succsessive differences between levels, which does not assume equal spacing of factor levels as does the polynomial contrasts (implicitly at least). This requires the contr.cum, which could be: contr.cum <- function (n, contrasts = TRUE) { if (is.numeric(n) && length(n) == 1) levs <- 1:n

Dear all, I do not seem to grasp how contrasts are set for ordered factors. Perhaps someone can elighten me? When I work with ordered factors, I would often like to be able to reduce the used polynomial to a simpler one (where possible). Thus, I would like to explicetly code the polynomial but ideally, the intial model (thus, the full polynomial) would be identical to one with an ordered factor.