similar to: Strange behaviour with the Homals Package

Dear R users and experts, I am using the package Homals to perform categorical PCA on some survey data. The documentation is a bit obscure, however, I followed the examples. X is a matrix with ordinal 1-5 Likert scale entries on 24 questions from about 2500 respondents. So I run X.nlpc = homals(X, rank = 1, level = "ordinal", ndim = 10) to get the first 10 pc's (the number of

Function homals() is not part of R. It appears to be in a contributed package `homals', but you did not say so. To quote the FAQ Bug reports on contributed packages should be sent first to the package maintainer, and only submitted to the R-bugs repository by package maintainers, mentioning the package in the subject line. I have closed this bug report on R-bugs. On Fri, 2 Jul

Hi All is there "homals" package prepared for win32? kind regards, Valery A.Khamenya --------------------------------------------------------------------------- Bioinformatics Department BioVisioN AG, Hannover

The homals package http://www.jstatsoft.org/v31/i04 will get a major programming overhaul. This will take some time, but what's a few years on a 40-year project. Suggestions from the audience are welcome. homals() has a core loop over the m variables in which 1. tapply is used to compute category quantification (centroids) 2. category quantifications are then adjusted to satisfy the

The homals package http://www.jstatsoft.org/v31/i04 will get a major programming overhaul. This will take some time, but what's a few years on a 40-year project. Suggestions from the audience are welcome. homals() has a core loop over the m variables in which 1. tapply is used to compute category quantification (centroids) 2. category quantifications are then adjusted to satisfy the

Hello The Homals package and its plot options are excellent. However, I am unable to manipulate the colour in the plots. In a call such as: plot(mc_analysis, plot.dim = c(1,3), plot.type = "jointplot", col = 1) this should be straightforward - but I can't seem to effect the plotted colours I have tried various combinations for "col" commands in other plot packages I

homals-0.9.0 is on CRAN -- by Jan de Leeuw and Patrick Mair This package implements the methods discussed in Gifi, Nonlinear Multivariate Analysis, Wiley, 1990. In the Gifi terminology it covers homals, princals, canals, morals, criminals, and overals. The R implementation fills several gaps in Gifi, adding multiple ordinal, numerical, and polynomial data transformations. Differences with

homals-0.9.0 is on CRAN -- by Jan de Leeuw and Patrick Mair This package implements the methods discussed in Gifi, Nonlinear Multivariate Analysis, Wiley, 1990. In the Gifi terminology it covers homals, princals, canals, morals, criminals, and overals. The R implementation fills several gaps in Gifi, adding multiple ordinal, numerical, and polynomial data transformations. Differences with

Dear guRus, I am doing a logistic regression using restricted cubic splines via rcs(). However, the fitted probabilities should be nondecreasing with increasing predictor. Example: predictor <- seq(1,20) y <- c(rep(0,9),rep(1,10),0) model <- glm(y~rcs(predictor,n.knots=3),family="binomial") print(1/(1+exp(-predict(model)))) The last expression should be a nondecreasing

Hey hey kids! I'm facing this rough problem: I have to modify a function in a R package (namely, homals) but I'm not able to do it from the R interface, since this function recalls other functions which looks like invisible... I downloaded the source package from CRAN, http://cran.r-project.org/src/contrib/homals_0.9-10.tar.gz http://cran.r-project.org/src/contrib/homals_0.9-10.tar.gz how

Hi, is there a function in R already doing what I try to do below: # Let 'x' be an array with *any* number of dimensions (>=1). x <- array(1:24, dim=c(2,2,3,2)) ... x <- array(1:24, dim=c(4,3,2)) i <- 2:3 ndim <- length(dim(x)) if (ndim == 1) y <- x[i] else if (ndim == 2) y <- x[i,] else if (ndim == 3) y <- x[i,,] else ... and so on. My current

Full_Name: Greg Kochanski Version: 2.2.1 OS: Debian Linux (testing) Submission from: (NULL) (212.159.16.190) This is really a feature request. When you do mosaicplot() on a data set where the probability of several nearby rows is small, then the labels for those rows are plotted overlapping each other. This situation can be improved by calling mosaicplot() with a large value of

I want to call dll from R but encounter problem in compiling the fortran program. First I try "Rcmd shlib prog.f", it failed and warning: make[1]: `libR.a' is up to date. make: *** No rule to make target `'prog.o', needed by `prog.a'. stop. Then I try to compile it by absoft fortran compiler, it works and produces prog.dll. But when this routine is called in R, it

Is it just my installation or bug in 1.9.0 ? The same thing works fine in 1.8.1 Best regards, Ryszard # R-1.9.0 library(pls.pcr) nr <- 8; ndim <- 2 x <- matrix(rnorm(nr*ndim), nrow=nr) y <- as.matrix(x[,1]) for (i in 2:ndim) y <- y + x[,i] y <- y + rnorm(length(y)) m <- pls(x,y,validation='CV') # Error in paste(ncomp, "LV's") : Argument

I was wondering if it is possible to do the following in a smarter way. I want get the mean value across the columns of a matrix, but I want to do this on subrows of the matrix, given by some vector(same length as the the number of rows). Something like nObs<- 6 nDim <- 4 m <- matrix(rnorm(nObs*nDim),ncol=nDim) fac<-sample(1:(nObs/2),nObs,rep=T) ##loop trough different

R-listers, In using the adapt function, I am getting the following warning: Ifail=2, lenwrk was too small. -- fix adapt() ! Check the returned relerr! in: adapt(ndim = 2, lower = lower.limit, upper = upper.limit, functn = pr.set, Would someone explain what the 'lenwrk' value indicates in order to help diagnose this issue. Also, what are the possible codes for Ifail, so I can set

Hi all I was wondering if someone could help me. I have to estimate some parameters, so I am using the function nlm. Inside this function I have to integrate, hence I am using the function adapt. I don't understand why it is giving the following warnings: At the beginning: Warning: a final empty element has been omitted the part of the args list of 'c' being evaluated was:

Hello, I am trying to call a FORTRAN subroutine from R. The Fortran code is @: http://lib.stat.cmu.edu/apstat/206 It performs a bivariate isotonic regression on a rectangular grid (m X n) matrix. I used the g77 compiler and successfully created a dll file and it also loads successfully from R. But somehow the programs fails to run properly. (I do get the correct result when I compile the

Hi, The man page for base::unique.matrix() and base::unique.array() says that MARGIN is expected to be a single integer. OTOH the code in charge of checking the user supplied MARGIN is: if (length(MARGIN) > ndim || any(MARGIN > ndim)) stop(gettextf("MARGIN = %d is invalid for dim = %d", MARGIN, dx), domain = NA) which doesn't really make sense. As

Hi all, I'm trying to do model reduction for logistic regression. I have 13 predictor (4 continuous variables and 9 binary variables). Using subject matter knowledge, I selected 4 important variables. Regarding the rest 9 variables, I tried to perform data reduction by principal component analysis (PCA). However, 8 of 9 variables were binary and only one continuous. I transformed the data by