similar to: Canonical variates

Dear R-users, I have a two - dimensional data set which needs to be clustered into groups: I'm searching for groups of points which show a positive correlation (in a twodimensional plot of the data set), but I do not have any knowledge about how many groups there might be. Do you know of a clustering algorithm in R (or in general) which can use a-priori information about the cluster's

Hi, I have been using some of the functions in r for classification purposes, chiefly lda, qda, knn and nnet. My problem is that the only one I can figure out how to represenent graphically is lda (using plot.lda). I have tried 'fooling' this function into accepting qda input for plotting but to no avail. I wonder if you have any suggestions? Thanks alot, Anne Marie Power Marine lab.

Hello! I need to classify a data set with 19 variables and 9 classes using Logistic Regression(on R). I know that when we have only 2 classes we can use glm() to estimate the coefficients of the model. But I don?t understand how can I do a classification task with Logistic Regression on a data set with 9 classes! Does anybody know how can I estimate these coefficients (of a model with 9

Hi, I am using the package PTAK and in particular the command Candpara to perform the Parafac factorizationor of a tensor. The results are not encouraging as I expected, I'm starting a phase of analysis to see if there are errors. I pose a question and I hope you can help me. The command to run the factorization is: ## CANDECOMP/PARAFAC results<- CANDPARA(data_matrix, dim=3)

I'm working on a rather interesting consulting problem with a client. A number of physical variables are measured on a number of cricket bowlers in the performance of a delivery. An example variable might be a directional component of angular momentum for a particular joint measured at a large number (101) of equally spaced timepoints. Each bowler generates a (fairly smooth) curve for

Hello! I need to test the adjustment of a (Negative) Exponential Distribution to a dataset. The parameter of the distribution is unknown. What is the appropriate test to do? I've tried the ks.test, although I think this isn't the appropriate one, as I don't know the parameter. Can anybody help me? Thanks in advance, Janete -- Janete da Silva Borges janeteborges at gmx.net Ab

PTAk_1.1-1 ( Principal Tensor Analysis on k modes) has been released on CRAN A multiway method to decompose a tensor (array) of any order, as a generalisation of SVD also supporting non-identity metrics and penalisations. 2-way SVD with these extensions is also available. The package includes also some other multiway methods: PCAn (Tucker-n) and

PTAk_1.1-1 ( Principal Tensor Analysis on k modes) has been released on CRAN A multiway method to decompose a tensor (array) of any order, as a generalisation of SVD also supporting non-identity metrics and penalisations. 2-way SVD with these extensions is also available. The package includes also some other multiway methods: PCAn (Tucker-n) and

Hello R-helpers, I need to generate standard variates normal to 'create' chi-squared variates. To make you more understand, (1)   a<-rnorm(3,0,1) *after do (1), I need to squared and summed the three values. My problem is, how am I going to continue the programming if I had to repeat the process for 15 times, which in the end I will get 15 values from the whole programme.Hope you can

I have been trying to use the rbeta() function to generate random beta variates with a particular mean and standard deviaiton for some simulation problems I am working on, however I am unsure of the relationship between the parameters shape1 and shape2 and the mean and standard deviation of the beta distribution. I am guessing that there is a formula for this, but I am unsure what it is. Can

Hi list, not really a problem to solve but a strange unwanted behavior (at least I do not like it): I'm loading library(multiv) to use sammon. Later I load library(MASS) which tells me that sammon is masked now, because there is also a MASS version of sammon. Since I believe that the MASS version is at least as good as the multiv version, I would also use the MASS version, but help(sammon)

Hi, I have run rpart to construct a regression tree. Is there any simple method to draw a nice picture of it, as it is usually done in books and paper to visualize the tree? Thank you, Christian *********************************************************************** Christian Hennig Fachbereich Mathematik-SPST/ZMS, Universitaet Hamburg hennig at math.uni-hamburg.de,

Hi, exist a special cluster-analysis algorithms which can work with NA's. a further "problem" is that i want cluster variables not cases to identify special variable-set's. Is it a common way turn the data.frame and use kmeans,because this works with NA's, or have anybody another method for finding "variable-sets" , with exception of factor analysis. thanks for

Hello list, I try to optimize some R code and it turns out that the function as.matrix takes a lot of time in my code. There is only one reason why I need as.matrix: My code should work for matrices with ncol=p and p should be allowed to be 1 or larger. Now I have a matrix x with dim(x)=(n,p), and I need to work with the submatrix y <-x[gv,], (gv a vector of n logicals) and to calculate

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

I've been asked about "Canonical Variance Analysis" (CVA). I don't see any reference to it searching the R site. Does it go by other names? Genstat describes it thus: Canonical variates analysis operates on a within-group sums of squares and products matrix, calculated from a set of variates and factor that specifies the grouping of units. It finds linear combinations of the

Hello list, I want to apply svm from library e1071, and I want to supply class weights. I do not really understand the help entry (and there is no example) class.weights: a named vector of weights for the different classes, used for asymetric class sizes. Not all factor levels have to be supplied (default weight: 1). All components have to be named. I have two

Dear group, I would like to document more than one function on one help page using the Rd language. I tried something like \value{ \code{foo1} returns a list with components \item{arg1}{Argument 1} \item{arg2}{Argument 2} \code{foo2} returns a list with components \item{arg3}{Argument 3} \item{arg4}{Argument 4} } Unfortunately this makes the text "\code{foo2} returns a list

Hi, I want to do something like this: par(mfrow=c(7,7)) for (i in 1:49) hist(RATDACOM[SUBJNO==i],breaks=0.5+(0:6), main="",xlab="",ylab="",xaxt="n",yaxt="n") (Don't think about what RATDACOM and SUBJNO are.) I get an error Error in plot.new() : Figure margins too large. 36 histograms with mfrow=c(6,6) work. But the 36 histograms

Hi list, is there any direct way to obtain confidence intervals for the regression slope from lm, predict.lm or the like? (If not, is there any reason? This is also missing in some other statistics softwares, and I thought this would be quite a standard application.) I know that it's easy to implement but it's for explanation to people who faint if they have to do their own programming...