similar to: Error in solve.default(V, W) with ProDenICA and the required sample size for stable result

Hi, I'm trying to use fastICA package but I only get an error message like following. > library(fastICA) Error in dyn.load(file, DLLpath = DLLpath, ...) : unable to load shared library '/usr/lib/R/library/fastICA/libs/fastICA.so': /usr/lib/R/library/fastICA/libs/fastICA.so: undefined symbol: sgesdd_ Error: package/namespace load failed for 'fastICA' Here's the

Hi all, Does anyone know how to get the independent components and loadings from an Independent Component Analysis (ICA), as well as principal components and loadings from a Pricipal Component analysis (PCA) using the fastICA package? Or perhaps if there's another way to do ICAs in R? Below is an example from the fastICA manual (http://cran.r-project.org/web/packages/fastICA/fastICA.pdf)

Hi everyone It looks like repeated runs of fastICA produce quite significantly different mixing matrices (not only in terms of sign and row order). I'm not a specialist, so would appreciate any advice on whether this should really be the case: > res3 = > fastICA(af[,2:20],4,alg.typ="parallel",fun="logcosh",alpha=1,method="C",row.norm=TRUE) colstandard >

Hi R useRs: How can I generate a random sample from Maxwell Distribution in C language? Since I know the function rmaxwell {VGAM}, I tried the function r2c {r2c}. But I still cannot find the source code. Thanks for your kind help. RegaRds, 2010-03-12 Randel [[alternative HTML version deleted]]

Hi, Looking at the fastICA library. I want to test separating out different sounds from a recorded wav file. But, I only have a SINGLE stream of data (one channel wav.) It appears as if the fastICA won't let me separate more sounds than I have columns in my data matrix. is this correct? If so, is there a work around that is commonly used? Thanks. -- Noah Silverman UCLA Department of

Hello, I think to use FastICA package for microarray data clusterization, but one question stops me: can I know how much variance explain each component (or all components together) ? I will be very thankful for the help. Thanks, Pavel [[alternative HTML version deleted]]

This is in a 64 bit CentOS 5.6 instance at Amazon AWS with R version 2.14.1 (2011-12-22). It happens on several packages: RMySQL, RODBC, FastICA. Many other packages install just fine. Here's an example error message: * installing *source* package 'RODBC' ... ** package 'RODBC' successfully unpacked and MD5 sums checked sh: ./configure: /bin/sh: bad interpreter:

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

Hello I have a time-comsuming program which need to simplify, I have tested the annotated program as follow: > #define function which will be call > calsta <- function(c, n=100000) + { + i <- seq(from=0, length=c) + logx <- lchoose(NT-n, CT-i) + lchoose(n, i) + logmax <- max(logx) + logmax + log(sum(exp(logx - logmax))) + } > CT=6000 #assignment to CT >

I have a column in a data frame that has a class of "Date" and a mode of "numeric". When I: max(df$Date) My output stays in Date format, i.e. "2006-01-03". However, when I run the following statment: tapply(df$Date, df$SomeFactor, max) my output looks like this: 9129 9493 9861 10226 10591 10956 11320 11687 12052 12417 The returned object is of

Hi Everyone, I want to get the commands history as a character vector instead of just displaying them, but the function history() just returns NULL. I checked the source code of 'history' and could not find a solution. Anybody has an idea? Thanks! P. S. My original problem is, when a user opens a graphics device like png() or pdf(), I want to know the file name used by this device. I

Hello I have 4 1000*1000 matrix A,B,C,D. I want to use the corresponding element of the 4 matrices. Using the "for loop" as follow: E<-o for (i in 1:1000) {for (j in 1:1000) { E<-fisher.test(matrix(c(A[i][j],B[i][j],C[i][j],D[i][j]),2))#call fisher.test for every element } } It is so time-consuming Need vectorization Yours sincerely ZhaoXing Department of

I have been playing around with the FastICA algorithm in R, and I have come across the following issue. I'm new to this topic, so I'm hoping someone will be able to shed some insight. We know that the variance of the Independent Components is constrained to be 1 and they are uncorrelated, so the sum of squares of the coefficients of the mixing matrix (a$A) will give the variance of the

Hi list members, ?"else" tells us In particular, you should not have a newline between '}' and 'else' to avoid a syntax error in entering a 'if ... else' construct at the keyboard or via 'source'. but there's no syntax error when you break the line between "}" and "else" in a function, e.g. f = function(x) {

Hi, I wonder whether there is any convenient function (or package) to extract tables from a HTML page? e.g. from http://www.google.com/finance/historical?q=SHE:002251 I know we can readLines('URL'), gsub('<td>...', '...', source), ... and at last get the numbers; I'm writing to ask whether someone has already contributed a more general function (with the package

Hello Suppose I have a matrix mat=(1:16,2) [,1] [,2] [,3] [,4] [1,] 1 5 9 13 [2,] 2 6 10 14 [3,] 3 7 11 15 [4,] 4 8 12 16 I just want to use the "apply" function to the matrix partially there is a vector end=c(2,3,1,3) #sum the 1st 2 numbers of the 1st column #sum the 1st 3 numbers of the 2nd column #sum the 1st number of the 3rd column

Hi, This is a problem about writing R documentation (R-exts 2.6). The command "\deqn" defined in "Rd.sty" is: \newcommand{\deqn}[2]{\[#1\]} which will put mathematical formulae in the "displaymath" environment; that means line breaks are not allowed (or will not be shown) in formulae, but sometimes we do need multiple lines of formulae. One solution is to write

Hello Suppose I have a matrix mat=(1:16,2) [,1] [,2] [,3] [,4] [1,] 1 5 9 13 [2,] 2 6 10 14 [3,] 3 7 11 15 [4,] 4 8 12 16 there is a vector end=c(2,3,1,3) #coerce the 1st 2 numbers of the 1st column to a vector [1] 1 2 #coerce the 1st 3 numbers of the 2nd column and append it to the previous vector [1] 1 2 5 6 7 #coerce the 1st number of

A colleague suggested that some R software was available for estimating "Independent Components Analysis" (ICA) (... signal separation). If so, I''d appreciate any pointers . Chris Barker Director of Statistical Research MEDTAP International, Inc. Redwood City, Ca www.medtap.com <http://www.medtap.com/> 650 632 4218

Hello I have a 100*100 matrix which is from a intensive computation, e.g. mat. Is there any method/function that return the max of every row and the subscript of maximum value simultaneously #define the function > returnfunction<-function(x){ + value<-apply(x,1,max) + index<-apply(x,1,which.max) + } > mat<-matrix(c(3,5,7,2,1,10,4,3,2),3)#initilize the matrix for test >